馬友平

(湖北民族學(xué)院，恩施，445000)

森林立地質(zhì)量是指具有相同立地條件的各個地段的綜合，是森林經(jīng)營管理中林地分類的重要依據(jù)。同一立地條件類型的林地具有相同的生產(chǎn)力并可采用相同的造林和育林技術(shù)措施。立地質(zhì)量在生態(tài)學(xué)上稱為生境，生境一詞首先由美國的Grinnel于1917年提出［1］。一定區(qū)域內(nèi)生境的好壞，一般都用生境質(zhì)量來表示，而在林學(xué)上多用立地質(zhì)量表示，可據(jù)一定時期內(nèi)林地生產(chǎn)林木的最高產(chǎn)量來測算。立地質(zhì)量在一定條件下決定了營林工作者在經(jīng)營林分中所面臨的經(jīng)營類型和范圍等問題及其可能性，是研究森林生長環(huán)境以及環(huán)境對于森林類型和生產(chǎn)力影響的一個重要手段［2-3］，是很多環(huán)境因子的總和，例如，土層厚度、坡度、坡向、樹種和其他因子等。開展森林立地分類研究，是因地制宜、適地適樹、合理利用林地、實現(xiàn)森林集約經(jīng)營的重要手段，可為科學(xué)造林、合理開發(fā)利用區(qū)域林地資源、充分挖掘林地生產(chǎn)潛力、實現(xiàn)森林生態(tài)系統(tǒng)的平衡、改善區(qū)域生態(tài)環(huán)境、科學(xué)制定營林規(guī)劃、實現(xiàn)林業(yè)經(jīng)營的可持續(xù)發(fā)展等提供基礎(chǔ)依據(jù)。

立地在測樹學(xué)上稱做“生境”，是森林或其他植被類型生存的空間及與之相關(guān)的自然因子的綜合，與生長于其上的樹種無關(guān)。但是立地質(zhì)量則是指在某一生境上既定森林或者其他植被類型的生產(chǎn)潛力，所以生境質(zhì)量與樹種相關(guān)聯(lián)，并有高低之分。一個既定的立地，對于不同的樹種來說，可能會得到不同的立地質(zhì)量［4］。張宋智等［5］進行了甘肅小隴山林區(qū)日本落葉松速生豐產(chǎn)林立地類型劃分及評價的研究，應(yīng)用多元回歸分析了小隴山林區(qū)日本落葉松幼林生長與立地因子的關(guān)系，篩選出影響日本落葉松幼林生長的主要立地因子是土壤腐殖質(zhì)層和有效土層厚度、坡度、海拔，坡位、石礫含量等因子被剔除;荀守華等［6］研究了日本落葉松生長量與立地因子的相互關(guān)系，以山東省膠東半島的嶗山林場和昆侖山林場的日本落葉松人工林為例，所選的5個因子中，坡位對日本落葉松生長量的影響不顯著，海拔、坡度、坡向、土層厚度4個因子影響顯著，海拔對日本落葉松生長量有影響，且呈正相關(guān)，說明日本落葉松生長量隨著海拔的升高而增大。分析其原因，一是因為日本落葉松原產(chǎn)地位于海拔1 000～2 500 m的日本富士山山區(qū)，它對高海拔山地有天然的生態(tài)適應(yīng)性。田志和等［7］在《日本落葉松育林學(xué)》一書中比較翔實地敘述了日本落葉松各生境因子間的關(guān)系。這些研究都從定量或定性的角度研究了不同生境因子與日本落葉松生長量的關(guān)系，但對日本落葉松立地質(zhì)量的定量研究僅限于生境因子的篩選，或僅考慮以林木的樹高為指標(biāo)來劃分，缺乏同時考慮多種因子，對每一個小的地塊進行精準(zhǔn)的量化研究，因此，本研究提出了運用灰色定權(quán)聚類法對湖北省建始縣長嶺崗林場日本落葉松人工林立地質(zhì)量的研究，充分考慮日本落葉松的生長特性和各生境因子，且實現(xiàn)了全程量化。

1 研究地概況

