環境因子對中條山油松人工林下物種分布的影響

　　摘要：【目的】量化環境要素對山西省中條山油松Pinustabuliformis人工林下物種分布影響，為進一步揭示中條山油松人工林下物種多樣性維持機制提供理論依據。【方法】基于中條山油松人工林樣地植物和土壤調查數據，分析環境因子對油松人工林下物種分布的影響�！窘Y果】環境因子對灌木、草本、藤本物種分布的解釋率(校正解釋率)分別為45.7%(34.5%)、37.0%(20.4%)、43.4%(30.6%);油松林下物種分布主要受非生物環境影響，海拔和土壤通氣性是影響該地區物種分布的關鍵環境因子;灌木和藤本分布主要受海拔影響，草本分布主要受土壤通氣性影響，且隨著海拔、土壤通氣性的升高油松人工林下物種多樣性下降。【結論】地形(海拔)和土壤(土壤通氣性)因子對區域群落生境的關鍵塑造作用所形成的環境篩是影響中條山油松人工林下物種分布的關鍵性因子。

　　關鍵詞：物種分布;物種多樣性;油松人工林;環境因子

　　物種多樣性的形成和維持機制是生物多樣性保護的基礎，長期以來，生態學家針對植物物種多樣性形成和維持的潛在機制先后提出了地質歷史過程、水熱動態假說、能量假說、寒冷忍耐假說及生境異質性假說等數十項假說[1]。森林群落的物種分布受多種因素影響，主要包括非生物環境因子、群落類型和空間距離等[2]，而林下物種又是人工林森林生態系統的重要組成部分，深入揭示環境因素對林下物種分布的影響，是解讀人工林森林生態系統物種多樣性的形成和維持的有效途徑之一。

　　油松Pinustabuliformis林是華北地區地帶性植被，也是華北地區頂級群落松櫟針闊混交林的重要組成部分。但由于人為活動干擾，該區域油松林以人工林為主，相比天然林，人工林存在病蟲害嚴重，生物多樣性、生產力、天然更新和生態穩定性下降等問題，森林生態系統健康面臨挑戰，系統原有的組成和結構被破壞，進而影響到生態過程和生態系統功能[3−4]。所以，開展環境因子對油松人工林下物種分布影響的研究，對于區域人工林生物多樣性保護和經營管理具有重要的意義。

　　大量研究[5]表明：地形、海拔、土壤、光照、植被類型等均對林下物種分布產生影響。對太岳山[6]與關帝山[7]油松林下物種分布研究表明：油松林下物種分布主要受非生物因子影響，受生物因子影響較小，且受地形因子的影響大于林分因子和土壤因子。以上2個研究地氣候類型為暖溫帶半干旱氣候區，而中條山處于暖溫帶半濕潤氣候區[8]，不同的氣候和土壤條件勢必會影響研究結果，對油松人工林經營管理帶來不確定性。基于此，本研究以山西中條山油松人工林為研究對象，基于植物和土壤調查數據，分析環境因子對影響油松人工林下物種分布的影響，為該地區油松人工林生物多樣性保護和經營管理提供數據支撐。

　　1 材料與方法

　　1.1 研究區概況

　　中條山位于山西省南部，是中國黃土高原向華北平原過渡的東南邊緣區，在氣候區劃上屬暖溫帶半濕潤季風氣候區。油松人工林分布區土壤為山地褐土及淋溶褐土，植被區劃上屬暖溫帶落葉闊葉林地帶，森林植被垂直分布明顯[9]。林下物種主要有連翹Forsythiasuspensa、土莊繡線菊Spiraeapubescens、披針薹草Carexlancifolia、早開堇菜Violaprionantha、茜草Rubiacordifolia和金線草R.membranace等。

　　1.2 樣地設置與植被調查

　　2018年6—8月，在對中條山油松人工林進行全面踏查的基礎上，根據海拔、林齡、密度、坡度、坡向、混交等選擇具有代表性的12個油松人工林樣地，每個樣地3個重復，共計36個喬木樣方。除2號樣地受條件限制，喬木樣方面積為15 m×15 m，其他各樣方面積為20 m×20 m，記錄樣方內出現的全部喬木種，并測定胸徑≥5 cm喬木樹種的名稱、胸徑、樹高。

