張崗崗，王得祥＊，張明霞，劉文楨，郭小龍

（1 西北農(nóng)林科技大學(xué) 林學(xué)院，陜西楊陵712100；2 甘肅小隴山森林生態(tài)系統(tǒng)國(guó)家定位觀測(cè)研究站，甘肅天水741020；3 青海省社會(huì)科學(xué)院，西寧810000；4 小隴山林業(yè)科學(xué)研究所，甘肅天水741020）

植物群落是由在時(shí)間和空間上彼此相關(guān)聯(lián)的共存物種組成，不同時(shí)空各物種之間的相互關(guān)系決定著群落的結(jié)構(gòu)特征和動(dòng)態(tài)［1］。研究種間關(guān)系不僅可以反映群落內(nèi)物種的數(shù)量特征，更重要的是能有效反映群落中各物種對(duì)環(huán)境因子的適應(yīng)程度及特定環(huán)境因子作用下種間相互依存和制約關(guān)系，對(duì)于正確認(rèn)識(shí)群落演替過(guò)程中種間替代和共存關(guān)系演變機(jī)制具有重要意義［2－5］。種間聯(lián)結(jié)是把成對(duì)物種的存在與否作為兩個(gè)物種出現(xiàn)的相似性的定性特征，種間相關(guān)則是一種定量關(guān)系，二者通常是由群落生境的差異影響了物種的分布和生長(zhǎng)而引起的，反映了不同物種在時(shí)間與空間上的相互關(guān)系［6－8］。基于種間關(guān)系的生態(tài)種組是從強(qiáng)聯(lián)結(jié)的種對(duì)中提煉出來(lái)，反映了種對(duì)在生態(tài)過(guò)程中的異質(zhì)性和物種生態(tài)需求的趨同性，也是研究種間關(guān)系的基礎(chǔ)［9］，對(duì)其進(jìn)行劃分有助于客觀地了解區(qū)域植物群落的結(jié)構(gòu)，并有助于探討不同物種之間、物種與環(huán)境間關(guān)系，為物種多樣性的保護(hù)、自然植被的恢復(fù)與重建提供理論依據(jù)［4－5，10］。國(guó)內(nèi)外學(xué)者對(duì)不同群落種間關(guān)系和生態(tài)種 組劃分進(jìn)行了大量研究［4－5，10－14］，但這些研究主要針對(duì)某一特定群落或群落某一演替階段優(yōu)勢(shì)種間關(guān)系進(jìn)行定性分析的，缺乏對(duì)不同演替階段群落種間關(guān)系的定性和定量雙重探討，在生態(tài)種對(duì)層面劃分生態(tài)種組的研究也鮮見(jiàn)報(bào)道。

在海拔800～2 300m 之間的秦嶺山地森林生態(tài)系統(tǒng)中，油松林（Pinustabulaeformisforest）、銳齒櫟林（Quercusalienavar.acuteserrataforest）及其混交林 （P.tabulaeformis－Q.alienavar.acuteserratamixed forest）是該區(qū)3 種典型的森林群落類型［15，16］。國(guó)內(nèi)學(xué)者關(guān)于秦嶺地區(qū)油松林、松櫟混交林和銳齒櫟林這3個(gè)群落是否處于同一演替序列中的不同階段還存在爭(zhēng)議，但部分學(xué)者認(rèn)為松類樹(shù)種是先鋒種，松櫟混交林是松類林受到干擾后整個(gè)群落處在被櫟類植物逐漸取代的初級(jí)演替階段，而櫟類的競(jìng)爭(zhēng)能力比較強(qiáng)，其取代油松的趨勢(shì)十分明顯［17，18］，而對(duì)這3個(gè)群落更新特征的研究也表明油松林群落終將為銳齒櫟林群落所取代［19－23］。最近，部分學(xué)者單獨(dú)對(duì)這3種群落進(jìn)行更新特征、生態(tài)位特征、種間聯(lián)結(jié)等方面研究［9，11，19－24］，張子良等［15］、張洪武等［16］則分別探討了這3種群落喬灌草層物種組成與環(huán)境因子、物種多樣性之間的關(guān)系，但針對(duì)秦嶺林區(qū)油松林→松櫟混交林→銳齒櫟林這一演替序列過(guò)程中種間關(guān)系動(dòng)態(tài)變化研究尚屬空白。秦嶺松櫟林演替過(guò)程中不同演替階段群落優(yōu)勢(shì)種的總體聯(lián)結(jié)性如何變化？主要種群種間聯(lián)結(jié)性的顯著程度及大小如何？群落演替過(guò)程中種間相關(guān)性如何？這些問(wèn)題和規(guī)律的深入研究和探討，有助于闡明不同演替階段物種之間替代和共存關(guān)系，并為推演秦嶺林區(qū)松櫟次生林演替趨勢(shì)奠定基礎(chǔ)。鑒于此，本研究采用空間替代時(shí)間的方法，選取秦嶺南坡典型的油松林、松櫟混交林、銳齒櫟林為研究對(duì)象，闡述了不同演替階段群落種間聯(lián)結(jié)性和相關(guān)性動(dòng)態(tài)變化規(guī)律，以期揭示優(yōu)勢(shì)種群種間關(guān)系和群落演替的相互作用，為秦嶺山地松櫟林演替過(guò)程中健康穩(wěn)定的群落關(guān)系構(gòu)建提供科學(xué)依據(jù)。

1 研究區(qū)概況

研究區(qū)域?yàn)榍貛X南坡中段主要包括秦嶺南坡至漢水以北的漢中市各（區(qū)）縣以及與之相鄰的鳳縣部分地區(qū)，即E105°30′～108°11′，N32°50′～33°53′。該區(qū)地形地貌極為復(fù)雜，全年雨量充沛，氣候溫暖濕潤(rùn)，差異明顯，海拔800 m 以下為北亞熱帶，800～2 400m為山地暖溫帶，2 400 m 以上為高山寒溫帶。本研究具體地點(diǎn)是以火地塘、旬陽(yáng)壩等為主的5個(gè)松櫟林典型分布區(qū)域。

秦嶺火地塘林區(qū)位于秦嶺南坡中山地帶中部（33°25′～33°29′N，108°25′～108°30′E），海拔800～2 500m，年均氣溫8～10 ℃，年均降雨量900～1 200mm，年日照時(shí)數(shù)1 327.5h，屬中國(guó)北亞熱帶和暖溫帶的過(guò)渡地帶。以山地棕壤、暗棕壤和山地草甸土為主，土層厚30～50cm，該區(qū)屬溫帶針闊混交林和寒帶針葉林，從低海拔到高海拔依次為：松櫟林亞帶，松樺林亞帶，云、冷杉針闊葉混交林亞帶。主要樹(shù)種有：銳齒櫟（Quercusalienavar.acuteserrata）、紅 樺（Betulaalbo－sinensis）、油 松（Pinus tabulaeformis）、華山松（Pinusarmandii）等。旬陽(yáng)壩林區(qū)位于秦嶺中段南麓的月河梁和平河梁之間（32°29′～33°13′N，103°58′～109°48′E），平均海拔1 300m，年均氣溫為10 ℃，年均降水量1 133mm，年日照時(shí)數(shù)1 638.3h，該區(qū)屬北亞熱帶溫暖濕潤(rùn)氣候區(qū)。土壤為礦礫質(zhì)壤黏土，呈微酸性，該區(qū)屬暖溫帶落葉闊葉林和針闊混交林向北亞熱帶常綠落葉混交林的過(guò)渡帶，林木主要有松、杉和落葉喬木等。

2 研究方法

2.1 野外調(diào)查

