鄒星晨,王欣苗,左亞凡,張澤鑫,賀康寧,*

1 北京林業(yè)大學(xué)水土保持學(xué)院,水土保持國(guó)家林業(yè)局重點(diǎn)實(shí)驗(yàn)室,北京 100083 2 北京市水土保持工程技術(shù)研究中心,北京 100083 3 水利部綜合事業(yè)局,北京 100053

新世紀(jì)以來,全球變暖、環(huán)境惡化和人類活動(dòng)的加劇導(dǎo)致了生物多樣性的銳減[1]。生物多樣性的保護(hù)已成為森林生態(tài)領(lǐng)域的研究熱點(diǎn)。林下植被作為森林生態(tài)系統(tǒng)一個(gè)重要組成成分,在維持生態(tài)系統(tǒng)物種多樣性方面有不可忽視的作用[2]。有研究表明:林下草本雖然生物量占比不高,但其物種多樣性卻能占據(jù)森林生態(tài)系統(tǒng)的90%[3]。影響草本多樣性的因素眾多,近年來,有關(guān)草本多樣性的研究多集中在分析氣候、林分結(jié)構(gòu)和土壤因子對(duì)草本多樣性的影響,并比較不同影響因素下草本多樣性的差異。這種僅僅分析單個(gè)因子對(duì)草本多樣性的研究,難以揭示影響因子之間的耦合關(guān)系,難以全面反映影響草本多樣性的機(jī)理。有研究表明:林下草本豐富度,不僅受到林分結(jié)構(gòu)等生物因子的影響,還與土壤肥力等非生物因子密切相關(guān)[4]。前人研究表明:樹高、胸徑、林分密度、郁閉度等林分結(jié)構(gòu)因子和土壤養(yǎng)分、水分會(huì)顯著影響林下草本多樣性[4-6]。探究影響林下草本多樣性的主要因素對(duì)維持森林生態(tài)系統(tǒng)穩(wěn)定具有重要意義。

青海云杉(Piceacrassifolia)是祁連山地區(qū)的主要森林類型,對(duì)于維系該地區(qū)的物種多樣性、改善區(qū)域小氣候、涵養(yǎng)水源和維持生態(tài)平衡等方面均有著至關(guān)重要的作用[7-8]。隨著天然林保護(hù)工程的實(shí)施,經(jīng)過20年的封禁治理,林區(qū)內(nèi)現(xiàn)存大量白樺(Betulaplatyphylla)、紅樺(Betulaalbosinensis)、山楊(Populusdavidiana)等天然林,形成了以青海云杉為頂級(jí)群落的完整演替序列,是研究青海云杉演替過程的良好場(chǎng)所[9]。以往對(duì)青海云杉的研究多集中在土壤質(zhì)量[10]、生態(tài)化學(xué)計(jì)量[11-12]、林分更新[13]、水源涵養(yǎng)[14]等方面,而關(guān)于青海云杉林下草本影響因素和演替過程中的草本多樣性變化的研究卻鮮有報(bào)道。

因此,本研究以闊葉純林(白樺純林、紅樺純林)→闊葉混交林(白樺×山楊、紅樺×山楊)→針闊混交林(白樺×云杉、紅樺×云杉)→常綠針葉林(青海云杉純林)完整的演替序列,以空間代替時(shí)間的方法,通過物種重要值、物種豐富度、Simpson指數(shù)、Shannon-Wiener指數(shù)和Pielou指數(shù)等多樣性指數(shù)對(duì)比分析,比較不同演替階段的草本多樣性特征,探究草本多樣性的影響因子,以期為祁連山地區(qū)物種多樣性保護(hù)和生態(tài)修復(fù)提供理論依據(jù)。

1 研究區(qū)概況

研究區(qū)位于青海省海東市互助縣東北部,地理坐標(biāo)為102°6′-102°43′E, 36°42′-37°60′N。該區(qū)處于青藏和黃土高原的過渡地帶,屬于溫帶大陸性季風(fēng)氣候,年平均氣溫5.8℃,年平均降水量477.4mm。林區(qū)內(nèi)生物資源豐富,主要喬木樹種有青海云杉、白樺、紅樺、山楊等;灌木主要有小檗(Berberisamurensis)、金露梅(Dasiphorafruticosa)、銀露梅(Dasiphoraglabra)等,草本主要有:野草莓(Fragariavesca)、黑麥草(Loliumperenne)、披堿草(Elymusdahuricus)、拳參(Bistortaofficinalis)。

