祝維 張西 賈黎明

(北京林業(yè)大學(xué)，北京，100083)

伏牛山地區(qū)栓皮櫟天然次生林地位指數(shù)表的編制1)

祝維 張西 賈黎明

(北京林業(yè)大學(xué)，北京，100083)

編制了伏牛山地區(qū)栓皮櫟天然次生林地位指數(shù)表，以期為當(dāng)?shù)厮ㄆ堤烊淮紊趾笃诘纳纸?jīng)營和天然林保護(hù)工作奠定基礎(chǔ)。從伏牛山南北坡88株平均優(yōu)勢木解析木數(shù)據(jù)獲得了619對樹高—年齡數(shù)據(jù)對，用SPSS軟件對優(yōu)勢木平均樹高與年齡的關(guān)系進(jìn)行擬合，從而獲得導(dǎo)向曲線方程。根據(jù)各齡階樹高變動系數(shù)(CH)及標(biāo)準(zhǔn)差(SH)的變化幅度確定基準(zhǔn)年齡，根據(jù)基準(zhǔn)年齡時(shí)的優(yōu)勢木樹高絕對變動幅度(ΔH)與指數(shù)級數(shù)量(n)的比值確定指數(shù)級距，采用標(biāo)準(zhǔn)差調(diào)整法導(dǎo)出地位指數(shù)曲線。從9個(gè)常用的數(shù)學(xué)模型中選定對數(shù)曲線式模型作為導(dǎo)向曲線；確定基準(zhǔn)年齡為45 a，指數(shù)級距為1 m；以15指數(shù)級樹高為基準(zhǔn)導(dǎo)出了9～18 m共10條地位指數(shù)曲線，最終通過上限排外法編制了伏牛山地區(qū)栓皮櫟天然次生林的地位指數(shù)表。經(jīng)χ2檢驗(yàn)和落點(diǎn)檢驗(yàn)后證明所編地位指數(shù)表精度較高，能夠客觀地評價(jià)伏牛山地區(qū)栓皮櫟天然次生林的立地質(zhì)量。

伏牛山地區(qū)；栓皮櫟；天然次生林；地位指數(shù)表；立地質(zhì)量評價(jià)

中國天然櫟類林面積廣闊，是我國最為重要的一類森林。栓皮櫟(Quercusvariabilis)是重要的櫟類樹種之一，具有極高的經(jīng)濟(jì)價(jià)值和生態(tài)效益。立地質(zhì)量直接影響林分生長和產(chǎn)量，因此立地質(zhì)量的確定是森林管理的重點(diǎn)[1]，地位指數(shù)表是森林經(jīng)營工作的基礎(chǔ)數(shù)表，用地位指數(shù)法對立地質(zhì)量進(jìn)行評價(jià)在我國應(yīng)用十分廣泛[2-3]。從1978年地位指數(shù)表引入中國以來，至今國內(nèi)已有62個(gè)樹種被編表，闊葉樹種有南方型楊樹(Populusdeltoides)[4]、栓皮櫟(Quercusvariabilis)[5-8]、巨桉(Eucalyptusgrandis)[9]等，針葉樹種有馬尾松(Pinusmassoniana)[10]、杉木(Cunninghamialanceolata)[11]、側(cè)柏(Platycladusorientalis)[12]、油松(Pinustabulaeformis)[13-14]等。在秦嶺、華北等地區(qū)已有關(guān)于栓皮櫟天然林的地位指數(shù)表[5-8]，雖伏牛山為秦嶺山脈向東最大的余脈，但它是中國天然栓皮櫟分布的地理中心之一。由于中國長期缺乏對栓皮櫟作用的正確認(rèn)識，導(dǎo)致伏牛山地區(qū)的栓皮櫟林至今仍沒有得到合理的撫育管理，結(jié)果林分內(nèi)密度過大，干形彎曲，材質(zhì)差。這樣的林分不僅損失巨大的經(jīng)濟(jì)價(jià)值，且對其潛在生態(tài)經(jīng)濟(jì)價(jià)值也影響重大。因此筆者使用在伏牛山南北坡獲得的619對樹高1 a齡數(shù)據(jù)對來編制地位指數(shù)表，利用該表可以科學(xué)、客觀地評價(jià)伏牛山地區(qū)栓皮櫟天然次生林的立地質(zhì)量，以期為后續(xù)的森林培育經(jīng)營方案提供基礎(chǔ)數(shù)表參考。

