林分密度對杉木人工林下物種多樣性和土壤養分的影響*

張勇強，李智超，厚凌宇，宋立國，楊洪國，孫啟武?

林分密度對杉木人工林下物種多樣性和土壤養分的影響*

張勇強，李智超，厚凌宇，宋立國，楊洪國，孫啟武?

（中國林業科學研究院林業研究所，林木遺傳育種國家重點實驗室，國家林業和草原局林木培育重點實驗室，北京 100091）

為研究不同林分密度對杉木人工林下物種多樣性和土壤養分的影響，以38 a生5種密度杉木人工林為研究對象，調查林下植被，測定土壤理化性質及酶活性，對物種多樣性指數和土壤指標進行單因素方差、相關性及主成分分析。結果顯示：杉木人工林下物種種類多達121種，灌木層以杜莖山（）為主，草本層以雙蓋蕨（）、黑足鱗毛蕨（）等蕨類為主。灌木層和草本層物種多樣性各指數大多在初植密度為5 000 hm–2時最高。不同土層間各種土壤養分變化趨勢基本一致，不同密度間土壤養分變化趨勢不同，更多的土壤養分在高密度或低密度林）分下達到最大。除了pH、有機碳以及纖維素酶外，其他土壤養分受林分密度變化響應均顯著（< 0.05）。0～20 cm土層的土壤養分與草本層物種多樣性關系更為密切，而20～40 cm土層的土壤養分與灌木層物種多樣性關系更為密切。pH、全氮、堿解氮、有效磷與灌草層多樣性指數的相關性最為密切。主成分分析綜合得分顯示，初植密度為6 667 hm–2時最高（1.17），第二為3 333 hm–2（0.93），其次為5 000 hm–2（0.28）、1 667 hm–2（0.12）和10 000 hm–2（–2.49），得分前兩名顯著大于其他得分。以上結果表明，初植密度在5 000 hm–2更有利于林下物種多樣性的穩定，但不利于土壤養分的累積；密度過低或過高皆不利于土壤理化性質和植物多樣性的發展，特別是密度過高時，對林地傷害巨大；基于土壤理化性質和植物多樣性的主成分分析結果出現“駝峰模式”，杉木人工林初植密度6 667 hm–2和3 333 hm–2更適合土壤理化性質和植物多樣性的發展。

林分密度；杉木人工林；林下植被；土壤養分；主成分分析

林下植被是森林生態系統的一個重要組成部分，林下植被在維護森林生態系統多樣性、保持林地生產力、促進森林演替發展和維持生態功能穩定性等方面具有重要作用[1-3]。植被和土壤是一個密不可分的整體[4]，土壤是植被生長的基質，植被的生長狀況與土壤密切相關，林下植被以及凋落物的種類和數量也會影響土壤養分循環。

林分生長主要受林齡、立地質量和林分密度3個因子的影響，其中密度是人工林經營最重要的可控因子之一[5]。研究表明，林分密度對林下植被多樣性及土壤質量影響顯著[6-12]。國內外學者[11-12]對不同密度的天然林及人工林的研究顯示，密度較低時，林木占有資源空間上升，對土壤養分的利用增多，同時種群競爭壓力減小，對林下植被的抑制作用更加明顯，導致土壤養分和林下物種多樣性降低；而密度過大時，大量單株對土壤養分過分消耗且不能及時歸還，高密度和高郁閉度導致林下水肥氣熱等條件變差，同樣不利于養分積累及林下植被生存。我國對密度林的研究主要集中于馬尾松（）、刺槐（）、柚木（）、柏木（）等樹種，且不同林齡或不同樹種林下植被及土壤養分對密度的響應各不相同，最適營林密度差異顯著[7-10]。杉木（）是我國南方最主要的用材樹種，第八次全國森林資源清查結果表明，杉木人工林面積約占我國人工林總面積的1/4[13]。基于杉木在我國重要的生態價值和經濟價值，眾多學者對杉木進行了廣泛且深入的研究，研究方向主要包括：間伐、連栽對杉木人工林的影響、不同密度或不同林齡下杉木林下的土壤養分狀況、生物量以及碳儲量、杉木人工林退化成因分析等[3，5，14-17]。林分密度是人工林培育的重要因素，通過改變密度實現人工林光照和林分溫濕度等條件的改善，可促進林下植被種類、數量、生物量和土壤理化性質的提升[18-19]。研究者們針對不同密度杉木人工林下植被發育狀況、生物量、徑材、理化性質等方面做了大量研究[20-24]，研究結果表明，杉木人工林林分密度對林下物種多樣性和林分生長有顯著影響[23-25]，且當人工林下植被達到一定覆蓋度和生物量時，對林地土壤肥力有明顯改善作用[18]。

江西省是我國杉木人工林主產區，其分布區屬南方丘陵山地屏障帶，是我國“兩屏三帶”生態安全戰略格局骨架的重要組成部分。目前，針對該地區杉木人工林下植物多樣性及土壤養分影響的相關研究較缺乏，許多研究者對于林分密度的研究范圍在650～5 000 hm–2[7-8，18，26]，且研究樣地多為頭年間伐，次年調查，密度條件改變的時間短，對植被及土壤養分的影響較小，同時難以避免間伐等人為干擾對林地的影響。本文以1980年造林，初植密度為1 667、3 333、5 000、6 666以及10 000 hm–2的五種杉木人工林為研究對象，調查林下植被群落，分析了不同密度杉木人工林下物種多樣性及土壤養分的變化規律，從而更清晰地評價密度對于已經生長發育至較穩定階段杉木人工林的地上和地下部分的影響，目的在于尋找合適的群落結構，使杉木人工林更加接近天然條件的森林林，為杉木人工林林分結構改造和經營管理提供理論依據。