2 研究方法

2.1 樣地設(shè)置與調(diào)查

2022年7月,在青海互助北山林場(chǎng),采用典型抽樣法,按照青海云杉演替序列,選擇闊葉純林(演替早期階段)、闊葉混交林(演替早中期階段)、針闊混交林(演替中后期階段)和青海云杉純林(頂級(jí)群落階段),各6個(gè)共24個(gè)20m×20m的樣地。記錄樣地的經(jīng)緯度、海拔、坡度和坡向等基本信息,樣地基本信息如表1所示。

表1 樣地基本信息表Table 1 Basic information table of sample plots

2.2 生物與非生物指標(biāo)測(cè)定

在各樣方內(nèi)對(duì)每株喬木進(jìn)行每木檢尺,測(cè)高儀測(cè)定每株樹高、胸徑尺測(cè)量胸徑,使用魚眼相機(jī)(WinSCANOPYTM 2010a) 獲取冠層照片, 運(yùn)用XLScanopy軟件處理照片,計(jì)算葉面積指數(shù)和郁閉度。統(tǒng)計(jì)樣地內(nèi)主要樹種的數(shù)量,推算各樣地的林分密度信息。采用五點(diǎn)取樣法,在各樣地取0-20cm的土層混合土樣帶回實(shí)驗(yàn)室,運(yùn)用環(huán)刀法測(cè)定土壤水分、采用重鉻酸鉀稀釋熱測(cè)定土壤有機(jī)質(zhì)[15]、靛酚藍(lán)比色法測(cè)定堿解氮、NaHCO3浸提-火焰光度法測(cè)定速效磷、醋酸銨浸提-火焰光度計(jì)法測(cè)定速效鉀、電位法測(cè)定土壤pH[16]。林分因子代表生物變量、土壤養(yǎng)分代表非生物變量,基礎(chǔ)數(shù)據(jù)如表2所示。

2.3 林下草本多樣性調(diào)查

按五點(diǎn)取樣法在各樣地布設(shè)5個(gè)1m×1m的草本調(diào)查樣方,參考中國(guó)在線植物志,記錄植物拉丁文、科、屬、種和生活習(xí)性等基本信息。統(tǒng)計(jì)計(jì)算各物種的數(shù)量、高度、蓋度,并計(jì)算各物種的重要值。采用草本的重要值(IV)、物種豐富度(S)、Shannon-Wiener多樣性指數(shù)(H)、Simpson多樣性指數(shù)(D)和Pielou均勻度指數(shù)(J)來分析林下植被的物種多樣性[8]。各指標(biāo)計(jì)算公式如下:

J=-∑Pi×lnPi/lnS;

式中:RA表示物種的相對(duì)多度、RC表示相對(duì)蓋度、RH表示物種的相對(duì)高度;Pi表示第i種的個(gè)數(shù)占樣方內(nèi)總物種的個(gè)數(shù);S代表樣方內(nèi)的總物種數(shù)。

2.4 數(shù)據(jù)分析

首先利用Excel 2019和SPSS 22.0對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,采用單因素方法分析(One-way ANOVA)中的LSD(最小二乘法)和鄧肯檢驗(yàn)對(duì)不同演替類型下的林下草本多樣性進(jìn)行差異性檢驗(yàn)。再用Spearman相關(guān)系數(shù)分析林下草本多樣性與影響因子的相關(guān)性,最后采用CANOCO 5.0中的RDA冗余分析來探討影響草本多樣性的主要影響因子。采用Origin 2021繪圖,Excel制表。

3 結(jié)果與分析

3.1 林下草本組成與重要值

林下草本組成如圖1所示,本研究共發(fā)現(xiàn)草本31科64屬81種,其中闊葉純林(Ⅰ)草本植物共21科33屬42種;闊葉混交林(Ⅱ)共22科40屬50種;針闊混交林(Ⅲ)共22科42屬51種;云杉純林(Ⅳ)共15科21屬23種。針闊混交林(Ⅲ)下的草本種類最豐富,而云杉純林(Ⅳ)下的草本種類最少。各林分類型下主要草本的重要值如表3所示:各林分類型之間的優(yōu)勢(shì)草種存在部分重疊,黑麥草、山尖子、野草莓在四種林分類型下均有分布。其中,黑麥草在林分Ⅰ、Ⅲ中的重要值最大,其重要值為:23%、24%;野草莓在林分Ⅱ的重要值最大,其值為26%;冷蕨在林分Ⅳ的重要值最大,其值為33%。