1 研究地概況

伏牛山位于河南省西部地區(qū)，山脈山勢自西北向東南逐漸變緩，海拔1 000 m左右，為秦嶺山脈向東最大的余脈。整個(gè)伏牛山區(qū)域約占河南省山地丘陵面積的40%[15]。山體主要由花崗巖組成，土壤隨著自然植被和氣候變化特點(diǎn)主要分為棕壤土和褐色土。由于伏牛山的地勢高差懸殊且正好處于南北分界線，導(dǎo)致其垂直氣候變化較為明顯且南北坡氣候也有所不同。南坡?lián)袊鴼夂騾^(qū)劃的指標(biāo)將其劃分為北亞熱帶，北坡可以將其劃分為暖溫帶[15]。伏牛山地區(qū)的植被以暖溫帶和亞熱帶天然闊葉林為主，同時(shí)兼有針闊混交林。栓皮櫟(Quercusvariabilis)是伏牛山地區(qū)海拔1 000 m以下的山地和丘陵的建群種，是該地區(qū)地帶性植被[16]。

2 研究方法

2.1 數(shù)據(jù)來源與立地類型劃分

2013年6月份將標(biāo)準(zhǔn)地設(shè)置在伏牛山北坡和南坡，北坡有洛陽市欒川縣的合峪鎮(zhèn)、廟子鎮(zhèn)，南坡有南陽市西峽縣米坪鎮(zhèn)、回車鎮(zhèn)和南召縣四棵樹鄉(xiāng)。標(biāo)準(zhǔn)地為長30 m、寬20 m的矩形，面積600 m2，共88塊標(biāo)準(zhǔn)地。對標(biāo)準(zhǔn)地內(nèi)的每一株達(dá)到起測徑階的樹木進(jìn)行每木檢尺，從5株優(yōu)勢木中選取1株平均優(yōu)勢木進(jìn)行樹干解析，求出其樹高、年齡及其各齡階樹高作為編表數(shù)據(jù)來源。

伏牛山區(qū)屬于北亞熱帶和暖溫帶的分界帶，因此其南北坡的氣候條件差異很大，有研究表明南坡極端低溫和均溫明顯高于北坡，而年降水量北坡則比南坡低14.41%[17]。因此本研究以南坡為陽坡，以北坡為陰坡作為劃分立地類型的依據(jù)之一；土壤厚度等土壤性狀因子是影響伏牛山南北坡喬木生長的重要因子[18]，所以土壤厚度將作為依據(jù)之二；由于栓皮櫟在伏牛山分布的垂直高度上存在差異，而海拔的不同水熱條件也會出現(xiàn)很大的變化，因此海拔作為依據(jù)之三。依據(jù)以上3個(gè)劃分立地類型的依據(jù)，我們將伏牛山栓皮櫟的立地類型劃分為5種。

本研究中設(shè)置的標(biāo)準(zhǔn)地覆蓋了伏牛山地區(qū)栓皮櫟次生林的所有齡組和有栓皮櫟純林分布的大部分立地類型(見表1、表2)。

表1 標(biāo)準(zhǔn)地年齡分布統(tǒng)計(jì)

2.2 編制天然林地位指數(shù)表的可行性分析