湖北省建始縣國營長嶺崗林場(包括長嶺崗場部和擦擦坡兩部分)地理位置為東經(jīng) 110°00'～110°04'，北緯 30°47'～30°50'，屬于云貴高原東緣，巫山山脈余脈，為長江、清江分水嶺，海拔1 500～1 920 m，以經(jīng)營日本落松為目標(biāo)，林下有少量天然植被馬尾松、化香、胡枝子、鹽膚木等。土壤類型為山地黃棕壤，土層厚度60～90 cm，pH值5.3～6.0，質(zhì)地為輕壤、中壤或重壤，有機質(zhì)含量高，缺磷。當(dāng)?shù)啬昶骄鶜鉁?1.70℃，極端最高溫度28.80℃，極端最低溫度 -17.20℃，大于10℃積溫2 244℃(5月22日至9月27日)，夏無酷暑，冬無嚴(yán)寒。春秋季4月30日至10月3日，冬季10月4日至翌年4月29日，沒有春夏季之分。無霜期4月6日至10月23日，最大積雪深度50 cm，最大凍土深度33 cm(1—2月份)，結(jié)冰期11月9日至翌年4月12日;年降水量1 884.3 mm(7—9月份最多)，年蒸發(fā)量933.3 mm，年平均相對濕度85%，年日照時數(shù) 1 533 h，年平均風(fēng)速 3.8 m·s-1。

2 材料與方法

為了對日本落葉松人工林的立地質(zhì)量進行量化研究，除了借用林場已有的3塊固定樣地資料外，設(shè)置了臨時標(biāo)準(zhǔn)地40塊，共計43塊標(biāo)準(zhǔn)地。據(jù)臨時標(biāo)準(zhǔn)地的位置和地況，面積為0.04～0.10 hm2，形狀為方形或矩形，在每塊標(biāo)準(zhǔn)地內(nèi)進行每木檢尺，計算林分的平均高和優(yōu)勢木高;同時記載標(biāo)準(zhǔn)地的坡度、坡位、海拔、坡向、土壤類型、林下植被、pH值，在具有代表性位置挖土壤剖面，測其土層厚度。用手持GPS72測定其坐標(biāo)位置。土層厚度的測定除了在43塊樣地中測定外，還采用隨機的方法在全場設(shè)置了土壤厚度測定點296個，即土壤厚度測定點總計為339個。收集林場1∶10 000BJ54地形圖6幅，為數(shù)字高程模型的建立做準(zhǔn)備;同時為確定各因子的權(quán)重，制定了專家調(diào)查表，以電子郵件或當(dāng)面發(fā)放的辦法填寫，共收回完全符合要求的專家調(diào)查表63份。

3 結(jié)果與分析

3.1 立地質(zhì)量劃分的主導(dǎo)因子確定

田志和、馬常耕等［7-8］都對日本落葉松與生境因子間的關(guān)系進行了論述，考慮到湖北省建始縣長嶺崗林場的實際情況，在調(diào)查的43塊樣地中選取了全部的處于數(shù)量成熟林［9］階段的標(biāo)準(zhǔn)地10塊，用樹高生長量作為因變量生長指標(biāo)，用海拔、土層厚度、坡向、pH值、坡度和坡位作為自變量指標(biāo);其中海拔、土層厚度、坡度和pH值為定量因子，坡向和坡位為定性因子。用分級方法對定性因子進行數(shù)量化處理［6］，坡向陽坡為1、半陽坡為2、半陰坡為3、陰坡為4;坡位上部為1、中部為2、下部與山腳為3.5，然后在Matlab7.1中進行逐步回歸分析，最終確定4個因子作為立地質(zhì)量的評定主導(dǎo)因子，分別是坡向、坡度、海拔、土壤厚度。

3.2 數(shù)字高程模型的建立與地形特征值的提取

數(shù)字高程模型(Digital elevation models，DEM)是地形地貌的一個非常直觀的表示方法，其數(shù)據(jù)來源于數(shù)字地形圖。為了建立湖北省建始縣長嶺崗林場的DEM模型，將6幅1∶10 000的地形圖在Able軟件公司的Raster2Vector 5.x(R2V)中進行矢量化［10-11］，然后在ArcGIS9.1中構(gòu)建了如圖1所示的長嶺崗林場的數(shù)字高程模型［12-14］。

