張春霞，馮自茂，李文鑫，楊培華，安向明，王明杰

(1西北農(nóng)林科技大學(xué) a國家林業(yè)和草原局油松工程技術(shù)研究中心，b林學(xué)院，陜西 楊凌712100；2延安市橋山國有林管理局，陜西 黃陵 727300)

立地質(zhì)量是影響林木生長發(fā)育的重要因素之一，評價不同立地條件下林木的生長對森林經(jīng)營和適地適樹造林具有非常重要的意義,根據(jù)不同立地類型條件下林木生長情況進行立地質(zhì)量評價是目前廣泛使用的評價方法之一[1-4]。不同立地條件下林木生長差別很大，比如生長較好地段的20年生油松優(yōu)勢木樹高可達8.6 m，而生長差的只有3 m左右[5]。由于立地類型影響林木的生長潛力，目前已有關(guān)于重要的林木如柚木、桉樹、油松、樺木、松類等立地指數(shù)表的編制、立地類型和立地質(zhì)量評價的研究報道[2-9]。但由于不同地區(qū)存在獨特的氣候因子、地形因子及其他因素的差異，比如溫度、降水、立地條件等，不同地區(qū)油松立地類型的劃分和立地質(zhì)量的評價結(jié)果存在較大差異[5,7,10-12]。對油松主產(chǎn)區(qū)的油松人工林進行立地條件評價和立地分類，并根據(jù)立地生產(chǎn)力確定不同的育林目標(biāo)和營林措施，依此指導(dǎo)荒山造林，是當(dāng)前亟待解決的一個問題。

油松(Pinustabulaeformis)是黃土高原地區(qū)的主要造林樹種，具有耐干旱、瘠薄，適應(yīng)性和抗逆性強等特點，已成為黃土高原人工植被的主要組成部分[13-15]。陜西省延安市黃陵地區(qū)是典型的黃土高原地形[16]，油松為黃陵地區(qū)主要的造林樹種[17]。本研究調(diào)查分析了陜西黃陵地區(qū)油松人工林林分的生長狀況，并對其立地質(zhì)量進行評價,以期為該地區(qū)油松人工林的經(jīng)營管理提供參考。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于陜西省延安市黃陵縣城以西45 km的雙龍林場，地處陜北黃土高原南部(108°45′～109°02′E、35°20′～35°50′N)。該區(qū)有溝壑、梁峁、殘原等特有的黃土地貌，黃土堆積較厚，為典型的黃土高原區(qū)域。該區(qū)海拔900～1 700 m，總體坡度較緩，平均坡度小于25°。年平均氣溫9.3 ℃，極端最高氣溫38 ℃，極端最低氣溫-23.7 ℃。年平均日照時數(shù)為2 528 h，初霜期多在10月上、中旬，終霜期多在4月下旬[18]。年平均降水量677.4 mm，屬暖溫帶半濕潤氣候區(qū)。土壤屬褐土地帶，多在暖溫帶大陸性半濕潤半干旱季風(fēng)氣候和森林條件下形成，成土母質(zhì)多為含碳酸鹽的黃土、沙巖和頁巖，土層深厚，一般在0.7～1.0 m，石灰淋溶和淀積作用明顯，黏化作用強烈，pH值7.0～8.0，腐殖質(zhì)層厚8～12 cm，土壤濕潤肥沃，有機質(zhì)含量15.262 g/kg，全氮含量1.06 g/kg，全磷含量1.06 g/kg。

1.2 研究方法

1.2.1 林地調(diào)查及取樣方法 采取典型樣地調(diào)查方法。2019年5月，在陜西黃陵雙龍林場設(shè)置52個20 m×20 m的標(biāo)準(zhǔn)樣地。樣地設(shè)置的具體原則為：(1)分別在林齡為20，30，40，50年的油松純林中，選取同一起源未進行間伐的區(qū)域設(shè)置樣地；(2)由于樣地內(nèi)未進行間伐，不同林分疏密度應(yīng)存在一定差異；(3)樣地選擇應(yīng)考慮不同的林齡、生長情況和立地條件，各種不同的立地類型應(yīng)涵蓋在所選樣地之中；(4)樣地為20 m×20 m，樣地內(nèi)至少有30株以上的油松。

1.2.2 林分因子調(diào)查 調(diào)查并記錄林分種源、起源、造林年度等；對樣地內(nèi)的每株油松進行每木檢尺，測定樣地內(nèi)油松的樹高、胸徑、枝下高、南北冠幅和東西冠幅。

