祝子梟 劉兆剛 董靈波 陳瑩 張凌宇 由珈齊

(東北林業(yè)大學(xué)，哈爾濱，150040)

森林天然更新，是利用林木自身繁殖和恢復(fù)能力，在林地或林跡地上形成新一代幼林的過程，是森林生態(tài)系統(tǒng)中資源的再生產(chǎn)[1]。天然更新，從種子產(chǎn)生、擴散、萌發(fā)、幼苗定居和建成到衰老枯倒，每一個階段都受到各種環(huán)境因子的影響[2]。影響更新因素的研究，主要集中在母樹與種源[3-5]、林分因子[6-9]、土壤因子[10-11]、氣候因子[12-13]、干擾因子[14-17]等幾方面；對土壤因子方面，并未考慮不同土壤厚度、理化性質(zhì)對更新的影響。更新苗在各個生長階段，根系扎根深度不同，直接影響林分的更新。按高度級對更新苗進行等級劃分，但是高度劃分并未考慮不同樹種之間的更新差異。國內(nèi)對更新的影響因子分析方法，主要有相關(guān)性分析、主成分分析、多元回歸統(tǒng)計分析、方差分析、灰色關(guān)聯(lián)分析等[18-21]。冗余分析是生態(tài)學(xué)領(lǐng)域一種常見的分析方法，在林分更新方面應(yīng)用的較少，本文引入冗余分析法探討林分更新的影響因子。

落葉松(Larixgmelinii)林是大興安嶺地區(qū)的頂級群落之一，是大興安嶺地區(qū)分布面積最大、分布最廣的基本森林類型[22-23]，在海拔300～1 200 m均有分布。天然林主要分布在海拔500 m以上的山區(qū)，由于人為破壞嚴(yán)重，現(xiàn)有的落葉松林多屬于次生林。本文研究多種因素對大興安嶺天然落葉松林更新的影響，了解森林更新與環(huán)境因素之間的關(guān)系，旨在為科學(xué)實施森林可持續(xù)經(jīng)營管理、實現(xiàn)人工林經(jīng)營過程的近自然化提高參考。

1 研究區(qū)概況

研究區(qū)位于大興安嶺地區(qū)翠崗林場(123°41′～125°25′E，51°21′～52°10′N)、新林林場(123°41′～125°25′E，51°20′～52°10′N)、壯志林場(123°29′～124°18′E，50°59′～51°24′N)。大興安嶺地區(qū)溫度變化較大，冬季寒冷漫長，達7個月之久，晝短夜長，日照少，平均氣溫-28 ℃；冬季降水量較少，僅占年降水量的10%。夏季時間6—8月份，日照時間較長，可達17 h，晝夜溫差大；夏季由于東南季風(fēng)氣候的影響，降水多，主要集中在7—9月份。由于降水量的差異，大興安嶺冬季干燥，而冬季氣候濕潤溫和，初霜、終霜時間都較早。主要林分類型為白樺林、落葉松林、落葉松與白樺為主的針闊混交林。主要樹種有落葉松(Larixgmelinii)、白樺(Betulaplatyphylla)、云杉(Piceaasperata)、樟子松(Pinussylvestrisvar.mongolicaLitv)、山楊(Populusdavidiana)等；灌木層樹種有篤斯越桔(Vacciniumuliginosum)、杜香(Ledumpalustre)、興安杜鵑(Rhododendrondauricum)；地表草本有蕨類(Adiantumspp.)、莎草(Cyperusspp.)、舞鶴草(Maianthemumdilatatum)、地榆(Sanguisorbaofficinalis)等。

2 研究方法

2.1 樣地設(shè)置

本文以2018、2019年，在大興安嶺地區(qū)翠崗林場、新林林場、壯志林場設(shè)定的57塊固定樣地作為原始數(shù)據(jù)(見表1)，其中：翠崗林場10塊0.06 hm2、10塊0.1 hm2固定樣地，新林林場18塊0.1 hm2固定樣地，壯志林場10塊0.1 hm2、9塊0.08 hm2固定樣地。調(diào)查并記錄喬木層(胸徑≥5 cm)每木樹種、狀態(tài)、胸徑、樹高、第一活枝高、第一死枝高、位置坐標(biāo)(x、y)、東南西北4個方向的冠幅等信息；調(diào)查并記錄更新層(胸徑<5.0 cm)每木樹種、狀態(tài)、地徑、胸徑、樹高、位置坐標(biāo)(x、y)、更新方式(實生和萌生)等。每個樣地的四角點，分別設(shè)置1個2 m×2 m的灌木調(diào)查樣方、1 m×1 m的草本調(diào)查樣方，調(diào)查每個小樣方中灌木和草本的種類、平均高、蓋度等信息，并記錄樣地的海拔、坡向、坡位、坡度、土壤類型等信息。

