李云琴，孟夢，胡光輝，周云，付玉嬪，邱瓊

(1.云南省林業(yè)和草原科學(xué)院，云南 昆明 650201；2.云南省森林植物培育與開發(fā)利用重點(diǎn)實(shí)驗(yàn)/國家林業(yè)和草原局云南珍稀瀕特森林植物保護(hù)和繁育重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650201；3.云南省高黎貢山生物多樣性重點(diǎn)實(shí)驗(yàn)室，云南 昆明 650201；4.高黎貢山森林生態(tài)系統(tǒng)云南省野外科學(xué)觀測研究站，云南 保山 678000)

森林是生態(tài)安全的根基，發(fā)揮著保障生態(tài)安全的主體責(zé)任，是我國綠色發(fā)展的重要載體，保持良好的生物多樣性是生態(tài)安全的基礎(chǔ)，是維系生態(tài)系統(tǒng)穩(wěn)定性和可持續(xù)性的關(guān)鍵[1-2]，而物種多樣性是生物多樣性的物種表現(xiàn)形式之一，是對(duì)生物多樣性較為直觀的認(rèn)識(shí)，是生物多樣性研究的核心內(nèi)容[3]。它是生態(tài)系統(tǒng)穩(wěn)定性的度量指標(biāo)，可以較為有效地反映群落發(fā)展階段、結(jié)構(gòu)類型、穩(wěn)定程度和功能特征的較有效的指標(biāo)，是備受關(guān)注的熱點(diǎn)與研究內(nèi)容[4]。

云南省森林面積較大，類型多樣，發(fā)展林下經(jīng)濟(jì)前景廣闊，已成為當(dāng)前山區(qū)致富增收的重要手段，促進(jìn)了區(qū)域經(jīng)濟(jì)的發(fā)展[5]。但在種植過程中，一味追求經(jīng)濟(jì)效益，出現(xiàn)林下土地過度開墾、生物多樣性減少、生態(tài)系統(tǒng)退化、水土流失等生態(tài)安全問題[6]，這與新時(shí)代綠色發(fā)展理念相悖，制約林下經(jīng)濟(jì)健康、持續(xù)發(fā)展。云南松(Pinusyunnanensis)林是云南省乃至我國西南地區(qū)最重要的森林類型之一，在云南省分布面積為167 572.55 km2，占云南省面積的 43.66%[7]，為開展林下種植提供了廣闊空間。前人對(duì)云南松開展了大量研究，如繁育[8]、遺傳變異[9]、種質(zhì)資源[10]等，但是關(guān)于云南松林下種植經(jīng)營模式對(duì)物種多樣性影響的研究尚屬空白。

木本森林蔬菜作為云南省特色優(yōu)勢資源，由于其風(fēng)味獨(dú)特、無污染、具有食療作用，深受廣大群眾喜愛與食用，市場開發(fā)前景巨大，經(jīng)濟(jì)效益顯著[11]。因此，充分利用云南省廣大的林下空間發(fā)展木本森林蔬菜種植，拓展林地空間利用效率，開發(fā)優(yōu)良種植資源，研發(fā)高效培育技術(shù)，推進(jìn)產(chǎn)業(yè)發(fā)展，將為云南省當(dāng)前綠色發(fā)展戰(zhàn)略提供新支撐。

本研究基于生態(tài)安全的前提，選擇在樹齡40 a左右的云南松林下種植市場認(rèn)知度高、口感好、營養(yǎng)豐富的木本森林蔬菜——大參(Macropanaxdispermus)，開展不同種植擾動(dòng)強(qiáng)度對(duì)云南松林內(nèi)植物物種組成和生物多樣性影響的研究，了解不同經(jīng)營水平下林內(nèi)物種多樣性的變化，評(píng)估林下種植對(duì)森林生態(tài)環(huán)境的干擾和影響，最終確定云南松林下生境種植森林蔬菜的開發(fā)承載量，為構(gòu)建一套云南松林下種植森林蔬菜開發(fā)控制體系，以期為解決當(dāng)前林下種植單一追求經(jīng)濟(jì)效益而過度開墾林下土地引起的生物多樣性減少等生態(tài)問題做參考。