圖1 DEM圖

嚴(yán)格地講，地表面任一點的坡度是指過該點的切平面與水平面的夾角。坡度表示了地面在該點的傾斜程度。坡面姿態(tài)(坡度及坡向)是指局部地表面在空間的傾斜程度和朝向。坡度表示了該局部地表坡面的傾斜程度，坡度大小直接影響著地表物質(zhì)流動與能量轉(zhuǎn)換的規(guī)模與強度，是制約生產(chǎn)力空間布局的重要因子［12-13］。

在3×3的DEM柵格窗口中，如果中心柵格無數(shù)據(jù)，則此柵格的坡度值也是無數(shù)據(jù);如果相鄰的任何柵格是無數(shù)據(jù)，那么它們被賦予中心柵格的值再被用來計算坡度值［12］。坡度值的范圍是0～90°。根據(jù)上述原理可在DEM的基礎(chǔ)上直接提取坡度信息(圖2)。坡度(P)級按國家森林資源連續(xù)清查技術(shù)規(guī)定通用標(biāo)準(zhǔn)劃分為6級：Ⅰ級為平坡，P＜5°;Ⅱ級為緩坡，5°≤P＜15°;Ⅲ級為斜坡，15°≤P＜25°;Ⅳ級為陡坡，25°≤P＜35°;Ⅴ級為急坡，35°≤P＜45°;Ⅵ級為險坡，P≥45°。

坡向的提取詳見表1及圖3、圖4。

表1 陰陽坡向劃分

圖2 坡度圖

圖3 坡向圖

圖4 陰陽坡圖

3.3 土層厚度的克立格

日本落葉松雖然是淺根性樹種，但隨著土層厚度不同卻有一定的可塑性。根據(jù)調(diào)查，生長在土層深厚地段的10年生日本落葉松，根主要分布在30～40 cm的土層中，20～30年生時則可深達70～80 cm。日本落葉松在遼寧東部山區(qū)，在土層厚度大于50 cm的土壤中生長良好;在厚度大于75 cm的壤土中則能速生豐產(chǎn);而在厚度小于25 cm的土壤中生長緩慢，有的形成小老樹。

為了應(yīng)用克立格法，對采樣的土層厚度值分別用Arc-GIS9.1自帶的Normal QQplot方法［12］和總體分布的卡方(Chi-square)檢驗法，得出同一結(jié)論總體服從正態(tài)分布。

克立格插值法(Kriging)又稱空間局部插值法，是以變異函數(shù)理論和結(jié)構(gòu)分析為基礎(chǔ)，在有限區(qū)域內(nèi)對區(qū)域化變量進行無偏最優(yōu)估計的一種方法，是地統(tǒng)計學(xué)的主要內(nèi)容之一。南非礦產(chǎn)工程師D.R.Krige(1951年)在尋找金礦時首次運用這種方法，法國著名統(tǒng)計學(xué)家G.Matheron隨后將該方法理論化、系統(tǒng)化，并命名為Kriging，即克里格方法［15-16］?？死锔袷且环N估計觀測樣點間內(nèi)插值的地統(tǒng)計學(xué)方法，它基于區(qū)域變量理論，當(dāng)獲得了某個變量的變異函數(shù)的模擬模型后，就可以利用樣點觀測值對研究區(qū)域上未取樣點的變量值進行估計。如果區(qū)域化變量滿足二階平穩(wěn)或本征假設(shè)，對點或塊段的估計可直接采用點克立格法(Puctual kriging)或者塊段克立格法(Block kriging)。這兩種方法是最基本的估計方法，也稱普通克立格法(Origing kriging，簡稱 OK 法)［14-16］，用普通克立格法得土層的厚度分布圖(圖5)。

圖5 基于OK法的土層厚度分布圖(土層厚度單位：cm)

3.4 基于層次分析法的各因子權(quán)重的確定