表1 樣地概況

土壤樣品采集：在每個樣地沿樣地對角線取中心及兩樣方角的土壤樣品，分別取A0層(枯枝落葉層，分解、半分解的有機物質(zhì)積累的層次，木本植被下的森林土壤最為明顯)、A1層(腐殖質(zhì)層，由于腐殖質(zhì)的累積，腐殖質(zhì)和礦質(zhì)養(yǎng)料含量豐富，且結(jié)合緊密，多呈良好的團粒結(jié)構(gòu)，土色較深)、A2層(淋溶層，由于雨水的淋洗作用，土體中易溶性鹽類及鐵、鋁水化物、腐殖質(zhì)膠體受到淋失，向下移動，使該層腐殖質(zhì)及養(yǎng)分含量減少，土色較淺)土壤，蔭干后置于密封袋中，帶回實驗室蔭干后參照《土壤農(nóng)業(yè)化學(xué)分析方法》測定其化學(xué)性質(zhì)，測量并記錄A0、A1、A2層的厚度[24-25]。

2.2 更新幼苗等級劃分

結(jié)合更新苗等級劃分標(biāo)準(zhǔn)(幼苗，株高(h)<30 cm；小苗，30 cm≤h<130 cm；幼樹，h≥130 cm或胸徑≤5 cm)，將更新苗按照株高劃分為3個等級，即幼苗、幼樹、小樹。按照落葉松、白樺、不分樹種3種情況，對每個樣地更新幼苗分等級計數(shù)，分別代表各個樣地主要樹種不同生長階段的更新情況，從而對天然落葉松次生林更新情況進行分析。其中，更新密度劃分為6個類別：落葉松幼苗、白樺幼苗、落葉松小苗、白樺小苗、落葉松幼樹、白樺幼樹。

2.3 環(huán)境因子的選取

本文選取25個環(huán)境因子(見表2)，以詮釋地形、林分條件、空間結(jié)構(gòu)、樹種組成、土壤環(huán)境等方面的環(huán)境變化，通過盡可能多的環(huán)境因子探討其對天然落葉松次生林更新的影響。將25個環(huán)境因子劃分為2類，即林分因子12個指標(biāo)、土壤因子13個指標(biāo)。土壤樣品分為枯枝落葉層、腐殖質(zhì)層、淋溶層3層，對腐殖質(zhì)層、淋溶層分別測量其化學(xué)性質(zhì)，包括pH值、有機質(zhì)、全氮、全磷、全鉀指標(biāo)。

2.4 非對稱約束排序

冗余分析(RDA)，是一種將回歸分析法與主成分分析法融合進行排序的分析方法，也是對多響應(yīng)變量回歸分析的延伸拓展。冗余分析是響應(yīng)變量矩陣與解釋變量矩陣之間，多元多重線性回歸擬合矩陣的主成分分析。冗余分析是生態(tài)學(xué)領(lǐng)域常用分析方法，Legender和Gallagher開發(fā)了一系列物種轉(zhuǎn)化數(shù)據(jù)，使物種數(shù)據(jù)適用于冗余分析，此后，冗余分析得到了更加廣泛的應(yīng)用[26]。

表2 冗余分析使用的生境因子

為了尋找簡約模型，解決解釋變量之間的共線性問題，本文采用解釋變量的前向選擇。變量的共線性程度可以用方差膨脹因子(FVI)度量，當(dāng)FVI大于20時，說明與其他解釋變量之間存在嚴(yán)重的共線性問題。

變差分解，即量化兩組或多組變量單獨及共同解釋的變差。1992年，Borcard et al.[27]首先提出變差分解的概念，并且發(fā)現(xiàn)其分解過程，這對于多元生態(tài)數(shù)據(jù)分析有重要的意義。Peres-Neto et al.[28]提出使用校正決定系數(shù)(R2)促進變差分解的使用。通過變差分解，將解釋變量分為3部分，即林分變量解釋部分、土壤變量解釋部分、共同解釋部分。