1 材料與方法

1.1 研究區概況及樣地選擇

研究區位于江西省新余市分宜縣西南部年珠林場（114°30′E～114°45′E，27°30′N～27°45′N）。該林場屬于低山丘陵地貌，氣候類型屬于亞熱帶季風氣候，全年平均溫度17.2°，年均降水量1 600 mm，年均無霜期270 d。雨量充沛，陽光充足，氣候溫和，無霜期長。母巖以千枚巖為主，土壤為黃、紅壤。地帶性植被為常綠闊葉林，喬木以杉木為主，林下層灌木主要以杜莖山（）、大青（cyrtophyllum）等為主，草本層主要以雙蓋蕨（）、黑足鱗毛蕨（）為主。

2018年10月，選擇1980年造五種不同密度（1 667 hm–2、3 333 hm–2、5 000 hm–2、6 667 hm–2、10 000 hm–2）杉木人工林作為實驗樣地，樣地海拔高度為190～210 m，坡向均為南坡，坡度約15°～20°。針對不同密度杉木人工林的研究樣地多通過人為間伐形成，有諸多間伐前后因素的干擾，而本次實驗的研究樣地為1980年造實驗樣地，該樣地在同一高度按水平分布分別造五種不同密度杉木林，每塊樣地造林面積約1 000 m2（33 m×33 m），并在相應樣地下方按照同樣造林密度設置兩組重復樣地，造林前三年進行撫育，后基本無人工干擾，為長期自然發展形成。樣地基本情況見表1。

表1 樣地基本概況

1.2 植被調查及土樣采集

在每個33 m×33 m樣地中間位置拉出20 m×20 m樣方，采用對角線法[9]，在每個樣方內選取5個2 m× 2 m的樣方調查灌木層，樣方內若有胸徑小于5 cm的小喬木記為灌木；同時選取5個1 m×1 m的小樣方調查草本層，記錄每種植物的種名、株數（叢數）、高度、蓋度。調查樣方內存活杉木數量，以計算存活率。

采用S形取樣法分別采集每個20 m×20 m樣地內5點的0～20 cm、20～40 cm土層土壤，同一土層土壤均勻混合，裝入土壤袋，風干過篩后進行土壤化學性質以及酶活性分析。在每塊20 m×20 m的樣地內用鋼制環刀（直徑5 cm，容積100 cm3）在不同土層各取3個原狀土，帶回實驗室做土壤物理性質分析。

1.3 土壤性質測定方法

土壤容重、含水率采用環刀法[28]測定。土壤pH采用電位法，用水浸提測定（水土比2.5︰1）；土壤有機質采用重鉻酸鉀氧化外加熱法測定；全氮用凱氏定氮法測定；堿解氮采用擴散法測定；全磷、有效磷采用鉬銻抗比色法測定；全鉀、速效鉀采用火焰光度計法測定[29]。土壤酶采用酶活性測定試劑盒（捷世康生物科技有限公司，青島）進行測定。

1.4 數據處理

林下灌草多樣性計算方法：

重要值（IV）：

IV=（相對多度+相對頻度+ 相對蓋度）/3（1）

豐富度指數：

辛普森（Simpson）多樣性指數：

香農-威納（Shannon-Wiener）多樣性指數：

式中，為物種總數；為物種的相對多度。

采用Microsoft Excel 2010、SPSS 20.0對數據進行單因素方差分析（Tukey’s HSD 檢驗，=0.05，=3）、相關性分析（皮爾森相關，< 0.05），用GraphPad PRISM?軟件繪圖。

2 結 果

2.1 不同密度杉木人工林下物種組成及重要值

重要值是以綜合數值來表示群落中不同種群的相對重要性，是研究群落結構的重要指標之一[30]。在所有樣方中，杉木人工林下出現的植物種類共56科98屬121種，其中，木本植物共39科72屬94種，草本植物17科26屬27種。分別對林下灌木層和草本層的重要值進行統計，結果表明：在灌木層中，杜莖山在各個密度林下均占有較大優勢，大青、珍珠蓮（）、毛冬青（）相對其他物種在各個密度林下稍占優勢；在草本層中，蕨類植物，尤其雙蓋蕨、黑足鱗毛蕨、邊緣鱗蓋蕨（）等占有較大優勢。在密度為1 667 hm–2時，菝葜（）和草珊瑚（）相對其他物種稍占優勢；在密度為3 333 hm–2時，五味子（）占有較大優勢，土茯苓（）、鉤藤（）稍占優勢；當密度達到10 000 hm–2時，灌木層出現玉葉金花（），還出現較之前密度林更多的木荷（）和刨花楠（），在草本層莎草（）成為占據較大優勢的草本。

2.2 不同密度杉木人工林下物種多樣性

將灌木和草本的豐富度指數、辛普森指數、香農-威納指數和均勻度（Peilou）指數分別按不同密度分組并進行單因素方差分析，結果表明，林分密度并不能顯著影響植物多樣性（表2），不同密度下灌木和草本的各指數雖各不相同，但組間差異不顯著（> 0.05）。