圖1 不同演替階段林下草本組成Fig.1 Understory herbaceous composition at different succession stages

表3 草本主要物種組成及其重要值Table 3 Main species composition of herbs and their important values

3.2 不同演替階段林下草本多樣性特征

不同演替階段草本物種多樣性如圖2所示,階段Ⅰ的物種豐富度均值為11.33,階段Ⅱ豐富度的均值為16.83,階段Ⅲ為19.29,階段Ⅳ為8;物種豐富度和物種多樣性指數(shù)在不同演替階段之間存在顯著差異,各指標(biāo)變化趨勢(shì)一致,隨著青海云杉演替正方向呈現(xiàn)先增加后減小的趨勢(shì),各多樣性指標(biāo)均在Ⅲ最大,在Ⅳ最小。

圖2 不同演替階段草本多樣性Fig.2 Herbaceous diversity at different succession stages不同字母表示差異顯著

3.3 草本多樣性與影響因素的相關(guān)分析

草本多樣性指標(biāo)與生物和非生物因素的相關(guān)分析如表4所示:生物因子中,郁閉度與草本豐富度呈顯著正相關(guān)(P<0.05);林分密度與Simpson指數(shù)和Pielou均勻度指數(shù)顯著負(fù)相關(guān);平均樹高與物種豐富度呈極顯著負(fù)相關(guān)(P<0.01),與Shannon-Wiener指數(shù)顯著負(fù)相關(guān);平均胸徑與物種豐富度、Shannon-Wiener指數(shù)極顯著負(fù)相關(guān),與Simpson指數(shù)顯著負(fù)相關(guān);LAI與4個(gè)多樣性指標(biāo)均顯著負(fù)相關(guān)。非生物因子中,土壤水分與Simpson指數(shù)、Shannon-Wiener指數(shù)、Pielou均勻度指數(shù)呈極顯著負(fù)相關(guān);堿解氮、速效磷與物種豐富度顯著正相關(guān);速效鉀與物種豐富度呈顯著正相關(guān),與Shannon-Wiener指數(shù)呈顯著正相關(guān);土壤有機(jī)質(zhì)與豐富度、Simpson指數(shù)和Shannon-Wiener指數(shù)顯著正相關(guān);土壤pH均大于7且與多樣性指標(biāo)無顯著相關(guān)性。

表4 草本多樣性與影響因素Spearman相關(guān)分析表Table 4 Spearman correlation analysis of herbaceous diversity and influencing factors

3.4 草本多樣性的驅(qū)動(dòng)因子分析

對(duì)各個(gè)環(huán)境因子與草本多樣性指數(shù)進(jìn)行RDA冗余分析(圖3)結(jié)果表明:第一主軸解釋變量為74.95%,第二主軸解釋變量為7.24%,兩軸共解釋草本多樣性與環(huán)境關(guān)系的82.19%,能較好的反映排序的大部分信息。實(shí)心箭頭代表草本多樣性指標(biāo),空心箭頭代表影響因子。影響因子的連線長(zhǎng)度代表該因子與草本多樣性指標(biāo)的影響程度,連線越長(zhǎng)表示該因子影響越大。箭頭之間夾角的余弦值表示相關(guān)系數(shù)[17]。其中,林分密度、LAI、胸徑、樹高與草本多樣性指標(biāo)S,D,H,J的夾角為鈍角且LAI連線最長(zhǎng),這表明:與草本多樣性負(fù)相關(guān)的指標(biāo)中,LAI具有較強(qiáng)的解釋能力。土壤水分和養(yǎng)分等指標(biāo)與草本多樣性指標(biāo)的夾角為銳角,土壤水分、速效鉀和土壤有機(jī)質(zhì)的連線較長(zhǎng),這表明:速效鉀、土壤有機(jī)質(zhì)和土壤水分對(duì)草本多樣性具有較強(qiáng)的解釋能力。采用蒙特卡洛置換檢驗(yàn)對(duì)11個(gè)影響因子進(jìn)行篩選排序(表5)結(jié)果表明:速效鉀變量解釋度為27%、土壤有機(jī)質(zhì)解釋度為28.5%、LAI變量解釋度為13.2%、土壤有機(jī)碳、速效鉀和LAI對(duì)草本多樣性指數(shù)具有顯著影響(P<0.05),是影響草本多樣性的主要因素。

