楊惠濱, 國慶喜

東北林業(yè)大學(xué)生態(tài)研究中心, 哈爾濱 150040

?

地形與競爭因子對紅松胸徑與年齡關(guān)系的影響

楊惠濱, 國慶喜*

東北林業(yè)大學(xué)生態(tài)研究中心, 哈爾濱 150040

森林年齡結(jié)構(gòu)一直是森林生態(tài)學(xué)研究的重要內(nèi)容之一,利用胸徑指標(biāo)間接推測樹木年齡是獲得森林年齡結(jié)構(gòu)的有效手段?；诤邶埥⌒∨d安嶺林區(qū)156株成熟紅松解析木數(shù)據(jù),分析成熟紅松個體胸徑與年齡之間關(guān)系,探討地形與競爭因子對紅松胸徑與年齡關(guān)系的影響。研究結(jié)果表明：對所搜集的紅松解析木數(shù)據(jù),海拔、坡度、坡向、競爭因子等因素都對紅松胸徑與年齡關(guān)系產(chǎn)生影響,其中坡向與坡位對紅松胸徑-年齡關(guān)系影響最為明顯,而競爭因子的影響較小并且在不同地形條件下表現(xiàn)出差異性。在南坡向-中坡位的紅松個體,其胸徑與年齡之間表現(xiàn)出了最佳的擬合效果。對地形與競爭因子的影響機(jī)理還有待于進(jìn)一步研究。

胸徑； 年齡； 競爭； 地形因子； 紅松

林分的年齡結(jié)構(gòu)一直是森林生態(tài)學(xué)研究的重要內(nèi)容之一。它是林木更新過程長短和更新速度快慢的反映[1- 2]。紅松(Pinuskoraiensis)作為東北東部地帶性頂極植被-闊葉紅松林的建群種,了解其年齡結(jié)構(gòu)對于研究其在氣候變化下生物群落演變、形成與維持機(jī)制,制訂合理的經(jīng)營方案,以及提高林分生產(chǎn)力等,都有一定的實際意義。但是由于空間分布復(fù)雜和發(fā)育階段各異等的影響,天然林中樹木年齡的測定難度很大[3],并容易對其造成傷害,因此對年齡的確定一般采用間接的推測方法。許多研究表明林木的直徑生長與年齡呈正相關(guān),并可建立一定的生長模型[4- 6]。通過建立樹齡與胸徑的回歸關(guān)系,能夠簡化測定方法,大大提高工作效率。在20世紀(jì)80年代,李克志[7]、孫玉軍[8]等利用伊春林區(qū)紅松解析木數(shù)據(jù),對紅松的胸徑、樹高生長過程以及各生長階段狀態(tài)作了系統(tǒng)的研究。同一時期,邵國凡[9]、蔣伊尹[10]則利用人工紅松林解析木數(shù)據(jù)建立了胸徑等生長的數(shù)學(xué)模型。近年來關(guān)于紅松胸徑與年齡關(guān)系的研究集中在長白山地區(qū)[1,3],并且考慮到了競爭因子在紅松胸徑與年齡關(guān)系中的影響[11- 13]。目前對林木胸徑和年齡的相互關(guān)系研究,均從單因子出發(fā),建立簡單的一元回歸統(tǒng)計模型。但由于數(shù)據(jù)難以獲得,迄今為止在這方面開展的研究很少,多數(shù)集中于紅松小徑木。從以往的研究和野外調(diào)查實踐都可以看出,相同徑級的林木,其年齡存在一定的差異,究其原因,除了本身的生物學(xué)特性,還有林木周圍的生長環(huán)境、周圍林木的相對空間位置以及相互之間的距離等因素[3]。但對這些影響因素的了解仍十分模糊,更缺乏影響因子的定量研究。有鑒于此,本文通過搜集黑龍江小興安嶺的樹齡在100a以上的成熟紅松解析木數(shù)據(jù)進(jìn)行分析,利用解析木數(shù)據(jù)中紅松生長的立地條件以及周圍鄰近木的信息,探討競爭、地形等有關(guān)因子對胸徑與年齡之間關(guān)系的影響,為闊葉紅松林年齡結(jié)構(gòu)的研究提供理論參考。

1 研究區(qū)概況