　　將2號樣地劃分為9個、其他樣地劃分為16個5 m×5 m的灌木樣方，并于每個灌木樣方中心設置1個1 m×1 m的草本樣方，灌木、草本樣方數量合計均為555個，分別記錄灌木和草本樣方內的灌木、草本植物的種類、數量，并記錄灌木樣方內藤本的種類、數量。共記錄178種植物，其中灌木47種，草本112種，藤本19種，分屬56科124屬，以36個喬木樣方的物種多度數據建立物種矩陣(Z47×36、Z112×36、Z19×36)。記錄每個樣地的經緯度、海拔、坡度、坡向，根據中條山森林二類調查數據獲取林齡和間伐時間，采用樹冠投影法測定大樣方郁閉度。

　　1.3 土壤樣品采集與測定

　　在每個喬木樣地內沿對角線隨機選擇3個點為土壤樣品采集點，在每個采樣點挖掘1個土壤剖面，按0~10 、10~20、20~40 cm劃分土壤層次。在每個土壤層次上，用100 cm3的環刀取樣。采用土壤袋采集各層土樣200 g，挑出侵入物質，帶回實驗室自然風干、研磨、過篩后進行化學分析和測定。

　　測定土壤質量含水量、最大持水量、毛管持水量、田間持水量、毛管孔隙、非毛管孔隙和總孔隙度等土壤物理性質，測定方法按照LY/T 1215-1999《森林土壤水分-物理性質的測定》[10]進行。pH采用pH酸度計法(HANNA pH 211 )測定;土壤有機質采用重鉻酸鉀氧化-外加熱法測定;全氮采用凱氏定氮法測定;速效氮采用堿解擴散法測定;全磷和速效磷采用等離子發射光譜法(IRIS Intrepid Ⅱ XSP)測定。每個樣地土壤理化性質以所有測定值的平均值表征。

　　1.4 物種與環境矩陣建立

　　生物因子數據選取喬木樣方林分平均胸徑、平均樹高、胸高斷面積、油松相對優勢度、株樹密度、林齡、疏密度、郁閉度、喬木公頃蓄積、喬木公頃生物量、喬木平均蓄積、喬木平均生物量，蓄積量采用一元材積表計算[11]，生物量采用方精云等[12]生物量/材積轉換系數計算。非生物因子數據包括土壤理化性質、坡度、坡向、海拔、間伐間隔。其中坡向為無序分類變量，按方位劃分為北、東北、東、東南、南、西南、西、西北，分析時按1~8依次賦值。將以上因子構成環境矩陣(Z36×29)。

　　1.5 數據處理

　　采用主成分分析(PCA)法對環境數據降維;采用除趨勢對應分析(DCA)法去除“弓形效應”;采用典范對應分析(CCA)對物種進行排序。數據整理采用Excel 2016完成;物種多樣性分析、作圖及環境因子PCA分析均采用R 4.0.3軟件中“vegan”程序包完成;DCA和CCA采用Canoco 5軟件完成。

　　2 結果與分析

　　2.1 環境因子主成分分析

　　由于環境因子較多，且因子間存在多重共線性，需降低數據維度。因土壤pH對物種豐富度及空間分布影響較大[13]，故將除土壤pH以外的28個環境因子進行主成分分析。

　　根據因子載荷矩陣，第1主成分主要反映土壤水分;第2主成分主要反映林分密度;第3主成分主要反映土壤養分;第4主成分主要反映土壤通氣性;第5主成分主要反映生物量;第6主成分主要反映林冠郁閉度;第7主成分為海拔。將以上7個主成分與土壤pH共同建立生境矩陣(Z36×8)。將土壤水分、養分、通氣性、pH、海拔歸為非生物環境因子，將林分密度、生物量、林冠郁閉度歸為生物環境因子�；�于主成分分析結果構建環境矩陣，分析環境因子對林下物種分布的影響。

　　2.2 環境因子對林下物種分布影響分析

　　為排除冗余變量共軛效應對排序結果的影響，在CCA排序前先使用前向選擇，并進行蒙特卡羅顯著性檢驗，綜合檢驗結果確定可以進入CCA排序分析的環境因子。