層次分析法(Analytia1 hierarchy process，簡稱AHP)是美國匹茲堡大學(xué)教授A.L.Saaty于20世紀(jì)70年代提出的一種系統(tǒng)分析方法［17］。是一種能將定性分析與定量分析相結(jié)合的系統(tǒng)分析方法。是分析多目標(biāo)、多準(zhǔn)則的復(fù)雜大系統(tǒng)的有力工具。

3.4.1 建立層次結(jié)構(gòu)

在確定了影響日本落葉松生長的4個主導(dǎo)因子后，可進一步將每個因子細(xì)分：坡向分為陽坡、半陽坡、陰坡、半陰坡和無坡向，其代碼分別為 B41、B42、B43、B44、B45。坡度劃分為6 級(圖2)，根據(jù)所掌握的日本落葉松的生長特性［7］，將Ⅰ級至Ⅲ級劃分為平緩斜坡，代碼B31;Ⅳ級為陡坡，代碼B32;Ⅴ、Ⅵ級為急險坡，代碼B33。土層厚度(H)分為H＜25 cm層，代碼B21;25 cm≤H＜50 cm層，代碼B22;50 cm≤H＜75 cm層，代碼B23;H≤75 cm層，代碼B24。海拔1520～1600 m，代碼B11;海拔1600～1720 m，代碼B12;海拔1 720～1 810 m，代碼B13;海拔1 810～19 200 m，代碼B14?？山⑷鐖D6所示的層次結(jié)構(gòu)圖。

圖6 立地因子評價層次結(jié)構(gòu)圖

3.4.2 構(gòu)造判斷矩陣

建立了層次分析結(jié)構(gòu)模型后，明確了上下層之間的隸屬關(guān)系，就可轉(zhuǎn)化為層次中的排序計算問題。在排序計算中，每一層次中的排序，又可簡化為將下一層次各因素(或方案)相對于上一層次每一準(zhǔn)則(或目標(biāo))的重要性進行兩兩比較，得到各因素相對重要性權(quán)重.對于定量數(shù)據(jù)，可用總和標(biāo)準(zhǔn)化得到它們各自相對重要性權(quán)重;對于定性數(shù)據(jù)，要對各因素的相對重要性作出判斷，并按Saaty創(chuàng)立的1～9級比率標(biāo)度將判斷定量化［17］，構(gòu)造出判斷矩陣。

經(jīng)計算4個因子中海拔的權(quán)重為0.061 6，土層厚度的權(quán)重為0.173 7，坡度的權(quán)重為0.291 0，坡向的權(quán)重為 0.473 7;陰坡排第1位，Ⅰ級至Ⅲ級坡度排第2位，半陰坡排第3位，海拔1 520～1 600 m 排第16位［18］。

3.5 灰色定權(quán)聚類

灰色聚類可分為灰色關(guān)聯(lián)聚類(主要用于同類因素的歸并，以使復(fù)雜系統(tǒng)簡化)和灰類白化權(quán)函數(shù)聚類(主要用于檢查觀測對象是否屬于事先設(shè)定的不同類別)。其中灰類白化權(quán)函數(shù)聚類又包含灰色變權(quán)聚類和灰色定權(quán)聚類兩種方法，灰色變權(quán)聚類適用于聚類指標(biāo)的意義、量綱皆相同的情形;而當(dāng)聚類指標(biāo)的意義、量綱不同且不同指標(biāo)的樣本值在數(shù)量上相差較懸殊時，宜采用灰色定權(quán)聚類方法［19-21］。鑒于劃分日本落葉松人工林立地質(zhì)量4個因子的特點，本研究宜采用灰色定權(quán)聚類方法進行立地質(zhì)量研究。

3.5.1 灰色定權(quán)聚類的基本原理

設(shè)xij(i=1，2，…，n;j=1，2，…，m)為對象i關(guān)于指標(biāo)j的樣本值.)(j=1，2，…，m;k=1，2，…，s)為j指標(biāo)k子類白化權(quán)函數(shù)。若j指標(biāo)關(guān)于k子類的權(quán)j=1，2，…，m;k=1，2，…，s)與k無關(guān)，記為 ηj(j=1，2…，m)，稱

為對象i屬于k灰類的灰色定權(quán)聚類系數(shù)。