對各指數對密度變化的響應進行回歸分析（圖1），結果表明，灌木和草本物種豐富度變化趨勢在不同密度杉木林下不一致，灌木豐富度指數隨著杉木林密度的增大而增大，在杉木密度最大時灌木豐富度指數達到最高；草本豐富度指數則隨杉木密度的增大呈現先增大后減小的趨勢，在杉木中等密度時草本豐富度指數最高，在杉木密度較大和密度較小時草本的豐富度指數均較低。灌木的辛普森指數雖在密度為5 000 hm–2時最高，但整體趨勢表現為隨密度的增大而減小；灌木的香農-威納指數在密度為5 000 hm–2時最高，但隨林分密度的變化趨勢并不明顯；灌木的均勻度指數同樣在密度為5 000 hm–2時最高，但整體趨勢表現為隨密度的增大而降低。草本層的四種指數均在5 000 hm–2左右時達到最大值，且各指數變化趨勢基本一致，均表現為隨著密度的增大先增大后減小的趨勢，從圖1中也不難看出，當密度過高時，各指數的降低趨勢更加明顯。

表2 不同密度杉木林下灌草多樣性

注：同列不同字母表示不同密度間差異顯著（< 0.05），= 豐富度指數，= 辛普森指數，= 香農威納指數，sw=均勻度指數。下同Note：Different letters in the same column indicate significant difference between plots different in density（< 0.05）.= richness index，= Simpson index，= Shannon index，sw= Pielou index. The same below

圖1 不同密度杉木人工林物種多樣性指數

2.3 不同密度杉木人工林土壤養分

在0～20 cm土層，各個林分密度間杉木人工林土壤全氮、堿解氮差異顯著（< 0.05）；在20～40 cm土層，不同林分密度間杉木人工林土壤全鉀和速效鉀含量差異均顯著；不同密度杉木林的土壤全磷和有效磷含量在各個土層差異均顯著。隨著杉木林密度的增大，在0～20 cm和20～40 cm土層，pH、全磷、速效鉀均表現出先降低后升高的變化趨勢，且均在密度為10 000 hm–2時達到最大；土壤有機碳、全氮、堿解氮以及全鉀的含量隨著密度的增大變化趨勢不明顯，有機碳和全氮在密度為3 333 hm–2時達到最大，全鉀和堿解氮在密度為6 667 hm–2時達到最大；有效磷隨著密度的增大呈現先增大后減小的變化趨勢，在密度為5 000 hm–2時達到最大（表3）。

表3 不同密度杉木林下土壤化學性質

注：相同土層中，同列不同字母表示不同密度間差異顯著（< 0.05），SL = 土層，SOM = 有機質，TN = 全氮，TP = 全磷，TK = 全鉀，AN = 堿解氮，AP = 有效磷，AK = 速效鉀。下同Note：Different letters in the same column mean significant difference between plots different in stand density and the same in soil layer（< 0.05），SL = Soil layer，SOM = Soil organic matter，TN = Total nitrogen，TP = Total phosphorus，TK = Total potassium，AN = Alkaline nitrogen，AP = Available phosphorus，AK =Available potassium. The same below

不同密度杉木林土壤含水率在各個土層差異均顯著；在0～20 cm土層，多酚氧化酶和纖維素酶在各個密度下差異均顯著；在20～40 cm土層，酸性磷酸酶在不同密度下差異顯著。在0～20 cm和20～40 cm土層，容重、多酚氧化酶以及酸性磷酸酶隨著密度的增大呈現出先下降，后上升，最后又下降的趨勢，且均在密度為5 000 hm–2時達到最大；含水率隨著密度的增大先升高后降低最后又升高，在密度為10 000 hm–2時達到最大；脲酶隨著密度的升高先增高后降低，在密度為1 667 hm–2時最大；纖維素酶隨著密度增大的變化規律不明顯，在密度3 333 hm–2時最大（表4）。

表4 不同密度杉木林土壤物理性質及酶活性

注：SBD = 容重，SM = 含水率，PPO = 多酚氧化酶，UE = 脲酶，CL = 纖維素酶，ACP = 酸性磷酸酶。下同Note：SBD = Soil bulk density，SM = Soil moisture，PPO = Polyphenol oxidase，UE = Urease，CL = Cellulase，ACP = Acid phosphatase. The same below

2.4 杉木人工林植物多樣性與土壤性質的關系

從表5可以看出，在0～20 cm土層，pH與草本層的香農-威納指數以及均勻度指數均呈顯著負相關；堿解氮與草本層均勻度指數呈顯著正相關；有效磷與草本層的辛普森指數和香農-威納指數呈極顯著正相關，與草本層的豐富度指數和均勻度指數呈顯著正相關；多酚氧化酶與草本層辛普森指數呈顯著正相關。0～20 cm土層的各個土壤養分指標與灌木層物種多樣性指數均無顯著相關性。

在20～40 cm土層，pH與草本層辛普森指數呈顯著負相關；全氮與草本層豐富度指數呈顯著負相關；纖維素酶與草本層辛普森指數呈顯著正相關。全氮與灌木層辛普森指數、香農-威納指數、均勻度指數均呈顯著負相關；堿解氮與灌木層辛普森指數、香農-威納指數、均勻度指數均呈極顯著負相關；土壤容重與灌木層豐富度指數呈顯著負相關；脲酶與灌木層香農-威納指數呈顯著負相關；酸性磷酸酶與灌木層均勻度指數呈顯著正相關。其他土壤養分指標與林下物種多樣性無顯著相關性。