采用空間代替時(shí)間的方法，在對(duì)油松林→松櫟混交林→銳齒櫟林這一演替序列分布較集中的火地塘、旬陽(yáng)壩等林區(qū)全面踏查的基礎(chǔ)上，選擇樣地時(shí)盡可能縮小各階段間的空間距離，并盡量保持立地因子一致、人為干擾較小的自然演替群落，每個(gè)群落類型（階段）采用典型取樣法設(shè)置重復(fù)樣地15個(gè)，樣地面積為20m×20m，在每個(gè)樣地四角和中心分別設(shè)置5個(gè)2m×2m 灌木樣方和1m×1m 草本樣方。對(duì)樣地內(nèi)所有喬木進(jìn)行每木檢尺、坐標(biāo)定位，記錄其物種名稱、胸徑、樹(shù)高、冠幅和生長(zhǎng)狀況等；灌草樣方記錄種類、株數(shù)、蓋度、分布和生長(zhǎng)狀況等；詳細(xì)記錄海拔、坡度、坡向、坡位等立地因子。45個(gè)調(diào)查樣地基本情況，即樣地編號(hào)、經(jīng)緯度、海拔、坡度、坡向、坡位、郁閉度、林齡、土壤等立地條件均已在參考文獻(xiàn)［15，16］中詳細(xì)列出，本文不再贅述。

2.2 數(shù)據(jù)統(tǒng)計(jì)與分析

分別計(jì)算秦嶺南坡油松林、松櫟混交林、銳齒櫟林各15個(gè)樣地喬木層物種重要值的平均值，列表統(tǒng)計(jì)重要值排序前20的種群（表1）在各樣方中出現(xiàn)與否的1、0二元數(shù)據(jù)矩陣，采用方差比率、χ2檢驗(yàn)、2×2聯(lián)列表、Jaccard指數(shù)、Spearman秩相關(guān)系數(shù)等分析方法，分別對(duì)松櫟林不同演替階段群落中主要種群間總體相關(guān)性、種對(duì)間的聯(lián)結(jié)性、種對(duì)間聯(lián)結(jié)強(qiáng)度和種間相關(guān)性進(jìn)行計(jì)算、統(tǒng)計(jì)和分析。具體方法參見(jiàn)文獻(xiàn)［3，6－8，25］。

表1 秦嶺南坡不同演替階段3種森林群落物種組成及重要值Table1 Importance values and species components of three forest communities at different stages of succession on the south－facing slopes in Qinling Mountains

3 結(jié)果與分析

3.1 群落不同演替階段優(yōu)勢(shì)種總體關(guān)聯(lián)性分析

群落中優(yōu)勢(shì)樹(shù)種間的總體聯(lián)結(jié)性反映了群落內(nèi)各樹(shù)種間關(guān)聯(lián)性的總體趨勢(shì)。不同演替階段優(yōu)勢(shì)種間的總體聯(lián)結(jié)性分析結(jié)果（表2）表明，油松林、松櫟混交林和銳齒櫟林優(yōu)勢(shì)種間的方差比率分別為1.81、1.97、2.45，均大于1，表明在獨(dú)立零假設(shè)條件下，總體聯(lián)結(jié)性表現(xiàn)為凈的正聯(lián)結(jié)；顯著性χ2檢驗(yàn)統(tǒng)計(jì)量W 分別為27.14、29.61、36.73，落在區(qū)間（7.261，27.996）外，說(shuō)明群落總體的正聯(lián)結(jié)程度均顯著，且隨演替進(jìn)程而加強(qiáng)，因此，不同演替階段的群落均表現(xiàn)為顯著正聯(lián)結(jié)，種間關(guān)系緊密聯(lián)系，各群落具有相對(duì)的穩(wěn)定性，也側(cè)面反映了主要種群所在的林區(qū)大環(huán)境與各群落生境需求的一致性，使其更易組成群落并保持其穩(wěn)定性。

3.2 群落不同演替階段主要種對(duì)聯(lián)結(jié)性分析

種間聯(lián)結(jié)的χ2檢驗(yàn)?zāi)鼙容^準(zhǔn)確地反映出種間聯(lián)結(jié)的顯著程度，JI值則能體現(xiàn)出χ2檢驗(yàn)不顯著種對(duì)聯(lián)結(jié)性的大小，其值越接近1，表明該種對(duì)的正聯(lián)結(jié)越緊密，其值越接近0，則表明負(fù)聯(lián)結(jié)越強(qiáng)。

表2 秦嶺南坡松櫟林群落不同演替階段群落優(yōu)勢(shì)種總體關(guān)聯(lián)性Table2 The overall interspecific connectivity of dominant species in P.tabulaeformis－Q.alienavar.acuteserrata mixed forest communities at different stages of succession on the south－facing slopes in Qinling Mountain

圖1 秦嶺南坡油松林群落主要種對(duì)間χ2 檢驗(yàn)和JI半矩陣圖Fig.1 The semi－matrix diagram aboutχ2 test and JI values of dominant species in P.tabulaeformis forest communities on the south－facing slopes in Qinling Mountain

在油松林群落中（圖1），主要種群χ2檢驗(yàn)中85對(duì)為正關(guān)聯(lián)，占總種對(duì)數(shù)的44.7%，其中極顯著正關(guān)聯(lián)的有1對(duì)，顯著正關(guān)聯(lián)的有6對(duì)；有103對(duì)為負(fù)關(guān)聯(lián)，占總對(duì)數(shù)的54.2%，均不顯著；有2對(duì)為無(wú)關(guān)聯(lián)，群落中正關(guān)聯(lián)的種對(duì)數(shù)明顯小于負(fù)關(guān)聯(lián)，顯著或極顯著的種對(duì)僅占4.7%，為銳齒櫟－漆樹(shù)、青榨槭－三椏烏藥、鵝耳櫪－三椏烏藥、梾木－栓皮櫟等，這些種對(duì)生態(tài)學(xué)特性相近，生境需求相似，種間關(guān)系協(xié)同進(jìn)化，趨于穩(wěn)定，其余大多數(shù)種對(duì)則處于不顯著水平，說(shuō)明油松群落種間聯(lián)結(jié)性較松散，這與總體聯(lián)結(jié)性分析結(jié)果一致。群落中JI＞0.75的有9對(duì)，0.75≥JI＞0.5的有28對(duì)，約占總種對(duì)數(shù)的14.7%；0.5≥JI＞0.25的有97對(duì)，約占總對(duì)數(shù)的51.1%；0.25≥JI≥0的有56對(duì)，約占總種對(duì)數(shù)的29.5%。進(jìn)一步說(shuō)明個(gè)別種對(duì)間對(duì)資源和空間存在一定的競(jìng)爭(zhēng)，但不激烈，多數(shù)種對(duì)種間關(guān)聯(lián)程度不緊密，這與χ2檢驗(yàn)結(jié)果相似。因此，油松林群落整體上處于較為穩(wěn)定的階段。