　　所有環境因子均可以進入CCA排序分析，且環境因子對灌木、草本、藤本物種的解釋率(校正解釋率)分別為45.7%(34.5%)、37.0%(20.4%)、43.4%(30.6%)，因此分析中很大一部分的物種信息可能會丟失。物種與環境變量CCA分析模型變量的方差膨脹因子分別為8.91、7.98、9.32，均小于10，說明模型可靠。由不同環境要素對油松人工林下物種分布的解釋率可知：油松林下物種分布主要受非生物環境影響。

　　由生境因子與物種分布的相關顯著性檢驗可知：灌木、草本、藤本分布與環境因子極顯著相關(P<0.001)，前兩軸的方差貢獻率分別為53.63%、49.12%、63.93%，排序的大部分信息被反映�；�于油松人工林下灌木、草本、藤本與環境因子的CCA二維排序結果可知：油松人工林下灌木、藤本分布主要受海拔影響，其次是土壤通氣性，而草本分布主要受土壤通氣性影響，其次是林冠郁閉度，即隨著海拔、土壤通氣性的升高油松人工林下物種多樣性下降。林下灌木物種分布除與土壤養分和生物量無顯著相關關系外，與其他環境因子均顯著相關，且灌木主要分布在海拔相對較低、土壤通氣條件較好的油松人工林下，而偶見種黃蘆木Berberisamurensis(S19)、美薔薇Rosabella(S28)主要分布在海拔較高的區域。

　　林下草本物種分布除與土壤水分無顯著相關關系外，與其他環境因子均顯著相關，且草本主要分布在土壤條件較好、林內郁閉度較低油松人工林下，而偶見種中華秋海棠Begoniagrandissubsp.sinensis(H22)、大油芒Spodiopogonsibiricus(H110)主要分布在土壤通氣條件較差的區域，小花燈心草Juncusarticulatus(H26)、獨活Heracleumhemsleyanum(H96)主要分布在土壤養分較好的區域，羊齒天門冬Asparagusfilicinus(H101)主要分布在郁閉度較高的區域。

　　林下藤本物種分布除與林分密度、土壤養分、生物量無顯著相關關系外，與其他環境因子均顯著相關，且藤本主要分布在低海拔、土壤通氣性良好、林分密度相對較低的油松人工林下，而偶見種粉綠鐵線蓮Clematisglauca(L10)主要分布在土壤通氣條件較好的高海拔區域。

　　3 討論與結論

　　物種分布受多種因素影響。一般認為大尺度下，氣候差異是物種多樣性變化的最主要原因[1，14−16]，其中寒冷地區非生物因子對物種分布影響較大，溫暖地區生物因子影響相對較大[17−18];隨著緯度升高，物種多樣性降低[19]。而在小尺度區域內，地形、土壤、生物等因素是物種多樣性變化的主要驅動因素[2，20−21]。

　　中條山油松人工林下物種主要受非生物環境因子影響，由非生物環境因子(海拔)和土壤理化性質形成的局域尺度上的環境異質性，是影響局域內植物物種分布的關鍵因子，即非生物因子形成的環境篩(environmental filtering)作用是影響該區域林下物種的擴散和生存[22−24]。海拔和土壤通氣性是影響油松林下物種分布的關鍵性因子，其中灌木和藤本分布主要受海拔影響，草本分布主要受土壤通氣性影響。研究[25]表明：海拔是影響植物物種分布的關鍵性因子之一，是反映溫度、濕度、光照等因素的綜合環境因素，對物種多樣性的影響較復雜。

　　一般認為海拔升高，水熱條件將變差，物種多樣性將下降;但在某些特殊環境中，海拔升高后水分條件會有所改善，人為干擾活動下降，物種多樣性將升高[26];當上升到水熱條件限制植物生存的海拔，物種多樣性又下降[27];此外，植物物種多樣性隨海拔升高還呈現其他更復雜變化[28]。本研究發現：中條山物種分布與海拔顯著負相關，即海拔升高，林下物種多樣性下降。KRAFT等[29]研究發現：生態策略趨同的物種主要受地形環境篩控制，而生態策略趨異的物種主要受生物因素(競爭或密度制約過程)的影響。