2.5 基于林下植被多樣性與土壤理化性質的綜合評分

五種不同密度杉木人工林下，基于林下植被多樣性與土壤理化性質共18個指標進行主成分分析，結果表明，第一主成分的方差貢獻率為39.47%，第二主成分的方差貢獻率為29.51%，第三主成分的方差貢獻率為21.34%，前三主成分的累積方差貢獻率為90.31%，足以代表所有指標的貢獻結果。因此，提取前三主成分進行分析，可以看出，第一主成分主要包含香農-威納指數、辛普森指數、有效磷、土壤容重、多酚氧化酶、均勻度指數、纖維素酶、堿解氮、酸性磷酸酶、全氮、全鉀和有機碳，第二主成分主要包括全氮、有機碳、堿解氮、纖維素酶、全鉀、全磷等，第三主成分主要包括酸性磷酸酶、均勻度指數、脲酶和含水率等。

根據因子載荷值和標準變量計算這三個主成分的權重，然后根據特征值計算每個林分密度的綜合得分，并對得分進行歸一化處理，得出不同密度杉木人工林基于物種多樣性和土壤理化性質的綜合得分（表6），由大到小依次為：初植密度6 667 hm–2（1.17）、3 333 hm–2（0.93）、5 000 hm–2（0.28）、1 667 hm–2（0.12）、10 000 hm–2（–2.49），說明初植密度10 000 hm–2或者1 667 hm–2均不利于杉木林下植被多樣性的建設和土壤理化性質的發展，特別是密度過大時，綜合得分極低，密度為5 000 hm–2時，對植物多樣性和土壤理化性質的綜合得分較6 667 hm–2或3 333 hm–2相差顯著。

表5 植物多樣性與土壤性質的相關系數

注：* 顯著相關（< 0.05），** 極顯著相關（< 0.01）。下同Note：* significant correlation（< 0.05），** extremely significant correlation（< 0.01）. The same below

表6 不同密度杉木人工林綜合評分

注：F1 = 第一主成分，F2 = 第二主成分，F3 = 第三主成分。Note：F1 = The first principal component，F2 = The second principal component，F3 = The third principal component.

3 討 論

3.1 密度對杉木人工林下物種組成及多樣性的影響

本研究地各個密度的杉木人工林下出現的植物種類共56科98屬121種，其中，木本植物39科72屬94種，草本植物17科26屬27種。本研究中杉木人工林下物種數量以及物種多樣性與其他研究地同齡杉木人工林的林下物種數以量及物種多樣性相比均較高[15，31]。在林下灌木層中，杜莖山在各個密度林下均占有較大優勢，大青、珍珠蓮、毛冬青相對其他物種稍占優勢；在草本層中，蕨類植物，尤其雙蓋蕨、黑足鱗毛蕨、邊緣鱗蓋蕨等占有較大優勢。說明這些物種對研究地杉木人工林林下環境的適應力很強。

本研究中，灌木的豐富度指數在林分密度為10 000 hm–2時達最大（表2），這與柚木、馬尾松等[7，8]研究結果不同，這一結果可能與杉木立木結構有關，單株杉木高大挺直，僅有林冠部分枝葉較多，較馬尾松或其他闊葉樹種枝葉稀少，太陽直射時仍有較高透光率，高密度種植增加了林下環境的復雜性，使林下灌木種類增多，但杉木林下灌木多樣性指數和均勻度指數均在5 000 hm–2時達最大（表2），與其他研究結果一致[8，28]。隨著密度的增大，林下草本豐富度指數、多樣性指數、均勻度指數均隨著密度的增大表現出先增大后減小的趨勢，其中，豐富度指數、辛普森指數、香農-威納指數均在密度為5 000 hm–2時達最大，均勻度指數在3 333 hm–2時最大（圖1）。灌木層和草本層物種多樣性大體呈現在5 000 hm–2時較高的趨勢，說明初植密度5 000 hm–2時的林內光照、水、熱等條件更加適合林下物種生長。在對不同密度的其他類型人工林的研究中，物種多樣性也有隨著密度增大而先增大后減小的趨勢，但是物種多樣性最高的林分密度與本結果相比均較低[15，16]。此外，除了在3 333 hm–2和6 666 hm–2密度下草本的均勻度指數高于灌木外，其他的各個物種多樣性指標灌木層均高于草本層（表2），在一定程度上說明灌木是較草本更適合當地環境特征的植被配置類型。

Alem等[32]發現人工林林分密度與林下植被物種多樣性無顯著關系，但是人工林下植被物種多樣性與到達植被層的光照有顯著相關性，這也說明了林分密度對林下植被的影響，主要是通過上層喬木對光照的截留，所到達植被的光照對林下植被生長發育產生影響，這一結論恰好解釋了本文中灌木豐富度隨林分密度增大而增大的結果（圖1）。本研究也發現林下灌草多樣性對于林分密度變化的響應不顯著（表2），這說明雖然林下物種多樣性在數值上各有高低，但密度對38 a自然生長的杉木人工林下植被的影響尚未達到顯著水平，這需要繼續監測植被與密度之間的關系，并且針對林分密度帶來的光照、溫濕度等改變展開更多研究。

3.2 密度對杉木人工林下土壤養分的影響