松櫟混交林群落中（圖2），主要種群χ2檢驗(yàn)中有81對(duì)為正關(guān)聯(lián)，占總對(duì)數(shù)的42.6%，其中顯著正關(guān)聯(lián)3對(duì)，分別為華山松－漆樹(shù)、漆樹(shù)－鵝耳櫪、野核桃－鵝耳櫪；有106對(duì)為負(fù)關(guān)聯(lián)，占總對(duì)數(shù)的55.6%，其中顯著負(fù)關(guān)聯(lián)的1對(duì)，為漆樹(shù)－燈臺(tái)樹(shù)；其余3對(duì)為無(wú)關(guān)聯(lián)，正負(fù)關(guān)聯(lián)比為0.76，絕大多數(shù)種對(duì)處于不顯著水平，說(shuō)明松櫟混交林群落中多數(shù)種對(duì)間存在微弱的相互競(jìng)爭(zhēng)或依賴關(guān)系。群落中JI＞0.75的僅有8對(duì)，占總種對(duì)數(shù)的4.2%；0.75≥JI＞0.5的有21對(duì)，約占總種對(duì)數(shù)的11.1%；0.5≥JI＞0.25的有113對(duì)，占總種對(duì)數(shù)的59.5%；0.25≥JI≥0的有48對(duì)，約占總種對(duì)數(shù)的25.3%，說(shuō)明群落中多數(shù)種對(duì)聯(lián)結(jié)程度較弱，群落結(jié)構(gòu)和狀態(tài)相對(duì)穩(wěn)定。值得注意的是，χ2檢驗(yàn)結(jié)果表明銳齒櫟、油松、華山松、漆樹(shù)之間正聯(lián)結(jié)關(guān)系并不顯著，但JI值遠(yuǎn)大于0.75，說(shuō)明這些樹(shù)種之間為顯著或極顯著正關(guān)聯(lián)，因此，為避免χ2檢驗(yàn)結(jié)果產(chǎn)生偏差，引入JI值測(cè)定種間關(guān)聯(lián)程度是必要的。

銳齒櫟林群落中（圖3），主要種群χ2檢驗(yàn)中有98對(duì)為正關(guān)聯(lián)，占總對(duì)數(shù)的51.6%，其中顯著正關(guān)聯(lián)2對(duì)，為鵝耳櫪－千金榆、建始槭－刺楸；有88對(duì)為負(fù)關(guān)聯(lián)，占總對(duì)數(shù)的55.6%，均不顯著；有4對(duì)為無(wú)關(guān)聯(lián)，正負(fù)關(guān)系比為1.11，檢驗(yàn)顯著率僅為1.05%，大部分種對(duì)間正或負(fù)關(guān)聯(lián)性不顯著，對(duì)資源和空間的競(jìng)爭(zhēng)較為緩和，種間關(guān)系趨于穩(wěn)定共存，群落總體處于比較穩(wěn)定的狀態(tài)。群落中JI＞0.75的種對(duì)僅有4 對(duì)，漆樹(shù)－華山松、梾木－鵝耳櫪也呈顯著正關(guān)聯(lián)；0.75≥JI＞0.5的有31對(duì)，占總對(duì)數(shù)的16.3%；0.5≥JI＞0.25的有118對(duì)，占總種對(duì)數(shù)的62.1%；0.25≥JI≥0的有37對(duì)，占總種對(duì)數(shù)的19.5%。說(shuō)明大多數(shù)種對(duì)間的關(guān)聯(lián)程度不高，物種分布具有一定的獨(dú)立性，群落穩(wěn)定性較高。

圖2 秦嶺南坡松櫟混交林群落主要種對(duì)間χ2 檢驗(yàn)和JI半矩陣圖Fig.2 The semi－matrix diagram aboutχ2 test and JI values of dominant species in P.tabulaeformis－Q.aliena var.acuteserrata mixed forest communities on the south－facing slopes in Qinling Mountain

圖3 秦嶺南坡銳齒櫟林群落主要種對(duì)間χ2 檢驗(yàn)和JI半矩陣圖Fig.3 The semi－matrix diagram aboutχ2 test and JI values of dominant species in Q.alienavar.acuteserrata forest communities on the south－facing slopes in Qinling Mountain

綜上所述，油松林→松櫟混交林→銳齒櫟林這一演替序列過(guò)程中，前期和中期群落物種間正聯(lián)結(jié)比例接近，分別為44.7%和42.6%，后期增加到51.6%；負(fù)聯(lián)結(jié)比例從54.2%略增大到55.6%后減小至46.3%；均有少量無(wú)關(guān)聯(lián)種對(duì)出現(xiàn)，說(shuō)明隨著演替的進(jìn)行，群落中主要種群對(duì)生態(tài)要求逐漸趨于相似；前期和中期群落中達(dá)到顯著或極顯著水平的正關(guān)聯(lián)種對(duì)占有一定比例，均多于后期，負(fù)關(guān)聯(lián)中達(dá)到顯著水平及以上的種對(duì)極少；JI≥0.5和0.25≥JI≥0的比例逐漸減少，0.5≥JI＞0.25逐漸增多，表明群落演替過(guò)程中種間關(guān)聯(lián)程度趨于弱度聯(lián)結(jié)，因此，群落在長(zhǎng)期的協(xié)同進(jìn)化過(guò)程中種間關(guān)系逐漸趨于穩(wěn)定共存。其次，χ2檢驗(yàn)中顯著正關(guān)聯(lián)種對(duì)JI值都比較高，顯著負(fù)關(guān)聯(lián)種對(duì)的JI值較小，即JI在反映種對(duì)之間關(guān)聯(lián)性和關(guān)聯(lián)程度方面與χ2檢驗(yàn)結(jié)果相一致，與郭垚鑫等［13］研究結(jié)果基本一致。

3.3 不同演替階段物種之間的聯(lián)結(jié)性分析

根據(jù)相關(guān)研究［26］及演替特性［27］，將各演替階段重要值排名前20的喬木劃分為4個(gè)樹(shù)種組，其中油松、華山松、山楊、冬瓜楊、糙皮樺為先鋒樹(shù)種，青榨槭、青麩楊、梾木、鵝耳櫪、刺楸、榛子、木姜子、三椏烏藥為前期過(guò)渡種，燈臺(tái)樹(shù)、化香樹(shù)、榛子、梾木、鵝耳櫪、千金榆、椴樹(shù)、建始槭、刺楸、榛子、木姜子、三椏烏藥為后期過(guò)渡種，銳齒櫟、漆樹(shù)、野核桃、鐵杉為頂級(jí)樹(shù)種。由演替過(guò)程中各樹(shù)種組組內(nèi)及組間種對(duì)聯(lián)結(jié)性統(tǒng)計(jì)結(jié)果（表3）可知，先鋒樹(shù)種之間正聯(lián)結(jié)和負(fù)聯(lián)結(jié)比例分別為33.3%和66.7%，均不顯著，華山松與油松之間競(jìng)爭(zhēng)較為激烈，JI值達(dá)0.67，總體處于不穩(wěn)定階段；先鋒種與前期過(guò)渡種之間正關(guān)聯(lián)種對(duì)占38.9%、負(fù)關(guān)聯(lián)種對(duì)占61.1%，正負(fù)關(guān)聯(lián)比增大，種間關(guān)系趨于緊密；隨著演替過(guò)程中的激烈競(jìng)爭(zhēng)，前期過(guò)渡種間、后期過(guò)渡種間及其前期過(guò)渡種和后期過(guò)渡種之間表現(xiàn)為較為有利的正聯(lián)結(jié)，分別占各自總種對(duì)數(shù)的59.2%、50.0%和53.4%，雖然正負(fù)聯(lián)結(jié)比均大于1，但不明顯，處于群落動(dòng)態(tài)演替過(guò)程中相對(duì)穩(wěn)定的階段；由于頂級(jí)樹(shù)種的侵入和競(jìng)爭(zhēng)，后期過(guò)渡種與頂級(jí)種之間正關(guān)聯(lián)比例明顯下降，負(fù)關(guān)聯(lián)比例高達(dá)63.4%，不同樹(shù)種之間對(duì)資源和空間實(shí)行再分配；頂級(jí)種間正聯(lián)結(jié)比例僅為22.2%，負(fù)關(guān)聯(lián)比例達(dá)77.8%，說(shuō)明這些頂級(jí)樹(shù)種之間競(jìng)爭(zhēng)較為激烈，尚未完全分化協(xié)調(diào)、演替至穩(wěn)定階段。整個(gè)演替過(guò)程中不同物種之間不斷競(jìng)爭(zhēng)資源和空間，物種組成及其種間關(guān)系處于動(dòng)態(tài)變化過(guò)程中，群落結(jié)構(gòu)也將不斷完善并趨于穩(wěn)定。