　　海拔是中條山油松人工林下灌木和藤本物種分布的重要影響因子之一，說明該地區林下物種生態策略存在趨同性。在局域范圍內，除海拔外，土壤是主導植物物種分布格局的另一個關鍵非生物環境因子[30−32]。本研究發現：土壤通氣性對林下灌木和藤本物種分布的影響僅次于海拔，但對草本物種分布的影響大于海拔，即土壤因素對草本的影響大于地形因素，而對灌木和藤本的影響小于地形因素。

　　該結論與前人研究結果相反，之前研究結果認為，在局域范圍內，地形因素是影響植物物種分布最關鍵的因素，而土壤因素次之[26]，其主要原因是，海拔升高導致溫度、土壤肥力、土層厚度下降，輻射增加[33−36]，即海拔是主導環境改變的關鍵因子。本研究的土壤通氣性是土壤非毛管孔隙度、坡度、油松相對優勢度、喬木平均生物量等4個環境因子的集合體，且土壤非毛管孔隙度是影響該綜合指標的關鍵環境因子，此外，本研究區域各海拔土壤厚度變化不大，且海拔越高林分密度越大，輻射對林下植物影響較小，所以當土壤通氣性增大時，土壤的保水能力下降，甚至消失，土壤孔隙被空氣占據，植物根系無法吸收水分與養分，進而導致植物物種無法生存，所以當其值增大時植物分布將受到限制。

　　光照也是影響物種分布的主要環境因子之一[37]。本研究發現：林冠郁閉度對林下灌草藤物種分布均存在顯著影響，說明光照是限制該區域林下物種多樣性的另一個關鍵環境因子。劉浩棟等[38]研究發現：冠層郁閉度坡向是解釋海南島霸王嶺陸均松Dacrydiumpetinatum天然群落物種分布的主要驅動因素。本研究中，林冠郁閉度是郁閉度、坡向、間伐間隔的等3個環境要素的綜合體，所以在今后的森林經營活動中，可以通過開林窗的方式改變林下光照環境，增加林下物種多樣性。油松人工林下物種分布主要受非生物環境因子影響，其中地形(海拔)和土壤(土壤通氣性)是影響中條山油松人工林下物種分布的關鍵環境因子。

　　由于山地生態系統的水熱主要受海拔影響，海拔對物種的分布的影響占主導地位，因此隨著海拔和土壤通氣性的升高，林下物種多樣性下降，即地形和土壤因子形成的非生物環境篩作用影響區域群落植物物種的擴散和生存，使物種廣泛分布于適宜其自身生存的環境中。本研究校正解釋率較低，說明分析中很大一部分的物種信息可能丟失。造成解釋率較低的原因可能是未考慮一些生物、非生物環境要素的結果，如非生物因素中的溫度、降水量、光照等，生物因素中的植物功能性狀、競爭、人類活動等[39]。所以為了增加環境因子對中條山油松人工林下物種多樣性的解釋率，今后的研究需增加相應的環境因子。

　　4 參考文獻

　　[1]王健銘, 崔盼杰, 鐘悅鳴, 等. 阿拉善高原植物區域物種豐富度格局及其環境解釋[J]. 北京林業大學學報, 2019,41(3):14 − 23.WANG Jianming, CUI Panjie, ZHONG Yueming,etal. Biogeographic patterns and environmental interpretation of plantregional species richness in Alxa Plateau of northern China[J].JBeijingForUniv, 2019,41(3): 14 − 23.

　　[2]翁昌露, 張田田, 巫東豪, 等. 古田山10種主要森林群落類型的α和β多樣性格局及影響因素[J].生物多樣性, 2019,27(1): 33 − 41.WENG Changlu, ZHANG Tiantian, WU Donghao,etal. Drivers and patterns of α- and β-diversity in ten main forestcommunity types in Gutianshan, eastern China[J].BiodiversitySci, 2019,27(1): 33 − 41.

　　[3]劉世榮, 代力民, 溫遠光, 等. 面向生態系統服務的森林生態系統經營: 現狀、挑戰與展望[J]. 生態學報, 2015,35(1):1 − 9.LIU Shirong, DAI Limin, WEN Yuanguang,etal. A review on forest ecosystem management towards ecosystem services:status, challenges, and future perspectives[J].ActaEcolSin, 2015,35(1): 1 − 9.

　　作者：裴順祥，法 蕾，杜滿義，辛學兵

轉載請注明來自發表學術論文網：http://www.keysida.com.cn/nylw/29258.html