根據(jù)林分樹高評定立地質(zhì)量的方法，其前提為研究對象必須是同齡純林或者相對同齡純林。然而天然次生林的年齡無法通過查閱造林記錄等方法確定其是否為同齡林，因此需要根據(jù)其他方法判斷其年齡。據(jù)文獻(xiàn)[2]的介紹，同齡純林的直徑結(jié)構(gòu)近似遵從正態(tài)分布。通過對林分直徑分布進(jìn)行擬合，根據(jù)擬合曲線的特征可以判斷其是否為同齡純林。根據(jù)上述方法，本研究從88塊標(biāo)準(zhǔn)地中隨機(jī)選擇一塊平均直徑為17.5 cm的標(biāo)準(zhǔn)地，運(yùn)用統(tǒng)計(jì)軟件SPSS18.0對該塊標(biāo)準(zhǔn)地內(nèi)樣木作直徑分布頻次直方圖(見圖1)。可知中等大小徑階(直徑在10.0～25.0 cm區(qū)間內(nèi))的林木數(shù)量占多數(shù),向兩端徑階的林木數(shù)量逐漸減少,在林分算數(shù)平均值(17.5 cm)時(shí)達(dá)到峰值。根據(jù)這些特征可以判斷該分布曲線近似于正態(tài)分布，與同齡純林直徑分布規(guī)律相吻合。運(yùn)用相同方法對其他標(biāo)準(zhǔn)地內(nèi)樣木作直徑分布頻次直方圖，得出的結(jié)果基本一致，可以認(rèn)為伏牛山地區(qū)栓皮櫟天然次生林為相對同齡林。

表2 標(biāo)準(zhǔn)地立地類型統(tǒng)計(jì)

圖1 直徑分布頻次圖

2.3 異常數(shù)據(jù)剔除

使用EXCEL2007，SPSS18.0等統(tǒng)計(jì)軟件對調(diào)查數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。首先將所有從平均優(yōu)勢木解析木獲得的樹高1 a齡數(shù)據(jù)對，按照5 a的階距將5～70 a劃分為14個(gè)齡階，然后統(tǒng)計(jì)出各齡階平均優(yōu)勢木解析木樣本數(shù)量及其平均樹高、標(biāo)準(zhǔn)差，再以每齡階平均樹高為基準(zhǔn)，使用3倍標(biāo)準(zhǔn)差法對該齡階內(nèi)異常的優(yōu)勢木樹高1 a齡數(shù)據(jù)對進(jìn)行剔除，最終得到619組數(shù)據(jù)對。

表3 解析木特征統(tǒng)計(jì)

2.4 地位指數(shù)曲線方程的選擇與擬合

根據(jù)標(biāo)準(zhǔn)地整理資料，使用統(tǒng)計(jì)軟件SPSS18.0，利用9個(gè)曲線方程來擬合樹高曲線，各方程表達(dá)式及計(jì)算結(jié)果見表4。比較9個(gè)樹高曲線方程，結(jié)合決定系數(shù)(R2)和殘差平方和(Q)最終確定對數(shù)曲線式lnH=a+blnA為導(dǎo)向曲線擬合公式，即

lnH=0.217 8+0.658 5lnA。

(1)

式中：H為優(yōu)勢木樹高；A為優(yōu)勢木年齡。R2=0.912，且殘差平方和在9個(gè)曲線中最小，故選擇其為導(dǎo)向曲線方程。

2.5 基準(zhǔn)年齡(A0)和地位指數(shù)級距(D)的計(jì)算

以標(biāo)準(zhǔn)地的619對樹高1 a齡數(shù)據(jù)樣本為研究對象，通過公式：ΔSH=SH(i+1)/SH(i)和ΔCH=CH(i+1)/CH(i)分別計(jì)算各齡階樹高標(biāo)準(zhǔn)差的變化幅度(ΔSH)和變異系數(shù)的變化幅度(ΔCH)。當(dāng)它們的變化幅度均接近于1時(shí)，說明該齡階樹高生長趨于穩(wěn)定，這時(shí)的年齡可以確定為基準(zhǔn)年齡[5,13]。根據(jù)計(jì)算的ΔSH和ΔCH的值繪制折線圖(圖2)?？梢钥闯銮?0 a樹高標(biāo)準(zhǔn)差(SH)與變異系數(shù)(CH)的變化幅度一直較大，隨著林木年齡增長逐漸減小，并在45齡級之后趨向穩(wěn)定，且伏牛山地區(qū)栓皮櫟林的成熟年齡為40～50 a。綜上，本研究將采用45 a作為基準(zhǔn)年齡。

表4 曲線方程表達(dá)式及其計(jì)算結(jié)果統(tǒng)計(jì)