3.4 群落不同演替階段主要種對(duì)相關(guān)性分析

Spearman秩相關(guān)系數(shù)是檢驗(yàn)2個(gè)物種線性關(guān)系的重要指標(biāo)，其處理的是定量數(shù)據(jù)，反映2個(gè)物種同時(shí)出現(xiàn)的可能程度［3］。不同演替階段群落優(yōu)勢(shì)種間Spearman秩相關(guān)性分析結(jié)果中（圖4～圖6），油松林群落正相關(guān)種對(duì)為88對(duì)，占總對(duì)數(shù)的46.3%，其中極顯著正相關(guān)8對(duì)，顯著正相關(guān)4對(duì)；負(fù)相關(guān)種對(duì)有100對(duì)，占總對(duì)數(shù)的52.6%，其中，極顯著負(fù)相關(guān)1對(duì)，顯著負(fù)相關(guān)3對(duì)；無(wú)關(guān)種對(duì)有2對(duì)；正負(fù)相關(guān)比為0.88，檢驗(yàn)顯著率為6.32%。松櫟混交林群落正相關(guān)種對(duì)為89對(duì)，占總對(duì)數(shù)的46.8%，其中極顯著正相關(guān)4 對(duì)，顯著正相關(guān)14 對(duì)；負(fù)相關(guān)種對(duì)100對(duì)，占總對(duì)數(shù)的52.6%，其中，極顯著負(fù)相關(guān)4對(duì)，顯著負(fù)相關(guān)3對(duì)，無(wú)關(guān)種對(duì)1對(duì)；正負(fù)相關(guān)比為0.89，檢驗(yàn)顯著率為13.16%。銳齒櫟林群落正相關(guān)種對(duì)為99對(duì)，占總對(duì)數(shù)的52.1%，其中極顯著正相關(guān)9對(duì)，顯著正相關(guān)5對(duì)；負(fù)相關(guān)種89對(duì)，占總對(duì)數(shù)的46.8%，其中，顯著負(fù)相關(guān)的有1對(duì)；無(wú)關(guān)種對(duì)有2對(duì)；正負(fù)相關(guān)比為1.11，檢驗(yàn)顯著率為7.89%。由此可知，與銳齒櫟林群落相比，油松林和松櫟混交林群落種間正負(fù)相關(guān)比接近，且均明顯小于銳齒櫟林，群落中多數(shù)物種間的關(guān)系不緊密，這與χ2檢驗(yàn)和JI半矩陣分析結(jié)果基本一致，這符合前人研究結(jié)果［3，28］；顯著率則先增大后減小，但均大于χ2檢驗(yàn)

結(jié)果，主要是因?yàn)槔忙?檢驗(yàn)來(lái)判斷種間關(guān)聯(lián)性會(huì)損失一定的信息量［28］，而Spearman秩相關(guān)系數(shù)檢驗(yàn)靈敏度較高，可以彌補(bǔ)χ2檢驗(yàn)的不足，這與之前研究結(jié)果相似［28，29］。

表3 秦嶺南坡松櫟林群落演替過(guò)程中先鋒種、過(guò)渡種以及頂級(jí)種間的種間聯(lián)結(jié)Table3 The interspecific associations of pioneer species，transitionary species and climax species in P.tabulaeformis－Q.alienavar.acuteserrata mixed forest communities at different stages of succession on the south－facing slopes in Qinling Mountain

圖4 秦嶺南坡油松林群落優(yōu)勢(shì)種間Spearman秩相關(guān)系數(shù)半矩陣圖Fig.4 The semi－matrix diagram about Spearman’s rank correlation coefficient of dominant species in P.tabulaeformis forest communities on the south－facing slopes in Qinling Mountain

圖5 秦嶺南坡松櫟混交林群落優(yōu)勢(shì)種間Spearman秩相關(guān)系數(shù)半矩陣圖Fig.5 The semi－matrix diagram about Spearman’s rank correlation coefficient of dominant species in P.tabulaeformis－Q.alienavar.acuteserrata mixed forest communities on the south－facing slopes in Qinling Mountain

圖6 秦嶺南坡銳齒櫟林群落優(yōu)勢(shì)種間Spearman秩相關(guān)系數(shù)半矩陣圖Fig.6 The semi－matrix diagram about Spearman’s rank correlation coefficient of dominant species in Q.alienavar.acuteserrataforest communities on the south－facing slopes in Qinling Mountain

3.5 群落不同演替階段生態(tài)種組劃分

為了更真實(shí)可靠地反映群落種間關(guān)聯(lián)以及群落演替方向，依據(jù)不同演替階段群落種間關(guān)聯(lián)和相關(guān)性計(jì)算結(jié)果，以顯著負(fù)相關(guān)性為界限，同一生態(tài)種組內(nèi)的物種有盡可能大的正相關(guān)性和關(guān)聯(lián)性，物種之間具有相似的生態(tài)要求和生活習(xí)性［2，29］，將不同演替階段群落生態(tài)種組進(jìn)行劃分。油松林群落第1生態(tài)種組為華山松、冬瓜楊、山楊、油松，其中油松和華山松為主要建群種，存在一定的競(jìng)爭(zhēng)，但生態(tài)習(xí)性相似，華山松與冬瓜楊、山楊則呈顯著正關(guān)聯(lián)和相關(guān)性；第2生態(tài)種組為刺楸、鹽膚木、梾木、鵝耳櫪、三椏烏藥、青榨槭、木姜子，處于群落的中下層，為主要伴生種，種間關(guān)聯(lián)性和相關(guān)性較顯著；第3生態(tài)種組為銳齒櫟、栓皮櫟、漆樹(shù)、糙皮樺、榛子、微毛櫻桃，多為陽(yáng)生樹(shù)種，分布廣，聯(lián)結(jié)不明顯，但相關(guān)性較顯著；第4生態(tài)種組為青麩楊、毛梾、毛栗，處于劣勢(shì)，彼此間相互獨(dú)立。松櫟混交林群落第1生態(tài)種組為銳齒櫟、油松、華山松、鐵杉，為群落主要優(yōu)勢(shì)樹(shù)種，兩兩之間多為正關(guān)聯(lián)且JI值極大；第2生態(tài)種組為野核桃、千金榆、鵝耳櫪、椴樹(shù)、青榨槭、華榛，數(shù)量較多，為主要伴生樹(shù)種，正相關(guān)性較顯著，多為不顯著正關(guān)聯(lián)；第3生態(tài)種組為化香、漆樹(shù)，二者之間呈不顯著正關(guān)聯(lián)和正相關(guān)；第4生態(tài)種組為燈臺(tái)樹(shù)、梾木、大果榆、刺楸，為林下次要伴生樹(shù)種，相關(guān)性較顯著，多與其它樹(shù)種呈不顯著負(fù)關(guān)聯(lián)；第5生態(tài)種組為四照花、建始槭、青麩楊、木姜子，數(shù)量較少，隨機(jī)分布，獨(dú)立性較強(qiáng)。銳齒櫟林群落第1生態(tài)種組為銳齒櫟、油松、漆樹(shù)、華山松，為群落主要優(yōu)勢(shì)種，多為正關(guān)聯(lián)，JI值極大；第2生態(tài)種組為梾木、椴樹(shù)、野核桃，為主要伴生種，競(jìng)爭(zhēng)上層資源和空間，多為顯著正相關(guān)和不顯著正關(guān)聯(lián)且JI值較大；第3生態(tài)種組為華榛、千金榆、燈臺(tái)樹(shù)、青檀、刺葉櫟、木姜子、鵝耳櫪、烏桕、三椏烏藥、刺楸、建始槭，為中下層優(yōu)勢(shì)種群，彼此間多為極顯著正相關(guān)和不顯著正關(guān)聯(lián)且JI值較??；第4生態(tài)種組為鐵杉、化香樹(shù)，數(shù)量較少，獨(dú)立性較強(qiáng)。