3 結(jié)果與分析

3.1 冗余分析

對物種數(shù)據(jù)進行趨勢對應(yīng)分析(DCA)，獲得物種數(shù)據(jù)4個排序軸，其中最長軸梯度為2.1109。如果最長排序軸小于3.0，排序分析方法選擇應(yīng)用線性模型的冗余分析更合理[29]。

為了確認(rèn)不同區(qū)域之間的是否存在顯著差異，運用組合聚類分析與排序分析的方法，在冗余分析排序圖基礎(chǔ)上添加聚類樹(見圖1)。3組聚類簇均包含3個區(qū)域的樣地，并且3組聚類簇的樣方點在排序坐標(biāo)上分布比較均勻，說明大興安嶺這3個區(qū)域無顯著性差異。

冗余分析結(jié)果顯示，25個環(huán)境因子對林分更新密度的貢獻率為77.16%(見表3)，利用蒙特卡羅(Monte Carlo)檢驗(置換次數(shù)為999)，全模型前兩個排序軸顯著(p<0.05)。前4軸累計貢獻率為97.33%，前2軸累計貢獻率為85.63%，因此，可以用冗余分析前2個排序軸反映大興安嶺天然落葉松林的更新狀況(見表4)。第1排序軸主要反映了林分平均胸徑、平均樹高、林分蓄積、草本蓋度、林分密度、枯落物厚度、腐殖質(zhì)層有機質(zhì)的影響；第2排序軸主要反映了海拔、林分蓄積、林分郁閉度、樹種組成、角尺度、腐殖質(zhì)層全磷質(zhì)量分?jǐn)?shù)的影響。

表3冗余分析前4個排序軸的累計貢獻率

模型典范特征值總和/%前4個排序軸的累計貢獻率/%第1排序軸第2排序軸第3排序軸第4排序軸全模型77.1668.6185.6394.1997.33簡約模型57.1084.0797.0499.4599.98

表4 環(huán)境因子與排序軸的相關(guān)關(guān)系

更新林分可以分為兩組：第一組為落葉松幼苗、白樺幼苗、落葉松小苗，3者聚集于坐標(biāo)原點附近，表明這些林分在整個區(qū)域分布均勻或者分布于環(huán)境梯度的中值區(qū)域；第二組為白樺小苗、落葉松幼樹、白樺幼樹，3者箭頭比較長，比第一組分布較散，表明這些林分受個別環(huán)境因子影響較大(見圖2)。

圖面中的數(shù)字1～57為樣地編號。ELE為海拔；DBH為平均胸徑；H為平均樹高；N為林分密度；M為林分蓄積量；CW為郁閉度；SCI為樹種組成指數(shù)；AS為角尺度；SR為大小比數(shù)：MA為混角度；SC為灌木蓋度；HC為草本蓋度；LTH為枯落物厚度；HTH為腐殖質(zhì)厚度；STH為淋溶層厚度；PH1為腐殖質(zhì)層pH；Oma1為腐殖質(zhì)層有機質(zhì)質(zhì)量分?jǐn)?shù)；TN1為腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)；TP1為腐殖質(zhì)層全磷質(zhì)量分?jǐn)?shù)；TK1為腐殖質(zhì)層全鉀質(zhì)量分?jǐn)?shù)；PH2為淋溶層pH；Oma2為淋溶層有機質(zhì)質(zhì)量分?jǐn)?shù)；TN2為淋溶層全氮質(zhì)量分?jǐn)?shù)；TP2為腐殖質(zhì)層全磷質(zhì)量分?jǐn)?shù)；TK2為淋溶層全鉀質(zhì)量分?jǐn)?shù)；ND1為落葉松幼苗；ND2為白樺幼苗；ND3為落葉松小苗；ND4為白樺小苗；ND5為落葉松幼樹；ND6為白樺幼樹。

圖2環(huán)境因子對林分更新密度影響在冗余分析中的排序

