彭科, 陳旭, 張宣, 房超, 楊彪

1. 彭州市規(guī)劃和自然資源局，四川 彭州 611930；

2. 成都暢新環(huán)?？萍加邢薰?，四川 成都 610047；

3. 西華師范大學，西南野生動植物資源保護教育部重點實驗室，四川 南充 637002

大熊貓（Ailuropodamelanoleuca）是我國珍稀瀕危特有種，被稱為國寶，同時也是全球生物多樣性保護的旗艦物種和傘護物種，對大熊貓的保護研究一直受到國內外學者的廣泛關注。大熊貓原本為大型的食肉動物，隨著自然環(huán)境的改變，大熊貓在數(shù)千年的演化與適應下逐漸轉變?yōu)橐詠喐呱叫街耦悶橹魇车闹彩承詣游颷1-2]。

據(jù)全國第四次大熊貓調查：我國野生大熊貓主要分布于秦嶺、岷山、邛崍山、大相嶺、小相嶺和涼山等六大山系，野生大熊貓取食竹種共有8屬38種，然而不同區(qū)域的優(yōu)勢竹種不同，大熊貓主要取食竹種亦有差異[1]。大熊貓主食竹是大熊貓繁衍和發(fā)展的基礎和保證，開展大熊貓主食竹的調查研究對于大熊貓的保護具有重要的意義[3]。以往研究已對大熊貓主食竹的多樣性[4]、分類與分布[5-6]、天然更新[7]、種群結構[8]、種群動態(tài)及其影響因素[9]、生物量分配與模擬[10]、生理生化特征[11]、營養(yǎng)成分[12-14]、粗蛋白、粗脂肪和氨基酸[15]及其季節(jié)動態(tài)[16]、適口性[14]、自然災害的影響[17]、對全球變化的響應（如CO2濃度升高和N沉降[3]）、全球變化情景下竹子分布區(qū)的可能變化[18]、病蟲害[19]、人為管理與干擾影響[13,20]等方面做了相關報道，但這些研究僅局限于六大山系的部分保護區(qū)，如王朗保護區(qū)[7]、臥龍保護區(qū)[8,9,12]、秦嶺保護區(qū)[11,16]、大風頂保護區(qū)[20]、栗子坪保護區(qū)[10]、鞍子河保護區(qū)[14]等，研究區(qū)域不平衡。已有研究指出大熊貓取食竹的種類存在明顯的區(qū)域差異[1,21]。加之20世紀80年代岷山山系大熊貓主食竹曾出現(xiàn)大量開花致死現(xiàn)象，2005年以來又陸續(xù)發(fā)現(xiàn)岷山山系大熊貓主食竹有新的開花現(xiàn)象。因此，有必要加緊積累不同大熊貓分布區(qū)的大熊貓主食竹研究數(shù)據(jù)從而為不同區(qū)域大熊貓的科學保護提供重要的數(shù)據(jù)支撐。

四川白水河國家級自然保護區(qū)位于四川盆地西北邊緣的彭州市龍門山和小魚洞鎮(zhèn)境內，是離成都市最近（僅50 km）的大熊貓自然保護區(qū)，多年來保護區(qū)雖開展了不少針對大熊貓的保護工作，但針對大熊貓的研究報道極其缺乏，僅見于李華東[22]對保護區(qū)內大熊貓生存環(huán)境評價及研究的報道，目前保護區(qū)內關于大熊貓主食竹的相關研究仍較為缺乏。本研究通過對該保護區(qū)大熊貓主食竹開展調查，旨在查清該區(qū)域大熊貓主食竹的種類、種群特征及其受環(huán)境因子（海拔和郁閉度）的影響，以期為不同大熊貓分布區(qū)的大熊貓主食竹研究積累基礎數(shù)據(jù)，同時也為該保護區(qū)大熊貓的保護提供數(shù)據(jù)支撐和理論依據(jù)。

1 研究區(qū)域與方法

1.1 研究區(qū)域概況