4 結(jié)論與討論

4.1 種間聯(lián)結(jié)與群落穩(wěn)定性

秦嶺南坡松櫟林群落不同演替階段群落總體關(guān)聯(lián)性均為顯著正關(guān)聯(lián)，各自具有相對(duì)的穩(wěn)定性，且總體關(guān)聯(lián)性隨著演替進(jìn)程而加強(qiáng)，穩(wěn)定性也逐漸增強(qiáng)，這與長(zhǎng)白山云冷杉針闊混交林演替過(guò)程中總體聯(lián)結(jié)性隨演替發(fā)展的變化規(guī)律一致［30］。前人對(duì)秦嶺林區(qū)松櫟林種間關(guān)系表明，秦嶺林區(qū)油松天然林大多數(shù)種群種間聯(lián)結(jié)很弱，具有相對(duì)獨(dú)立分布的特征［9］，華山松天然次生林優(yōu)勢(shì)種群間聯(lián)結(jié)較松散，物種之間存在一定的獨(dú)立性，趨于形成一個(gè)較為穩(wěn)定的群落［11－28］，銳齒櫟林絕大多數(shù)種對(duì)間的獨(dú)立性相對(duì)較強(qiáng)，群落處于較為穩(wěn)定的次生林演替階段［12－13］，這些不同演替階段特定群落類型總體均表現(xiàn)出一定的松散性，但獨(dú)立程度較低，本研究結(jié)果與之較為一致。這種種間聯(lián)結(jié)的松散性可能主要決定于物種的生物學(xué)和生態(tài)學(xué)特性［4，7－8］，同時(shí)受群落演替的時(shí)間和空間［11］及其生態(tài)位分化的影響，此外，外界干擾和病蟲害等也會(huì)影響群落結(jié)構(gòu)的穩(wěn)定性。

4.2 種間聯(lián)結(jié)與群落演替的關(guān)系

一般來(lái)說(shuō)，隨著植被群落演替的進(jìn)展，群落結(jié)構(gòu)及種類組成將逐漸趨于完善和穩(wěn)定，種間關(guān)系也將同步趨于正關(guān)聯(lián)，以求得多物種間的穩(wěn)定共存［25，31］。秦嶺南坡松櫟林演替過(guò)程中主要植物種間聯(lián)結(jié)性隨演替進(jìn)展所表現(xiàn)出的特點(diǎn)也遵循這一規(guī)律，與油松林和松林混交林群落相比，銳齒櫟林群落物種間正聯(lián)結(jié)比例較高，這與物種對(duì)群落的適應(yīng)性和群落對(duì)物種的選擇和淘汰密切相關(guān)［11，30］，但檢驗(yàn)顯著率和JI＞0.75或0.25≥JI≥0的比例較低，趨于相對(duì)獨(dú)立分布，穩(wěn)定性明顯高于演替前期和中期。Spearman秩相關(guān)性研究結(jié)果進(jìn)一步證實(shí)油松林和松櫟混交林群落中多數(shù)物種間的關(guān)系不緊密，銳齒櫟群落種間關(guān)系則達(dá)到某種協(xié)調(diào)，生態(tài)位彼此分化表現(xiàn)出明顯的正相關(guān)特性，表明群落在長(zhǎng)期的協(xié)同進(jìn)化過(guò)程中種間關(guān)系趨于穩(wěn)定共存。前人對(duì)中國(guó)森林演替模式的研究表明，次生演替一般遵循針葉林到 針闊混交林再到闊葉林的規(guī)律［7，32，33］，就種間關(guān)系來(lái)說(shuō)，本研究結(jié)果支持這一演替規(guī)律，但僅從種間聯(lián)結(jié)和相關(guān)性來(lái)反映秦嶺南坡松櫟林群落演替動(dòng)態(tài)是不夠全面的，今后還需結(jié)合群落結(jié)構(gòu)、更新特征等進(jìn)一步分析考證。

對(duì)油松林→松櫟混交林→銳齒櫟林這一演替序列整個(gè)過(guò)程中的先鋒樹(shù)種、前期過(guò)渡種、后期過(guò)渡種和頂級(jí)種種間聯(lián)結(jié)性研究發(fā)現(xiàn)，種對(duì)之間的正負(fù)關(guān)聯(lián)比先增大后減小，近似呈單峰分布，即油松、華山松等先鋒樹(shù)種在群落中均占優(yōu)勢(shì)，且生態(tài)習(xí)性和生境需求相似，表現(xiàn)出對(duì)資源和空間的激烈競(jìng)爭(zhēng)，正負(fù)聯(lián)結(jié)種對(duì)分別占總種對(duì)數(shù)的33.3%和66.7%，尚未形成穩(wěn)定的群落關(guān)系，必然導(dǎo)致林木間出現(xiàn)自疏和他疏，形成林隙，林下充足的光照、水熱、養(yǎng)分和生存空間有利于青榨槭、鵝耳櫪、木姜子、千金榆、椴樹(shù)等闊葉樹(shù)種的侵入和定居，重新形成激烈的種間競(jìng)爭(zhēng)格局，經(jīng)過(guò)長(zhǎng)期的生態(tài)位分化，種間關(guān)系趨于松散，因此，這些過(guò)渡種間（包括前期過(guò)渡種間、后期過(guò)渡種間及其相互之間）多表現(xiàn)為有利的正聯(lián)結(jié)，正負(fù)聯(lián)結(jié)比均大于1，但不明顯，處于群落動(dòng)態(tài)演替過(guò)程中相對(duì)穩(wěn)定的階段；而適應(yīng)性強(qiáng)且生長(zhǎng)迅速的頂級(jí)樹(shù)種銳齒櫟、漆樹(shù)、野核桃等闊葉樹(shù)種則與后期過(guò)渡種間的負(fù)關(guān)聯(lián)比例達(dá)63.4%，經(jīng)過(guò)不斷的競(jìng)爭(zhēng)和淘汰逐漸占據(jù)上層空間，相同的空間分布和生長(zhǎng)資源競(jìng)爭(zhēng)使得這些樹(shù)種之間競(jìng)爭(zhēng)較為激烈，多形成顯著的負(fù)聯(lián)結(jié)，正負(fù)關(guān)聯(lián)比僅為0.29，而鵝耳櫪、千金榆等主要伴生樹(shù)種則在林分中下層占有一定的優(yōu)勢(shì)，與上層頂級(jí)種間多形成較為廣泛的不顯著正關(guān)聯(lián)，正是這些不同層次和不同生活型物種間對(duì)資源和空間需求的差異，使得物種組成及其種間關(guān)系處于動(dòng)態(tài)變化過(guò)程中，群落結(jié)構(gòu)也將不斷完善并趨于穩(wěn)定。馬映棟等［13］對(duì)銳齒櫟林種群空間分布格局和關(guān)聯(lián)性研究也證實(shí)了這一點(diǎn)，林冠上層銳齒櫟分布穩(wěn)定，主要為隨機(jī)分布，林間層多均勻分布，林下層聚集分布明顯，膀胱果（Staphyleaholocarpa）與銳齒櫟在幼樹(shù)期因相互競(jìng)爭(zhēng)而大部分尺度上互為負(fù)相關(guān)，其生長(zhǎng)發(fā)育到一定階段相互促生且小部分尺度上相互獨(dú)立；樹(shù)種銳齒櫟與青榨槭、膀胱果與青榨械在整體上保持無(wú)關(guān)聯(lián)，逐漸發(fā)展為小范圍內(nèi)互為相關(guān)性。群落中種間關(guān)系是復(fù)雜的，是多種因子耦合的結(jié)果［34］，僅僅依靠同一時(shí)空尺度下的調(diào)查資料和種對(duì)聯(lián)結(jié)統(tǒng)計(jì)結(jié)果分析，尚不能完全揭示松櫟林演替過(guò)程中的種間關(guān)系，更需要對(duì)演替過(guò)程中多物種的聯(lián)結(jié)關(guān)系進(jìn)行長(zhǎng)期定位研究，同時(shí)考慮環(huán)境因子和人為因素的影響，才能更好揭示群落種間相互關(guān)系和演替趨勢(shì)。

