古玉, 陳明華, 曾進科, 馬曉平

1.四川農(nóng)業(yè)大學 a.生命科學學院； b.動物醫(yī)學院，四川 成都 611130；

2.四川美姑大風頂國家級自然保護區(qū)管理局，四川 美姑 616450

竹類植物屬禾本科（Poaceae）竹亞科（Bambusoideae），有68屬1 500余種，廣布大洋洲、非洲、美洲和亞洲等熱帶、亞熱帶地區(qū)。中國有32（37）屬450~500種，主要分布在長江流域及以南各省，四川分布160余種[1]。竹子占大熊貓食譜的99%，經(jīng)數(shù)千年演化適應(yīng)，大熊貓形成了以亞高山竹類為主要食物的覓食策略，竹葉、莖稈、竹枝、竹筍均會取食。并非所有竹種都為大熊貓所喜愛，具匯集資料，有13屬、79種、1變種、3栽培品種計83種及種下分類群作為分布于四川、甘肅、陜西等省的野生大熊貓主食竹加以利用[2]。大熊貓主食竹的地理分布是決定大熊貓地理分布的重要影響因素之一。主食竹的生長情況決定大熊貓的遷徙和繁衍。

美姑大風頂國家級自然保護區(qū)位于四川省涼山彝族自治州美姑縣東北部，地處橫斷山脈中段，四川盆地與青藏高原東南緣間，黃茅埂山脈頂峰大風頂以西，生物學地理位置上處于古北界青藏區(qū)東南緣。保護區(qū)屬深切割中山地貌類型，地勢由西南向東北傾斜，海拔1 356~4 042 m，垂直高差2 686 m[3]。該區(qū)是大熊貓涼山山系種群的集中分布區(qū)和腹心區(qū)，其連接起馬邊大風頂、越西申果莊、甘洛馬鞍山、峨邊黑竹溝自然保護區(qū)，構(gòu)成了涼山山系大熊貓保護網(wǎng)絡(luò)，在該山系的大熊貓及生物多樣性保護中具有不可替代的地位[4]。據(jù)四川省第四次大熊貓調(diào)查報告顯示，本保護區(qū)大熊貓棲息地40 531.58 hm2，潛在棲息地701.19 hm2[5]。

保護區(qū)于1978年批準成立時位于美姑縣瓦侯鄉(xiāng)境內(nèi), 東臨羅姑坡、 羅姑列諧、 覺羅豁、 大風頂, 山嶺線與馬邊大風頂國家級自然保護區(qū)相連, 南以覽波加界 、 牙斯普山脊為界, 西面以此補聶山脊與瓦侯鄉(xiāng)龍窩村接壤, 北面以瓦候河為界，介于東經(jīng)103°50′~103°20′，北緯28°36′~28°45′，面積21 864 hm2，主要分布熊竹、冷箭竹、大風頂玉山竹、馬邊玉山竹、八月竹，海拔范圍1 240~4 042 m[3]。2005年將蘇洛鄉(xiāng)、依果覺鄉(xiāng)和炳途鄉(xiāng)部分區(qū)域劃入保護區(qū)范圍，擴建后面積為50 655 hm2，地理坐標為東經(jīng)102°52′~103°20′，北緯28°30′~28°50′；擴建區(qū)域內(nèi)大面積分布短錐玉山竹、熊竹和冷箭竹，海拔范圍為1 988~3 765 m。日常巡護和專項調(diào)查均發(fā)現(xiàn)2018年前大熊貓在擴建前區(qū)域內(nèi)活動頻繁，但2018年起未見大熊貓活動痕跡，大熊貓遷徙出擴建前區(qū)域。

為科學合理保護大熊貓，掌握本保護區(qū)竹子種類和分布情況，并從食物和環(huán)境上調(diào)查大熊貓遷徙出擴建前區(qū)域的原因，于2020年開展了該區(qū)竹子資源及擴建前后區(qū)域共有竹種生長情況專項調(diào)查，掌握了本保護區(qū)竹子種類、分布和大熊貓食物喜好等情況，并從食物因素分析大熊貓遷徙出原活動區(qū)域原因。