1 材料與方法

1.1 研究區(qū)概況

試驗(yàn)地位于昆明樹木園的云南松人工林，地理位置為25°8′N、102°45′E，海拔2 000 m左右。占地面積52.4 hm2。冬無嚴(yán)寒，夏無酷暑，夏季溫暖的半濕潤北亞熱帶高原季風(fēng)氣候類型。年均氣溫14.8 ℃，絕對(duì)最低氣溫-7.6 ℃，最熱月(7月)平均氣溫21 ℃，≥10 ℃的積溫4 807.6 ℃；年降水量700～1 100 mm，6—10月為雨季，降水量占全年降水量的85%，干濕季節(jié)分明，年均相對(duì)濕度72%。土壤為第四紀(jì)紅土層發(fā)育的山地紅壤和玄武巖發(fā)育的山地紅壤，pH5～6[12]。

林分內(nèi)上層僅為云南松，林相整齊，為同齡林，林下灌木零星分布有蒙桑(Morusmongolica)，刺槐(Robiniapseudoacacia)，灌木狀喀西茄 (Solanumaculeatissimum)等，草本層以菊科(Asteraceae)、禾本科(Poaceae)、蝶形花科(Papilionaceae)等植物為主。

1.2 研究方法

于2020年7月，在云南松林下選用無干擾自然經(jīng)營、30%擾動(dòng)(近自然經(jīng)營)和50%擾動(dòng)(中強(qiáng)度)和100%擾動(dòng)(高強(qiáng)度)等4種種植模式種植購自紅河哈尼族彝族自治州金平縣的大參實(shí)生苗，苗平均高40 cm，每種模式面積100 m2。

(1)100%擾動(dòng)，即林下灌木、草本層全部伐除，再以橫坡，寬窄行方式雙行種植大參。株距0.4 m，窄行距0.4 m，寬行距1.6 m。后期管理過程中，每年5月、10月進(jìn)行中耕除草。(2)50%擾動(dòng)，即林下以橫坡向2 m帶寬間隔伐除灌木、草本層，在伐除帶中，株行距均為0.4 m，雙行種植大參。(3)30%擾動(dòng)，即林下以橫坡向2.2 m帶寬伐除灌木、草本層，在伐除帶中，株行距均為0.4 m，雙行種植大參；伐除帶下方設(shè)4.4 m原生植被保留帶，下方再以橫坡向2.2 m帶寬伐除灌木、草本層，在伐除帶中，株行距均為0.4 m，雙行種植大參，以此類推。(4)無干擾，即林下灌木、草本層不擾動(dòng)，保留原始狀態(tài)。后期管理過程中，每年5月、10月，僅種植帶中進(jìn)行中耕除草，原生植物保留帶中則不進(jìn)行擾動(dòng)。

2021年10月，在4 種種植模式的大參人工林中分別設(shè)置1塊規(guī)格 10 m×10 m樣地，調(diào)查喬木的株數(shù)、胸徑、高度等，同時(shí)記載樣地的海拔、坡向、坡度等環(huán)境因子，分別在每1塊樣地下設(shè)置4塊規(guī)格5 m×5 m的小樣方，進(jìn)行灌木、草本調(diào)查，分別記錄物種名、高度、蓋度等指標(biāo)，為確保調(diào)查準(zhǔn)確性，樣地的立地條件盡量保持一致[13]。

表1 云南松林樣方的自然概況Tab.1 Sample situation

1.3 多樣性指數(shù)及計(jì)算方法