諸多研究者在對不同類型人工林的研究中均指出，林分密度對土壤容重、全氮等諸多土壤理化性質有顯著影響[8]。本研究中，不同密度杉木人工林對林下土壤含水率、全磷以及有效磷的影響顯著，對于0～20 cm土層的全氮、堿解氮、多酚氧化酶以及纖維素酶影響顯著，對于20～40 cm土層的全鉀、速效鉀以及酸性磷酸酶影響顯著（表3，表4）。

在不同的土層，同一種土壤養分含量隨林分密度的變化趨勢基本一致。對不同密度杉木林下的土壤理化性質分析表明，密度為1 667 hm–2、3 333 hm–2時有利于杉木林生長發育至成熟林齡內土壤有機碳、全氮、纖維素酶以及脲酶的積累（表4），這與胡小燕等[33]的研究結果一致；中密度（5 000 hm–2）有利于有效磷、多酚氧化酶、酸性磷酸酶的積累及土壤容重的改善；較高密度（6 667～10 000 hm–2）有利于全磷、全鉀、速效鉀、堿解氮的積累及pH、含水率的改善（表3）。對不同密度馬尾松林下土壤的研究[7]表明，土壤養分含量在密度為1 136 hm–2時取得最大值；而對不同密度柚木人工林下土壤的研究[8]表明，在密度為650～1 450 hm–2的范圍內，土壤理化性質隨著密度的增大先有所改善之后又退化。而依據本研究，杉木初植密度在5 000 hm–2以上更有利于本研究地杉木人工林土壤肥力的維持。

3.3 杉木人工林下物種多樣性與土壤養分的相關關系

根據不同土層土壤養分與林下物種多樣性的相關關系發現，0～20 cm土層的土壤養分與草本層物種多樣性關系更為密切，而20～40 cm土層的土壤養分與灌木層物種多樣性關系更為密切（表5）。這是由于草本植物覆蓋土壤表層且根系較淺，直接為土壤表層輸送養分，對表層土影響較大，而對較深層土壤影響很小[34]。灌木層的凋落物和根系分泌物會對土壤表層和較深層均有影響，且從本研究中可以發現，灌木層根系對較深層土壤養分的影響較灌木層凋落物對表層土壤養分的影響更大，因為灌木層植物多樣性指數與20～40 cm土層的土壤因子相關性更加密切。土壤各養分指標中，pH、全氮、堿解氮、有效磷與灌草層多樣性指數的相關性最為密切，而有機碳卻與灌草多樣性均不顯著（表5），可能是因為杉木林下氮、磷含量較低，植物多樣性受其限制更明顯[35-36]。

3.4 林分密度對杉木人工林經營的影響

主成分分析結果顯示，初值密度在6 667 hm–2時和3 333 hm–2時，土壤理化性質和植物多樣性分析的綜合得分明顯高于其他密度（表6），與其他研究中密度影響土壤理化性質和植物多樣性的單峰模式不同[1，19]，本研究出現了駝峰模式，表明綜合土壤和植物兩種因素，可能存在兩個差別較大的密度級別，均能夠較好保持土壤性質并增加生態系統中的植物多樣性，改善生態環境，這更有利于多目標培養杉木人工林，實現經營的可持續。例如，想實現大徑級用材林的需求，可采用初植密度3 333 hm–2；實現杉木量產而不要求材徑，可采用初植密度6 667 hm–2；若想實現生態經營，也可采用初植密度6 667 hm–2，可帶來更多的生態效益和經濟效益。本研究發現，低密度種植和高密度種植同樣不利于土壤養分的積累和植物多樣性的發展，低密度浪費營林空間，土壤養分被某些林下植被過度消耗，高密度營林樹種過度消耗地力，影響養分循環，同時不利于林下植被發育，對土地損耗巨大。依據不同密度水平下杉木的存活生長狀況、土壤養分含量以及林下物種多樣性這三個方面，本研究建議將初植密度6 667 hm–2作為杉木人工林造林的參考，在此初植密度之下可保持林下物種多樣性以及土壤養分的較高水平，且杉木的現存活率達到77%以上，現密度約5144 hm–2。在實際的營林過程中，要根據對杉木人工林的具體需求而采取如抽針補闊、補植鄉土樹種等不同的經營方式，以保持更加適宜的林分密度，促進杉木林群落結構向近自然群落方向發展。

4 結 論

各個密度下的杉木人工林下物種種類總數多達121種，灌木層以杜莖山為主，草本層以雙蓋蕨、黑足鱗毛蕨等蕨類為主。灌木層和草本層物種多樣性各指數大多在初植密度為5 000 hm–2時最高。不同土層間各個土壤養分變化趨勢基本一致，不同密度間土壤養分變化趨勢多樣，大多在較高密度下（≥5 000 hm–2）達到最大。除了pH、有機碳以及纖維素酶外，其他土壤養分受林分密度變化影響顯著。0～20 cm土層的土壤養分與草本層物種多樣性關系更為密切，而20～40 cm土層的土壤養分與灌木層物種多樣性關系更為密切。土壤各養分指標中，pH、全氮、堿解氮、有效磷與灌草層多樣性指數的相關性最為密切。綜合土壤理化性質和植被多樣性的主成分分析結果顯示，初植密度6 667 hm–2和3 333 hm–2較初植密度為5 000 hm–2更有利于二者的可持續發展。