灰色定權(quán)聚類可按下列步驟進行：第一步給出j指標(biāo)k子類白化權(quán)函數(shù)(.)(j=1，2，…，m;k=1，2，…，s);第二步確定各指標(biāo)的聚類權(quán)，根據(jù)白化權(quán)函數(shù)(.)(j=1，2，…，m;k=1，2，…，s)、聚類權(quán) ηj(j=1，2，…，m)，以及對象i關(guān)于指標(biāo)j的樣本值xij(i=1，2，…，n;j=1，2，…，m)算出灰色定權(quán)聚類系數(shù)第三步，若，則斷定對象i屬于灰類k*。

3.5.2 白化權(quán)函數(shù)(.)的確定

在進行立地質(zhì)量研究時，將立地質(zhì)量劃分為好、中、差3類，即s=3，在n個樣本單元中，每個樣本單元選取4個聚類指標(biāo)，海拔、土層厚度、坡度、坡向，即m=4，可得出j指標(biāo)k子類的白化權(quán)函數(shù))(j=1，2，3，4;k=1，2，3)。

當(dāng)j=1，指標(biāo)為海拔時，［1 720，1 810，-，-］［1 520，1 600，1 720，1 810］［-，-，1520，1600］所對應(yīng)的函數(shù)分別為：

當(dāng)j=2，指標(biāo)為土層厚度時［25，75，-，- ］［25，50，- ，100］［-，-，25，75］所對應(yīng)的函數(shù)分別為：

當(dāng)j=3，指標(biāo)為坡度時［-，-，25，35［0，25，35，45］25，45，-，-］所對應(yīng)的函數(shù)分別為：

當(dāng)j=4，指標(biāo)為坡向時，所對應(yīng)的函數(shù)分別為：

由式(1)～式(13)計算出各因子的白化權(quán)函數(shù)值，綜合層次分析法計算出的各因子的權(quán)重，結(jié)合圖1～圖4，在Matlab軟件中計算出各柵格點屬于日本落葉松人工林立地質(zhì)量好、中、差的灰色定權(quán)聚類系數(shù)，按照其判別原則生成長嶺崗林場的日本落葉松人工林立地質(zhì)量圖(圖7)。

經(jīng)統(tǒng)計得到日本落葉松人工林立地質(zhì)量好的面積為354.60 hm2(占總面積的39.26%);中等的面積為429.75 hm2(占總面積的47.58%);差的面積為118.95 hm2(占總面積的13.17%)。

4 結(jié)論與討論

在確定4個主導(dǎo)因子的白化權(quán)函數(shù)時，坡度、海拔、土壤厚度均可找到典型的白化權(quán)表達函數(shù)，而坡向卻找不到。原因有二：一是因坡向值的提取是在ArcGIS9.1軟件中進行的，無坡向被賦值-1，其他值則為0°～360°;二是據(jù)已有的研究資料表明，日本落葉松雖是陽性樹種，但不耐上方庇蔭，在陰坡、半陰坡的生長比陽坡、半陽坡的要好［7-8］，考慮到坡向值取值的連續(xù)性和人們的習(xí)慣，找不到現(xiàn)有的典型白化權(quán)函數(shù)和以前相似的研究。

圖7 日本落葉松人工林立地質(zhì)量圖

在用339個土壤剖面數(shù)據(jù)對土層厚度估計時，因其采樣數(shù)據(jù)服從正態(tài)分布，故采用普通克立格法;若不服從正態(tài)分布，需進行數(shù)據(jù)變換，使其服從正態(tài)分布后方可使用普通克立格法，否則結(jié)果會出現(xiàn)較大的偏差。

雖然日本落葉松喜微酸性土壤，不耐鹽堿，但全場的pH值為5.3～6.5，均為微酸性土壤，因此該因子在全場內(nèi)對日本落葉松的立地質(zhì)量劃分影響不顯著;若林場內(nèi)有堿性土壤或強酸性土壤，該因子一定會被選中，參與立地質(zhì)量的評定。