1 方法

1.1 野外調(diào)查

1.1.1 竹種分布和面積調(diào)查

對保護區(qū)50 655 hm2全域設(shè)置樣線，抽樣比例達0.37 km·100 hm-2，樣線覆蓋所有區(qū)域和森林類型。樣線上每100 m海拔差設(shè)置1個樣方；如100 m海拔間樣線距離大于500 m，則在每間隔500 m處設(shè)置1個樣方。各種竹類出現(xiàn)處和結(jié)束處設(shè)置樣方，且需沿著該處海拔高度向樣線兩側(cè)各500 m處設(shè)置1個樣方。一條樣線上的同種竹類采集一式3份標本。

1.1.2 大熊貓遷出和現(xiàn)活動區(qū)域內(nèi)熊竹生長情況調(diào)查

設(shè)置1 m×1 m樣方調(diào)查各齡級竹子數(shù)量、死亡株數(shù)、各生境因子、干擾情況，每樣方中隨機選擇10株多年生竹測量株高、竹竿基部直徑。

1.2 竹種鑒定

對所有采集標本的竹筍，竹鞭，1年生幼竹和桿籜，多年生竹的主莖、枝條和葉進行形態(tài)觀測，以《中國植物志》和《四川竹類植物志》的分類系統(tǒng)進行鑒定[6,7]。

1.3 竹類分布和生長情況與大熊貓活動范圍間關(guān)系

根據(jù)同一海拔高度竹類分布基本一致的原則，對調(diào)查樣方中竹子種類沿海拔高度向樣線兩側(cè)延伸，采用1∶50 000地形圖在ARGIS軟件中繪制各竹種分布區(qū)域、計算分布面積。同時根據(jù)日常紅外相機監(jiān)測和2018年大熊貓野外調(diào)查糞便采集項目獲得的大熊貓活動痕跡，在竹類分布圖上標注出2018年前后大熊貓活動痕跡點，以分析大熊貓取食竹子偏好。并對大熊貓遷出區(qū)域和現(xiàn)活動區(qū)域內(nèi)熊竹的生長情況做單因素ANOVA方差分析，對喬木層郁閉度、干擾因素和生境因子進行χ2檢驗，差異顯著性檢驗標準為p＜0.05，以期分析竹類生長情況及環(huán)境與大熊貓活動范圍間關(guān)系。

2 結(jié)果

2.1 竹子種類鑒定

根據(jù)鑒定結(jié)果，本保護區(qū)主要分布6種竹種。分別為熊竹（Yushania ailuropodinaT. P. Yi）、短錐玉山竹（Yushania brevipaniculata（Hand.-Mazz.）Yi）、大風頂玉山竹（Yushania dafengdingensisT.P.Yi）、馬邊玉山竹（Yushania mabianensisT. P. Yi）、冷箭竹（Arundinaria faberiRendle）、八月竹（Chimonobambusa szechuanensis（Rendle） Keng f.）。

2.2 竹類資源分布與大熊貓活動范圍比較

2.2.1 竹類資源水平分布與大熊貓活動空間關(guān)系

保護區(qū)內(nèi)竹子分布在東經(jīng)102.884 2°~103.332 0°、北緯28.517 0°~28.828 3°間。

2018年后大熊貓活動痕跡點位分布在東經(jīng)102.887 7°~103.138 8°、北 緯28.609 5~28.810 0°間，主要在欽波惹、竹魯依付、椅子埡口、椅子河壩、阿米都爾、毛洪覺、馬加耶依達、維核洛和濫龍等區(qū)域。2018年調(diào)查巡護中首次發(fā)現(xiàn)大熊貓在毛洪覺區(qū)域活動，2019年持續(xù)觀察到活動，但2020和2021年未觀察到活動痕跡。區(qū)域內(nèi)分布熊竹、短錐玉山竹和冷箭竹。大熊貓對竹子的利用面積分別為短錐玉山竹4 517 hm2，冷箭竹2 871 hm2，熊竹2 049 hm2（見表1）。

