馬占霞,張 玲,甘建民,武珊珊,蘇建美,尚宇梅,李鑫鑫

(1.滇西應(yīng)用技術(shù)大學(xué)普洱茶學(xué)院,云南 普洱 665000;2.中國科學(xué)院西雙版納熱帶植物園熱帶森林生態(tài)學(xué)重點(diǎn)實(shí)驗(yàn)室,云南 勐侖 666303)

位于印度-緬甸生物多樣性熱點(diǎn)區(qū)域的西雙版納,是我國熱帶森林植被面積最大、類型最多樣、保存最完整、生物多樣性最豐富的地區(qū),也是具有全球意義的生物多樣性關(guān)鍵地區(qū)之一[1]。由于處于熱帶生物區(qū)系向亞熱帶生物區(qū)系過渡的生物地理群落交錯(cuò)帶上,西雙版納地區(qū)生物多樣性較為豐富,但生態(tài)系統(tǒng)也較為脆弱[2-3]。近幾十年來,西雙版納地區(qū)大面積的熱帶雨林逐漸被單一種植的橡膠林替代,橡膠林面積由1976年占全區(qū)面積的1.3%擴(kuò)展到2007年的11.3%[4],2017年西雙版納傣族自治州橡膠種植面積已達(dá)50.04萬hm2[5-6]。橡膠種植帶來了一系列生態(tài)問題[7-8],其中最直接的生態(tài)效應(yīng)是生物多樣性降低[9-10]。

土地利用變化造成生物多樣性喪失的原因主要是:使原有生境消失,從而導(dǎo)致依賴這些生境的物種隨之滅絕;改變了原有生境的空間格局,影響物種的擴(kuò)散、遷移和建群,從而縮小了物種的分布范圍,降低個(gè)體和種群的存活機(jī)會(huì)[11-12]。WATT等[13]研究發(fā)現(xiàn),在間伐熱帶天然林林冠上,螞蟻物種保留在87%以上,是皆伐林地的20倍。熱帶森林中甲蟲的物種豐富度比橡膠種植園高,且隨著森林破碎化程度的增高而減少[14]。不同節(jié)肢動(dòng)物類群對(duì)土地利用方式的響應(yīng)不同,BENG等[15]對(duì)西雙版納茶園和橡膠林中節(jié)肢動(dòng)物的調(diào)查結(jié)果顯示,鞘翅目和雙翅目與土地利用方式相關(guān),而鱗翅目和蜚蠊目對(duì)土地利用沒有響應(yīng)。隨著西雙版納地區(qū)土地利用過程中干擾強(qiáng)度的加大及橡膠林的擴(kuò)張,甲蟲、螞蟻等節(jié)肢動(dòng)物物種多樣性和群落組成均受到影響[16-19]。目前,西雙版納土地利用方式對(duì)生物多樣性影響的研究主要側(cè)重于節(jié)肢動(dòng)物的物種多樣性和不同土地利用方式下節(jié)肢動(dòng)物多樣性的比較分析,然而在單一橡膠林種植背景下,森林鑲嵌格局在不同景觀尺度上對(duì)節(jié)肢動(dòng)物物種分布格局的影響研究還鮮見報(bào)道。綜上所述,探究具有指示環(huán)境變化功能的節(jié)肢動(dòng)物群落對(duì)橡膠種植的響應(yīng),可為西雙版納生物多樣性保護(hù)提供數(shù)據(jù)支撐。