研究選用α 物種多樣性測度方法中較為常用Simpson指數(shù)(D)、Shannon-Wiener指數(shù)(H′)、Pielou均勻度指數(shù)(Jsw、Jsv)來描述物種多樣性[14-15]。公式如下。

喬木、灌木重要值Pi=(相對(duì)蓋度+相對(duì)頻度+相對(duì)多度)/3

草本重要值Pi=(相對(duì)蓋度+相對(duì)頻度)/2

H′=-∑PilnPi

D=1-∑Pi2

式中：S為樣本的植物種總數(shù)，即豐富度指數(shù)；Pi為種i的重要值。

1.4 數(shù)據(jù)處理和分析

用Excel 2010軟件進(jìn)行試驗(yàn)數(shù)據(jù)整理、分析，Origin 2021進(jìn)行作圖，并用SPSS 27.0軟件進(jìn)行相關(guān)性分析。

2 結(jié)果與分析

2.1 物種構(gòu)成

2.1.1 不同種植模式下植物物種組成

4種模式中共調(diào)查記錄植物101種，包括46科85屬(包括云南松和大參)?？品N數(shù)量大小依次為菊科14種，薔薇科(Rosaceae)7種，禾本科7種，蝶形花科5種，唇形花科(Lamiaceae)4種，茜草科(Rubiaceae)5種，蓼科(Polygonaceae)，桑科(Moraceae)分別各3種，其余科均為1～2種。樣方中的物種以草本為主。

無干擾模式記錄到植物22科37種，其中喬木8種，灌木4種，草本25種。草本中菊科5種，禾本科5種，薔薇科4種，茜草科3種。高強(qiáng)度100%擾動(dòng)模式調(diào)查記錄到33科56種，物種數(shù)較無干擾模式高56.77%。其中喬木5種，灌木10種，草本41種，草本菊科11種，禾本科6種，薔薇科3種，蝶形花科3種。50%擾動(dòng)模式記錄到26科53種，物種數(shù)較無干擾模式高43.24%。其中喬木12種，灌木9種，草本32種，草本中草本菊科11種，薔薇科6種，禾本科5種，蝶形花科3種。30%擾動(dòng)的模式記錄到37科66種，物種數(shù)較無干擾模式高81.08%。其中喬木9種，灌木14種，草本43種，草本中菊科9種，禾本科7種，薔薇科6種，蝶形花科3種。不同程度的種植模式下物種數(shù)均比無干擾模式多，其中以30%擾動(dòng)模式物種數(shù)最多。

不同擾動(dòng)模式物種科種數(shù)量均比無干擾模式有所增加(圖1)，植物物種的科種數(shù)目均以30%擾動(dòng)最多，表現(xiàn)規(guī)律為30%擾動(dòng) >100%擾動(dòng)>50%擾動(dòng)>無干擾。

圖1 不同擾動(dòng)模式下植物構(gòu)成Fig.1 Plant composition of different disturbance modes

圖2 物種多樣性變化情況

2.1.2 擾動(dòng)模式下物種重要值

不同干擾模式下植物種的重要值見表2。4種模式下喬木重要值最大的均為云南松，灌木為桑科的蒙桑，蝶形花科的刺槐，茄科(Solanaceae)的灌木狀喀西茄等，草本層主要以菊科紫莖澤蘭(Ageratinaadenophora)、小蓬草(Erigeroncanadensis)、野茼蒿(Crassocephalumcrepidioides)，禾本科的矛葉藎草(Arthraxonprionodes)、細(xì)柄草(Capillipediumparviflorum)，蝶形花科的兩型豆(Amphicarpaeaedgeworthii)等為主。

表2 不同擾動(dòng)模式下主要物種的重要值Tab.2 Important values of main species of different disturbance modes

2.2 不同擾動(dòng)模式下物種多樣性分析