四川白水河國家級自然保護區(qū)（下文簡稱“保護區(qū)”）位于龍門山褶皺帶的中南段、橫斷山東部，地理位置北緯31°10’—31°29’、東經103°41’—103°57’。該區(qū)域年平均降水量1 500 mm左右，主要集中在7、8、9這3個月；年平均氣溫為15.6 ℃，年平均無霜期277 d，日照僅711 h。保護區(qū)最低海拔1 481 m（鍋框巖）、最高海拔4 814 m（太子城），海拔高差跨度極大、地形地貌復雜，溫暖濕潤，降水豐富，云霧多日照少，且富垂直變化的氣候條件，形成了水熱條件各異的多種生態(tài)環(huán)境。土壤具有明顯的垂直地帶性，自下而上是山地黃棕壤、暗棕壤、棕色針葉林土、亞高山草甸土和高山草甸土。植被類型完整豐富，有亞熱帶山地常綠、落葉闊葉混交林、落葉闊葉林、針闊混交林、高山山地硬葉常綠闊葉林、暗針葉林、高山灌叢和高山草甸等。

1.2 野外調查

野外調查采用樣線法和樣方法進行。在收集保護區(qū)本底調查等歷史資料的基礎上共設置25條調查樣線（見圖1），每條樣線兩側50 m內布設不少于4個20 m×20 m的樣方進行調查。2018年5—6月，在白水河保護區(qū)及周邊共設置植物樣方70個，竹子樣方197個。由于第一次調查未能發(fā)現(xiàn)所有竹種的竹筍，同年9—10月在同一區(qū)域設置植物樣方36個，竹子樣方123個以補充調查。另外，調查期間在保護區(qū)邊界外圍發(fā)現(xiàn)有與保護區(qū)內同種的“小徑竹”分布，故做樣方調查以增加低海拔區(qū)域的樣方，所有調查樣方最低海拔1 240 m，最高海拔3 200 m。對樣方中的竹子實測其克隆分株數(shù)（筍、1年生克隆植株和超過1年生克隆植株）、基徑、高度，采用GPS和羅盤儀測量海拔、坡度和坡向等環(huán)境因子，采用“多點人工目測”方式測定郁閉度，同時記錄竹子的病蟲、開花、大熊貓和其他動物取食等情況。

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

以本研究獲得的所有“小徑竹種類”來組織研究區(qū)域的全部竹種名錄。在做環(huán)境因子相關分析時，由于刺黑竹（Chimonobambusa neopurpurea）、箬 葉 竹（Indocalamus longiauritus）和 柺 棍 竹（Fargasia robusta）的分布海拔范圍較窄，其調查得到的樣方數(shù)達不到統(tǒng)計的最小樣本數(shù)（各僅有6個小樣方），故僅對調查樣本量較多的白夾竹（Phyllostachys bissetii）、短錐玉山竹（Yushania brevipaniculata）、冷箭竹（Bashania fangiana）和細枝箭竹（Fargesia stenoclada）4種竹種做環(huán)境要素（海拔和郁閉度）與其種群特征之間的關系分析，以初步評估研究區(qū)域種群分布較廣的竹種受環(huán)境因子的變化情況。

圖1 調查樣線分布圖Fig.1Distribution of survey sample lines in Baishuihe National Nature Reserve, Sichuan

所有數(shù)據(jù)的統(tǒng)計檢驗分析在Origin Pro 8.0軟件上完成，數(shù)據(jù)作圖在 Excel 2016軟件上完成。

2 結果和分析

2.1 種類與分布

本研究共調查到竹類6屬7種，分別為剛竹屬（Phyllostachys）的白夾竹、方竹屬（Chimonobambusa）的刺黑竹、玉山竹屬（Yushania）的短錐玉山竹、箬竹屬（Indocalamus）的箬葉竹、巴山木竹屬（Bashania）的冷箭竹以及箭竹屬（Fargesia）的拐棍竹和細枝箭竹。7種竹子海拔分布范圍主要在3 000 m以下，郁閉度在0.3～0.9之間。本研究中僅發(fā)現(xiàn)短錐玉山竹和細枝箭竹2種竹子的莖葉被大熊貓取食后留下食痕，7種竹子均未發(fā)現(xiàn)被其他動物取食的痕跡，未發(fā)現(xiàn)有開花和病蟲害現(xiàn)象（見表1）。