筆者探究熱帶森林轉(zhuǎn)化為橡膠林后,橡膠林和森林鑲嵌景觀格局對(duì)節(jié)肢動(dòng)物群落組成和物種多樣性的影響,分析影響節(jié)肢動(dòng)物分布格局的環(huán)境因素,從而探討森林破碎化產(chǎn)生的生態(tài)學(xué)效應(yīng),為保護(hù)區(qū)規(guī)劃提供數(shù)據(jù)支持。采用野外調(diào)查和室內(nèi)分析相結(jié)合方法,探討西雙版納納板河流域橡膠林和森林鑲嵌格局對(duì)節(jié)肢動(dòng)物群落結(jié)構(gòu)和生物多樣性的影響,擬回答以下主要科學(xué)問題:(1)森林和橡膠林節(jié)肢動(dòng)物群落組成和多樣性是否有差異?(2)橡膠林景觀背景下森林鑲嵌格局導(dǎo)致森林和橡膠林節(jié)肢動(dòng)物群落結(jié)構(gòu)發(fā)生變化的主導(dǎo)因素是什么?

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)域位于云南省西雙版納傣族自治州景洪市與勐?？h接壤處的納板河流域國家級(jí)自然保護(hù)區(qū)及其周邊森林和橡膠林,總面積為270.1 km2,海拔為800～2 304 m,全區(qū)地勢(shì)西北高東南低,中部為納板河流域谷地;氣候類型屬北熱帶濕潤氣候區(qū),年降水量為1 100～1 600 mm;保護(hù)區(qū)按小流域生物圈保護(hù)區(qū)理念建設(shè),劃分核心區(qū)、緩沖區(qū)和試驗(yàn)區(qū)3大功能區(qū)。在緩沖區(qū)和試驗(yàn)區(qū)遍布數(shù)個(gè)村寨及大面積橡膠林,其中,零星散布著被橡膠林隔斷的殘余森林片段。

1.2 取樣方法

通過實(shí)地情況調(diào)查并結(jié)合使用IKONOS衛(wèi)星地圖[20],選取9個(gè)片段森林及其相毗鄰的9個(gè)橡膠林為取樣樣地(圖1、表1)。節(jié)肢動(dòng)物的獲取采用馬氏網(wǎng)(長×寬×高為3.5 m×2.0 m×1.5 m)的方法[21]。在所選取的每個(gè)樣地內(nèi)分別設(shè)置4個(gè)馬氏網(wǎng),且兩網(wǎng)間距大于100 m,距離林地邊緣大于100 m,以避免邊緣效應(yīng)對(duì)實(shí)驗(yàn)結(jié)果的影響。定期收集樣品,包括旱季(2014年3—5月)和雨季末(2014年9—11月),每15 d清空1次收集瓶,旱季和雨季各收集6次節(jié)肢動(dòng)物樣品。所收集樣品均保存在φ=95%的酒精中帶回實(shí)驗(yàn)室進(jìn)行分子生物學(xué)實(shí)驗(yàn),利用metabarcoding技術(shù)鑒定至可操縱的分類單元(operational taxonomic units,OTUs)。

表1 納板河研究樣地基本信息

1.3 樣地環(huán)境因子的測(cè)定

調(diào)查過程中,用手持GPS儀(GARMIN 62S)記錄相關(guān)樣點(diǎn)經(jīng)緯度、海拔、坡度和坡向等樣點(diǎn)信息。每個(gè)取樣點(diǎn)的土地利用和森林覆蓋情況,使用分辨率圖像(來自DigitalGlobe衛(wèi)星的ArcGIS 10.5 Baseman圖像)進(jìn)行繪制。在景觀尺度上,以樣地為中心,分別以500、600、700、800、900和1 000 m為半徑畫圈,得到6個(gè)尺度的覆蓋區(qū),并從該區(qū)內(nèi)提取所需要的森林覆蓋率指數(shù)。對(duì)不同景觀尺度的森林覆蓋率和節(jié)肢動(dòng)物物種多樣性做logistic回歸分析,結(jié)果顯示,在以500 m為半徑的區(qū)域內(nèi),節(jié)肢動(dòng)物物種多樣性和森林覆蓋率的相關(guān)性尤為顯著(R2=0.75,P=0.001),同時(shí),鞘翅目(R2=0.48,P=0.017)、雙翅目(R2=0.62,P=0.004)和膜翅目(R2=0.50,P=0.025)與森林覆蓋率的相關(guān)性也在該區(qū)域內(nèi)尤為顯著。因此,最終確定景觀尺度為以500 m為半徑的圓。在樣方尺度上,變量在半徑為50 m的覆蓋區(qū)內(nèi)進(jìn)行評(píng)估。選擇不同的空間尺度,以覆蓋具有不同覓食距離的一些節(jié)肢動(dòng)物的“效應(yīng)尺度”。景觀尺度環(huán)境變量包括天然森林覆蓋率、取樣點(diǎn)到森林距離。橡膠林內(nèi)的環(huán)境變量包括膠樹年齡、橡膠樹間距、取樣點(diǎn)到片段森林邊緣距離、海拔、坡度和坡向。森林內(nèi)的環(huán)境變量包括森林面積、樹種多樣性、取樣點(diǎn)到片段森林邊緣距離(表1～2)。所有景觀變量采用ArcGIS 10.5計(jì)算。