4.3 種間聯(lián)結(jié)與物種多樣性關(guān)系

物種多樣性是反映植物群落組成結(jié)構(gòu)和穩(wěn)定性的重要指標(biāo)，種間關(guān)系是群落形成、演化和穩(wěn)定性的重要數(shù)量特征［35］，二者之間存在一定的聯(lián)系。對(duì)油松林→松櫟混交林→銳齒櫟林這一演替序列整個(gè)過(guò)程中的先鋒樹(shù)種、前期過(guò)渡種、后期過(guò)渡種和頂級(jí)種種間聯(lián)結(jié)性研究發(fā)現(xiàn)，種對(duì)之間的正負(fù)關(guān)聯(lián)比先增大后減小，近似呈單峰分布，也即先鋒種定居后，青榨槭、鵝耳櫪等闊葉樹(shù)種的侵入促使種間競(jìng)爭(zhēng)更加激烈，而過(guò)渡種間關(guān)系趨于松散，頂級(jí)樹(shù)種銳齒櫟、漆樹(shù)、野核桃等競(jìng)爭(zhēng)也較為激烈。而秦嶺南坡松櫟林群落演替過(guò)程中雖然群落間物種豐富度無(wú)顯著差異，但群落總體物種多樣性隨演替進(jìn)展總體上呈現(xiàn)明顯的單峰型變化，在松櫟林階段物種多樣性達(dá)到最大，與中期物種多樣性假說(shuō)一致［16］。這與沙棘人工林群落物種多樣性及種間關(guān)聯(lián)變化特征較為一致，沙棘人工林群落各層及總體物種多樣性指數(shù)基本呈“S”形波動(dòng)變化的趨勢(shì)，種間關(guān)聯(lián)正負(fù)關(guān)聯(lián)比呈拋物線形變化，說(shuō)明群落發(fā)育到13a前后呈現(xiàn)相對(duì)穩(wěn)定的態(tài)勢(shì)［35］。這可能是因?yàn)槿郝浣?jīng)過(guò)闊葉樹(shù)種入侵過(guò)程的激烈競(jìng)爭(zhēng)與淘汰，種間關(guān)系則達(dá)到某種協(xié)調(diào)，生態(tài)位彼此分化，群落組成和結(jié)構(gòu)達(dá)到了穩(wěn)定共存的階段。

4.4 種間關(guān)系與種群生態(tài)位關(guān)系

種間聯(lián)結(jié)與其生態(tài)位重疊之間是緊密聯(lián)系的。油松林中銳齒櫟－油松、漆樹(shù)－華山松、銳齒櫟－漆樹(shù)、鹽膚木－栓皮櫟、銳齒櫟－鵝耳櫪χ2檢驗(yàn)值較大，呈顯著正關(guān)聯(lián)或一定的正關(guān)聯(lián)，其生態(tài)位重疊值相對(duì)較大［9，24］，分 別 為0.589、0.561、0.231、0.189、0.176，而鵝耳櫪－梾木、栓皮櫟－鵝耳櫪呈顯著負(fù)關(guān)聯(lián)，其生態(tài)位重疊值［9］僅為0.033、0。銳齒櫟林中銳齒櫟－鵝耳櫪、華山松－油松、華山松－鵝耳櫪、漆樹(shù)－油松、漆樹(shù)－三椏烏藥、油松－鵝耳櫪、三椏烏藥－鵝耳櫪等呈一定程度的正關(guān)聯(lián)，其生態(tài)位重疊值［36］分別為0.881、0.892、0.956、1.00、0.971、1.00、0.836；銳齒櫟－椴樹(shù)、銳齒櫟－華山松、銳齒櫟－漆樹(shù)、鵝耳櫪－漆樹(shù)、椴樹(shù)－華山松等呈一定的負(fù)關(guān)聯(lián)，其生態(tài)位重疊值［37］僅為0.441、0.502、0.301、0.467、0.317，由此可見(jiàn)，優(yōu)勢(shì)種群間的正聯(lián)結(jié)性越強(qiáng)，其生態(tài)位重疊值就越大；反之，種間的負(fù)聯(lián)結(jié)性越強(qiáng)，其生態(tài)位重疊值就越小，這與前人的研究結(jié)果一致［38，39］。

4.5 種間聯(lián)結(jié)與CCA 排序關(guān)系

通過(guò)對(duì)比種間聯(lián)結(jié)性及其種對(duì)在物種與環(huán)境因子CCA 排序圖［15］中的分布可以發(fā)現(xiàn)，銳齒櫟林中的銳齒櫟與鵝耳櫪在CCA 排序圖中的距離較近，而與椴樹(shù)、華山松、漆樹(shù)相距較遠(yuǎn)；鵝耳櫪與銳齒櫟相距較近，而與漆樹(shù)相距較遠(yuǎn)；漆樹(shù)與油松相距較近，而與銳齒櫟、鵝耳櫪相距較遠(yuǎn)。松櫟混交林中的漆樹(shù)與華山松、鵝耳櫪相距較近，而與油松、鐵杉相距較遠(yuǎn)。由此可見(jiàn)，這些種對(duì)在物種與環(huán)境因子CCA 排序圖中分布距離遠(yuǎn)近很大程度上取決于其正負(fù)關(guān)聯(lián)性，也就是說(shuō)，對(duì)于同一物種，若與其它物種間存在較為顯著的正聯(lián)結(jié)性，則正聯(lián)結(jié)性越顯著，物種在CCA 排序圖上的距離越近；若與其它物種間存在較為顯著的負(fù)聯(lián)結(jié)性，則負(fù)聯(lián)結(jié)性越顯著，物種在CCA 排序圖上的距離越遠(yuǎn)。這與鄭超超等［39］對(duì)浙江江山公益林種間關(guān)系的研究結(jié)果較為一致。這可能是因?yàn)檎P(guān)聯(lián)種對(duì)能夠充分利用周圍資源空間，生態(tài)位重疊程度較高，導(dǎo)致其在CCA 排序圖中集中分布，而負(fù)關(guān)聯(lián)種對(duì)偏好于某些環(huán)境因子，整體上對(duì)環(huán)境資源的利用能力較差，因而分散分布于CCA 排序圖中。此外，油松林關(guān)聯(lián)種對(duì)間分布距離與其他2種群落類型有較大差異，這可能是因?yàn)橛退闪值姆植寂c銳齒櫟林和松林混交林在CCA 排序圖中的分布相距較遠(yuǎn)，即其群落在結(jié)構(gòu)組成和生境等方面差異較大［15］。CCA 排序同時(shí)也表明，坡度是影響喬木層物種分布的主要環(huán)境因子，其次是海拔、坡度等［15］，正聯(lián)結(jié)的物種生境特征比較一致，負(fù)聯(lián)結(jié)的物種生境特征剛好相反。物種之間的關(guān)系也不完全受環(huán)境條件制約，還需進(jìn)一步結(jié)合生態(tài)位理論、CCA 協(xié)變量矩陣偏典范對(duì)應(yīng)分析等探討。

