小隴山國家自然保護區(qū)毛栗群落生態(tài)位及空間分布格局研究

王振乾++巨天珍++劉洋洋+傅晨玲++任藝彬++高鳳如+張宋智

摘要：在野外調(diào)查的基礎(chǔ)上，對小隴山國家自然保護區(qū)內(nèi)毛栗（Castanea mollissima Bl.）群落中主要種群的Shannon-Wiener多樣性指數(shù)及Levins生態(tài)位寬度、生態(tài)位重疊及空間分布格局各項指標(biāo)進行了分析。結(jié)果表明，毛栗的生態(tài)位寬度最大，為0.81，楓楊（Pterocarya stenoptera C.DC.）的生態(tài)位寬度值最小，為0.50；毛栗與其他種的生態(tài)位重疊值在0.33～0.72。毛栗的生態(tài)位重疊值比較大，群落中由于利用相同資源位而產(chǎn)生的競爭相對較大。6個樣地中毛栗群落的分布格局均表現(xiàn)為聚集分布，樣地3中的毛栗群落聚集強度最低，樣地4中的毛栗群落聚集強度最高。

關(guān)鍵詞：毛栗（Castanea mollissima Bl.）；群落；生態(tài)位；空間分布格局；小隴山國家自然保護區(qū)

中圖分類號：S664.2+54.93 文獻標(biāo)識碼：A 文章編號：0439-8114（2017）10-1878-04

DOI：10.14088/j.cnki.issn0439-8114.2017.10.020

Study on Niche and Spatial Distribution Pattern of Castanea mollissima Community in Xiaolongshan National Nature Reserve

WANG Zhen-qian1，JU Tian-zhen1，LIU Yang-yang1，F(xiàn)U Chen-ling1，REN Yi-bin1，GAO Feng-ru1，ZHANG Song-zhi2

（1. College of Geography and Environment， Northwest Normal University， Lanzhou 730070， China；

2. Xiaolongshan Forestry Science Institute， Tianshui 741000， Gansu， China）

Abstract： On the basis of field investigation， the index of Shannon-Wiener diversity indicators and Levins niche breadth and niche overlap and spatial distribution pattern of the dominant populations in Castanea mollissima Bl. community in Xiaolongshan National Nature Reserve were analyzed. The results showed that the widest niche width of C. mollissima was 0.81， while the niche breadth of Pterocarya stenoptera C.DC. was the smallest（0.50）； The ecological niche overlap of C. mollissima and other species was 0.33 to 0.72. The niche overlap value of C. mollissima community was relatively high， and the competition is relatively large due to the use of same resources. The six samples of C. mollissima communities all showed aggregated distribution pattern； the aggregation intensity of sample 3 C. mollissima community was the lowest， while of sample 4 was the highest.

Key words： Castanea mollissima Bl.； community； niche； spatial distribution pattern； Xiaolongshan National Nature Reserve