全模型雖然能較大程度反映落葉松林更新苗木間的顯著差異，但是環(huán)境因子較多，大部分環(huán)境因子對更新狀態(tài)的影響很小，而且不利于后續(xù)分析。因此，采用解釋變量的前向選擇法進行變量篩選，尋找冗余分析簡約模型。簡約冗余分析結(jié)果顯示，枯落物厚度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度、草本層蓋度、林分平均5個因子對林分更新密度影響最大，因此這5個因子進入簡約模型(見圖3)。這5個環(huán)境因子(方差膨脹因子均小于2.5)對林分更新密度影響的累計貢獻率為57.1%，占25個環(huán)境因子對林分更新密度貢獻率的74%，前2個排序軸的累計貢獻率為97.04%(見表3)。

3.2 簡約模型的林分更新密度等值線

為了更加直觀的分析環(huán)境因子對林分更新密度的影響，提取冗余分析簡約模型中的變量，即林分平均高、林分郁閉度、枯落物厚度、草本蓋度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)5個環(huán)境因子。作出不同樹種、不同等級更新苗密度的冗余分析等值線圖(見圖4～圖9)，這5個環(huán)境因子對林分更新密度影響的貢獻率為57.1%，占全部環(huán)境因子對林分更新密度影響貢獻率的74%。

由圖4可見：落葉松幼苗更新密度等值線形成1個半封閉圓弧，涵蓋了枯落物厚度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度3個因子箭頭經(jīng)過區(qū)域，因而落葉松幼苗更新密度與這3個因子相關(guān)性不明顯。在垂直于落葉松幼苗更新密度等值線上，隨著平均高、草本蓋度增大，更新層密度逐漸減小。等值線密度整體變化很小，更新密度多為10株/hm2，說明落葉松更新幼苗稀少。

圖面中的數(shù)字1～57為樣地編號。LTH為枯落物厚度；TN1為腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)；CW為郁閉度；HC為草本蓋度；H為平均樹高；ND1為落葉松幼苗；ND2為白樺幼苗；ND3為落葉松小苗；ND4為白樺小苗；ND5為落葉松幼樹；ND6為白樺幼樹。

圖35個環(huán)境因子對林分更新密度影響在冗余分析中的排序

圖面中的數(shù)字1～57為樣地編號。LTH為枯落物厚度；TN1為腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)；CW為郁閉度；HC為草本蓋度；H為平均樹高。

圖4落葉松幼苗更新密度等值線

由圖5可見：白樺幼苗更新密度等值線，在坐標(biāo)原點呈橢圓，外圍呈半封閉圓弧，而橢圓等值線為零值，說明草本蓋度、林分平均高變化對白樺幼苗更新密度影響不顯著。在垂直于白樺幼苗更新密度等值線上，隨著枯落物厚度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度的增大，白樺更新幼苗密度增大。除了44、45、55、56、57號樣地外，其余樣地基本無白樺幼苗更新，說明翠崗、新林地區(qū)樣地白樺幼苗更新不良；壯志地區(qū)部分樣地雖然有白樺更新幼苗，但其更新密度也較低。

圖面中的數(shù)字1～57為樣地編號。LTH為枯落物厚度；TN1為腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)；CW為郁閉度；HC為草本蓋度；H為平均樹高。

圖5白樺幼苗更新密度等值線

由圖6可見：在垂直于落葉松小苗更新密度等值線上，更新密度從左至右逐漸增大；即更新密度，與枯落物厚度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度3個因子呈正相關(guān)，與草本蓋度、林分平均高呈負(fù)相關(guān)。

圖面中的數(shù)字1～57為樣地編號。LTH為枯落物厚度；TN1為腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)；CW為郁閉度；HC為草本蓋度；H為平均樹高。

圖6落葉松小苗更新密度等值線

由圖7可見：在垂直于白樺小苗更新密度等值線上，隨著草本蓋度、林分平均樹高的增大，白樺小苗更新密度減小。隨著枯落物厚度、郁閉度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)的增大，更新密度增大。從左至右，密度等值線逐漸變密，白樺小苗更新密度差異較大，大部分樣地白樺小苗更新不良，缺乏更新。

由圖8可見：沿著草本蓋度、林分平均高箭頭方向形成一層層半閉合圓?。谎刂捷S方向等值線密度逐漸變??；在垂直于落葉松幼樹更新密度等值線上，隨著草本蓋度、林分平均高的增加，落葉松幼樹更新密度先增加后減小。更新密度與枯落物厚度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度3個變量呈負(fù)相關(guān)。