4.6 種間聯(lián)結(jié)對(duì)林分經(jīng)營(yíng)的指導(dǎo)

秦嶺南坡油松林、銳齒櫟林及其混交林穩(wěn)定性程度直接關(guān)系到秦嶺生物多樣性維持、水碳平衡、氣候調(diào)節(jié)、水土保持以及水源涵養(yǎng)等生態(tài)服務(wù)功能的發(fā)揮。闡釋植物種群種間關(guān)系旨在揭示種對(duì)間競(jìng)爭(zhēng)的結(jié)果或群落的現(xiàn)狀，這有助于對(duì)地帶性植被進(jìn)行改造和恢復(fù)，有效縮短演替進(jìn)程。油松林中銳齒櫟－漆樹(shù)、鹽膚木－梾木，松櫟林中的銳齒櫟－油松、銳齒櫟－漆樹(shù)等呈顯著正關(guān)聯(lián)的種對(duì)應(yīng)作為生態(tài)目標(biāo)樹(shù)或經(jīng)營(yíng)目標(biāo)樹(shù)予以保留；相反，油松林中的油松－三椏烏藥、漆樹(shù)－冬瓜楊，松櫟混交林中的燈臺(tái)樹(shù)－漆樹(shù)、銳齒櫟－化香樹(shù)等種對(duì)呈不顯著負(fù)關(guān)聯(lián)且JI值較大或顯著負(fù)關(guān)聯(lián)的種對(duì)應(yīng)選擇性的伐除；此外，還可依據(jù)種間關(guān)系的變動(dòng)趨勢(shì)確定撫育間伐時(shí)植株的去留［30］和采伐階段，例如華山松－梾木種間聯(lián)結(jié)性演替中期呈不顯著正聯(lián)結(jié)、前期和后期均為不顯著負(fù)聯(lián)結(jié)，應(yīng)選擇松櫟混交林向銳齒櫟林過(guò)渡階段調(diào)整華山松和梾木的種間關(guān)系；生態(tài)種組的劃分也為篩選群落內(nèi)環(huán)境監(jiān)測(cè)的指示種［28］，滿足林區(qū)營(yíng)林和重建需要提供了保障，通過(guò)這些有效的經(jīng)營(yíng)措施促進(jìn)林分向健康穩(wěn)定的地帶性頂級(jí)群落發(fā)展，縮短秦嶺山地松櫟林演替進(jìn)程。

［1］ ARMAS C，PUGNAIRE F I.Plant interactions govern population dynamics in semi－arid plant community［J］.JournalofEcology，2005，93（5）：978－989.

［2］ LIU J H（劉玨宏），GAO H（高 慧），etal.Comparative analysis of interspecific association within theStipagrandis－S.kryloviicommunity in typical steppe of Inner Mongolia，China［J］.ChineseJournalofPlantEcology（植物生態(tài)學(xué)報(bào)），2010，34（9）：1 016－1 024（in Chinese）.

［3］ GUO Y X（郭垚鑫），HU Y N（胡有寧），KONG L T（孔令童），etal.Interspecific association and correlation of main plant species in subalpine meadow of Zibai Mountain［J］.ActaEcologicaSinica（生態(tài)學(xué)雜志）2011，30（8）：1 775－1 780（in Chinese）.

［4］ SU S J，LIU J F，HE Z S，etal.Ecological species groups and interspecific association of dominant tree species in Daiyun Mountain National Nature Reserve［J］.JournalofMountainScience，2015，12（3）：637－646.

［5］ LI Y D，XU H，etal.Division of ecological species groups and functional groups based on interspecific association——a case study of the tree layer in the tropical lowland rainforest of Jianfenling in Hainan Island，China［J］.FrontiersofForestryinChina，2008，3（4）：407－715.

［6］ GREIG－SMITH P.Quantitative Plant Ecology（3rd edition）［M］.BlackwellScientificPublications，1983.

［7］ PENG SH L（彭少麟），ZHOU H CH（周厚誠(chéng)），GUO SH C（郭少聰），etal.Studies on the changes in interspecific association of zonal vegetation in Dinhushan［J］.ActabotanicaSinica（植物學(xué)報(bào)），1999，41（11）：1 239－1 244（in Chinese）.

［8］ WANG B S（王伯蓀），PENG SH L（彭少麟）.Study on the measuring techniques of interspecific association of lower－subtropical evergreenbroadleaved forests I.The exploration and the revision on the measuring formulas of interspecific association［J］.ActaPhytoecologicaet GeobotanicaSinica（植物生態(tài)學(xué)與地植物學(xué)叢刊），1985，9（4）：274－285（in Chinese）.

［9］ LIN Y Y（藺雨陽(yáng)），WANG D X（王得祥），LEI R D（雷瑞德），etal.Classification of ecological twin species and plant functional type for naturalPinustabulaeformisForests in middle part of Qinling Mountains［J］.JournalofNortheastForestryUniversity（東北林業(yè)大學(xué)學(xué)報(bào)），2009，37（11）：46－51（in Chinese）.

［10］ XU J SH（許金石），CHEN Y（陳 煜），WANG G X（王國(guó)勛），etal.Interspecies association of dominant woody plant communities in Qinling forest region，Shaanxi［J］.ActaBot.Boreal.－Occident.Sin.（西北植物學(xué)報(bào)），2014，34（7）：1 467－1 475（in Chinese）.

［11］ ZHANG M X（張明霞），WANG D X（王得祥），KANG B（康 冰），etal.Interspecific associations of dominant plant populations in secondary forest ofPinusarmandiiin Qinling Mountains［J］.ScientiaSilvaeSinicae（林業(yè)科學(xué)），2015，51（1）：12－21（in Chinese）.

［12］ YANG H Y（楊海裕），ZHANG S ZH（張宋智），etal.Interspecific associations among tree species in natural deciduous broadleaved forests in western Qinling，Northwest China［J］.ChineseJournalofEcology（生態(tài)學(xué)雜志），2012，31（10）：2 513－2 520（in Chinese）.

［13］ MA Y D（馬映棟），LIU W ZH（劉文楨），etal.Change analysis of spatial pattern and correlation forQuercusalienavar.acuteserrataforest population in Xiaolongshan Mountains［J］.ActaBot.Boreal.－Occident.Sin.（西北植物學(xué)報(bào)），2014，34（9）：1 878－1 886（in Chinese）.

［14］ LIU Y Y，LI F R，JIN G Z.Spatial patterns and associations of four species in an old－growth temperate forest［J］.JournalofPlantInteractions，2014，9（1）：745－753.

［15］ ZHANG Z L（張子良），YU F（于 飛），WANG D X（王得祥）.Effects of environmental factors on the species composition of three typical secondary forests in Qinling Mountains［J］.JournalofNorthwestA＆FUniversity（Nat.Sci.Ed.）（西北農(nóng)林科技大學(xué)學(xué)報(bào)·自然科學(xué)版），2014，42（8）：69－78（in Chinese）.