［1］ Zou Y，He K N，Zhao C，et al. Relationship between plantation density and the plants diversity under forest ofalpine region. Acta Botanica Boreali-Occidentalia Sinica，2013，33（12）：2543—2549. [鄒揚，賀康寧，趙暢，等. 高寒區青海云杉人工林密度與林下植物多樣性的關系. 西北植物學報，2013，33（12）：2543—2549.]

［2］ Lin K M，Huang B L. Studies on β-diversity index of undergrowth plant in Chinese fir plantation. Chinese Biodiversity，2001，9（2）：157—161. [林開敏，黃寶龍. 杉木人工林林下植物物種β多樣性的研究. 生物多樣性，2001，9（2）：157—161

［3］ Sheng W T，Yang C D，Fan S H. Variation of soil properties of Chinese fir plantation. Forest Research，2003，16（4）：377—385. [盛煒彤，楊承棟，范少輝. 杉木人工林的土壤性質變化. 林業科學研究，2003，16（4）：377—385.]

［4］ Pohl M，Alig D，Korner C，et al. Higher plant diversity enhances soil stability in disturbed alpine ecosystems. Plant and Soil，2009，324（1/2）：91—102.

［5］ Che S H，Zhang J G. Comparative study on stand density measures based on self-thinning theory for Chinese fir plantation. Bulletin of Botanical Research，2012，32（3）：343—347. [車少輝，張建國. 基于自稀疏理論的杉木人工林密度指標研究. 植物研究，2012，32（3）：343—347.]

［6］ Ahmad B，Wang Y H，Hao J，et al. Optimizing stand structure for trade-offs between overstory timber production and understory plant diversity：A case-study of a larch plantation in Northwest China. Land Degradation & Development，2018，29（9）：2998—3008.

［7］ Sun Q H，Wu X，Wang M Z，et al. Effects of stand density on understory species diversity and soil physicochemical properties ofplantation. Chinese Journal of Applied Ecology，2018，29（3）：732—738. [孫千惠，吳霞，王媚臻，等. 林分密度對馬尾松林林下物種多樣性和土壤理化性質的影響. 應用生態學報，2018，29（3）：732—738.]

［8］ Zhou S P，Liang K N，Du J，et al. Effects of stand density on understory species diversity and soil physicochemical properties ofplantation. Bulletin of Botanical Research，2017（2）：42—52. [周樹平，梁坤南，杜健，等. 不同密度柚木人工林林下植被及土壤理化性質的研究. 植物研究，2017（2）：42—52.]

［9］ Mao Z H，Zhu J J，Liu Z G，et al. Effects of thinning on species diversity and composition of understory herbs in a larch plantation. Chinese Journal of Ecology，2006，25（10）：1201—1207. [毛志宏，朱教君，劉足根，等. 間伐對落葉松人工林內草本植物多樣性及其組成的影響. 生態學雜志，2006，25（10）：1201—1207.]

［10］ Jin Y Q，Bao W K. Relationships of the understory biomass with stand structure if the Sichuan cypress plantation forests across Sichuan Basin，China. Acta Ecologica Sinica，2014，34（20）：5849—5859. [金艷強，包維楷. 四川柏木人工林林下植被生物量與林分結構的關系. 生態學報，2014，34（20）：5849—5859.]

［11］ Xu X L，Sun Y J，Zhou H，et al. Effects of thinning intensity on understory growth and soil properties in Chinese fir plantation. Scientia Silvae Sinicae，2019，55（3）：1—12. [徐雪蕾，孫玉軍，周華，等. 間伐強度對杉木人工林林下植被和土壤性質的影響. 林業科學，2019，55（3）：1—12.]

［12］ Tu Y H. Dynamic characteristics of undergrowth vegetations inwith different management density. Journal of Northwest Forestry University，2005（4）：52—55. [涂育合. 杉木不同經營密度的林下植被變化. 西北林學院學報，2005（4）：52—55.]

［13］ Li Y Y，Zheng J Y，Shao M A. Comparison of the attributes of natural forests and plantations in Ziwuling Mountain. Acta Botanica Boreali-Occidentalia Sinica，2005，25（12）：2447—2456. [李裕元，鄭紀勇，邵明安. 子午嶺天然林與人工林群落特征比較研究. 西北植物學報，2005，25（12）：2447—2456.]

［14］ Sun D J，Wen Y G，Luo Y H，et al .Effect of close-to-nature management on species diversity in aplantation. Forest Research，2015，28（2）：202—208. [孫冬婧，溫遠光，羅應華，等. 近自然化改造對杉木人工林物種多樣性的影響. 林業科學研究，2015，28（2）：202—208.]

［15］ Chao L，Hong T，Li J，et al. Analyses on diameter class structure and species diversity of understory of artificial forest ofwith different forest ages in mid-subtropical region. Journal of Plant Resources and Environment，2015（2）：88—96. [巢林，洪滔，李鍵，等. 中亞熱帶不同林齡杉木人工林徑級結構與林下物種多樣性分析. 植物資源與環境學報，2015（2）：88—96.]

［16］ Ma R F，Liu J F，Zhang G S，et al. Soil ecological stoichiometric and microbial community characteristics offorest and their relevance analysis. Journal of Plant Resources and Environment，2015（1）：19—27. [馬瑞豐，劉金福，張廣帥，等. 格氏栲林土壤生態化學計量和微生物群落特征及其關聯性分析. 植物資源與環境學報，2015（1）：19—27.]

［17］ Chen Y，Lü G H，Li Y. Soil microbial functional diversity of rhizosphere and non-rhizosphere of three dominant herbaceous plants in the Dushanzi District. Acta Ecologica Sinica，2018，38（9）：3110—3117. [陳悅，呂光輝，李巖. 獨山子區優勢草本植物根際與非根際土壤微生物功能多樣性. 生態學報，2018，38（9）：3110—3117.]

［18］ Sheng W T. A long-term study on development and succession of undergrowth vegetations in Chinese fir（）plantations with different density. Forest Research，2001，14（5）：463—471. [盛煒彤. 不同密度杉木人工林林下植被發育與演替的定位研究. 林業科學研究，2001，14（5）：463—471.]

［19］ He Y L，Fu M Y. Review of studies on understorey of plantations. Forest Research，2002，15（6）：727—733. [何藝玲，傅懋毅. 人工林林下植被的研究現狀. 林業科學研究，2002，15（6）：727—733.]

［20］ Huang X S，Wu C Z，Hong W，et al. The relationship between stand density and biomass of two rotation Chinese fir plantations. Journal of Fujian College of Forestry，2011，31（2）：102—105. [黃賢松，吳承禎，洪偉，等. 2種杉木人工林密度與立木生物量的研究.福建林學院學報，2011，31（2）：102—105.]

［21］ Ye G F，Tu Y H，Lin R R，et al. On big-diameter-oriented cultivation techniques ofof different density measures. Journal of Beihua University（Natural Science），2005，6（6）：544—549. [葉功富，涂育合，林瑞榮，等. 杉木人工林不同密度管理定向培育大徑材. 北華大學學報（自然科學版），2005，6（6）：544—549.]

［22］ Cao X Y，Li J P，Yang J，et al. Stoichiometric characterization of soil C，N，and P of different age-group Chinese fir plantation. Soils，2019，51（2）：290—296. [曹小玉，李際平，楊靜，等. 不同齡組杉木林土壤碳、氮、磷的生態化學計量特征. 土壤，2019，51（2）：290—296.]

［23］ Wang H L，Liu J，Jiang P K，et al. Advancement in researches on effect of forest management on soil N2O emission in forest. Acta Pedologica Sinica，2018，55（1）：18—30. [王會來，劉娟，姜培坤，等. 營林措施對森林土壤N2O排放影響的研究進展. 土壤學報，2018，55（1）：18—30.]

［24］ Wu X Y，Wang W F，Yuan T T，et al. Carbon storage distribution characteristics of Chinese fir plantation in different age group and density. Forest Inventory and Planning，2016，41（4）：26—31. [武曉玉，王偉峰，袁婷婷，等. 不同齡組和密度杉木人工林碳儲量分配特征. 林業調查規劃，2016，41（4）：26—31.]

［25］ Liu Y B. Relationship between undergrowth species diversity and stand density of 29-year-old Chinese fir plantations. Journal of Fujian College of Forestry，2005，25（1）：27—30. [劉玉寶. 29年生杉木林下植物多樣性與密度的關系.福建林學院學報，2005，25（1）：27—30.]

［26］ Xiang C W. Effects of site quality and planting density on wood assortment rate in Chinese fir plantation. Forest Research，2015，28（5）：654—659. [相聰偉. 杉木人工林材種結構的立地及密度效應研究. 林業科學研究，2015，28（5）：654—659.]

［27］ Jia Y Y，He Z M，Zhou L L，et al. Effects of planting densities on the growth and space utilization of youngplantation. Chinese Journal of Ecology，2016，35（5）：1177—1181. [賈亞運，何宗明，周麗麗，等. 造林密度對杉木幼林生長及空間利用的影響. 生態學雜志，2016，35（5）：1177—1181.]

［28］ Yang Y Y，Chen Q B，Li J Q，et al. Vegetation diversity and soil biological fertility ofin central Yunnan plateau，China. Soils，2017，49（1）：90—96. [楊媛媛，陳奇伯，黎建強，等. 滇中高原華山松植物多樣性與土壤生物肥力特征. 土壤，2017，49（1）：90—96.]

［29］ State Forestry Administration. Forest soil analysis methods. Beijing：China Standard Press，2000. [國家林業局. 森林土壤分析方法. 北京：中國標準出版社，2000.]

［30］ Zhang H F. Species quantity change and ecosystem dynamics in the lower reaches of Tarim River. Chinese Journal of Ecology，2004，23（4）：21—24. [張宏鋒. 塔里木河下游植物群落的物種數量變化與生態系統動態研究. 生態學雜志，2004，23（4）：21—24.]

［31］ Lin K M，Yu X T，Huang B L，et al.Dynamical characterusitics of undergrowth plant diversity in Chinese fir plantations. Chinese Journal of Applied and Environmental Biology，2001，7（1）：13—19. [林開敏，俞新妥，黃寶龍，等. 杉木人工林林下植物物種多樣性的動態特征. 應用與環境生物學報，2001，7（1）：13—19.]

［32］ Alem S，Pavlis J，Urban J，et al. Pure and mixed plantations ofand：their growth interactions and effect on diversity and density of undergrowth woody plants in relation to light. Open Journal of Forestry，2015，5（4）：375—386.

［33］ Hu X Y，Duan A G，Zhang J G，et al. Effect of stand density on soil nutrient of Chinese fir mature plantations in south asia subtropical zone. Forest Research，2018，31（3）：15—23. [胡小燕，段愛國，張建國，等. 南亞熱帶杉木人工成熟林密度對土壤養分效應研究. 林業科學研究，2018，31（3）：15—23.]

［34］ Yu Y，Zou L，Sun T T，et al. Relationship between vegetation diversity and soil functional diversity in native pine forests. Pratacultural Science，2014，31（7）：1241—1247. [于洋，鄒莉，孫婷婷，等. 天然紅松林植被多樣性與土壤功能多樣性的關系. 草業科學，2014，31（7）：1241—1247.]

［35］ Andrew S M，Totland O，Moe S R. Spatial variation in plant species richness and diversity along human disturbance and environmental gradients in a tropical wetland. Wetlands Ecology and Management，2015，23（3）：395—404.

［36］ Chen C H，Tian D L，Fang X，et al. Correlativity between undergrowth vegetation species diversity/biomass and soil fertility of four types planted forest in suburbs. Journal of Soil & Water Conservation，2010，24（6）：213—217. [陳彩虹，田大倫，方晰，等. 城郊4種人工林林下植被物種多樣性、生物量與土壤養分相關性. 水土保持學報，2010，24（6）：213—217.]

Effects of Stand Density on Understory Species Diversity and Soil Nutrients in Chinese Fir Plantation

ZHANG Yongqiang, LI Zhichao, HOU Lingyu, SONG Liguo, YANG Hongguo, SUN Qiwu?

(Research Institute of Forestry, Chinese Academy of Forestry; State Key Laboratory of Tree Genetics and Breeding; Key Laboratory of Tree Breeding and Cultivation, State Forestry and Grassland Administration, Beijing 100091, China)

To study effects of stand density on understory species diversity and soil nutrients in Chinese fir plantations, a field survey was carried out in a Chinese fir plantation, 38 years old, varying in stand density in Dagang Mountain, Jiangxi Province.Five plots of woodlands different in stand density were delineated for investigation of understory species therein, and soil samples were collected from each plot for analysis of soil physical and chemical properties and soil enzyme activities. And One-way ANOVA analysis and correlation analysis and principal component analysis of the obtained data were performed for determination of relationships of understory species diversity and soil properties with stand density. Besides. principal component analysis and factor analysis were conducted to calculate scores of various index factors of the plantation, in an attempt to screen out a suitable density for growth of Chinese fir trees in plantations.Results show that there were 121 species of understory plants found under the Chinese fir trees in the plantation. The shrub layer was dominated by, while the herb layer was by ferns, such asandThe two layers both were the highest in diversity in the plot, 5 000 trees hm–2in initial density. The variation of soil nutrients did not differ much with soil layer, but did with stand density. All the soil nutrients, except pH, organic carbon and celllase, were the highest in the plots either the highest or the lowest in stand density and obviously varied with stand density. Soil nutrients in the 0～20 cm soil layer were more closely related to herbaceous diversity, while soil nutrients in the 20～40 cm soil layer were more to shrub diversity. The factors of pH, total nitrogen, alkalytic nitrogen, available phosphorus were the most closely related to the species diversity index of the understory. The principal component analysis shows that the plot 6 667 trees hm–2in initial stand density was the highest in comprehensive score (1.17), and followed by the plots 3 333 trees hm–2(0.93), 5 000 trees hm–2(0.28), 1 667 trees hm–2(0.12) and 10 000 trees hm–2(–2.49). obviously the first two were significantly higher than the other three in score.The research on plant diversity index and soil physicochemical properties reveals that the initial density of 5 000 trees hm–2is good to stability of understory species diversity, but not so to accumulation of soil nutrients. The principal component analysis shows that stand density, either too high or too low, is detrimental to the development of soil physical and chemical properties and plant diversity, especially when the density is too high, its damage to forest land is enormous, exhibiting a “hump model”, which indicates the initial densities of 6 667 trees hm–2and 3 333 trees hm–2are more suitable to the development of soil physical-chemical properties and plant diversity. The “Hump model” shows that taking into comprehensive account, the co-existence of different stand densities in a plantation may better maintain soil properties, increase plant diversity in the ecosystem and improve the ecological environment, thus contributing to multi-purpose cultivation and sustainable development of of Chinese fir plantations.

Stand density; Chinese fir plantation; Understory species; Soil nutrients; Principal component analysis

S714.2

A

10.11766/trxb201904080633

張勇強，李智超，厚凌宇，宋立國，楊洪國，孫啟武. 林分密度對杉木人工林下物種多樣性和土壤養分的影響[J]. 土壤學報，2020，57（1）：239–250.

ZHANG Yongqiang，LI Zhichao，HOU Lingyu，SONG Liguo，YANG Hongguo，SUN Qiwu. Effects of Stand Density on Understory Species Diversity and Soil Nutrients in Chinese Fir Plantation [J]. Acta Pedologica Sinica，2020，57（1）：239–250.

* 國家重點研發專項項目（2017YFC0505605）資助 Supported by the National Key Research and Development Program of China（No.2017YFC0505605）

，E-mail：soil2010@qq.com

張勇強（1994—），男，山東人，碩士研究生，主要從事森林土壤研究。E-mail：zyq199445@163.com

2019–04–08；

2019–05–31；

2019–07–05

（責任編輯：陳榮府）