由圖9可見：在垂直于白樺幼樹更新密度等值線上，隨著草本蓋度、林分平均高的增大，白樺幼樹更新密度減小。隨著枯落物厚度的增大，白樺幼樹更新密度增大。腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度箭頭與密度等值線基本保持平行，故腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度變化對白樺幼樹更新密度影響不顯著。

圖面中的數(shù)字1～57為樣地編號。LTH為枯落物厚度；TN1為腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)；CW為郁閉度；HC為草本蓋度；H為平均樹高。

圖7白樺小苗更新密度等值線

圖面中的數(shù)字1～57為樣地編號。LTH為枯落物厚度；TN1為腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)；CW為郁閉度；HC為草本蓋度；H為平均樹高。

圖8落葉松幼樹更新密度等值線

圖面中的數(shù)字1～57為樣地編號。LTH為枯落物厚度；TN1為腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)；CW為郁閉度；HC為草本蓋度；H為平均樹高。

圖9白樺幼樹更新密度等值線

結(jié)合圖4、圖5分析，落葉松和白樺幼苗均存在更新不良的問題，絕大部分樣地出現(xiàn)無更新幼苗的情況。結(jié)合圖6、圖7發(fā)現(xiàn)，白樺小苗更新密度整體好于落葉松小苗的更新密度；落葉松小苗更新密度各樣地間差異較小，而白樺小苗更新密度等值線疏密變化較大，即各樣地間更新差異較大。結(jié)合圖8、圖9分析，落葉松和白樺幼樹更新密度各樣地間相差不大，二者均存在更新不良的問題。

整體看，林分更新密度與林分平均高、草本蓋度呈負(fù)相關(guān)，與枯落物厚度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度3個因子呈正相關(guān)。其中：落葉松幼樹與其他苗木更新密度稍有差異，落葉松幼樹更新密度與林分平均高、草本蓋度相關(guān)性不呈線性增長趨勢，呈先增大后減小的趨勢；與枯落物厚度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度呈負(fù)相關(guān)。

3.3 環(huán)境因子對林分更新密度的影響

為了量化林分因子、土壤因子兩組因子單獨及共同影響的差異，本文采用變差分解的方法，分別對林分因子、土壤因子進行解釋變量前向選擇。其中，林分因子簡約組合為林分平均高、草本蓋度、郁閉度，土壤因子簡約組合為枯落物厚度、淋溶層厚度、淋溶層全磷質(zhì)量分?jǐn)?shù)、腐殖質(zhì)層厚度。林分因子單獨對林分更新密度影響的貢獻，占環(huán)境因子對林分更新密度影響總貢獻率的72.28%；而土壤因子對林分更新密度影響貢獻率較低，只占環(huán)境因子對林分更新密度影響總貢獻率的13.14%；兩組因子共同對林分更新密度影響的貢獻率，占環(huán)境因子對林分更新密度影響總貢獻率的14.58%(見圖10)。總體看，林分因子對林分更新密度的影響，大于土壤因子對林分更新密度的影響。

4 討論

張樹梓等[30]研究表明，Ⅲ齡級林分天然更新苗在早期生長中，受土壤枯落物因子限制比較大，隨著生長繼續(xù)，林分結(jié)構(gòu)因子逐漸代替土壤枯落物因子成為主要的更新限制因子，而這種現(xiàn)象在Ⅳ齡級林分中體現(xiàn)不明顯。王笑梅等[31]研究發(fā)現(xiàn)，當(dāng)冬青重要值較高時，林分結(jié)構(gòu)因子是主要影響因子，土壤因子影響較?。划?dāng)冬青重要值較低時，土壤因子的影響大于林分結(jié)構(gòu)因子。而本文的變差分解結(jié)果表明，林分因子對林分更新密度的影響大于土壤因子。根據(jù)林分更新密度等值線圖(見圖4、圖5)發(fā)現(xiàn)，大興安嶺地區(qū)更新幼苗低萌發(fā)量，導(dǎo)致在排序時，無法體現(xiàn)土壤因子對更新苗在前期生長的限制影響；而更新小苗、幼樹的密度與更新幼苗密度相比，明顯增大；因此結(jié)果顯示，林分因子對天然更新的影響，大于土壤因子。此外，本文只選取了土壤因子6個指標(biāo)，即全氮、全鉀、全磷、pH值、有機質(zhì)、土壤厚度，還不能全面反映林分的養(yǎng)分狀況和物理性質(zhì)。實際林分中，對更新產(chǎn)生影響的土壤因子還有很多，在以后的研究中，應(yīng)對土壤的速效鉀、有效磷、孔隙度、凋落物持水量等方面進一步研究，從而全面揭示土壤因子對更新的影響。