表1 四川白水河國家級自然保護區(qū)內的竹子種類及概況Tab.1Bamboo species and general information in Baishuihe National Nature Reserve, Sichuan

2.2 種群特征

2.2.1 竹筍與幼竹

竹筍是評價竹子種群更新的重要指標之一。本研究期間調查樣方中未發(fā)現(xiàn)刺黑竹和箬葉竹的竹筍。其原因一方面可能與刺黑竹種群較少導致發(fā)筍情況較少有關，據(jù)保護區(qū)以往調查資料得知雖該區(qū)域刺黑竹的發(fā)筍期在農歷8月中旬之后，但本次調查獲得的刺黑竹的樣方數(shù)較少，刺黑竹種群數(shù)量不多從而導致發(fā)筍情況較少可能是本次調查難以見到其發(fā)筍情況的重要原因；另一方面可能與箬葉竹新筍識別難度較大有關，箬葉竹的新克隆分株是一邊生長，一邊發(fā)側枝，野外調查時極易忽略是“新筍”。而白夾竹、短錐玉山竹、柺棍竹、冷箭竹和細枝箭竹均有竹筍出現(xiàn)，但其竹筍樣方的百分比均小于50%（見表2）。

表27 種竹子的竹筍、1年生克隆分株、死亡成竹Tab.2Bamboo shoots, 1-year-old clonal ramets and dead bamboo of 7 bamboo species

1年生克隆分株是竹子種群更新擴展的新生力量，也是經過自然淘汰后保留的最具生命力的克隆分株。雖本研究在調查樣方中7種竹子均有出現(xiàn)“自然死亡的成熟克隆分株”，其中出現(xiàn)箬葉竹死亡的樣方占比最高（100%），然而除細枝箭竹（47.1%）外，其余6種竹子出現(xiàn)1年生克隆分株的樣方比均超過65%，而刺黑竹和箬葉竹出現(xiàn)1年生克隆分株的樣方比達到100%（見表2），表明調查區(qū)內竹子更新良好。

西南山地竹子是以克隆繁殖為主的多年生植物，1年生克隆分株比是衡量竹子種群更新和種群穩(wěn)定性的重要指標。研究期間發(fā)現(xiàn)1年生克隆分株比最高的是冷箭竹，其均值達25.9%；其次是白夾竹，達16.3%；再次是柺棍竹、短錐玉山竹和刺黑竹，分別為15.0%、13.5和11.7%；而箬葉竹和細枝箭竹的1年生克隆分株比的均值最低，在7%～9%之間（見表2）。

2.2.2 克隆分株密度

本研究發(fā)現(xiàn)，樣方內7種竹子的1年生克隆分株數(shù)量為0～12株。1年生克隆分株數(shù)量最多的是冷箭竹，最多達12株；其次是短錐玉山竹，最多達11株；而其余竹種的1年生克隆分株數(shù)量較低。1年生克隆分株密度最大的是冷箭竹，均值達6株·m?2；其次是箬葉竹，均值達4株·m?2；密度最小的是細致箭竹，為1.3株·m?2。

樣方內7種竹子的1年生以上克隆分株數(shù)量為2—103株。1年生以上克隆分株數(shù)量最多的是短錐玉山竹，最多達103株。1年生以上克隆分株密度最大的是箬葉竹，為44.8株·m?2；其次是刺黑竹，為33.4株·m?2；密度最小的是細枝箭竹，為11.1株·m?2。

樣方內7種竹子的數(shù)量為3～103株。其中，密度最大的是箬葉竹，為48.8株·m?2；其次是刺黑竹，為36.6株·m?2；密度最低的是細枝箭竹，為14.5株·m?2（見表3）。

2.2.3 基徑與稈高