小興安嶺位于黑龍江省東北部,地理位置大約在46°28′—49°21′N, 127°42′—130°14′E之間。該區(qū)域?qū)俦睖貛Т箨懠撅L(fēng)氣候區(qū),四季分明,最冷月1月平均氣溫-20—-25℃,最低溫度可達(dá)-40℃,最熱月7月平均氣溫19—21℃左右,最高溫度可達(dá)31℃,極端最高氣溫為35℃,全年平均氣溫-1—1℃；年降水量為500—700 mm；土壤類型為典型暗棕壤。小興安嶺在中國植被區(qū)劃上屬于溫帶針闊混交林區(qū)域內(nèi)的北部地區(qū),原生地帶性植被是紅松為主的溫帶闊葉紅松林,是我國主要林區(qū)之一,擁有豐富的珍貴用材樹種,木材蓄積量大。木本植物種類多樣,喬木樹種主要包括紅松 (Pinuskoraiensis)、紫椴(Tiliaamurensis)、楓樺(Betulacostata)、白樺(Betulaplatyphylla)、臭冷杉(Abiesnephrolepis)、興安落葉松(Larixgmelini)、蒙古櫟(Quercusmongolica)、裂葉榆(Ulmuslaciniata)、紅皮云杉(Piceakoraiensis)、魚鱗云杉(Piceajezoensis)、水曲柳(Fraxinusmandshurica)、黃菠蘿(Phellodendronamurense)、胡桃楸(Juglansmandshurica)、色木槭(Acermono)和大青楊(Populusussuriensis)等,灌木主要包括暴馬丁香(Syringareticulata)、毛榛子(Corylusmandshurica)、瘤枝衛(wèi)矛(Euonymuspauciflorus)、和東北鼠李(Rhamnusdavurica)等。

2 數(shù)據(jù)來源及處理

數(shù)據(jù)來自黑龍江省小興安嶺林區(qū)各地紅松解析木數(shù)據(jù)資料且全部為林分平均木,搜集的地區(qū)其南北緯度范圍在46°56′—48°06′N之間,北自豐林國家級自然保護(hù)區(qū),南至朗鄉(xiāng)、鐵力林業(yè)局其主要分布如下：五營、紅星、湯旺河、朗鄉(xiāng)、美溪、翠巒、南岔、鐵力、新清、烏伊嶺等林業(yè)局。整理每株紅松解析木的胸徑、年齡、生長狀況、周圍鄰近木的胸徑以及與解析木的距離,海拔、坡度、坡向、坡位等地形因子的數(shù)據(jù)信息。所有的解析木都為天然林紅松針闊葉混交林,個體海拔主要分布在300—700 m之間共計156株,周圍鄰近木共計566株,其中紅松占總株數(shù)的44.9%,云冷杉占29.8%,其他闊葉樹占25.3%。解析木與周圍鄰近木生長狀況及解析木胸徑與年齡分布散點(diǎn)圖分別見表1和圖1。

表1 解析木與周圍鄰近木狀況一覽表

圖1 解析木胸徑與年齡關(guān)系散點(diǎn)圖 Fig.1 Scatter diagram of the relatoinship between the DBH and age of analysis wood

2.1 簡單競爭指數(shù)

林木間的競爭對林木生長的影響是很難量化的[14],學(xué)者們通常采用競爭指標(biāo)來度量林木所遭受的競爭強(qiáng)弱,競爭指標(biāo)可分為與距離有關(guān)的競爭指標(biāo)和與距離無關(guān)的競爭指標(biāo)[15],一般認(rèn)為與距離有關(guān)的競爭指標(biāo)能更好的反應(yīng)競爭關(guān)系。為探討每株解析木受到的競爭影響,采用了Hegyi簡單競爭指數(shù)[16],它包括了樹木大小、距離和空間分布等多方面的信息,并且簡單直觀,測算容易。Hegyi簡單競爭指數(shù)是我國研究林木種群競爭關(guān)系最多的與距離有關(guān)的競爭指數(shù)。表達(dá)式為[17]:

(1)

式中,CI為對象木的簡單競爭指數(shù),對象木指計算競爭指標(biāo)時所針對的樹木；D為對象木的胸徑；Dj為競爭木i的胸徑；DISTi為對象木與競爭木i的距離；n為競爭木的數(shù)量。

關(guān)于競爭木的距離范圍界定,國內(nèi)有采用處于對象木5、6 m或8 m鄰體半徑范圍內(nèi)的林木作為競爭木[18- 19],也有采用與對象木最近的4 株競爭木計算競爭指數(shù)[20]。考慮到競爭木樹冠半徑一般為2—4 m,因此為了包含更多的競爭信息,本文將鄰體半徑設(shè)置為8 m,以盡量選取更多競爭木。參考有關(guān)研究報道[12],如果鄰體半徑內(nèi)林木胸徑小于對象木胸徑,則不作為競爭木。