1.2.3 環(huán)境因子調(diào)查 調(diào)查和記錄樣地地形環(huán)境因子，根據(jù)樣地所處的山形位置記錄坡位和坡向，采用 Vertex Ⅳ 測定坡度大小，羅盤儀確定方位角，GPS儀確定樣地的經(jīng)緯度和海拔高度。

1.3 內(nèi)業(yè)處理及數(shù)據(jù)分析

1.3.1 油松生長指標(biāo)分析 對不同林齡的油松進行樹高、胸徑、單株材積等生長指標(biāo)分析，根據(jù)每株油松的樹高和胸徑計算單株材積，計算公式為V=0.102 5πD2(H+3)[19]，其中V為單株材積，D為胸徑，H為樹高。

1.3.2 因子類目劃分 根據(jù)立地因子的變化幅度及其對林木生長發(fā)育的影響程度，并結(jié)合相關(guān)林業(yè)調(diào)查標(biāo)準(zhǔn)，對坡度、坡位、坡向和海拔等立地因子進行分級處理。由于林齡和林分密度也影響油松生長，所以對林齡和林分密度也進行了類目劃分，具體因子類目劃分及賦值見表1。

1.3.3 數(shù)量化理論Ⅰ模型建立 利用數(shù)量化理論Ⅰ，將林齡、林分密度、海拔、坡向、坡位、坡度等作為類目，以各因子不同的分類類型(表1)作為對應(yīng)的類別[20]。用每一類目對應(yīng)的觀察值為0或者1，建立[0,1]反應(yīng)矩陣表。基于52個油松樣地調(diào)查數(shù)據(jù)，以每塊標(biāo)準(zhǔn)樣地內(nèi)油松的樹高平均值為因變量，林齡、林分密度、海拔、坡向、坡位、坡度等因子為自變量，構(gòu)建數(shù)量化理論Ⅰ模型，并對油松人工林立地質(zhì)量進行評價。

表1 立地因子類目劃分及賦值

2 結(jié)果與分析

2.1 黃陵油松人工林油松生長指標(biāo)分析

由表2可知，20年生油松的平均樹高為6.70 m，平均胸徑為11.47 cm，平均單株材積為0.045 m3；30年生油松樹高、胸徑和單株材積平均值分別較20年生油松增加了3.31 m、3.26 cm和0.065 m3。在林齡30～40年，平均樹高增加了5.46 m，胸徑增加了5.16 cm，單株材積增加了0.140 m3。由此可以看出，在林齡30～40年，油松樹高、胸徑、單株材積較林齡20～30年增加明顯，而在林齡40～50年時平均樹高增加較小，僅為1.27 m，胸徑增加了3.02 cm，單株材積增加了0.110 m3，說明油松在40～50年林齡段的生長速度較慢。盡管不同林齡段油松生長指標(biāo)每10年增加的幅度不同，但是油松樹高、胸徑和單株材積均隨著林齡的增長而增加。

表2 不同林齡油松平均樹高、胸徑和單株材積

2.2 黃陵油松人工林原始數(shù)據(jù)分類及各因子方差分析

應(yīng)用數(shù)量化理論Ⅰ，根據(jù)52個標(biāo)準(zhǔn)樣地采集的數(shù)據(jù)建立油松平均樹高與林齡、林分密度、海拔、坡度、坡位和坡向之間的關(guān)系模型,并運用Forstat軟件對各立地因子及水平進行數(shù)量化綜合系數(shù)檢驗和方差分析統(tǒng)計[20]，結(jié)果見表3和表4。由表3可知，林齡對油松平均樹高的影響最大，這說明林齡對油松樹高的生長起著決定性作用。從表3還可知，模型經(jīng)F檢驗達到極顯著水平(P<0.001)，可在實際生產(chǎn)中應(yīng)用。

2.3 數(shù)量化理論Ⅰ模型的建立

根據(jù)樣地調(diào)查情況，結(jié)合立地因子的類目劃分(表1)，按照數(shù)量化理論Ⅰ模型計算各類目與油松樹高平均值的回歸系數(shù)及有關(guān)參數(shù)[20],以樹高(Y)為因變量，以林齡(X1)、林分密度(X2)、海拔(X3)、坡向(X4)、坡度(X5)和坡位(X6)為自變量構(gòu)建數(shù)量化理論模型，得到回歸模型為：

Y=2.47+X11+3.81X12+9.83X13+10.39X14+0.53X21+X22+1.66X31+1.97X32+X33+X41+0.98X42+0.90X43+1.21X51+X52+0.066X61+X62+0.11X63。