本研究發(fā)現(xiàn)7種竹子中，克隆分株基徑最大的是刺黑竹，均值為8.4 mm；其次是白夾竹，為7.5 mm；最小的是冷箭竹，為4.6 mm。調查到最粗壯的克隆分株是細枝箭竹，基徑達15.6 mm；其次是短錐玉山竹，基徑達11.9 mm。另外，克隆分株最高大的是白夾竹，均值為189.0 cm；其次是刺黑竹，為180.7 cm；最矮小的是冷箭竹，為110.9 cm。調查到最高的克隆分株是白夾竹，達297.0 cm；其次是短錐玉山竹，達288.0 cm（見表3）。

表3 七種竹子的密度、基徑和高度Tab.3Density, basal diameter and height of 7 bamboo species

2.3 竹子與環(huán)境要素的關系

2.3.1 海拔的影響

統(tǒng)計分析發(fā)現(xiàn)，海拔變化僅對冷箭竹的密度有顯著影響（P<0.05），而對白夾竹、短錐玉山竹以及細枝箭竹的密度均無顯著影響（P>0.05），同時海拔變化對以上4種竹種的1年生克隆分株比、基徑以及高度均無顯著影響（冷箭竹：P>0.05，見圖2；白夾竹：P>0.05，見圖3；短錐玉山竹：R2<0.1，見圖4；細枝箭竹：P>0.05，見圖5）。

2.3.2 郁閉度的影響

統(tǒng)計分析發(fā)現(xiàn)，郁閉度變化對細枝箭竹的密度、基徑和高度均有極顯著影響（P<0.01），而對1年生克隆分株比無顯著影響（P>0.05，見圖6）；郁閉度變化對冷箭竹和白夾竹的密度有顯著影響（P<0.05），但對其1年生克隆分株比、基徑和高度無顯著影響（P>0.05；見圖7，圖8）；而郁閉度變化對短錐玉山竹的密度、1年生克隆分株比、基徑和高度均無顯著影響（R2<0.1，見圖9）。

3 討論

竹類是大熊貓的食物資源，是大熊貓生存和繁衍的基礎，不同區(qū)域竹子資源及種類的多樣性直接影響該分布區(qū)大熊貓的食譜情況。與其他大熊貓保護區(qū)相比，如四川王朗自然保護區(qū)（5種）[23]、臥龍保護區(qū)（6種）[8,24]等，四川白水河保護區(qū)的竹子種類相對較多，高達7種，且7種竹類的不同克隆分株年齡階段的無性系分株都有分布，各類竹種的年齡結構相對合理，除刺黑竹、箬葉竹和柺棍竹的種群數(shù)量相對較少、分布區(qū)域相對集中外，其余4種竹類資源相對豐富，能夠為該區(qū)域大熊貓?zhí)峁┹^為充足的食物資源。調查區(qū)域的7種竹類海拔分布多在3 000 m以下，竹類分布海拔略低于其他保護區(qū)。

已有研究指出，在大熊貓分布區(qū)不同大熊貓主食竹的種群數(shù)量、生長情況以及分布面積往往存在差異[4,24]。本研究取得了類似的研究結果，如研究發(fā)現(xiàn)該區(qū)域短錐玉山竹、白夾竹、刺黑竹增長率較高，而冷箭竹增長率相對較低；短錐玉山竹在數(shù)量和分布面積上占有絕對優(yōu)勢等。本研究發(fā)現(xiàn)，雖白水河保護區(qū)竹子種類較多，但大熊貓取食的種類較少（僅發(fā)現(xiàn)短錐玉山竹和細枝箭竹有取食現(xiàn)象），這與全國第四次大熊貓調查的結果一致[1]，即同一區(qū)域，大熊貓一般僅取食少數(shù)幾種竹子，這可能與大熊貓的取食偏好有關。已有研究發(fā)現(xiàn)不同區(qū)域大熊貓取食竹種存在差異[1,24,25]，并且相同區(qū)域的大熊貓食性也會存在季節(jié)性變化[26]。