2018年前，除上述區(qū)域外，波么古資和大風頂區(qū)域大熊貓活動頻繁，活動范圍為東經(jīng)102.886 1°~103.323 9°、北緯28.612 0°~28.811 7°。大熊貓對竹子的利用面積分別為熊竹5 563 hm2、短錐玉山竹4 152 hm2，冷箭竹2 504 hm2，大風頂玉山竹2003 hm2，馬邊玉山竹335 hm2。

表1 2018年前、后大熊貓活動區(qū)域各竹種面積Tab. 1 Areas of Bamboo Species in giant pandas activity area before and after 2018

對比可知，該保護區(qū)還有較大面積竹林沒有大熊貓活動，2018年后大熊貓在美姑大風頂?shù)幕顒酉蛭髌?，緯度上則無太大變化，對竹類利用面積減小5 120 hm2。遷出區(qū)域主要生長大風頂玉山竹、馬邊玉山竹、八月竹、熊竹和冷箭竹。

2.2.2 竹類資源垂直分布與大熊貓活動空間關(guān)系

根據(jù)調(diào)查，美姑大風頂竹子分布的海拔范圍為1 356~3 770.86 m。熊竹主要分布于大風頂、維核洛區(qū)域，海拔1 356~2 820.75 m連續(xù)分布，面積達16 988 hm2。冷箭竹在大風頂、維核洛、濫龍、椅子河壩等區(qū)域分布，面積13 383 hm2，海拔2 844.95~3 770.86 m間呈帶狀分布。短錐玉山竹分布于濫龍和椅子河壩區(qū)域，面積11 114 hm2，海拔2 454.14~3 296.74 m成片分布。大風頂玉山竹分布于大風頂區(qū)域，面積較小，僅有2003 hm2, 海拔2 816~3 098 m呈帶狀分布。馬邊玉山竹緊鄰大風頂玉山竹分布，同屬于大風頂區(qū)域，海拔2 562~2 816 m間呈狹長條帶分布，有335 hm2。八月竹呈點狀區(qū)塊分布于大風頂、維核洛區(qū)域，海拔2 070~2 170 m間，僅有263 hm2。（見圖1）

比較可知，大熊貓在2018年前和2018年后活動區(qū)域的海拔范圍基本一致，均為2 246~3 616 m，分布著除八月竹以外的其他5種竹子?？梢姾０尾皇怯绊懘笮茇垙拇箫L頂區(qū)域遷出的原因。

2.2.3 大熊貓對竹種食性偏好

圖1 竹類分布區(qū)域圖Fig. 1 Distribution map of bamboo species

根據(jù)歷年來大熊貓活動區(qū)域可知，大熊貓不在八月竹林區(qū)生活，其余5種竹林中均活動頻繁。本保護區(qū)大熊貓喜食除八月竹外的其他5種竹種，喜食竹面積達43 823 hm2，而現(xiàn)在大熊貓僅利用了其中的14 557 hm2。據(jù)魏輔文[8]測算臥龍大熊貓的巢域面積為3.9~6.4 km2，由此可測算理論上美姑大風頂保護區(qū)可容納68~112只大熊貓，現(xiàn)保護區(qū)內(nèi)僅有大熊貓23只，完全可滿足更多大熊貓生活。（見見見見見見見見見見見見見見見圖1）

2.3 大熊貓遷出和現(xiàn)活動區(qū)域內(nèi)熊竹生長及生境分析

熊竹為大熊貓遷出區(qū)域和2018年后活動區(qū)域共有竹種，故我們對熊竹生長情況、環(huán)境因子和干擾因素進行統(tǒng)計分析。

2.3.1 熊竹林區(qū)生境及干擾情況