2.2 大樹斷面積競爭指標(biāo)的改進(jìn)

大樹斷面積(BAL)為大于對象木胸徑的所有林木的斷面積之和[11,21]。目前其在國際上應(yīng)用廣泛,是最好的競爭指標(biāo)之一[13,22- 23],國內(nèi)則研究較少,是一種與距離無關(guān)的競爭指標(biāo)。Busing認(rèn)為,某株林木周圍鄰體對其遮蔽和擁擠效應(yīng)的大小,均與它們之間距離的遠(yuǎn)近有關(guān),距離越遠(yuǎn)的鄰體對這株林木生長的影響越??；反之則影響越大[24]。因此采用林分大樹斷面積作為競爭指標(biāo)時,可以只計算對象木周圍鄰體樹木的大樹斷面積。有的研究根據(jù)樹木周圍鄰體半徑內(nèi)其他林木的生物量來計算該樹木所受到的競爭[25]。本文據(jù)此對BAL加以改進(jìn),采用鄰體半徑內(nèi)所有大于對象木胸徑的林木胸高斷面積的距離加權(quán)之和作為對象木所受到的競爭壓力WBAL。

(2)

PAi=[T′HEi-sin(THEi)]/π

(3)

圖2 中心樹木與其鄰體間距離與夾角THE之間的關(guān)系 Fig.2 The relationship between the distance and intersection angle of the center tree and its neighbous夾角THE=∠AOB;OA、OE、OB等于鄰體半徑RADIUS;OC、OC′為鄰體與中心樹木之間的距離,這一角度可以由下式計算,式中DISTi為一樹木與其鄰體之間的距離；RADUS為鄰體半徑

式中,Ai為鄰體半徑內(nèi)第i株大于對象木胸徑的林木個體的胸高斷面積；PAi為對象木與其第i株鄰體之間的距離加權(quán)系數(shù)取值范圍為[0,1]。當(dāng)鄰體與對象木距離等于0時,PA=1,即鄰體對于對象木的遮蔽和擁擠效應(yīng)達(dá)到最大；當(dāng)鄰體距離等于鄰體半徑時,PA=0,鄰體對對象木生長無影響；加權(quán)系數(shù)PA的計算公式是根據(jù)圖2中角度與距離的關(guān)系所提出,從圖2中可以看出,鄰體C和C′與中心樹木O的距離越近,其夾角AOB越大,距離越遠(yuǎn),夾角越小。當(dāng)鄰體距離等于鄰體半徑時,夾角等于0；當(dāng)鄰體與中心樹木在一起時,夾角等于180度按弧度制則等于π。THEi為一樹木與其第i個鄰體之間的夾角,單位為弧度。

THEi=2×arccos(DISTi/RADIUS)

目前,國內(nèi)外常用的競爭指標(biāo)包括與距離有關(guān)的競爭指標(biāo)和與距離無關(guān)的競爭指標(biāo)兩種,每個競爭指標(biāo)都有其各自的優(yōu)勢,大樹斷面積(BAL)競爭指標(biāo)的改進(jìn)使兩類計算方式相結(jié)合,能更加有效的描述對象木所受的競爭壓力。

2.3 回歸方程的建立

應(yīng)用SPSS 19.0軟件,以林木年齡為因變量,以胸徑、競爭指標(biāo)CI或WBAL、地形等因子為自變量建立回歸方程。采用逐步回歸分析法選取方程中所需的自變量。以F概率作為引入和剔除變量的判據(jù),當(dāng)一個變量的雙側(cè)檢驗的顯著性水平≤0.05 時,該變量被引入方程；當(dāng)顯著性水平≥0.10時,該變量被剔除。為消除共線性,方差膨脹因子(VIF)大于5 的自變量也被排除在外[26]。

3 結(jié)果與分析

3.1 紅松胸徑與年齡的關(guān)系