表2 景觀尺度和樣方尺度的解釋變量及其描述

1.4 分子實(shí)驗(yàn)流程

節(jié)肢動(dòng)物樣品參照YU等[3]方法進(jìn)行挑撿,即從每一只節(jié)肢動(dòng)物上獲取組織材料,不同節(jié)肢動(dòng)物取樣位置不同,如較小的節(jié)肢動(dòng)物取整只,較大的取前足;如鞘翅目、鱗翅目和蜘蛛目等一般取前足,大小如蚊子大小。用德國QIAGEN公司提供的DNA提取試劑盒(QIAGEN DNeasy Blood & Tissue)對(duì)混合樣品提取總DNA。利用美國invitrogenQubit 2.0熒光定量儀和QubitdsDNA HS分析試劑盒對(duì)基因組DNA精確定量,以確定PCR反應(yīng)加入的DNA量。PCR選取的引物和反應(yīng)體系參照YU等[3]的實(shí)驗(yàn)。PCR產(chǎn)物定量混合后制備文庫,用美國Illumina公司提供的建庫試劑盒,按照建庫說明書建立文庫,并用北京索萊寶科技有限公司提供的2×SYBRGreen qpcr 試劑盒對(duì)文庫進(jìn)行絕對(duì)定量。等濃度混合文庫后,按照說明書,使用MiSeq ReagentKit v3(Illumina,nc.,2015)進(jìn)行測(cè)序。用Soapbarcode[22]流程對(duì)所得序列進(jìn)行分析,包括降噪、嵌合體識(shí)別、OTU挑選和分類歸屬鑒定等。

1.5 統(tǒng)計(jì)分析

利用配對(duì)F檢驗(yàn)分析片段森林和橡膠林中節(jié)肢動(dòng)物物種多樣性差異(R軟件包vegan)。對(duì)不同土地利用方式節(jié)肢動(dòng)物物種組成進(jìn)行維恩圖(Venn diagram)分析(R軟件包VennDiagram),并對(duì)節(jié)肢動(dòng)物物種組成進(jìn)行稀疏曲線分析和制圖(R軟件包vegan和elasticnet包)。為了探究森林和橡膠林中節(jié)肢動(dòng)物群落組成,對(duì)節(jié)肢動(dòng)物群落組成結(jié)構(gòu)差異性進(jìn)行分析:用bray相異性指數(shù)度量森林之間和橡膠林之間節(jié)肢動(dòng)物群落的beta多樣性(R軟件包vegan),并做非度量多維尺度法(nonmetric multidimensional scaling,NMDS)和差異性檢驗(yàn)(R軟件包vegan和cluster),用脅強(qiáng)系數(shù)(stress,S)衡量NMDS分析結(jié)果的優(yōu)劣。通常認(rèn)為當(dāng)S<0.2時(shí),可用NMDS的二維點(diǎn)圖表示,其圖形有一定的解釋意義;當(dāng)S<0.1時(shí),可認(rèn)為是一個(gè)好的排序;當(dāng)S<0.05時(shí),認(rèn)為具有很好的代表性。為了確定不同景觀尺度下每個(gè)變量對(duì)響應(yīng)變量的影響,在R軟件包中挖掘1組預(yù)測(cè)變量,對(duì)候選模型(廣義線性模型,GLM)進(jìn)行方差分析,選出最優(yōu)模型(R軟件包lme4)。然后對(duì)最佳候選模型進(jìn)行多模型推理(R軟件包MuMIn),采用赤池信息準(zhǔn)則(AIC，CAI)比較不同模型的擬合優(yōu)度,其中,AIC的截止閾值(ΔwAIC,ΔCwAI)<2[23],并對(duì)最佳模型集進(jìn)行平均化。

2 結(jié)果與分析

2.1 橡膠林和森林林下節(jié)肢動(dòng)物群落組成和alpha多樣性