自“生態(tài)位”提出之后，國內(nèi)外學(xué)者們紛紛展開對生態(tài)位的研究，國內(nèi)從20世紀(jì)80年代開始陸續(xù)進行了有關(guān)生態(tài)位方面的理論探索[1]。在此之后，又有不少學(xué)者提出了不同的生態(tài)位定義，但目前在應(yīng)用生態(tài)位理論范疇內(nèi)，多認(rèn)為物種在多維環(huán)境梯度空間的位置就是生態(tài)位，這在森林資源保護與利用、群落演替、群落結(jié)構(gòu)、生物多樣性、物種進化及瀕危物種評價等方面有著廣泛的應(yīng)用前景[2]。研究種群的核心內(nèi)容就是種群的分布格局及年齡結(jié)構(gòu)，研究種群分布格局的目的，一方面是對種群的水平結(jié)構(gòu)進行定量描述，另一方面可以揭示格局的成因，闡述種群及群落的動態(tài)變化[3]。毛栗（Castanea mollissima Bl.）是殼斗科（Fagaceae）栗屬（Castanea Mill.）植物，原產(chǎn)于中國，多生長于海拔370～2 800 m的山地；現(xiàn)已由人工廣泛栽培，主要分布在華北地區(qū)的燕山、太行山及鄰近地區(qū)，也包括北京、河南、山東、陜西、甘肅、江蘇?。ㄊ校┑鹊亍８拭C省的小隴山是國家級自然保護區(qū)，毛栗也廣泛分布于此。關(guān)于小隴山整個林區(qū)植物優(yōu)勢種群生態(tài)位和空間分布格局的研究比較多，如巨天珍等[4，5]對該保護區(qū)的油松（Pinus tabuliformis Carrière）、華山松（P. armandii Franch.）種群的生態(tài)位進行了研究，曹春等[6]、馬映棟等[7]對銳齒櫟（Quercus aliena var. acuteserrata Maxim.）種群的生態(tài)位和種群空間格局及關(guān)聯(lián)性變化進行過分析；高洋等[8]對白皮松（P. bungeana Zucc.）種群的生態(tài)位進行過比較；劉小林等[9]、郭垚鑫等[10]對小隴山紅樺（Betula albosinensis Burk.）進行過種群分布格局的調(diào)查，鄭智等[11]驗證了小隴山維管植物豐富度和種域海拔梯度格局及Rapoport法則。而關(guān)于小隴山毛栗種群生態(tài)位及分布格局兩大特征相結(jié)合的研究未曾報道，鑒于此，課題組以前人探索工作為基礎(chǔ)，對小隴山國家自然保護區(qū)毛栗種群生態(tài)位及分布格局做了進一步探討，這不僅有助于保護區(qū)管理人員了解毛栗種群數(shù)量動態(tài)及發(fā)展趨勢，同時對毛栗種群的可持續(xù)經(jīng)營、撫育及保護等都具有指導(dǎo)作用。

1 研究區(qū)概況

小隴山國家自然保護區(qū)地處北緯34°0′-34°40′，東經(jīng)105°30′-106°30′。屬秦嶺西部山地的一部分，區(qū)內(nèi)地形復(fù)雜，相對高差約500～1 000 m。由于處在溫帶南緣與北亞熱帶的過渡地帶，主要屬于暖溫帶-中帶半濕潤大陸性季風(fēng)類型，氣候比較溫暖濕潤，大多數(shù)地域的極端最高氣溫為39.2 ℃，極端最低氣溫-23.2 ℃，年平均氣溫約7～12 ℃，>10 ℃年有效積溫在2 444～3 825 ℃。該區(qū)域降水適中，主要集中于7-9月，年降水量在460～800 mm，春、秋雨占年降雨量的比例大，約 70%～80%，且降水年際間變化大，年蒸發(fā)量989～1 658 mm，相對濕度為68%～78%。區(qū)內(nèi)的地帶性土壤在秦嶺以南為黃褐土，以北為灰褐土。植被種類眾多，有苔蘚、蕨類、裸子、被子植物224科945屬2 700多種，其中草本植物1 900多種，木本植物800多種（不包括引進栽培種），溫帶植物占優(yōu)勢，熱帶、亞熱帶植物比較少[12]。毛栗群落里生長有毛栗、銳齒櫟、華山松、楓楊（Pterocarya stenoptera C.DC.）、油松、栓皮櫟（Quercus variabilis Bl.）等喬木植物。

2 研究方法

2.1 樣地調(diào)查

在小隴山保護區(qū)內(nèi)的黨川林場、百花林場等四大林場內(nèi)，分別設(shè)置6個面積為20 m×20 m的樣地，每個樣地內(nèi)選取5個4 m×4 m的小樣方進行植被調(diào)查，測定每塊樣地的海拔、坡向、坡位、坡度、土壤條件和群落類型等因子，采用每木調(diào)查法記錄樣地內(nèi)所有個體的種名、樹高、胸徑、冠幅、枝下高（起測徑階≥4 cm）和各樹種的密度、蓋度等。

2.2 生態(tài)位指標(biāo)

2.2.1 生態(tài)位寬度 生態(tài)位寬度采用Levins指數(shù)，公式為：

式（1）中，B（L）i為物種i的Levins生態(tài)位寬度指數(shù)；Pij為物種i對第j資源梯度級的利用占它對全部資源利用的百分率，Pij=nij/Ni，而Ni=∑nij，nij為物種i在資源梯度級j的數(shù)量特征值，如蓋度、重要值、密度等，本試驗中為種i在第j樣方的重要值； r為樣方數(shù)。此外，生態(tài)位寬度也可采用Shannon-Wiener多樣性指數(shù)計測，公式為：