胸徑隨年齡的生長一般呈S型曲線[9],林木一般在幼年時期胸徑生長較慢,而在中期,尤其是在150a以后生長量大大提高[5],直到生長后期又將呈現(xiàn)出生長緩慢趨勢。反過來說,本文利用胸徑信息估算紅松年齡大小,考慮到生物學(xué)意義,其總的擬合曲線應(yīng)呈現(xiàn)出快—慢—快的變化趨勢。又考慮到本文所選取的數(shù)據(jù)均為年齡在100a以上的近成熟紅松,變化趨勢應(yīng)表現(xiàn)為先慢后快,所以用指數(shù)模型進(jìn)行擬合更符合其生物學(xué)意義。并采用線性回歸模型作為對比分析,紅松個體胸徑與年齡的相關(guān)性分析結(jié)果表明Pearson相關(guān)系數(shù)為0.562,表明對于所搜集的所有解析木,二者之間線性相關(guān)程度并不高。隨著年齡增加,紅松的胸徑隨之增大,二者必然存在著緊密的相關(guān)性。然而由于各種影響因素的存在,如競爭、地形、氣候、甚至紅松自身的生長節(jié)律等都會降低胸徑與年齡之間的相關(guān)程度。

采用指數(shù)、線性兩類模型,以年齡為因變量,胸徑為自變量進(jìn)行擬合(表2)。從擬合結(jié)果看出,在胸徑與年齡變量均取對數(shù)的前提下,指數(shù)模型決定系數(shù)達(dá)到0.507。即對于所搜集的解析木資料來說,胸徑變量通過指數(shù)模型對年齡的解釋程度可以達(dá)到50.7%。而蔣伊尹對黑龍江省東部林區(qū)紅松人工林胸徑生長所做的分析表明,胸徑生長曲線用冪函數(shù)模型擬合比較理想,模型的相關(guān)系數(shù)基本都在0.99以上[10]。胡云云等[3]對吉林省汪清林業(yè)局紅松所做的胸徑生長方程,其冪函數(shù)模型決定系數(shù)也達(dá)到了0.97以上,但所涉及的紅松個體均為小徑木。盡管上述報道是以胸徑為因變量,以年齡為自變量建立的回歸方程,但仍然說明對于本研究所涉及的天然林紅松個體,其胸徑與年齡關(guān)系更加復(fù)雜,受到的影響也更多。

表2 紅松解析木年齡與胸徑關(guān)系方程

A: 年齡, D: 胸徑

那么是否由于所搜集的解析木數(shù)據(jù)分布空間較大而影響了胸徑年齡回歸方程的擬合效果？比如研究地區(qū)南北相差約1個緯度,年均溫相差1.7℃,可能會對紅松胸徑與年齡關(guān)系產(chǎn)生影響。我們挑選包含解析木數(shù)量最多(共計63個)的伊春地區(qū)五營林業(yè)局,對其解析木進(jìn)行胸徑與年齡(均取對數(shù))關(guān)系的曲線擬合,其結(jié)果以指數(shù)模型擬合的決定系數(shù)只有0.402。說明即使消除了緯度變化因子,紅松胸徑-年齡之間的擬合效果也并沒有得到改善。

3.2 競爭因子對紅松胸徑與年齡關(guān)系的影響

林木在生長過程中會受到種內(nèi)或者種間競爭的影響。趙俊卉等人在對長白山主要針葉樹種胸徑和樹高隨年齡的變異系數(shù)與競爭因子的關(guān)系研究表明,競爭是引起胸徑變動的主要因子,推斷如果將競爭因子添加到傳統(tǒng)的生長方程將會大幅度提高模型的精度[13]。因此本研究分別將CI、WBAL等競爭指標(biāo)加入胸徑-年齡回歸方程,探討競爭因子對紅松胸徑與年齡關(guān)系的影響。以年齡為因變量,以胸徑、競爭因子為自變量,采用多元逐步回歸模型進(jìn)行年齡的預(yù)測(表3)。

表3 模型參數(shù)與擬合統(tǒng)計量

從表3可以看出,不同的競爭指標(biāo)對紅松胸徑-年齡回歸關(guān)系的影響是不同的。引入簡單競爭指數(shù)CI之后,回歸模型的修正決定系數(shù)增加為0.542；而引入WBAL指標(biāo)的結(jié)果則使胸徑-年齡擬合效果有所改善達(dá)到0.583。但是這種改善效果也是很有限的,沒有大幅度提高擬合精度,至少對于所搜集的解析木數(shù)據(jù)是這樣。因此進(jìn)一步探討其他因子的影響是十分必要的。

3.3 地形因子對紅松胸徑與年齡關(guān)系的影響