在片段森林中采集節(jié)肢動(dòng)物共14 451個(gè)個(gè)體,隸屬于27個(gè)類群。其中,優(yōu)勢(shì)類群包括彈尾目(Collembola)、鱗翅目(Lepidoptera)、蜱螨目(Acarina)和雙翅目(Diptera),4者個(gè)體數(shù)之和占總個(gè)體數(shù)的70.22%;常見類群包括膜翅目(Hymenoptera)、鞘翅目(Coleoptera)、嚙目(Psocoptera)、蜘蛛目(Araneae)、同翅目(Homoptera)、等翅目(Isoptera)、纓翅目(Physopoda)和直翅目(Orthoptera),占總個(gè)體數(shù)的26.42%;其余稀有類群有15目,占總個(gè)體數(shù)的3.36%。橡膠林中采集到節(jié)肢動(dòng)物10 754個(gè)個(gè)體,隸屬于24個(gè)類群。其中,優(yōu)勢(shì)類群包括蜱螨目、鱗翅目、雙翅目和膜翅目,占總個(gè)體數(shù)的71%;常見類群包括彈尾目、鞘翅目、嚙目、直翅目和蜘蛛目,占總個(gè)體數(shù)的24.69%;其余稀有類群有15目,占總個(gè)體數(shù)的4.31%。橡膠林和森林中節(jié)肢動(dòng)物群落物種組成的稀疏曲線見圖2,森林中節(jié)肢動(dòng)物alpha多樣性(OTUs)明顯高于橡膠林。在森林和橡膠林中共得到397個(gè)OTUs,其中,片段森林中有347個(gè),橡膠林中有136個(gè)。片段森林和橡膠林中共有OTUs為86個(gè),片段森林節(jié)肢動(dòng)物特有OTUs為261個(gè),占物種總數(shù)的65.74%,橡膠林中節(jié)肢動(dòng)物特有OTUs為50個(gè),占物種總數(shù)的12.59%。森林節(jié)肢動(dòng)物群落組成中,鞘翅目(F=12.01,P=0.003)、雙翅目(F=12.43,P=0.003)、鱗翅目(F=7.365,P=0.015)和膜翅目(F=20.98,P<0.001)均顯著高于橡膠林(圖3)。

2.2 橡膠林和森林林下節(jié)肢動(dòng)物群落beta多樣性

單因素方差分析結(jié)果表明,林型對(duì)節(jié)肢動(dòng)物的beta多樣性(bray相異性指數(shù))有顯著影響(F=4.47,P=0.038),森林之間、橡膠林之間beta多樣性平均值分別為0.760和0.706。NMDS分析結(jié)果(圖4)表明,片段森林和橡膠林中節(jié)肢動(dòng)物物種組成存在差異。脅強(qiáng)系數(shù)(stress)為0.19,表明前兩個(gè)維度解釋了高百分比的排序。

2.3 森林和橡膠林鑲嵌景觀對(duì)林下節(jié)肢動(dòng)物多樣性影響的因素分析

用廣義線性模型對(duì)節(jié)肢動(dòng)物物種多樣性和不同尺度下的環(huán)境解釋變量做回歸分析,在所得到的模型中選擇ΔCwAI<2的模型,所選模型中的環(huán)境解釋變量為在該尺度下對(duì)節(jié)肢動(dòng)物多樣性有影響的環(huán)境因子。在景觀尺度(表3)上,最優(yōu)模型1環(huán)境解釋變量為森林覆蓋率和到森林距離;優(yōu)化模型2環(huán)境解釋變量為森林覆蓋率。模型預(yù)測(cè)結(jié)果表明,節(jié)肢動(dòng)物物種多樣性與森林覆蓋率的增加有關(guān)。在樣方尺度上,橡膠林中的多模型分析(表3)表明模型6最優(yōu),其環(huán)境解釋變量包括膠樹年齡和到森林邊緣距離。對(duì)于片段森林中節(jié)肢動(dòng)物物種(表3)而言,模型8最優(yōu),其環(huán)境解釋變量包括片段森林面積、樹種多樣性、到森林邊緣距離、坡向和海拔。

表3 影響節(jié)肢動(dòng)物物種多樣性(OTUs)的優(yōu)化模型

3 討論