不同擾動(dòng)模式對(duì)草本層物種多樣性影響較大，對(duì)灌木影響較小，而對(duì)喬木幾乎沒有影響[16]。不同擾動(dòng)模式下喬木灌木草本的Shannon-Wiener 指數(shù)均是先升高后下降，喬木層Shannon-Wiener 指數(shù)在50%擾動(dòng)模式下最高，其次是30%擾動(dòng)模式。灌木和草本層的Shannon-Wiener 指數(shù)在擾動(dòng)為30%模式下最高，其次是50%擾動(dòng)模式。Simpson多樣指數(shù)也是隨著擾動(dòng)強(qiáng)度升高呈現(xiàn)先升高后下降的趨勢，喬木、灌木、草本都是以30%模式下值最高，其次是50%模式。喬木層的Pielou均勻度指數(shù)Jsw以無干擾模式最高，其次是30%擾動(dòng)模式，Jsv以30%擾動(dòng)模式最高，其次是無干擾模式；灌木和草本的Jsw和Jsv也和Shannon-Wiener 指數(shù)、Simpson多樣指數(shù)一樣呈現(xiàn)先升高后下降的趨勢，并且以30%擾動(dòng)模式最高，是因?yàn)镻ielou均勻度指數(shù)是根據(jù)物種數(shù)、 Shannon-Wiener 指數(shù)以及Simpson多樣指數(shù)計(jì)算得到的。

根據(jù)多樣性指數(shù)可以看出，林分?jǐn)_動(dòng)強(qiáng)度要適中，過度的人為干擾或者無干擾均對(duì)物種在林分內(nèi)的競爭產(chǎn)生影響，進(jìn)而影響了林分植物的物種多樣性。適度的擾動(dòng)優(yōu)化了林下環(huán)境，開闊了林下空間，改善了光照條件，增加了光合利用效率，使得物種可以充分利用立地條件，有利于一些原本沒有生存機(jī)會(huì)的植物萌發(fā)生長，增加物種多樣性。而無干擾模式下林下草本植物處于林相底層，不利于與上層的云南松及其他喬灌木競爭，由于得不到生長所需的光照、水分、溫度、養(yǎng)分等，致使一些草本植物滅絕，因此，無干擾模式的林下草本植物的生物多樣性指數(shù)尤其是Shannon-Wiener 指數(shù)、Simpson多樣指數(shù)顯著小于30%和50%擾動(dòng)模式，無干擾模式不利于林下植被的生長更新，而30%(近自然)擾動(dòng)模式更有利于植被的生長[12]。

2.3 相關(guān)性分析

對(duì)擾動(dòng)強(qiáng)度和喬木、灌木、草本的物種多樣性指數(shù)進(jìn)行相關(guān)性分析(表3)，由表3可知，喬木層Simpson多樣指數(shù)和擾動(dòng)強(qiáng)度呈顯著負(fù)相關(guān)(P<0.05)，均勻度指數(shù)Jsw和Jsv均和擾動(dòng)強(qiáng)度呈極顯著負(fù)相關(guān)(P<0.01)；灌木層和草本層僅均勻度指數(shù)Jsw與擾動(dòng)強(qiáng)度呈極顯著負(fù)相關(guān)(P<0.01)。

表3 擾動(dòng)強(qiáng)度和物種多樣性指數(shù)的相關(guān)性Tab.3 Correlation between disturbance intensity and species diversity indices

對(duì)不同多樣性指數(shù)和喬木、灌木、草本層分層多樣性指數(shù)之間進(jìn)行相關(guān)性分析(表4)，結(jié)果發(fā)現(xiàn)，喬木層Shannon-Wiener 指數(shù)和Simpson多樣指數(shù)呈極顯著正相關(guān)(P<0.01)，Pielou的均勻度指數(shù)Jsv與Simpson多樣指數(shù)和Pielou的均勻度指數(shù)Jsw呈極顯著正相關(guān)(P<0.01)。灌木層Shannon-Wiener 指數(shù)與Simpson多樣指數(shù)、Pielou的均勻度指數(shù)Jsv呈極顯著正相關(guān)(P<0.01)，Pielou的均勻度指數(shù)Jsv與Simpson多樣指數(shù)和Pielou的均勻度指數(shù)Jsw呈極顯著正相關(guān)(P<0.01)。草本層Shannon-Wiener 指數(shù)與Simpson多樣指數(shù)、Pielou的均勻度指數(shù)Jsv、Pielou的均勻度指數(shù)Jsw兩兩之間均呈極顯著正相關(guān)(P<0.01)。