在前面所做的多元線性回歸模型基礎(chǔ)上,進(jìn)一步加入坡度、海拔等地形因子。結(jié)果表明坡度與海拔因子均被剔除,說明對于所搜集的全部解析木,坡度與海拔因子沒有顯示出影響。本文擬將坡向與坡位因子分別加以限定,繼續(xù)開展相關(guān)探討。

3.3.1 坡向?qū)t松胸徑與年齡關(guān)系的影響

參照有關(guān)文獻(xiàn)[27],將方位角0°—180°范圍內(nèi)的山坡定義為東坡向；90°—270°為南坡向；180°—360°為西坡向；270°—90°為北坡向。各坡向范圍之間存在一定的重疊。

各坡向紅松解析木胸徑-年齡最佳回歸方程擬合效果依次為南坡>西坡>東坡>北坡。在各個坡向中,處于南坡的回歸方程以指數(shù)模型擬合效果最佳,達(dá)到了0.594,而以北坡的擬合效果最差。從這一現(xiàn)象來看,不同坡向?qū)t松胸徑-年齡之間關(guān)系的影響存在差異,尤其在北坡受到的影響最大。

表4 各坡向胸徑與年齡最佳擬合方程

分別用各個坡向的數(shù)據(jù)建立胸徑與年齡、競爭指標(biāo)、坡度、海拔等之間的多元回歸方程(表5),并與表4中相對應(yīng)坡向的單因子回歸方程決定系數(shù)進(jìn)行比較。

由表5可知,采用東坡向數(shù)據(jù)建立的方程,競爭因子、海拔等變量皆被剔除,說明對分布于東坡的紅松解析木,在其胸徑-年齡關(guān)系中加進(jìn)入競爭因子和海拔變量對提高擬合優(yōu)度沒有貢獻(xiàn)。在南坡,基于表4中胸徑-年齡最佳線性擬合曲線(指數(shù)模型)效果(R2=0.594),加入WBAL與海拔指標(biāo)的回歸方程,其擬合效果有所改善(修正R2=0.673)。在各坡向中,加入競爭指標(biāo)之后對提高回歸方程擬合優(yōu)度最明顯的是北坡與西坡,例如加入WBAL指標(biāo),決定系數(shù)分別提高了0.355和0.269。說明競爭因子對北坡與西坡的紅松解析木胸徑-年齡關(guān)系有較大的影響。

表 5 各坡向多元線性回歸方程

3.3.2 坡位對紅松胸徑與年齡關(guān)系的影響

從表6可以看出分布于中坡位的紅松解析木胸徑-年齡回歸方程其修正系數(shù)達(dá)到了0.717。而上坡位與下坡位的擬合效果較差。

表6 各坡位胸徑與年齡最佳擬合方程

分別用各個坡位的數(shù)據(jù)建立年齡與胸徑、競爭指標(biāo)、坡度、海拔等之間的多元回歸方程(表7)。并與表6中相對應(yīng)坡位的單因子回歸方程決定系數(shù)進(jìn)行比較。

在各坡位,采用CI與WBAL競爭指標(biāo)的回歸方程,其擬合優(yōu)度都得到了改善,其中以采用WBAL指標(biāo)的效果更好一些。在中坡位,加入WBAL與海拔因子之后,回歸方程修正決定系數(shù)達(dá)到了0.788,是擬合效果最好的一組方程。而加入WBAL對改善下坡位的胸徑-年齡回歸方程效果最為明顯,這也說明下坡位的紅松解析木的胸徑生長受到競爭因子的影響最大。

表7 各坡位多元線性回歸方程

本文將分布于南坡向與中坡位的解析木篩選出來,共計39株,所做的胸徑-年齡回歸方程以指數(shù)函數(shù)模型擬合效果最佳,方程式為lnA=2.115e0.239lnD,決定系數(shù)達(dá)到0.809,進(jìn)一步建立年齡與胸徑、競爭指標(biāo)、坡度、海拔等之間的多元回歸方程(表8)。結(jié)果表明采用WBAL、胸徑變量對年齡的擬合效果有所改善；進(jìn)一步引入海拔因子,修正決定系數(shù)提高到了0.858。

表8 南向中坡位多元回歸方程

4 結(jié)論與討論