圖3 林下草本多樣性影響因子RDA分析結(jié)果圖Fig.3 RDA analysis results of understory herb diversityS:物種豐富度 Species richness index;D:Simpson多樣性指數(shù) Simpson index; H:Shannon-Wienerd多樣性指數(shù) Shannon-Wiener index; J: Pielou均勻度指數(shù) Pielou evenness index;CD:郁閉度 Canopy density;SD:林分密度 Stand density; MTH:平均樹高 Mean tree height; DBH:平均胸徑 Mean diameter at breast height; LAI:葉面積指數(shù) Leaf area index; SWC:土壤含水量 Soil water content; AN:堿解氮 Alkali-hydrolyzable nitrogen;AP:速效磷 Available phosphorus; AK:速效鉀Available potassium; SOM:土壤有機(jī)質(zhì) Soil organic matter; pH:土壤酸堿度 Soil pH

表5 生物非生物指標(biāo)解釋度與蒙特卡洛檢驗(yàn)表Table 5 Biometric Interpretation and Monte Carlo Test

4 討論

4.1 群落演替與草本多樣性的關(guān)系

本文采用了空間代替時(shí)間的方法研究不同演替階段下的草本多樣性差異,此方法雖簡(jiǎn)單可行,但缺少時(shí)間維度的長(zhǎng)期觀測(cè)。今后可進(jìn)行多年連續(xù)檢測(cè)觀察不同演替階段的草本多樣性的變化。研究表明:祁連山東部地區(qū)森林在不同演替階段的林下草本組成和多樣性存在顯著差異。研究發(fā)現(xiàn):黑麥草、野草莓是不同演替階段的優(yōu)勢(shì)草種。這可能是黑麥草因其根系發(fā)達(dá)、分蘗能力強(qiáng)、繁殖力高,能適應(yīng)各種林下環(huán)境[18];而野草莓通過匍匐莖能占據(jù)更多的生境空間,獲取更多的養(yǎng)分和資源[19]。不同演替階段下,林下草本優(yōu)勢(shì)種以陽(yáng)生草本為主,而陰生草本隨著演替的進(jìn)展其重要值在不斷增加(表3)。在青海云杉純林階段,冷蕨等陰生草本占據(jù)絕對(duì)優(yōu)勢(shì),這主要是所調(diào)查的青海云杉純林密度較大(表2),林內(nèi)形成了偏陰生的環(huán)境且調(diào)查的樣方多在陰坡(表1),光照時(shí)間短,水分較好,適合耐陰喜濕的草本生長(zhǎng)[20]。在針闊混交林階段,草本多樣性指標(biāo)達(dá)到峰值,這可能是一方面:針闊混交增加了生境異質(zhì)性,為不同生活型草本的生長(zhǎng)提供條件[21],另一方面,針闊混交林通過增加枯落物種類,促進(jìn)林內(nèi)養(yǎng)分循環(huán),改良土壤養(yǎng)分狀況[22],進(jìn)而促進(jìn)了林下草本的生長(zhǎng)。本研究發(fā)現(xiàn),在闊葉純林向青海云杉的演替序列中,林下草本多樣性呈先增加而后減少的規(guī)律?？赡苁且环矫媪址盅萏娴角嗪Ｔ粕技兞蛛A段時(shí),林內(nèi)陰暗,氣溫低,濕度大,林下草本難以生長(zhǎng)[23-24];另一方面青海云杉作為淺根性樹種,其根系通過化感作用,產(chǎn)生能夠抑制種子萌發(fā)所需的關(guān)鍵酶類, 抑制了林下草本的生長(zhǎng)[25]。

4.2 生物因素對(duì)草本多樣性的影響