表4 多性指數(shù)之間的相關(guān)性Tab.4 Correlation between diversity indices

2.4 生活型譜

生活型是植物對(duì)于外界環(huán)境條件長期適應(yīng)而在外貌上反映出來的植物類型[17]。生活型的形成是不同的植物類型對(duì)同一環(huán)境條件下在形態(tài)、生理及適應(yīng)方式等方面產(chǎn)生趨同適應(yīng)的結(jié)果。通過生活型譜可以了解植物物種隨著生境變化而改變的一些外貌特征[17]。

不同模式下植物生活型譜見圖3，無干擾模式下植物生活型以高位芽植物最多，有14 種，占總數(shù)的37.84%，其次是地面芽植物有12種，占總數(shù)的32.43%。30%擾動(dòng)以高位芽植物最多，有25種，占總數(shù)的37.31%，其次是地面芽植物，有23種，占總數(shù)的34.33%。50%擾動(dòng)模式下，高位芽植物最多，有20種，占總數(shù)的37.74%，其次是地面芽植物16種，占總數(shù)的30.19%。100%擾動(dòng)模式下以地面芽植物最多，有21種，占總數(shù)的37.5%，如求米草(Oplismenusundulatifolius)、皺葉狗尾草(Setariaplicata)等。其次是一年生植物，有13種，占總數(shù)的23.21%。

圖3 植物生活型譜Fig.3 Plant lifeform spectrum

高位芽植物主要是喬木和灌木，如云南油杉(Keteleeriaevelyniana)、云南松等，在林地中本來就存在的，且受擾動(dòng)影響較小。除去高位芽植物，4個(gè)擾動(dòng)模式植物生活型譜都以地面芽植物為主，一年生植物次之。且在不同擾動(dòng)強(qiáng)度下，二者所占比重均比無干擾模式上升。

3 討論與結(jié)論

森林經(jīng)營是維護(hù)森林系統(tǒng)生態(tài)平衡的關(guān)鍵技術(shù)措施，也是維持林業(yè)可持續(xù)健康發(fā)展的重要組成部分。目前已有學(xué)者開展了關(guān)于不同種植經(jīng)營模式對(duì)林下植物多樣性影響的研究。龔固堂等[18]研究表明經(jīng)營措施促使林內(nèi)灌木層和草本層的生長以及增加物種數(shù)量，還有研究[19]表明，林下灌木和草本層的物種數(shù)隨著經(jīng)營強(qiáng)度增加而顯著增加。本研究也得出一致結(jié)果：不同程度的擾動(dòng)模式下物種數(shù)均比無干擾模式多，其中以30%擾動(dòng)模式物種數(shù)最多，即干擾后，灌木層和草本層以及總體物種豐富度均明顯增加，可能是擾動(dòng)導(dǎo)致地面小生境發(fā)生變化，使得土壤種子庫中一些在競爭中處于弱勢或沒有機(jī)會(huì)萌發(fā)生長的物種得到生長空間，因此增加了物種數(shù)。