注：H為優(yōu)勢木各齡階平均樹高(m)；A為優(yōu)勢木各齡階平均年齡(a)；a、b、c、d為相關(guān)參數(shù)；R2為決定系數(shù)；Q為殘差平方和。

圖2 樹高標(biāo)準(zhǔn)差變化幅度與變異系數(shù)變化幅度

本研究中栓皮櫟在基準(zhǔn)年齡45 a時(shí)的優(yōu)勢木樹高范圍為8.8～17.9 m，即ΔH=9.1。通常情況下指數(shù)級距(D)為1～4 m，指數(shù)級的數(shù)量(n)在10個(gè)左右為宜[2]，本研究以n=10個(gè)為標(biāo)準(zhǔn)，通過公式D=ΔH/n確定D為1 m，最終得到9、10、11、12、13、14、15、16、17和18 m共10個(gè)指數(shù)級。

3 結(jié)果與分析

3.1 地位指數(shù)表的編制

以導(dǎo)向曲線為基礎(chǔ)，根據(jù)指數(shù)級距(D)和基準(zhǔn)年齡(A0)時(shí)樹高值，通過標(biāo)準(zhǔn)差調(diào)整法形成地位曲線簇[2]。將各齡階年齡(A)代入導(dǎo)向曲線方程(1)，可以計(jì)算得出各齡階導(dǎo)向曲線樹高(Hik)。根據(jù)各齡階樹高標(biāo)準(zhǔn)差(SH)與各齡階平均年齡(Ai)的值，擬合得到方程表達(dá)式

lnSH=-1.160+1.157lnA。

(2)

再將各齡階年齡A代入方程(2)中，可以計(jì)算得到各齡階的樹高標(biāo)準(zhǔn)差理論值(SAi)。本研究中，基準(zhǔn)年齡45 a時(shí)，導(dǎo)向曲線的理論樹高值(15.25 m)不等于地位指數(shù)級數(shù)值(15 m)，因此要進(jìn)行相應(yīng)的調(diào)整。首先計(jì)算基準(zhǔn)年齡時(shí)導(dǎo)向曲線樹高值與基準(zhǔn)年齡時(shí)第15指數(shù)級樹高值的差值d=Hoj-Hok=0.25m，然后計(jì)算調(diào)整系數(shù)Kj=((Hoj-Hok)/SA0)≈0.077，最后將調(diào)整系數(shù)(Kj)與各齡階樹高標(biāo)準(zhǔn)差理論值(SAi)相乘，即得到各齡階調(diào)整值(Kj·SAi)。將各齡階導(dǎo)向曲線樹高值(Hik)與各齡階調(diào)整值相減即可得到以15指數(shù)級為基礎(chǔ)各齡階調(diào)整后的樹高值(Ho)。具體計(jì)算結(jié)果見表5。

表5 各齡階的樹高調(diào)整值統(tǒng)計(jì)

注：Hik為第i齡階的導(dǎo)向曲線樹高(m)；SAi為第i齡階樹高標(biāo)準(zhǔn)差理論值；kj為調(diào)整系數(shù)；Ho為15指數(shù)級各齡階樹高值(m)。

以15指數(shù)級的導(dǎo)向曲線為準(zhǔn)，按照指數(shù)級距1 m逐齡階導(dǎo)算出各地位指數(shù)曲線上的樹高值，最終通過上限排外法編制出伏牛山地區(qū)的栓皮櫟林地位指數(shù)表(表6)，根據(jù)各齡階各地位指數(shù)級的樹高值繪制的地位指數(shù)曲線簇見圖3。

3.2 地位指數(shù)表的適用性檢驗(yàn)

必須對編制好的地位指數(shù)表進(jìn)行檢驗(yàn)，在保證其精度和適用性后，方能投入到實(shí)際工作中使用。本研究將以常用的方法落點(diǎn)檢驗(yàn)和擬合顯著性檢驗(yàn)對所編制的伏牛山地區(qū)栓皮櫟天然次生林進(jìn)行檢驗(yàn)。

3.2.1 落點(diǎn)檢驗(yàn)