式（2）中，B（sw）i為物種i的生態(tài)位寬度指數(shù)；Pij是物種i對第j資源梯度級的利用占它利用全部資源的百分率；s為種群數(shù)；r為資源位數(shù)。其中：

式（3）中，nij為物種i在第j資源位的重要值，Yi為物種i所在利用全部資源位的重要值之和，生態(tài)位寬度B（sw）i值域[0，1]。

2.2.2 生態(tài)位重疊 生態(tài)位重疊是指一定資源序列上，2個物種利用同等級資源而相互重疊的情況，可用對稱a法（Pianka公式）測算，公式為：

式（4）中，NOik為生態(tài)位重疊值，Pij和Pkj為種i和種k在資源j上的優(yōu)勢度（為物種重要值）。

2.3 徑級劃分

由于喬木種群樹齡比較難以測定，通常采用徑級結(jié)構(gòu)代替年齡結(jié)構(gòu)來進行研究[13-15]。試驗在已有的統(tǒng)計數(shù)據(jù)基礎(chǔ)之上，根據(jù)毛栗個體胸徑的大小分為5個等級，具體為：幼苗（Ⅰ級），胸徑<5.0 cm；幼樹（Ⅱ級），5.0 cm≤胸徑<11.0 cm；小樹（Ⅲ級），11.0 cm≤胸徑<17.0 cm；中樹（Ⅳ級），17.0 cm≤胸徑<23.0 cm；大樹（Ⅴ級），胸徑≥23.0 cm。對于齡級個數(shù)較少或缺失的植株不進行統(tǒng)計。

2.4 空間分布格局分析方法

種群分布格局研究主要選取以下基礎(chǔ)指標(biāo)對毛栗種群的各個聚集強度及分布格局指數(shù)進行測算[16，17]，從而綜合分析毛栗種群的空間分布格局。

2.4.1 方差、均值比 此乃基礎(chǔ)指標(biāo)，公式為：

C=S2/X （5）

式（5）中，S2為種群多度的方差，也是種群多度均值。當(dāng)C<1時，為均勻分布；C=1時，為隨機分布；C>1時，為聚集分布。由于實測與預(yù)期結(jié)果存在一定的偏離，因此用t檢驗來進行確定。t檢驗以為標(biāo)準(zhǔn)差，由此t=，N為所采用的樣方數(shù)。

2.4.2 負二項分布參數(shù)（k） 該指標(biāo)用于判別種群中植株的聚集程度，公式為：

k=X2/（S2-X） （6）

式（6）中，當(dāng)k<0時，為均勻分布；k>0時，為聚集分布，k值越接近于0，則聚集度越大；當(dāng)k>8時，種群接近泊松分布，即隨機分布。

2.4.3 聚集度指數(shù)（I） 指植物在空間分布的集中程度，公式為：

式（7）中，當(dāng)I<0時，為均勻分布；I=0時，為隨機分布；I>0時，為聚集分布。

2.4.4 格林指數(shù)（GI） 公式為：

式（8）中，GI在0（隨機）和1（最大集聚）之間變化，當(dāng)GI<0時，為均勻分布；GI=0時，為隨機分布；GI>0時，為聚集分布[18]。

2.4.5 Cassie指數(shù) 公式為：

式（9）中，當(dāng)Ca<0時，為均勻分布；Ca=0時，為隨機分布；Ca>0時，為聚集分布。

2.4.6 聚塊性指數(shù)（PAI） 該指標(biāo)不受空樣地的影響，公式為：

PAI=m*/X

m*=X+[S2/X-1]=X+1 （10）

式（10）中，當(dāng)PAI<1時，為均勻分布；PAI=1時，為隨機分布；PAI>1時，為聚集分布。

3 結(jié)果與分析

3.1 生態(tài)位寬度分析