3.1 森林和橡膠林對(duì)節(jié)肢動(dòng)物物種多樣性和beta多樣性的影響

在西雙版納納板河自然保護(hù)區(qū),橡膠林種植導(dǎo)致的橡膠林和森林鑲嵌性格局中,節(jié)肢動(dòng)物群落物種多樣性和群落結(jié)構(gòu)因土地利用方式的改變而變化。森林中節(jié)肢動(dòng)物的多度和物種豐富度均顯著高于與之相毗鄰的橡膠林。這一結(jié)果與其他關(guān)于土地利用變化對(duì)節(jié)肢動(dòng)物影響的研究結(jié)果一致,表明將自然林轉(zhuǎn)變?yōu)槿斯ち謺?huì)降低節(jié)肢動(dòng)物物種多樣性。森林雖然在結(jié)構(gòu)上已經(jīng)被破壞,但是其相對(duì)于環(huán)境單一、異質(zhì)性低的橡膠林,擁有更豐富的食物資源和更適宜的環(huán)境,利于對(duì)環(huán)境要求高的生物類群生存[24]。MENG等[17]對(duì)納板河流域甲蟲的調(diào)查結(jié)果顯示,橡膠林中甲蟲多度和物種多樣性顯著低于自然林;也有研究表明橡膠林步行蟲物種豐富度比森林低,原因可能是活動(dòng)能力強(qiáng)的昆蟲會(huì)選擇較適宜的生境[25]。

對(duì)節(jié)肢動(dòng)物beta多樣性的研究結(jié)果表明森林中節(jié)肢動(dòng)物與橡膠林中具有顯著差異,森林中的物種更替速率快于橡膠林。盡管由于取樣量的問題,在森林和橡膠林中節(jié)肢動(dòng)物取樣量均未達(dá)到可以對(duì)其進(jìn)行全面調(diào)查的程度,但是通過采樣量外推后,兩種生境中節(jié)肢動(dòng)物物種豐富度依然表現(xiàn)為森林顯著高于橡膠林,這說明對(duì)于節(jié)肢動(dòng)物群落分布,森林生境更加異質(zhì)化。HOUSTON等[26]研究發(fā)現(xiàn),白蟻群落組成受植被類型的影響,植被類型相似性越高,其中的白蟻群落組成相似性就越高。LIU等[19]對(duì)西雙版納地區(qū)螞蟻的研究結(jié)果也表明,橡膠林環(huán)境的同質(zhì)性導(dǎo)致橡膠林中螞蟻beta多樣性較低,即螞蟻的群落組成比較相似。在納板河地區(qū),各個(gè)片段森林因地理因素和植被類型不同,生境異質(zhì)性較高且存有豐富的節(jié)肢動(dòng)物物種多樣性。另一方面,在自然林向橡膠林轉(zhuǎn)變的過程中,環(huán)境受到很大破壞,較為單一的橡膠林生境只能保存對(duì)生境質(zhì)量要求不高的類群,所以橡膠林中節(jié)肢動(dòng)物物種的beta多樣性較低。

值得注意的是,即使橡膠林中節(jié)肢動(dòng)物物種數(shù)量和物種多樣性均小于森林,但橡膠林之間節(jié)肢動(dòng)物物種組成也存在差異,當(dāng)采樣量加大時(shí),橡膠林中也會(huì)發(fā)現(xiàn)更多的節(jié)肢動(dòng)物物種。這可能是由不同橡膠林樣地的不同管理方式所致,如殺蟲劑的使用與否,施肥量的不同以及刈割的頻率與程度不同等,使不同的橡膠林樣地所受干擾程度有所差異。另外,也可能是由于橡膠林的林齡不同,納板河流域大多數(shù)橡膠林都是低齡林,隨著橡膠林齡的增長,棲息地的生存環(huán)境質(zhì)量會(huì)下降,這對(duì)某些昆蟲物種極為不利,尤其對(duì)專性昆蟲和稀有種的影響更大。

3.2 橡膠林和森林鑲嵌格局中影響節(jié)肢動(dòng)物多樣性的主要環(huán)境變量