［16］ ZHANG H W（張洪武），YU F（于 飛），WANG D X（王得祥），etal.Changes of species composition and diversity in the process of community succession of pine oak forest on the south－facing slopes in Qinling Mountains［J］.ActaBot.Boreal.－Occident.Sin.（西北植物學(xué)報(bào)），2014，34（1）：169－176（in Chinese）.

［17］ BARBéRO M，LOISEL R，QUéZEL P，etal.Pines of the Mediterranean Basin［M］／／RICHARDSON D M.Ecology and Biogeography of Pinus，Cambridge University Press，Cambridge.1998：153－170.

［18］ GRACIA M，RETANA J，ROIG P.Mid－term successional patterns after fire of mixed pine－oak forests in NE Spain［J］.ActaOeclolgica，2002，（23）：105－411.

［19］ KANG B（康 冰），LIU SH R（劉世榮），WANG D X（王得祥），etal.Regeneration characteristics of woody plant seedlings in typical secondary forests in Qinling Mountains［J］.ChineseJournalofAppliedEcology（應(yīng)用生態(tài)學(xué)報(bào)），2011，22（12）：3 123－3 130（in Chinese）.

［20］ YU F，WANG D X，SHI X X，etal.Seed dispersal by small rodents favors oak over pine regeneration in the pine－oak forests of the Qinling Mountain，China［J］.ScandinavianJournalofForestResearch，2013，28（6）：540－549.

［21］ YU F，WANG D X，SHI X X，etal.Effects of environmental factors on tree seedling regeneration in a pine－oak mixed forest in the Qinling Mountains，China［J］.JournalofMountainScience，2013，10（5）：845－853.

［22］ YU F，SHI X X，WANG D X，etal.Effects of insect infestation onQuercusalienavar.acuteserrataacorn dispersal in the Qinling Mountains，China［J］.NewForests，2015，（46）：51－61.

［23］ YU F，SHI X X，WANG D X，etal.Seed predation patterns favor the regeneration of dominant species in forest gaps compared with the understory in an oak－pine mixed forest［J］.ActaTheriol.，2014，（59）：495－502.

［24］ CHAI Z ZH（柴宗政），WANG D X（王得祥），etal.Niche characteristics of main plant populations in naturalPinustabulaeformiscommunities in Qinling Mountains，Northwest China［J］.ChineseJournalofEcology（生態(tài)學(xué)雜志），2012，31（8）：1 917－1 923.

［25］ DU D L（杜道林），LIU Y CH（劉玉成），LI R（李 睿）.Studies on the interspecific association of dominant species in a subtropicalCatanopsisfargesiiforest of Jinyun Mountain，China［J］.ActaPhytoecologicaSinica（植物生態(tài)學(xué)報(bào)），1995，19（2）：149－157（in Chinese）.

［26］ LEI R D（雷瑞德），PENG H（彭 鴻），CHEN C G（陳存根）.Types and phytoenosis of natural secondary forests at Huoditang Forest Region［J］.JournalofNorthwestForestryCollege（西北林學(xué)院學(xué)報(bào)），1996，11（S1）：43－52（in Chinese）.

［27］ YANG Q（楊 琦），ZHOU J（周 婧），TAO CH（陶 楚），etal.Comparative study on interspecific association in tropical evergreenmonsoonelfinforest at two successional stages in Tongguling of Hainan Island［J］.ChineseAgriculturalScienceBulletin（中國(guó)農(nóng)學(xué)通報(bào)），2014，30（22）：8－15（in Chinese）.

［28］ JU T ZH（巨天珍），KANG L D（康麗丹），WANG L F（王立峰），etal.Interspecific association analysis of dominant species in shrub layer ofPinusarmandiicommunity in Xiaolongshan Mountain of Gansu Province［J］.ActaBot.Boreal.－Occident.Sin.（西 北 植 物 學(xué) 報(bào)），2010，30（6）：829－834（in Chinese）.

［29］ WANG L（王 琳），ZHANG J T（張金屯）.Interspecific association and correlation of dominant species of Lishan Mountain meadow in Shanxi Province［J］.ActaBot.Boreal.－Occident.Sin.（西北植物學(xué)報(bào)），2004，24（8）：1 435－1 440（in Chinese）.

［30］ GONG ZH W（龔直文），KANG X G（亢新剛），GU L（顧 麗），etal.Interspecific association among arbor species in two succession stages of spruce－fir conifer and broadleaved mixed forest in Changbai Mountains，northeastern China［J］.JournalofBeijingForestryUniversity（北京林業(yè)大學(xué)學(xué)報(bào)），2011，33（5）：29－33（in Chinese）.

［31］ ZHOU X Y（周先葉），WANG B S（王伯蓀），etal.An analysis of interspecific associations in secondary succession forest communities in Heishiding Natural Reserve，Guangdong Province［J］.ActaPhytoecologicaSinica（植物生態(tài)學(xué)報(bào)），2000，24（3）：332－339（in Chinese）.

［32］ LEI N F（雷濘菲），SU ZH X（蘇智先），etal.A comparative study on life－form spectra of evergreen broad－leaved forest in different successions in Jinyun Mountain［J］.ChineseJournalofAppliedEcology（應(yīng)用生態(tài)學(xué)報(bào)），2002，13（3）：267－270（in Chinese）.

［33］ PETRITAN I C，MARZANO R，PETRITAN A M，etal.Overstory succession in a mixedQuercuspetraea－Fagussylvaticaold growth forest revealed through the spatial pattern of competition and mortality［J］.ForestEcologyandManagement，2014，326：9－17.

［34］ TAN Y B（譚一波），ZHAN CH A（詹潮安），YANG H D（楊海東），etal.Inter－specific associations among main tree species inMachilus chinensiscommunities in Nan’ao Island，Guangdong Province［J］.JournalofCentralSouthUniversityofForestry＆Technology（中南林業(yè)科技大學(xué)學(xué)報(bào)），2012，32（11）：92－99（in Chinese）.

［35］ GUO L J（郭連金），ZHANG W H（張文輝），LIU G B（劉國(guó)彬）.Species diversity and interspecific association in development sequence ofHippophaerhamnodiesplantations in loess hilly region［J］.ChineseJournalofAppliedEcology（應(yīng)用生態(tài)學(xué)報(bào)），2007，18（1）：9－15（in Chinese）.

［36］ ZHAO Y H（趙永華），LEI R D（雷瑞德），HE X Y（何興元），etal.Niche characteristics of plant populations inQuercusalienavar.acuteserratastands in Qinling Mountains［J］.ChineseJournalofAppliedEcology（應(yīng)用生態(tài)學(xué)報(bào)），2004，15（6）：913－918（in Chinese）.

［37］ CAO CH（曹 春），JU T ZH（巨天珍），etal.Study on niche of main population ofQuercusalienavar.acuteserrataforest on to the Mt.Xiaolongshan in Gansu Province［J］.JournalofGansuAgriculturalUniversity（甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào)），2006，41（3）：65－69（in Chinese）.

［38］ SHI Z M（史作民），CHENG R M（程瑞梅），LIU SH R（劉世榮）.Niche characteristics of plant populations in deciduous broad－leaved forest in Baotianman［J］.ChineseJournalofAppliedEcology（應(yīng)用生態(tài)學(xué)報(bào)），1999，10（3）：265－269（in Chinese）.

［39］ ZHENG CH CH（鄭超超），YI L T（伊力塔），etal.Interspecific relationship and canonical correspondence analysis of the dominant species in ecological service forest of Jiangshan City in Zhejiang Procince［J］.ActaEcologicaSinica（生態(tài)學(xué)報(bào)），2015，35（22）：1－11（in Chinese）.