影響紅松個體生長的因素是多方面的。立地條件、生長節(jié)律、周圍林木競爭等等,都會造成個體之間差異,導(dǎo)致胸徑-年齡關(guān)系復(fù)雜化,增加依靠胸徑推斷年齡的困難。本文利用所搜集的小興安嶺紅松解析木資料,對影響成熟紅松個體胸徑-年齡關(guān)系的因素作了初步探討。針對所搜集的全部解析木估計其胸徑與年齡的最佳擬合曲線。結(jié)果表明在胸徑與年齡變量均取對數(shù)的前提下,以指數(shù)模型進(jìn)行擬合,其決定系數(shù)也只有0.507。侯向陽認(rèn)為,紅松林中林木年輪徑向生長受其周圍其他林木的空間競爭影響是巨大的。這種生存空間競爭對林木生長的影響往往超過任何其他因素的影響[28]。但是在引入競爭因子后,并沒有像預(yù)計的那樣大幅度提高對年齡的預(yù)測精度,盡管采用距離加權(quán)的大樹斷面積指標(biāo)(WBAL)比簡單競爭指數(shù)(CI)更好地體現(xiàn)了競爭因子的影響。李克志對小興安嶺紅松生長的分析表明,年齡在100a以后,林木已擺脫被壓狀態(tài)[7]。孫玉軍將小興安嶺紅松的生長劃分為4個階段,在90a以下為樹木生長緩慢階段,100—130a為生長加速階段,在此階段,競爭減弱,胸徑生長明顯加快[8]。而本研究所涉及的紅松解析木個體,年齡都在100a以上,胸徑基本都大于20 cm,所受到的競爭應(yīng)該是比較小的。以往有研究表明,紅松胸徑小于25 cm時,所承受的競爭壓力較大；胸徑大于25 cm以后,競爭指數(shù)變化不大,并維持較低的水平[13]。

為了探討地形因子對成熟紅松胸徑-年齡關(guān)系的影響,分別按照坡向與坡位對紅松解析木數(shù)據(jù)進(jìn)行分類,建立胸徑-年齡擬合方程。結(jié)果表明南坡向或中坡位的紅松解析木,其胸徑-年齡單因子回歸方程的擬合優(yōu)度大幅提高,而南坡向與中坡位的地形因子組合使紅松胸徑-年齡單因子回歸方程的擬合效果達(dá)到了最優(yōu),引入競爭因子與海拔之后進(jìn)一步改善了擬合效果。在其他坡向或坡位建立的方程雖然擬合效果較差,但在引入競爭因子之后,則明顯提高了對年齡的預(yù)測精度。由以上模型得出紅松胸徑生長是十分復(fù)雜的過程,在不同坡向、坡位其生長規(guī)律有所不同,在南坡向、中坡位其生長規(guī)律較為一致,這可能與其充足光照和相對較高的溫度有關(guān)。而在西、北坡向其生長差異較大,但是引入競爭指標(biāo)后其擬合精度雖然不是很理想但是卻有了較大的提高,可能由于某些生長因子缺乏如光、熱量,而不能滿足所有個體的生長需求從而凸顯競爭優(yōu)勢的重要性,使得量化的競爭指標(biāo)對其擬合優(yōu)度有了大為促進(jìn)作用。

通過本文分析可知,對于所搜集的小興安嶺地區(qū)紅松解析木個體,坡向、坡位、海拔,以及競爭指標(biāo)等都影響著紅松胸徑與年齡的關(guān)系,其中對紅松胸徑-年齡關(guān)系影響最明顯的是地形因子中的坡向與坡位,在南坡向與中坡位的紅松個體,其胸徑與年齡之間表現(xiàn)出了最佳的擬合效果。而加入競爭因子卻沒有大幅度提高紅松胸徑-年齡回歸方程的擬合優(yōu)度,并且競爭因子作用的發(fā)揮可能受到地形因子尤其是坡向與坡位的制約, 坡向、坡位能夠控制太陽輻射和降水的空間再分配,營造局部小氣候[29],但二者又是通過哪些具體機(jī)制來影響紅松胸徑與年齡的關(guān)系還有待進(jìn)一步研究。并且現(xiàn)有的分析方法只能體現(xiàn)地上部分的競爭,并不能體現(xiàn)地下部分根系的競爭情況和林木對光、溫度的利用情況[30]。

[1] 代力民, 孫偉中, 鄧紅兵, 肖以華. 長白山北坡椴樹闊葉紅松林群落主要樹種的年齡結(jié)構(gòu)研究. 林業(yè)科學(xué), 2002, 38(3): 73- 77.

[2] Zhang Y, Drobyshev I, Gao L S, Zhao X H, Bergeron Y. Disturbance and regeneration dynamics of a mixed Korean pine dominated forest on Changbai Mountain, North-Eastern China. Dendrochronologia, 2014, 32(1): 21- 31.