毛栗群落主要種群的生態(tài)位寬度測定情況見表1，分析比較見圖1。從表1、圖1可見，毛栗群落主要種群生態(tài)位寬度在采用不同公式計算后，結(jié)果基本相似。B（sw）i、B（L）i按照數(shù)值的大小依次為毛栗、銳齒櫟、華山松、楓楊、油松、栓皮櫟，說明毛栗生態(tài)位寬度最寬，它在每個資源位中都有出現(xiàn)，表明毛栗在保護區(qū)里分布廣、數(shù)量多、利用資源也較為充分；其次是銳齒櫟，該種群雖然在保護區(qū)分布也廣，但在建群種林中依然是處于次優(yōu)勢種的地位；其余種群也是毛栗的伴生種群，但生態(tài)位寬度值相對較窄，比如華山松、楓楊、油松、栓皮櫟，說明這幾種植物對資源的利用具有較大的局限性，并且對生態(tài)因子有一定的依賴，或者也有可能是適合毛栗的生境并不一定適合這些種群的生長。因此，各個種群的生態(tài)位寬度排序表明了它們在毛栗群落中的地位與作用的大小。

3.2 生態(tài)位重疊分析

毛栗群落各種群的生態(tài)位重疊值測定結(jié)果見表2。從表2可知，華山松-油松的生態(tài)位重疊值最大，為0.84，說明它們利用生境的模式最相似。華山松-栓皮櫟、油松-栓皮櫟、毛栗-銳齒櫟、銳齒櫟-楓楊之間也具有較大的生態(tài)位重疊，重疊值分別為0.78、0.76、0.72、0.72。由于六者均為喬木樹種，所以它們在同一群落中會產(chǎn)生一定程度的競爭。同時，從表2還可看出，生態(tài)位寬度越寬的物種，與其他種的生態(tài)位重疊不一定大。說明生態(tài)位寬度只表示一個種群在資源軸上的空間占據(jù)能力、而生態(tài)位重疊表示各對植物種群之間的競爭或者共存關(guān)系[19]。重疊值越大說明其生態(tài)學(xué)特性越相似，利用資源的能力比較相近；或者說其在同一空間內(nèi)對生態(tài)的需求互補性越強。

3.3 毛栗種群空間格局分析

通過上述指標(biāo)對小隴山毛栗種群的6個樣地個體數(shù)進行測算，結(jié)果見表3，從表3可見，各個樣地毛栗種群的分布格局參數(shù)C均大于0，且通過了t檢驗，顯著性較高。I、GI、Ca也均大于0，聚塊性指數(shù)PAI均明顯大于1，聚集強度參數(shù)k值也大于0。以上分析表明，各樣地毛栗種群的分布格局均表現(xiàn)為聚集分布，其中樣地3的k值最大，表明該樣地的毛栗種群聚集強度最低；其次為2號樣地。4號樣地的k值最小，說明樣地4的毛栗種群聚集強度最高。這說明由于種群生存策略及適應(yīng)機制的不同，會引起其聚集強度也會不同[20]。

4 小結(jié)

物種生態(tài)位寬度表現(xiàn)的是它們利用環(huán)境資源能力的程度，也表征了物種的生態(tài)適應(yīng)能力和分布幅度。從小隴山國家自然保護區(qū)的喬木層種群的生態(tài)位寬度來看，毛栗、銳齒櫟和華山松種群的生態(tài)位寬度數(shù)值遠大于其他種群，為保護區(qū)主要的森林群落建群種，其對地域性的環(huán)境資源利用具有明顯優(yōu)勢。但由于三者生物學(xué)及生態(tài)學(xué)特性存在差異，且經(jīng)過長期環(huán)境的選擇與生存競爭，形成了各自明顯的環(huán)境生態(tài)位特征。

生態(tài)位重疊體現(xiàn)了物種對相同等級資源的利用程度以及空間配置關(guān)系，反映了物種之間對資源利用的相似程度和競爭關(guān)系。小隴山國家自然保護區(qū)的華山松-栓皮櫟、油松-栓皮櫟、毛栗-銳齒櫟、銳齒櫟-楓楊是生態(tài)位重疊值相對較大（0.78、0.76、0.72、0.72）的種對，它們的生態(tài)學(xué)特性與適應(yīng)方式相似，生態(tài)位分化的程度相對低，它們必然對有限的資源進行激烈的競爭。

小隴山國家自然保護區(qū)的6個樣地中毛栗種群的分布格局均表現(xiàn)為聚集分布，其中樣地3的毛栗群落聚集強度最低，樣地4的毛栗群落聚集強度最高。

參考文獻：