共調(diào)查到遷出區(qū)域內(nèi)熊竹樣方27個，在常綠闊葉林、落葉闊葉林區(qū)、針闊混交林區(qū)和針葉林區(qū)各分布1、9、13和4個樣方。建群樹種主要為冷杉、云杉、鐵杉、樺木、楓楊、水青樹、連香樹等。2018年后活動區(qū)域內(nèi)49個熊竹樣方，其中常綠闊葉林區(qū)1個、落葉闊葉林區(qū)18個、針闊混交林區(qū)25個、針葉林區(qū)5個。建群樹種主要為冷杉、云杉、鐵杉、落葉松、樺木、楓楊、水青樹、漆樹、楊樹、杜鵑等?？梢姡苤裨趦蓚€區(qū)域內(nèi)生境一致，χ2檢驗亦顯示無顯著差異。

遷出區(qū)域內(nèi)熊竹分布在0.2~0.7郁閉度的喬木層下，其中尤以0.5郁閉度為主，占該區(qū)域內(nèi)熊竹樣方數(shù)的40.74%，0.2、0.4和0.6郁閉度喬木層下分布樣方數(shù)分別占總樣方比例的10.20~18.52%，0.3和0.7郁閉度喬木層下僅分布少量樣方。現(xiàn)活動區(qū)域內(nèi)熊竹林喬木層郁閉度范圍為0.2~0.6，0.2~0.4郁閉度下分布的樣方數(shù)比較均勻，比例為7.41%~12.24%；5郁閉度下樣方數(shù)最多，占比為57.14%；在0.6郁閉度的喬木層下僅分布2.04%的熊竹樣方?？梢妰煞N區(qū)域內(nèi)熊竹林的喬木層郁閉度都以0.2~0.5為主，0.7以上的高郁閉度下熊竹難以生長。χ2檢驗結(jié)果顯示熊竹生長范圍內(nèi)兩個林區(qū)郁閉度不存在顯著差異。

χ2檢驗顯示兩個林區(qū)范圍內(nèi)熊竹均有大于40%的樣方受到放牧或割竹等明顯人為干擾，其無顯著差異。生境因子中，兩個熊竹林區(qū)在山谷分布較少，坡位無顯著差異；但地貌和坡度顯示出極顯著差異，大熊貓現(xiàn)活動區(qū)域熊竹生長在16~35°緩坡和中坡更多，主要呈現(xiàn)中山地貌；坡向亦呈現(xiàn)顯著差異，生長在東南坡的竹子樣方明顯增加。（見表2）

2.3.2 熊竹生長情況

從表3統(tǒng)計結(jié)果看，大熊貓遷出和現(xiàn)活動區(qū)域內(nèi)熊竹的竹高、地徑、死亡竹數(shù)、一年生活竹數(shù)、多年生活竹數(shù)、樣方內(nèi)竹子株數(shù)均無顯著差異。熊竹高度約2.7 m，基部竹竿直徑1.1 cm左右。1 m2單位面積內(nèi)約有竹子30株，一年生活竹約5株，多年生活竹約18株，死亡竹約6株。

表2 大熊貓遷出和現(xiàn)活動區(qū)域內(nèi)熊竹林生境分析Tab. 2 Habitat analysis of Y. ailuropodina in the reserve where giant pandas moved out and current living

表3 大熊貓遷出和現(xiàn)活動區(qū)域內(nèi)熊竹生長狀況Tab. 3 Growth condition of Y. ailuropodina in the reserve where giant pandas moved out and current living

圖2 大熊貓活動痕跡點與竹類分布關(guān)系Fig. 2 Relationship between giant pandas activities area and bamboo species distribution

但竹林蓋度和當年生筍數(shù)有顯著差異。遷出區(qū)域內(nèi)當年發(fā)筍1±2棵，顯著少于現(xiàn)活動區(qū)域內(nèi)的3±2棵；竹林蓋度為0.54±0.1，顯著高于現(xiàn)活動區(qū)域內(nèi)的0.47±0.15。

3 討論