林分更新與環(huán)境因子的關(guān)系比較復(fù)雜，不同的研究結(jié)果不盡相同。灌草植被的地上、地下競爭，改變了林下光照、土壤濕度、溫度以及凋落物的分布[9]。林下層對幼苗出現(xiàn)和定居的干預(yù)，不僅直接與其進行資源競爭，也間接影響種子捕食者行為，從而對幼苗定居與生長產(chǎn)生影響[32-33]。不同灌草的地上、地下競爭，會影響更新苗的光合作用與生物量的積累[34]。雜草的地上競爭加劇更新苗的空間競爭壓力，并限制更新苗的光合作用；地下競爭又爭奪土壤養(yǎng)分與水分，導(dǎo)致更新苗的光補償點提高，凈光合速率下降限制幼苗的生長發(fā)育。即隨著林下草本蓋度增加，更新密度呈下降趨勢。

林地凋落物對種子萌發(fā)、幼苗建立和生長有雙重作用[35]。有研究表明，枯落物層對種子存在自毒作用[36]。即使有充足的種源也很難形成更新，種子散落在枯落物層，林內(nèi)濕度造成種子發(fā)霉，失去活力而無法萌發(fā)?？萋湮锔采w的種子存活率明顯高于裸露微生境，且與枯落物厚度呈正相關(guān)[37]。枯落物對種子有保護作用，防止動物覓食。此外，枯落物層導(dǎo)熱不良，防止土壤升溫過快而造成水分蒸騰，枯落物分解后還能為土壤提供營養(yǎng)元素。白樺種子和落葉松通常停留在枯落物層[38-39]，這也驗證了本文更新密度與枯落物呈正相關(guān)的結(jié)論。

土壤氮素質(zhì)量分?jǐn)?shù)是幼苗生長的限制因子之一[40]。宋新章[41]研究發(fā)現(xiàn)，全氮質(zhì)量分?jǐn)?shù)對幼苗的密度有顯著的決定性作用，全磷質(zhì)量分?jǐn)?shù)對幼苗的徑生長有較大的促進作用，全鉀質(zhì)量分?jǐn)?shù)對幼苗的高生長有很強的促進作用。而本文也證明了更新密度與腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)呈正相關(guān)。

在郁閉度小于0.7時，白皮松林下天然更新密度較大；當(dāng)郁閉度大于0.8時，天然更新密度小，而且成活率降低[42]。興安落葉松與白樺皆為陽性樹種，當(dāng)林分郁閉度過高時，林下光照不足，將限制白樺與落葉松更新，大部分更新未進入主林層之前便逐漸枯死。但當(dāng)郁閉度過小時，林下光照過強，土壤溫度過高，水分流失較大，也不利于林分更新。而大興安嶺大部分林分郁閉度偏低，因而結(jié)果顯示，林分更新密度與郁閉度呈正相關(guān)。

葉面積指數(shù)與林分平均高呈顯著負(fù)相關(guān)[43]。葉面積指數(shù)越大，林分郁閉程度越差，即林分平均高與郁閉度呈負(fù)相關(guān)，更新密度與林分平均高呈負(fù)相關(guān)，本研究結(jié)果與其他學(xué)者研究結(jié)果一致[44]。

5 結(jié)論

大興安林地區(qū)，主要影響林分更新的環(huán)境因子是：林分平均高、林分郁閉度、枯落物厚度、草本蓋度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)，這5個環(huán)境因子對林分更新密度影響的貢獻率為57.1%。從變量分解結(jié)果看，林分因子對林分更新密度的影響大于土壤因子。

大興安嶺地區(qū)主要樹種更新不良，缺乏幼苗更新。3個區(qū)域的更新無明顯差異，落葉松、白樺樹種之間更新密度差異明顯，幼苗、小苗、幼樹3個等級更新密度差異明顯。整體看，大興安嶺地區(qū)主要樹種更新密度，與林分平均高、草本蓋度呈負(fù)相關(guān)，與枯落物厚度、腐殖質(zhì)層全氮質(zhì)量分?jǐn)?shù)、郁閉度呈正相關(guān)。