用灰色定權(quán)聚類法雖已取得了日本落葉松人工林立地質(zhì)量的量化劃分，但在實際應(yīng)用中林業(yè)上常常以一個小班為經(jīng)營單位，用該法可能會產(chǎn)生在同一個小班內(nèi)，不同像元的立地質(zhì)量可能不一致，這時可以立地質(zhì)量類型所占的面積百分比最大者為該小班的立地質(zhì)量歸屬，即以最大者命名。

［1］張洪亮，李芝喜，王人潮，等.基于GIS的貝葉斯統(tǒng)計推理技術(shù)在印度野牛生境概率評價中的應(yīng)用［J］.遙感學(xué)報，2000，4(1)：66-70.

［2］鄧送求，閆家鋒，關(guān)慶偉，等.基于聚類分析的風(fēng)景林立地類型劃分［J］.南京林業(yè)大學(xué)學(xué)報：自然科學(xué)版，2009，33(3)：73-77.

［3］《中國森林立地分類》編寫組.中國森林立地分類［M］.北京：中國林業(yè)出版社，1989：23-97.

［4］孟憲宇.測樹學(xué)［M］.北京：中國林業(yè)出版社，1996：97-128.

［5］張宋智，張會軍，蔣明，等.小隴山林區(qū)日本落葉松速生豐產(chǎn)林立地類型劃分及評價的研究［J］.甘肅林業(yè)科技，1999，24(2)：1-6.

［6］荀守華，曹漢玉，吳德軍.日本落葉松生長量與立地因子的相關(guān)分析［J］.山東師大學(xué)報：自然科學(xué)版，1995，10(3)：264-266.

［7］田志和，董鍵，王喜武，等.日本落葉松育林學(xué)［M］.北京：北京農(nóng)業(yè)大學(xué)出版社，1995：1-67.

［8］馬常耕，王建華.我國發(fā)展日本落葉松區(qū)域的探討［J］.林業(yè)科技通訊，1990(4)：26-28.

［9］沈作奎，魯勝平.日本落葉松人工林撫育間伐起始期的研究［J］.湖北民族學(xué)院學(xué)報：自然科學(xué)版，2003，21(4)：44-46.

［10］韓晨，衛(wèi)海燕.基于R2V、ArcView的地形圖屏幕矢量化和三維顯示的方法［J］.山西師范大學(xué)學(xué)報：自然科學(xué)版，2006，20(1)：92-95.

［11］羅云啟，曾琨，羅毅.數(shù)字化地理信息系統(tǒng)建設(shè)與Maplnfo高級應(yīng)用［M］.北京：清華大學(xué)出版社，2003：164-179.

［12］李志林，朱慶.數(shù)字高程模型［M］.武漢：武漢測繪科技大學(xué)出版社，2005：16-78.

［13］湯國安，劉學(xué)軍，閭國年.數(shù)字高程模型及地學(xué)分析的原理與方法［M］.北京：科學(xué)出版社，2005：234-325.

［14］湯國安，楊昕.ArcGIS地理信息系統(tǒng)空間分析實驗教程［M］.北京：科學(xué)出版社，2006：363-455.

［15］候景儒，尹鎮(zhèn)南，李維明，等.實用地質(zhì)統(tǒng)計學(xué)［M］.北京：地質(zhì)出版社，1998：31-72.

［16］王政權(quán).地統(tǒng)計學(xué)及在生態(tài)學(xué)中的應(yīng)用［M］.北京：科學(xué)出版社，1999：1-102.

［17］朱建軍.層次分析法的若干問題研究及應(yīng)用［D］.東北大學(xué)機械工程與自動化學(xué)院，2005：1-69.

［18］馬友平.日本落葉松人工林林分結(jié)構(gòu)與生長量研究［D］.北京：北京林業(yè)大學(xué)資源與環(huán)境學(xué)院，2007：45-78.

［19］劉思峰，郭天榜，黨耀國.灰色系統(tǒng)理論及其應(yīng)用［M］.北京：科學(xué)出版社，1999：102-155.

［20］易德生，郭萍.灰色理論與方法：提要·題解·程序·應(yīng)用［M］.北京：石油工業(yè)出版社，1992：128-162.

［21］趙憲文.林業(yè)遙感定量估測［M］.北京：中國林業(yè)出版社，1997：97-113.