美姑大風頂國家級自然保護區(qū)擴建前區(qū)域歷來是大熊貓頻繁活動區(qū)域，但近年來卻少有大熊貓活動痕跡，大熊貓遷出此區(qū)域。從竹子分布和生長情況看，此處的熊竹、馬邊玉山竹、大風頂玉山竹均是大熊貓喜食竹種。統(tǒng)計分析表明，熊竹的生長勢在大熊貓遷出區(qū)域和現(xiàn)活動區(qū)域內(nèi)基本一致，但發(fā)筍率和竹林蓋度有顯著差異。有研究指出，濕度和降雨量主要通過控制竹子發(fā)筍和生長來影響大熊貓的遷徙[9]；大熊貓取食和巴山木竹發(fā)筍空間上呈現(xiàn)正相關(guān)[10]。通過對佛坪自然保護區(qū)大熊貓活動與巴山木竹間關(guān)系的研究表明，春季大熊貓更偏好選擇竹子密度和蓋度略低、幼竹與竹筍比例高的竹林活動[11]。張澤鈞[12]通過對全國第三次和第四次大熊貓調(diào)查數(shù)據(jù)分析都發(fā)現(xiàn)竹類資源狀況在大熊貓選擇微生境時具有頭等重要地位，大熊貓偏好選擇高度高，蓋度大的竹林作為棲息區(qū)域。竹子密度過大導致其營養(yǎng)質(zhì)量較差且大熊貓在其中穿行困難，而過稀則能量收入少；攝食竹筍的效率高則單位收益耗能相對較少，故竹子密度和發(fā)筍率是大熊貓對微生境的選擇因素。擴建前保護區(qū)內(nèi)熊竹發(fā)筍率顯著低于擴建后區(qū)域，而竹林蓋度則顯著高于擴建后區(qū)域內(nèi)熊竹，這可能是大熊貓在此區(qū)域活動減少的因素。因此，通過適當?shù)氖璺ゴ胧?，砍掉部?年以上的老竹子，降低竹林密度，從而降低蓋度，增大通光量，加速林下竹筍的發(fā)筍率，提高枝葉的生物量。同時降低竹林密度有利于大熊貓穿行，降低其能力消耗，以促進大熊貓在遷出范圍內(nèi)的活動恢復。

熊竹林區(qū)的生境和受干擾程度在大熊貓遷出和現(xiàn)活動區(qū)域內(nèi)基本一致，但地貌、坡向和坡度存在顯著差異?？梢娙祟惢顒?，如放牧、割竹等，不是導致大熊貓遷徙出原保護區(qū)的根本原因。魏輔文等[13]對馬邊大風頂自然保護區(qū)大熊貓對生境選擇研究發(fā)現(xiàn)，其喜歡在坡度小于30°的東南緩坡活動，這與我們的調(diào)查結(jié)果一致。研究表明動物選擇生境時會權(quán)衡代價與收益，最大能量收入是其選擇的主要驅(qū)動力[14]。美姑大風頂保護區(qū)大熊貓遷出區(qū)域內(nèi)熊竹部分分布在高山地區(qū)，而現(xiàn)生活區(qū)域熊竹則全部分布在中山地貌區(qū)，大熊貓生活于中山地貌耗能更少。北半球山脈的東南坡植被較西北坡茂盛，能提供的單位面積能量更多，也符合動物生境選擇原則。但周世強[15]等的研究指出，坡度、坡向、坡位、地貌與水源等是大熊貓選擇的因子或偏好的生境結(jié)構(gòu)，而不可能成為大熊貓地理分布的限制條件，因為野生大熊貓可以通過移動來找尋其較為理想的生活區(qū)域?？梢姷孛?、坡向和坡度間的差異，是導致大熊貓遷移出原保護區(qū)的因素之一，但不應(yīng)該是決定因素，還應(yīng)結(jié)合食物條件、隱蔽條件、氣候條件、共生動物等各因素綜合分析。

致謝四川美姑大風頂國家級自然保護區(qū)管理局提供項目資金，四川涼山州飛播林區(qū)保護站站長方奎，四川美姑大風頂國家級自然區(qū)管理局局長曲比石布、辦公室主任龔億華、保護科科長李良及沙馬牛布等巡護人員的指導和大力支持，在此一并感謝！