[1] 劉金福，洪 偉.格氏栲群落生態(tài)學(xué)研究——格氏栲林主要種群生態(tài)位的研究[J].生態(tài)學(xué)報，1999，19（3）：59-64.

[2] 康永祥，岳軍偉，雷瑞德，等.陜北黃龍山遼東櫟群落優(yōu)勢種群生態(tài)位研究[J].西北植物學(xué)報，2008，28（3）：574-581.

[3] 韓 路，席琳喬，王家強，等.塔里木河上游灰胡楊種群生活史特征與空間分布格局[J].生態(tài)學(xué)報，2013，33（19）：6181-6190.

[4] 巨天珍，張彥文，葛建團，等.小隴山國家級自然保護區(qū)油松群落生態(tài)位研究[J].干旱區(qū)資源與環(huán)境，2013，27（4）：171-175.

[5] 巨天珍，王立峰，康麗丹，等.小隴山自然保護區(qū)華山松群落8個主要種群的生態(tài)位研究[J].生態(tài)與農(nóng)村環(huán)境學(xué)報，2010，26（5）：461-465.

[6] 曹 春，巨天珍，王勤花.甘肅小隴山銳齒櫟林主要種群生態(tài)位的研究[J].甘肅農(nóng)業(yè)大學(xué)學(xué)報，2006，41（3）：65-69.

[7] 馬映棟，劉文楨，趙中華，等.小隴山銳齒櫟林種群空間格局及關(guān)聯(lián)性變化分析[J].西北植物學(xué)報，2014，34（9）：1878-1886.

[8] 高 洋，巨天珍，葛建團，等.國家級自然保護區(qū)小隴山白皮松優(yōu)勢種群生態(tài)位研究[J].福建林業(yè)科技，2010，37（1）：12-16.

[9] 劉小林，郭小龍，李安民，等.秦嶺西段小隴山林區(qū)紅樺天然種群分布格局研究[J].農(nóng)業(yè)科技與信息，2007（11）：48-50.

[10] 郭垚鑫，康 冰，李 剛，等.小隴山紅樺次生林物種組成與立木的點格局分析[J].應(yīng)用生態(tài)學(xué)報，2011，22（10）：2574-2580.

[11] 鄭 智，龔大潔，孫呈祥，等.秦嶺小隴山保護區(qū)維管植物豐富度和種域海拔梯度格局及其對Rapoport法則驗證[J].應(yīng)用生態(tài)學(xué)報，2014，25（9）：2477-2485.

[12] 索安寧.天水小隴山銳齒櫟群落生物多樣性研究[D].蘭州：西北師范大學(xué)，2003.

[13] 張毓?jié)?，李吉玫，常順利，?天山中部天山云杉種群空間分布格局及其與地形因子的關(guān)系[J].應(yīng)用生態(tài)學(xué)報，2011，22（11）：2799-2806.

[14] 吳承禎，洪 偉，謝金壽，等.珍稀瀕危植物長苞鐵杉種群生命表分析[J].應(yīng)用生態(tài)學(xué)報，2000，11（3）：333-336.

[15] 江 洪. 云杉種群生態(tài)學(xué)[M].北京：林業(yè)出版社，1992.

[16] NAVAS M L，GOULARD M. Spatial pattern of a clonal perennial weed Rubia peregrine（Rubuaceae） in vineyards of Southern France[J].J Appl Ecol，1991（28）：1118-1129.

[17] 張 峰，上官鐵梁.山西翅果油樹群落優(yōu)勢種群分布格局研究[J].植物生態(tài)學(xué)報，2000，24（5）：590-594.

[18] 張 瑾，陳文業(yè)，張繼強，等.甘肅敦煌西湖荒漠濕地生態(tài)系統(tǒng)優(yōu)勢植物種群分布格局及種間關(guān)聯(lián)性[J].中國沙漠，2013， 33（2）：349-357.

[19] 李軍玲，張金屯.太行山中段植物群落優(yōu)勢種生態(tài)位研究[J].植物研究，2006，26（2）：2156-2162.

[20] 李先琨，黃玉清，蘇宗明.元寶山南方紅豆杉種群分布格局及動態(tài)[J].應(yīng)用生態(tài)學(xué)報，2000，11（2）：169-172.