Spearman相關(guān)分析(表4)和RDA冗余分析(圖3)的結(jié)果表明:生物因子中,LAI是顯著影響林下草本多樣性的主要林分因子,LAI與林下草本多樣性呈顯著負(fù)相關(guān)。這可能是LAI代表冠層葉片生長(zhǎng)狀況,反映喬木對(duì)光的獲取能力,喬木冠層通過葉片對(duì)光截獲,抑制林下低矮草本對(duì)光的獲取[26],而光作為植物生長(zhǎng)必要的環(huán)境因子之一[27],光的缺乏會(huì)影響林下草本生長(zhǎng)而導(dǎo)致林下草本多樣性偏低。林分密度是進(jìn)行林分調(diào)控的重要因子,林分密度可通過控制林下光照條件和種間競(jìng)爭(zhēng)影響林下草本的生長(zhǎng)狀況[28-29]。本研究得出:林分密度的增加會(huì)導(dǎo)致物種多樣性的降低,這與部分學(xué)者得出的結(jié)論相同[30-31]。此外,樹高和胸徑能反映喬木生長(zhǎng)狀況[32],是反映林冠結(jié)構(gòu)的重要指標(biāo)。喬木可通過林冠結(jié)構(gòu)改變林內(nèi)光照強(qiáng)度、林分濕度的等林內(nèi)環(huán)境,影響林下草本植物的生長(zhǎng)[33-34]。本研究中:樹高和胸徑與草本多樣性負(fù)相關(guān),這可能是喬木的胸徑越大,樹高越高,其冠幅也越大,喬木對(duì)土壤養(yǎng)分、水分的吸收能力越強(qiáng),導(dǎo)致林下草本土壤水分和養(yǎng)分的缺乏,草本多樣性降低[35-36]。

4.3 非生物因素對(duì)草本多樣性的影響

非生物因子中,速效鉀和土壤有機(jī)質(zhì)是顯著影響林下草本多樣性的主要土壤養(yǎng)分。本研究發(fā)現(xiàn)速效鉀與物種豐富度和Shannon-Wiener指數(shù)顯著正相關(guān),這可能與當(dāng)?shù)氐臍夂蚝统赏聊笌r有關(guān)[37]。土壤中的速效鉀主要取決于土壤中微生物對(duì)礦物分解的活性,雖然祁連山地區(qū)土壤發(fā)育于玄武巖、火山巖等富含鉀元素的礦物,但由于研究區(qū)普遍海拔較高,土壤微生物活性較低,使得土壤中的速效鉀含量不高,成為限制該地區(qū)草本多樣性的主要環(huán)境因子[38]。本研究得出:土壤有機(jī)質(zhì)含量越高,林下草本多樣性越好,可能是土壤有機(jī)質(zhì)作為土壤養(yǎng)分的載體[39-40],通過礦化作用分解成能被林下草本吸收的礦物養(yǎng)分[41]。因此,有機(jī)質(zhì)含量越高,能被植被吸收利用的養(yǎng)分就越多,草本多樣性就越高[42]。土壤水分與草本多樣性正相關(guān),這是因?yàn)橥寥浪肿鳛橄拗浦脖簧L(zhǎng)的重要環(huán)境因子[43],水分的缺乏會(huì)抑制植物對(duì)養(yǎng)分的吸收[44]。氮、磷元素是組成植物蛋白質(zhì)的主要元素,由于在土壤中供應(yīng)受限,通常被認(rèn)為是限制植物生長(zhǎng)的重要因素[45]。楊學(xué)亭通過研究祁連山不同草地類型與養(yǎng)分的關(guān)系得出,氮元素與草本多樣性正相關(guān),磷元素是影響草本多樣的關(guān)鍵因素[46],這與本文研究結(jié)果部分一致。然而,在本研究中卻未體現(xiàn)磷元素的重要性,這可能是研究區(qū)磷元素含量較少,未能體現(xiàn)差異性(表2)。

5 結(jié)論

在闊葉純林→闊葉混交林→針闊混交林→云杉純林的演替序列過程中,草本多樣性總體呈先增大而后減少的規(guī)律,并在針闊混交林達(dá)到最大,云杉純林最小。LAI是影響草本多樣性的主導(dǎo)生物因子,速效鉀,土壤有機(jī)質(zhì)是影響草本多樣性的主要非生物因子。今后,在草本多樣性低下的青海云杉林進(jìn)行林分改造時(shí),可以適當(dāng)?shù)恼{(diào)整林分密度和補(bǔ)植樺樹以營(yíng)造針闊混交林。在對(duì)祁連山東部地區(qū)的森林群落多樣性研究過程中,可考慮加入角尺度、競(jìng)爭(zhēng)指數(shù)、光照條件等其他環(huán)境因子的作用,以期為祁連山地區(qū)森林生態(tài)恢復(fù)和林分改造提供理論支持。