[3] 胡云云, 亢新剛, 趙俊卉. 長白山地區(qū)天然林林木年齡與胸徑的變動關(guān)系. 東北林業(yè)大學(xué)學(xué)報, 2009, 37(11): 38- 42.

[4] 徐德應(yīng). 林分和林木生長預(yù)測模型. 世界林業(yè)研究, 1988, (3): 42- 49.

[5] 鄧?yán)っ? 邵彬, 李飛. 長白山北坡云冷杉林胸徑、樹高結(jié)構(gòu)及其生長規(guī)律的分析. 資源科學(xué), 1999, 21(1): 77- 84.

[6] Wimberly M C, Bare B B. Distance-dependent and distance-independent models of Douglas-fir and western hemlock basal area growth following silvicultural treatment. Forest Ecology and Management, 1996, 89(1/3): 1- 11.

[7] 李克志. 紅松生長過程的研究. 林業(yè)科學(xué), 1983, 19(2): 126- 136.

[8] 孫玉軍, 劉吉春, 陳霖生, 張鵬, 劉傳照, 戴興安, 李?？? 天然紅松單木生長階段的多元定量化研究. 東北林業(yè)大學(xué)學(xué)報, 1993, 21(6): 1- 5.

[9] 邵國凡. 紅松人工林單木生長模型的研究. 東北林業(yè)大學(xué)學(xué)報, 1985, 13(3): 38- 46.

[10] 蔣伊尹. 紅松人工林生長與生長模型. 東北林學(xué)院學(xué)報, 1985, 13(2): 5- 15.

[11] 黃新峰, 亢新剛, 孫玲, 馮啟祥, 姚景春. 紅松單木斷面積生長模型. 西北林學(xué)院學(xué)報, 2011, 26(3): 143- 146.

[12] 劉紅潤, 李鳳日. 紅松天然林種內(nèi)和種間競爭關(guān)系的研究. 植物研究, 2010, 30(4): 479- 484.

[13] 趙俊卉, 亢新剛, 張慧東, 劉燕. 長白山主要針葉樹種胸徑和樹高變異系數(shù)與競爭因子的關(guān)系應(yīng)用生態(tài)學(xué)報, 2009, 20(8): 1832- 1837.

[14] Holmes M J, Reed D D. Competition indices for mixed species Northern Hardwoods. Forest Science, 1991, 37(5): 1338- 1349.

[15] Avery T E, Burkhart H E. Forest Measurements. New York: McGraw-Hill, 1983: 1- 331.

[16] Radtke P J, Westfall J A, Burkhart H E. Conditioning a distance-dependent competition index to indicate the onset of inter-tree competition. Forest Ecology and Management, 2003, 175(1/3): 17- 30.

[17] 沈琛琛, 雷相東, 王褔有, 馬武, 沈劍波. 金蒼林場蒙古櫟天然中齡林競爭關(guān)系研究. 林業(yè)科學(xué)研究, 2012, 25(3): 339- 345

[18] 張思玉, 鄭世群. 筆架山常綠闊葉林優(yōu)勢種群種內(nèi)種間競爭的數(shù)量研究. 林業(yè)科學(xué), 2001, 7(1): 185- 188.

[19] 鄒春靜, 王慶禮, 韓士杰. 長白山暗針葉林建群種競爭關(guān)系的研究. 應(yīng)用與環(huán)境生物學(xué)報, 2001, 7(2): 101- 105.

[20] 李建軍, 李際平, 劉素青, 趙春燕, 張宏偉, 馮湘蘭. 基于Hegyi改進(jìn)模型的紅樹林空間結(jié)構(gòu)競爭分析. 中南林業(yè)科技大學(xué)學(xué)報, 2010, 30(12): 23- 27.

[21] Wykoff W R. A basal area increment model for individual conifers in the northern Rocky Mountains. Forest Science, 1990, 36(4): 1077- 1104.

[22] Stage A R, Ledermann T. Effects of competitor spacing in a new class of individual-tree indices of competition: semi-distance-independent indices computed for Bitterlich versus fixed-area plots. Canadian Journal of Forest Research, 2008, 38(4): 890- 898.

[23] Roberts S D, Harrington C A. Individual tree growth response to variable-density thinning in coastal Pacific Northwest forests. Forest Ecology and Management, 2008, 255(7): 2771- 2781.

[24] Busing R T. A spatial model of forest dynamics. Vegetatio, 1991, 92: 167- 179.