張象君等[20]研究發(fā)現(xiàn)近自然改造促進(jìn)林下草本層植物的生長發(fā)育，提高了物種多樣性；林同龍[21]以杉木(Cunninghamialanceolata)人工林為研究對(duì)象，發(fā)現(xiàn)近自然經(jīng)營模式的林內(nèi)物種數(shù)、喬灌草多樣性指數(shù)等高于常規(guī)經(jīng)營模式。李瑞霞等[19]對(duì)馬尾松(Pinusmassoniana)人工林開展不同強(qiáng)度經(jīng)營模式研究發(fā)現(xiàn)間伐顯著影響了馬尾松人工林林下植被多樣性，其中中度間伐經(jīng)營林下物種多樣性最高。馮穎等[22]對(duì)華北落葉松(Larixgmeliniivar.principis-rupprechtii)人工林采用無干擾、常規(guī)和近自然經(jīng)營3種經(jīng)營方式，發(fā)現(xiàn)在近熟林階段，近自然經(jīng)營模式的物種數(shù)，多樣性指數(shù)最高，而無干擾經(jīng)營模式阻礙了林下植被的生長。趙蘇亞等[23]研究表明撫育間伐提高了杉木林下植被多樣性。本研究綜合各處理措施分析，30%擾動(dòng)(近自然經(jīng)營)模式下，多樣性指數(shù)較高，表明該強(qiáng)度的擾動(dòng)模式為試驗(yàn)區(qū)云南松林下種植大參的最好模式，與上述研究結(jié)果一致。適度地種植經(jīng)營模式，拔除林下一些灌木和草本，增加了到達(dá)地面的光照，進(jìn)而促進(jìn)地表枯落物的分解，改變植物生長的環(huán)境條件，環(huán)境異質(zhì)性增加，環(huán)境復(fù)雜性就越高，生物群落復(fù)雜性和生物多樣性就越高[24-26]。100%擾動(dòng)模式下，物種多樣性偏低。主要是因?yàn)楦邚?qiáng)度的全面清除林下植被，極大程度上破壞了生態(tài)環(huán)境，物種對(duì)環(huán)境的適應(yīng)具有滯后效應(yīng)，需要一定時(shí)間恢復(fù)，從而降低了植物多樣性，與涂丹丹等[4]研究結(jié)果一致，即高強(qiáng)度采伐導(dǎo)致林下植被物種多樣性降低。而無干擾模式隨著林分密度越來越大，透光度大大降低，下層植被由于得不到生長需要的養(yǎng)分、水分、光照等條件，使很多林下植物尤其是喜光植物無法生長，從而導(dǎo)致林下生物多樣性降低。因此，無干擾模式植物的生物多樣性也偏低。本研究結(jié)果與 “中度干擾假說”相吻合，即弱中度干擾后植物物種數(shù)量增加，高強(qiáng)度干擾后植物物種減少[27]。

生活型是植物對(duì)所在環(huán)境綜合條件長期適應(yīng)而表現(xiàn)在外貌上的類型。它是不同的植物體與相同環(huán)境之間長期趨同適應(yīng)的結(jié)果[28-29]。本研究除了100%擾動(dòng)模式外，生活型均以高位芽植物最多，除去高位芽植物，各擾動(dòng)模式均以地面芽為主，一年生植物次之，且比重均比無干擾模式上升。擾動(dòng)改變了植物生境的異質(zhì)性，從而改變了植物的生活型譜[16]。

本研究僅是在云南松林下種植大參1年后，不同擾動(dòng)模式對(duì)林下生物多樣性變化情況的短期研究結(jié)果。為了更加科學(xué)合理的評(píng)價(jià)種植擾動(dòng)強(qiáng)度對(duì)物種多樣性的影響，將持續(xù)對(duì)該林分開展監(jiān)測，并與森林撫育相結(jié)合，研究其長期穩(wěn)定的變化規(guī)律，為該林型不同種植擾動(dòng)強(qiáng)度對(duì)森林的影響效果給出合理的評(píng)價(jià)，為該地區(qū)不同森林蔬菜林下種植給出合理建議。

致謝：本研究在植被調(diào)查過程中,得到云南省林業(yè)和草原科學(xué)院司馬永康研究員熱情幫助,在此表示衷心感謝!