環(huán)境因子諸如溫度、濕度、光照強度等對竹林密度、高度、基徑、竹筍的高生長、出筍量、無性系根莖等具有顯著影響[27]。而海拔因子包含了諸多環(huán)境因子，因此竹林各項指標往往會隨著海拔變化而變化[27,28]。本研究表明，海拔變化對冷箭竹的密度影響顯著，海拔2 400～3 200 m的冷箭竹明顯高于低海拔區(qū)域情況，說明2 400～3 200 m為該區(qū)域冷箭竹生長更新的最適海拔，天然更新良好。然而本研究發(fā)現(xiàn)海拔對白夾竹、短錐玉山竹以及細枝箭竹的密度無顯著影響，同時海拔對以上4種竹種的1年生克隆分株比、基徑以及高度也無顯著影響，說明以上4種竹種在海拔分布上相對均勻（冷箭竹：2 150～3 200 m；白夾竹：1 240～1 760 m；短錐玉山竹：1 540～2 800 m；細枝箭竹：1 750～2 500 m）。喬木作為森林群落的最高層，其郁閉度的大小會對林下環(huán)境（如光照強度、溫度、濕度等）產生影響，從而影響竹子的生長和天然更新。本研究發(fā)現(xiàn)，郁閉度變化對細枝箭竹的密度、基徑和高度以及對冷箭竹和白夾竹的密度影響顯著，說明不同喬木郁閉度下竹子生長和更新往往存在差異。這與李波等對岷山北部箭竹的研究結果[29]以及賈炅和武吉華[7]對王朗自然保護區(qū)缺苞箭竹的研究結果類似。已有研究指出，郁閉度過高的過分蔭蔽環(huán)境會導致竹子缺乏光照而影響根系的發(fā)育，使竹子幼苗的根系生長脆弱，根系競爭力差，進而影響幼苗對水分和養(yǎng)分的吸收，造成幼苗死亡[30]。而在郁閉度很低的開闊地段，因為林下地面的光照充足，這些區(qū)域林下往往會灌木雜草叢生，從而不利于竹子的生長和天然更新[7,31]。

圖2 海拔對冷箭竹的影響Fig.2Effects of altitude on Bashania fangiana

圖3 海拔對白夾竹的影響Fig.3Effects of altitude on Phyllostachys bissetii

圖4 海拔對短錐玉山竹的影響Fig.4Effects of altitude on Yushania brevipaniculata

圖5 海拔對細枝箭竹的影響Fig.5Effects of altitude on Fargesia stenoclada

圖6 郁閉度對細枝箭竹的影響Fig.6Effects of canopy density on Fargesia stenoclada

圖7 郁閉度對冷箭竹的影響Fig.7Effects of canopy density on Bashania fangiana

圖8 郁閉度對白夾竹的影響Fig.8Effects of canopy density on Phyllostachys bissetii

圖9 郁閉度對短錐玉山竹的影響Fig.9Effects of canopy density on Yushania brevipaniculata

已有研究發(fā)現(xiàn)，相同竹種的植株基徑和株高存在顯著的正相關性[7]。本研究發(fā)現(xiàn)短錐玉山竹基徑-株高變化系數(shù)較高，說明在相同粗壯程度下，短錐玉山竹能長更高以獲得更多陽光，進而易于在取食竹群落中不斷積累優(yōu)勢，成為優(yōu)勢竹種。

本研究發(fā)現(xiàn)，盡管全光照下細枝箭竹有較高的密度、基徑和高度，但由于全光照不利于大熊貓的活動，因而在管理中建議傾向于維持在0.4左右的郁閉度為宜。另外，研究發(fā)現(xiàn)該區(qū)域大部分竹類偏好郁閉度較低的生境，因此在管理中可考慮對茂密森林（尤其是棲息地內的人工純林）進行適當透光疏伐或以林窗形式疏伐而形成林窗，使棲息地的喬木層平均郁閉度保持在一個相對較為適宜的范圍內（以0.4～0.7為宜）。

竹子具有周期性開花死亡的特點，大熊貓在較為單一竹種組成的生境中，為避免將來可能會面臨由竹類周期性開花死亡造成的食物相對短缺問題，可以適當進行不同竹種的鑲嵌種植以提高竹類的多樣性，同時也可以種植不同竹鞭年齡的竹子以形成異齡竹林。此外，在實際管理中還應該與科研單位建立合作關系，探索大熊貓取食較多的短錐玉山竹管理模式。