在88塊有解析木的標(biāo)準(zhǔn)地中，統(tǒng)計(jì)出每一塊樣地的優(yōu)勢木樹高均值，將得到的88株優(yōu)勢木高均值繪在地位指數(shù)曲線簇上，有84株優(yōu)勢木樹高落在所編的地位指數(shù)曲線內(nèi)(如圖4)，精度達(dá)到95.45%。結(jié)果表明，所編的地位指數(shù)表具有良好地反映伏牛山地區(qū)栓皮櫟天然次生林立地質(zhì)量的能力。

表6 伏牛山地區(qū)栓皮櫟天然林地位指數(shù)

圖3 伏牛山地區(qū)栓皮櫟天然林地位指數(shù)曲線簇

3.2.2 擬合顯著性檢驗(yàn)

為了判斷解析木樹高生長趨勢與地位指數(shù)曲線趨勢的差異是否顯著，需要進(jìn)行擬合顯著性檢驗(yàn)。從地位指數(shù)級j=12、15、17中各隨機(jī)抽取1株50年生以上的解析木，將處于3個(gè)不同地位指數(shù)級(j=12、15、17)的3株解析木的10個(gè)齡階(i=5、10、15、…、50)的樹高值按公式Hoj=Hok-SA0(Hik-Hij)/SAi轉(zhuǎn)化為基準(zhǔn)年齡(45 a)時(shí)第j指數(shù)級的樹高，即地位指數(shù)。式中：Hoj為基準(zhǔn)年齡時(shí)第j指數(shù)級的樹高值(即地位指數(shù))；Hok為基準(zhǔn)年齡時(shí)導(dǎo)向曲線樹高；Hik為第i齡階導(dǎo)向曲線樹高；Hij為第i齡階第j指數(shù)級的樹高；SA0為基準(zhǔn)年齡所在齡階樹高標(biāo)準(zhǔn)差理論值；SAi為第i齡階樹高標(biāo)準(zhǔn)差理論值。本研究中SA0=3.25，Hok=15.25。抽取的3個(gè)指數(shù)級的3株解析木的各齡階樹高值Hij和基準(zhǔn)年齡時(shí)第j指數(shù)級樹高值Hoj的具體計(jì)算結(jié)果見表7。

圖4 伏牛山地區(qū)栓皮櫟天然林地位指數(shù)表的落點(diǎn)檢驗(yàn)

表7 伏牛山地區(qū)栓皮櫟天然林地位指數(shù)表的顯著性檢驗(yàn)

注：Hij為第i齡階第j指數(shù)級的樹高(m)；Hoj為基準(zhǔn)年齡時(shí)第j指數(shù)級的樹高值(即地位指數(shù))(m)。

4 結(jié)論與討論

本研究利用從88塊伏牛山地區(qū)栓皮櫟天然次生林標(biāo)準(zhǔn)地的調(diào)查資料和88株平均優(yōu)勢木解析得到619對樹高1 a齡數(shù)據(jù)對研制伏牛山地區(qū)栓皮櫟天然次生林地位指數(shù)表。經(jīng)過統(tǒng)計(jì)分析，從9個(gè)常用的數(shù)學(xué)模型中選定lnH=0.217 8+0.658 5lnA作為導(dǎo)向曲線，選定基準(zhǔn)年齡為45年，指數(shù)級距為1 m。運(yùn)用標(biāo)準(zhǔn)差調(diào)整法，導(dǎo)出了9～18 m共10條地位指數(shù)曲線，最終編制了伏牛山地區(qū)栓皮櫟天然次生林的地位指數(shù)表。運(yùn)用χ2檢驗(yàn)和落點(diǎn)檢驗(yàn)對該地位指數(shù)表進(jìn)行檢驗(yàn)，結(jié)果表明所編地位指數(shù)表精度較高，可用來對伏牛山地區(qū)栓皮櫟天然次生林的立地質(zhì)量進(jìn)行科學(xué)的評價(jià)。

在海拔較低或人為干擾嚴(yán)重的地區(qū)，形成了部分多代萌生林，萌生林與實(shí)生林的生長過程存在差異，為提高所編制的地位指數(shù)表對立地質(zhì)量反映的準(zhǔn)確性，本研究所布置的樣地均在實(shí)生林內(nèi)。由于伏牛山地區(qū)的萌生林面積很小，因此本研究所編制的地位指數(shù)表能用來評價(jià)伏牛山地區(qū)的大部分栓皮櫟林立地質(zhì)量。