在景觀尺度上,研究顯示節(jié)肢動(dòng)物物種多樣性與森林覆蓋率顯著相關(guān)。當(dāng)森林覆蓋率大于50%時(shí),節(jié)肢動(dòng)物物種多樣性大幅度增加?？赡艿脑蚴窃谏指采w率大于50%的區(qū)域內(nèi),景觀受人類活動(dòng)的干擾較小,基質(zhì)連通性高,且植物物種資源豐富,可以容納更多的生物種群。

在樣方尺度上,地形因子坡向和海拔對(duì)森林節(jié)肢動(dòng)物群落組成有不同程度的影響。海拔、坡向和坡度等地形因子影響光照、溫濕度和養(yǎng)分等的空間分布,對(duì)生物多樣性也會(huì)產(chǎn)生影響[27-28]。其中,海拔影響著森林土壤、植被、小氣候以及樹木的生長和分布等,對(duì)昆蟲的生長、發(fā)育、存活和生殖力產(chǎn)生影響,進(jìn)而影響昆蟲群落的組成、結(jié)構(gòu)和功能[29];因橡膠的種植有人為因素的影響,多數(shù)橡膠林坡向向陽,所以筆者研究所調(diào)查的森林和橡膠林之間存在一定的海拔、坡度和坡向差異,給分析節(jié)肢動(dòng)物群落的分布格局帶來一定影響。

在橡膠林中樣方尺度上,橡膠林林齡和到森林邊緣距離對(duì)節(jié)肢動(dòng)物物種多樣性產(chǎn)生顯著影響。林齡小于20 a的橡膠林(低齡林)具有開放式林冠,郁閉度較低,而林齡在20～40 a(高齡林)的橡膠林具有較密閉的林冠層,郁閉度相對(duì)較高。研究發(fā)現(xiàn),低齡橡膠林中食蚜蠅的物種豐富度提高,而野蜂物種多樣性下降[17];步甲在低齡橡膠林中的個(gè)體數(shù)量雖與森林中相似,但在高齡橡膠林中的數(shù)量下降,這些都說明隨著膠樹年齡的變化,節(jié)肢動(dòng)物的群落組成也會(huì)有差別[30]。在森林中樣方尺度上,片段森林面積、樹種多樣性和樣點(diǎn)到森林邊緣距離是影響節(jié)肢動(dòng)物多樣性的主要環(huán)境變量。有關(guān)種-面積曲線的研究結(jié)果表明,與小片段島嶼相比,大片段島嶼可容納更多的物種數(shù),這是因?yàn)楦竺娣e的生境具有更高的環(huán)境異質(zhì)性,因此生物類群的數(shù)量和組成均顯著優(yōu)于小片段生境,較高的生境異質(zhì)性在提高物種多樣性的同時(shí),也為更多的特有種生存提供了可能[31-32]。小面積片段森林受生存資源不足、邊緣效應(yīng)即樣點(diǎn)到森林邊緣距離和樹種多樣性較低等因素的影響,生存環(huán)境穩(wěn)定性減弱,可能易導(dǎo)致物種間的競(jìng)爭(zhēng)加強(qiáng),促使動(dòng)物群落中優(yōu)勢(shì)物種過度繁殖,從而可能導(dǎo)致小片段森林中節(jié)肢動(dòng)物優(yōu)勢(shì)種群所占比例比大片段森林要高,但物種豐富度降低。

4 結(jié)論

(1)以橡膠林為背景的土地利用景觀格局中,橡膠種植選址應(yīng)充分考慮周邊森林的樹種多樣性,適當(dāng)保留向陽坡向的森林片段,提高森林與橡膠林鑲嵌格局的異質(zhì)性,將有利于節(jié)肢動(dòng)物多樣性的保護(hù)。

(2)在500 m研究尺度上,當(dāng)森林覆蓋率大于50%時(shí),節(jié)肢動(dòng)物物種多樣性顯著增加,即在該尺度上,森林覆蓋率至少要維持在50%以上,才能更好地保護(hù)節(jié)肢動(dòng)物多樣性。

(3)隨著膠齡的增加,開展多層次、多物種的人工橡膠林群落構(gòu)建模式,包括橡膠與茶、咖啡、可可、大葉千斤拔,甚至多物種混合的復(fù)合混農(nóng)林體系,豐富樹種多樣性將有利于保護(hù)鳥類、土壤動(dòng)物和昆蟲等群落的物種多樣性。