[25] 葛劍平. 森林生態(tài)學(xué)建模與仿真. 哈爾濱: 東北林業(yè)大學(xué)出版社, 1996.

[26] 雷相東, 張則路, 陳曉光. 長白落葉松等幾個樹種冠幅預(yù)測模型的研究. 北京林業(yè)大學(xué)學(xué)報, 2006, 28(6): 75- 79.

[27] 童躍偉, 項文化, 王正文, Durka W, Fischer M. 地形、鄰株植物及自身大小對紅楠幼樹生長與存活的影響. 生物多樣性, 2013, 21(3): 269- 277.

[28] 侯向陽, 韓進(jìn)軒, 陽含熙. 長白山紅松闊葉林林冠木競爭生長及林冠空隙動態(tài)研究. 生態(tài)學(xué)報, 2000, 20(1): 68- 72.

[29] Cantón Y, Del Barrio G, Solé-Benet A, Lázaro R. Topographic controls on the spatial distribution of ground cover in the Tabernas badlands of SE Spain. Catena, 2004, 55(3): 341- 365.

[30] Bugmann H. A review of forest gap models. Climatic Change, 2001, 51(3/4): 259- 305.

Influence of topography and competitive factors on the relationship between DBH and age of Korean pine

YANG Huibin, GUO Qingxi*

CenterforEcologicalResearch,NortheastForestryUniversity,Harbin150040,China

The ability to accurately age a forest has always been an important component of forest ecology research. Traditionally, the diameter of a tree at breast height (DBH) has been used to estimate tree age and subsequently identify the age structure of a forest. Unfortunately, the DBH index is not an accurate estimator of tree age for natural Korean pines (Pinuskoraiensis) because of the complexity of spatial distribution and different developmental stages of this species. Previous studies have examined the relationship between DBH and age of natural Korean pines using a simple regression analysis. In our study, the wood data of 156 mature Korean pines in the Lesser Khingan Mountains of the Heilongjiang province was collected to analyze the relationship between DBH and age of Korean pine, as well as the influences of topographic and competition factors. In our wood data analysis, we preformed simple regression matching with an exponential model to better suit the biological significance and distributional characteristics of age for natural Korean pines. In addition, we adopted the widely used Hegyi simple competition index (CI) and improved Basal-Area-in-Larger Trees (WBAL) as competition indexes. We imported the competition indexes into the (DBH-age) regression equation and analyzed them by different slopes and slope positions using multiple regression analysis to understand how topography and competitive factors influence the relationship between DBH and age of natural Korean pines. Our results showed that the Pearson correlation coefficient of DBH and age was 0.562, and for our wood data collected and analyzed, the linear correlation between the DBH and age was not high. When we fit our wood data with exponential models the determination coefficient R2was only 0.507. However, when we imported the competition indexes for multiple regression analysis, the adjusted determination coefficient of the regression model increased to 0.542; when we imported the WBAL index, the fitting effect of DBH-age improved to 0.583. We conducted fitting analysis using different slopes and slope positions and found that natural Korean pines in the middle slope and southward slopes showed the best fitting effects between the DBH and age. However, after importing the competition factors, we significantly improved our precision in estimating age, especially with pines in the northward and westward slopes. For example, after WBAL was imported, the determination coefficient for DBH and age increased by 0.355 and 0.269, respectively. Our study found that altitude, slope, aspect, and competition factors could all influence the relationship between DBH and age of natural Korean pines. Among them the aspect and slope position had the strongest influence; however, the competition factors showed limited influence, depending on topography. However, the underlying mechanism of effects of topography and competitive factors on the relationship between DBH and age of Korean pine is complicated. This article has only carried on the preliminary exploration. More investigations should be conducted in the future to further improve fitness of the relationship between DBH and age of Korean pine.

DBH; forest age; competition factors; topographic factors; Korean pine

東北溫帶針闊混交林區(qū)生物多樣性調(diào)查(2012FY112000)

2015- 03- 09;

日期:2016- 01- 22

10.5846/stxb201503090457

*通訊作者Corresponding author.E-mail: gqx@nefu.edu.cn

楊惠濱, 國慶喜.地形與競爭因子對紅松胸徑與年齡關(guān)系的影響.生態(tài)學(xué)報,2016,36(20):6487- 6495.

Yang H B, Guo Q X.Influence of topography and competitive factors on the relationship between DBH and age of Korean pine.Acta Ecologica Sinica,2016,36(20):6487- 6495.