利用樣地平均優(yōu)勢木解析數(shù)據(jù)編制地位指數(shù)曲線，能準(zhǔn)確描述林分在各年齡階段的優(yōu)勢高生長情況和反映樹高生長與林地生產(chǎn)力之間的實(shí)際關(guān)系，因?yàn)榻馕瞿舅跇拥氐牡匚恢笖?shù)是實(shí)測值[5]。

將本研究與華北地區(qū)、秦嶺地區(qū)地位指數(shù)表進(jìn)行對比[5-6]，可以發(fā)現(xiàn)3個(gè)地區(qū)栓皮櫟林立地質(zhì)量是：秦嶺地區(qū)最好，伏牛山地區(qū)次之，華北地區(qū)相對最差；其中秦嶺與伏牛山地區(qū)的差距較小，說明栓皮櫟與秦嶺的立地條件有更高的相適度。

該地區(qū)的成過熟林相對較少，且伏牛山地區(qū)的成過熟林受到當(dāng)?shù)叵嚓P(guān)部門的保護(hù)，因此該年齡層能夠獲得的解析木數(shù)據(jù)較少，數(shù)據(jù)多集中在中幼齡林，為完善伏牛山地區(qū)栓皮櫟天然次生林地位指數(shù)表，在后續(xù)的工作中還應(yīng)再加入成過熟林的解析木數(shù)據(jù)。

在全球氣候變化的大背景下，伏牛山地區(qū)近年多次發(fā)生嚴(yán)重干旱等等各種因素均有可能對栓皮櫟樹高生長造成影響。國外已有加入氣象的相關(guān)因素的文獻(xiàn)[19-20]，后期研究可加入相關(guān)因素，以對伏牛山地區(qū)立地質(zhì)量進(jìn)行更科學(xué)的評價(jià)。

在伏牛山地區(qū)，人為干擾是影響林分生長的重要且可控因素，因此，必須要采取合理的森林撫育措施，對森林進(jìn)行適當(dāng)?shù)娜斯じ深A(yù)，讓森林按照人們預(yù)想的方向發(fā)展，使林地發(fā)揮最大的生產(chǎn)潛力。在后期的森林經(jīng)營工作中，建議根據(jù)近自然森林經(jīng)營模式對伏牛山地區(qū)的栓皮櫟天然次生林進(jìn)行撫育改造[21]。

[2] 孟憲宇.測樹學(xué)[M].3版.北京:中國林業(yè)出版社，2006.

[3] 詹昭寧.森林生產(chǎn)力的評定方法[M].北京:中國林業(yè)出版社，1982.

[4] 肖君，方升佐，徐錫增.南方型楊樹人工林立地指數(shù)表的編制[J].福建農(nóng)林大學(xué)學(xué)報(bào)(自然科學(xué)版)，2002，20(4)：604-609.

[5] 鄭聰慧，賈黎明，段劼，等.華北地區(qū)栓皮櫟天然次生林地位指數(shù)表的編制[J].林業(yè)科學(xué)，2013，49(2)：79-85.

[6] 張瑜，賈黎明，鄭聰慧,等.秦嶺地區(qū)栓皮櫟天然次生林地位指數(shù)表的編制[J].林業(yè)科學(xué)，2014，50(4)：47-54.

[7] 侯箕，門海元.中條山天然栓皮櫟林地位指數(shù)表的編制[J].山西林業(yè)科技，1993(2)：1-5.

[8] 馬增旺，畢君，劉德章，等.太行山天然次生林主要樹種地位指數(shù)表的編制[J].山西林業(yè)科技，1995(3)：13-16.

[9] 黃從德，胡庭興，賴家明.四川巨桉短周期工業(yè)原料人工林地位指數(shù)表的編制[J].四川農(nóng)業(yè)大學(xué)學(xué)報(bào)，2006，35(6)：347-350.