式中:Y代表樹高平均值;X11、X12、…、X63代表6個因子劃分的17個變量，其中X1代表林齡(X11、X12、X13、X14分別代表20，30，40和50年)，X2代表林分密度(X21和X22分別代表低密度和高密度)，X3代表海拔(X31、X32、X33分別代表低海拔、中海拔、高海拔)，X4代表坡向(X41、X42、X43分別代表陽坡、半陰坡、陰坡)，X5代表坡度(X51和X52分別代表緩坡和陡坡)，X6代表坡位(X61、X62、X63分別代表上部、中部、下部)，系數(shù)表示各個變量的得分值。模型的復(fù)相關(guān)系數(shù)R=0.956，達到極顯著水平，說明油松平均樹高與上述6個因子高度相關(guān)，所建回歸模型是可靠的。

表3 黃陵油松人工林原始數(shù)據(jù)分類

表4 黃陵油松人工林各立地因子的方差分析

由回歸分析結(jié)果(表5)及表2油松生長指標(biāo)測定結(jié)果可知，隨著林齡的增加油松樹高呈增加趨勢;根據(jù)得分值可以看出,對油松生長影響最大的因子為林齡，林齡與油松平均樹高的相關(guān)性達極顯著水平(P<0.001)。立地因子得分排序為：海拔>坡向>坡度>林分密度>坡位。將回歸模型的6個因子(林齡、林分密度、海拔、坡向、坡度、坡位)與油松平均樹高生長模型進行適用性檢驗，結(jié)果顯示6個因子與樹高平均值有極其緊密的相關(guān)性(P<0.001)，可用該方程來評價黃陵地區(qū)油松人工林的立地質(zhì)量。

表5 黃陵油松人工林各因子數(shù)量化回歸分析結(jié)果

2.4 黃陵油松人工林立地因子和樣地的得分值

根據(jù)數(shù)量化理論模型的得分范圍、方差比和偏相關(guān)系數(shù)判斷自變量的影響大小，數(shù)值大的作用程度大。由表5可知，在海拔>900～≤1 100 m時，隨著海拔的增加得分值增大，而海拔高于1 100 m時得分值降低，說明油松適宜在中低海拔生長。林分密度較低時油松生長較好，適當(dāng)?shù)氖璺ス芾韺Υ龠M油松人工林的生長很有必要。油松在陰坡和半陰坡的得分值高于陽坡，雖然油松是陽性樹種，但是往往在陰坡的生長狀況好于陽坡，這可能和陰坡的保水保肥能力高于陽坡有關(guān)。在緩坡的得分值高于陡坡，在下坡位和上坡位的得分值高于中坡位。綜合來看，黃陵地區(qū)油松適宜的生長環(huán)境為中低海拔、陰坡半陰坡、緩坡、上坡位或者下坡位。根據(jù)海拔、坡向、坡位、坡度等立地因子的得分值，計算出52個標(biāo)準(zhǔn)樣地的得分值見表6。

表6 黃陵地區(qū)油松人工林樣地的得分值

表6(續(xù)) Continued table 6

根據(jù)樣地得分值劃分立地質(zhì)量評價等級，將立地因子代數(shù)和的最大值和最小值均分為三等分，根據(jù)需要將樣地分為Ⅰ(優(yōu),得分值≥2.847)、Ⅱ(中，1.423<得分值<2.847)、Ⅲ(差,得分值≤1.423)3個等級。表6表明，52個調(diào)查樣地中有4個樣地的得分值低于1.423，評價結(jié)果為差，其他樣地評價為中等或者優(yōu)，即在所調(diào)查的樣地中，有92.3%的立地質(zhì)量為中等或者優(yōu)，說明黃陵地區(qū)油松人工林總體生長條件良好。

2.5 立地分類結(jié)果與立地質(zhì)量評價

根據(jù)海拔、坡度、坡位、坡向?qū)S陵地區(qū)油松人工林進行了立地類型劃分。先依據(jù)海拔高度劃為立地類型區(qū)，再根據(jù)坡度劃分立地類型組，最后依據(jù)坡位和坡向劃分立地類型，結(jié)果見表7。

表7 黃陵地區(qū)油松人工林立地分類與立地質(zhì)量評價

表7(續(xù)) Continued table 7

表7表明，52個油松人工林樣地可劃分為36種立地類型。利用得分值對其進行立地質(zhì)量評價，表7顯示，在所劃分的36種立地類型中，有75%立地類型的立地質(zhì)量處在中等及中等以上水平，且黃陵地區(qū)中低海拔(>900～≤1 100 m)生長的油松都處于中等以上水平。