[10] 遲健，李桂英，陳家明,等.浙江省馬尾松人工林多形地位指數(shù)表及林分生長過程表的編制[J].林業(yè)科學(xué)研究，1996，9(1):68-74.

[11] 沈俊強(qiáng)，蘇上真，鄭孫登.福建省杉木立地指數(shù)表編制[J].林業(yè)勘察設(shè)計(jì)，1996(1)：46-51.

[12] 馬豐豐，賈黎明，段劼，等.北京山區(qū)側(cè)柏人工林立地指數(shù)表的編制[J].北京林業(yè)大學(xué)學(xué)報(bào)，2008，30(6)：78-82.

[13] 李佩萍，武建林.晉中東部山區(qū)油松人工林地位指數(shù)表的編制[J].山西林業(yè)科技，1999(4)：19-27.

[14] 段劼，馬履一，賈黎明，等.北京低山地區(qū)油松人工林立地指數(shù)表的編制及應(yīng)用[J].林業(yè)科學(xué)，2009，45(3)：7-12.

[15] 宋朝樞.伏牛山自然保護(hù)區(qū)科學(xué)考察集[M].北京：中國林業(yè)出版社，1994.

[16] 羅偉祥.中國栓皮櫟[M].北京：中國林業(yè)出版社，2009.

[17] 陳志超，郝成元，劉昌華，等.伏牛山阻隔作用下水熱因子、土地覆被與NPP的時(shí)空特征[J].生態(tài)環(huán)境學(xué)報(bào)，2013，22(5)：761-766.

[18] 郝成元，周見.伏牛山南、北坡植物物種多樣性及其與主要生態(tài)因子的相關(guān)性[J].植物資源與環(huán)境學(xué)報(bào)，2013，22(3)：38-44.

[19] SABATIA C O， BURKHART H E. Predicting site index of plantation loblolly pine from biophysical variables[J]. Forest Ecology and Management，2014，326：142-156.

[21] 陸元昌，雷相東，洪玲霞，等.近自然森林經(jīng)理計(jì)劃技術(shù)體系研究[J].西南林學(xué)院學(xué)報(bào)，2010，30(1)：1-5.

EstablishmentofSiteIndexTableforQuercusvariabilisNaturalSecondaryForestinFuniushanMountainousArea

Zhu Wei, Zhang Xi, Jia Liming

(Beijing Forestry University, Beijing 100083, P. R. China)Journal of Northeast Forestry University，2017,45(12)：32-37.

Funiu mountainous area;Quercusvariabilis; Natural secondary forest; Site index table; Stand site evaluation

1)北京林業(yè)大學(xué)青年教師科學(xué)研究中長期項(xiàng)目(2015ZCQLX02)；北京市教育委員會科學(xué)研究與科研基地建設(shè)項(xiàng)目(SYSBL2009)。

祝維，女，1994年12月生，北京林業(yè)大學(xué)林學(xué)院，碩士研究生。E-mail：18811375082@163.com。

賈黎明，北京林業(yè)大學(xué)林學(xué)院，教授。 E-mail：jlm@bjfu.edu.cn。

2017年4月7日。

戴芳天。

S757.1

In order to lay the foundation for the later forest management and natural forest protection program, we established the site index table forQuercusvariabilisnatural secondary forest in the Funiu Mountain area. In this study, 619 height-age data from 88 analytic and average dominant individuals which grow in the southern and northern slopes of the Funiu Mountain were used to construct oriented curve with the software of SPSS. The datum age was ensured according to the varied range of variation coefficient and standard deviation of tree height in various ages. The exponential distance was ensured according to the ratio of absolute the varied range of dominant tree height at the datum age to the number of exponential. The site index curves were developed by using the method of standard deviation adjustment. Logarithmic curve model was selected as oriented curve from nine mathematical models. The datum age was 45 a, and the exponential distance was 1 m. Ten site index curves (9-18 m) were developed on the base of the tree height in the fifteenth site index. Eventually, the site index table forQ.variabilisnatural secondary forest in the Funiu Mountain area was established. The table was demonstrated to be of high accuracy tested by theχ2test and falling point test. This site index table is accurate in evaluating the site quality ofQ.variabilisnatural secondary forest in the Funiu Mountain area.