3 討論與結(jié)論

海拔、坡向、坡度、坡位、林分密度等對油松樹高的生長都有一定的影響[21]。本研究采用數(shù)量化理論分析表明，海拔得分值所占比例最大，說明海拔對黃陵地區(qū)油松樹高的生長影響最大，這與段高輝等[7]的研究結(jié)果相一致。根據(jù)數(shù)量化理論Ⅰ模型構(gòu)建了黃陵地區(qū)油松優(yōu)勢木平均樹高預(yù)測方程，模型檢驗結(jié)果表明該方程可用于黃陵地區(qū)油松人工林的立地質(zhì)量評價，為該林區(qū)油松林的種植、撫育、管理等提供了理論基礎(chǔ)。

立地質(zhì)量評價結(jié)果顯示，黃陵地區(qū)油松人工林立地質(zhì)量評價得分為0～4.27，整體得分值較其他研究[7]偏低。這可能是由于本研究調(diào)查的是不同林齡階段的樣地，林齡顯著影響油松樹高的生長[22]，導(dǎo)致林齡的得分值最高，其他立地因子的得分值相對林齡偏低，所以樣地總體得分值較其他根據(jù)同齡林進行的數(shù)量化理論分析的得分值偏低。

采用數(shù)量化理論Ⅰ對黃陵地區(qū)油松人工林52個樣地的立地質(zhì)量進行了評價，結(jié)果表明有92.3%的油松標(biāo)準(zhǔn)樣地的立地質(zhì)量評價為“中等”或“優(yōu)”，說明黃陵地區(qū)油松人工林整體種植區(qū)域分布比較合理，且生長情況良好。立地質(zhì)量評價結(jié)果為“差”的有4個樣地，說明此區(qū)域不適宜栽植油松，可考慮栽植其他更適宜的樹種來提升土地的利用效率。從整體評價結(jié)果來看，立地評價結(jié)果較好的種植區(qū)域仍有較大的增長空間，比如林分密度對油松人工林平均樹高生長有很大的影響[23]，在林分經(jīng)營管理中應(yīng)該考慮以合理的種植密度來提升林分質(zhì)量。數(shù)量化立地質(zhì)量評價表可以解決黃陵地區(qū)油松林的統(tǒng)一評價問題，為實現(xiàn)黃陵地區(qū)油松人工林的合理經(jīng)營管理提供了一定的理論依據(jù)。

油松生長的影響因素很多，比如海拔、坡向、坡位、坡度，外界因素還有林分密度、土壤養(yǎng)分狀況等，另外還有降雨、光照等其他氣候因素[21,24-26]。由于本研究立地因子和樣本數(shù)量的限制，后期在研究區(qū)范圍和立地因子方面應(yīng)加大選擇力度，以增強油松人工林立地分類系統(tǒng)的準(zhǔn)確性和針對性。選擇合適的立地因子對立地質(zhì)量進行評價研究非常重要，立地因子不僅要反映林地生產(chǎn)潛力，也要便于獲取、計算與運用。本研究以油松平均樹高作為因變量，以實際調(diào)查的海拔、坡度、坡向、坡位等因子作為自變量來建立立地質(zhì)量評價模型，經(jīng)分析檢驗，表明該模型能真實反映林區(qū)立地的生產(chǎn)潛力。

已有研究表明，松類林木在林分密度適宜的情況下以油松的生長最好[23,27-28]。本研究結(jié)果也表明，林分密度對油松樹高生長有一定的影響，林分密度過大不利于油松的生長，合理的定植油松苗和疏伐管理能夠促進油松的生長發(fā)育。油松是喜陽樹種，有研究表明，黃龍山林區(qū)陽坡的油松生長狀況要優(yōu)于陰坡[7]。而本研究和其他一些研究[10,24]均發(fā)現(xiàn)陰坡和半陰坡油松的生長狀況優(yōu)于陽坡，這可能與不同研究選取的樣地位置不同有關(guān)。黃土高原地區(qū)比較干旱，陰坡土壤的保水保肥能力較強，從而導(dǎo)致陰坡油松的生長狀況較陽坡好，以后的研究需要綜合考慮地理環(huán)境因子、立地因子、土壤施肥保濕條件、氣候因子等多重因素。本研究數(shù)量化理論分析建模評價結(jié)果表明，黃陵地區(qū)油松最適宜在半陰坡、中海拔、下坡位且地勢較平緩的地方生長，此研究結(jié)果為黃陵地區(qū)油松人工林的營造和管理奠定了理論基礎(chǔ)。