晏婷婷,冉江洪,趙晨皓,鐘 雪,梁敏儀

四川大學(xué)生命科學(xué)學(xué)院生物資源與生態(tài)環(huán)境教育部重點實驗室,成都 610064

?

氣候變化對邛崍山系大熊貓主食竹和棲息地分布的影響

晏婷婷,冉江洪*,趙晨皓,鐘 雪,梁敏儀

四川大學(xué)生命科學(xué)學(xué)院生物資源與生態(tài)環(huán)境教育部重點實驗室,成都 610064

氣候變化對生物多樣性的影響,特別是珍稀瀕危物種的影響是當(dāng)前的研究熱點。全球氣候變化對大熊貓的影響一直受到廣泛關(guān)注。根據(jù)野外調(diào)查的大熊貓活動痕跡點、竹類分布點和主食竹擴(kuò)散距離數(shù)據(jù),采用Maxent模型,利用植被、地形、氣候等因素,在RCP8.5下分析了2050年和2070年邛崍山系大熊貓主食竹分布及棲息地變化趨勢。結(jié)果顯示：(1)未來大熊貓適宜生境及主食竹氣候適宜區(qū)面積均有所減少,到2070年分別減少37.2%和4.7%；(2)未來主食竹分布范圍總體向高海拔擴(kuò)展,但面積持續(xù)減少,到2070年分布面積比當(dāng)前減少8.3%；(3)大熊貓棲息地未來有向高海拔擴(kuò)張的趨勢,在低海拔地區(qū)退縮明顯,到2070年較當(dāng)前減少27.2%；但到2070年大熊貓棲息地面積加上非棲息地有主食竹分布的面積,較現(xiàn)有大熊貓棲息地面積大1.5%；(4)受氣候變化影響較嚴(yán)重的區(qū)域是邛崍山系南部以及低海拔地區(qū),其余區(qū)域所受影響相對較?。?5)未來需要加強對受氣候變化影響嚴(yán)重區(qū)域的監(jiān)測與保護(hù),特別是邛崍山系中部的大熊貓集中分布區(qū)。

氣候變化；大熊貓；邛崍山系；竹子；保護(hù)

大量陸地和海洋觀測證據(jù)表明,地球上許多自然系統(tǒng)正在受到區(qū)域氣候變化的影響,特別是溫度升高的影響,氣候變暖可能會對全球生態(tài)系統(tǒng)和生物多樣性造成極大的影響[1]。大熊貓(Ailuropodamelanoleuca)是全球珍稀瀕危野生動物,是生物多樣性保護(hù)的旗艦物種[2- 3]。在全球氣候變暖的大背景下,大熊貓棲息地氣候總體向干暖趨近[4],有研究預(yù)測氣候變化將導(dǎo)致岷山山系67%的大熊貓棲息地向高海拔、高緯度遷移[5]。并且適宜生境面積減小[6- 7]。秦嶺山系2080年大熊貓適宜生境可能喪失60%,在現(xiàn)有的大熊貓保護(hù)區(qū)內(nèi)僅有15%適合大熊貓棲息,且破碎化嚴(yán)重[8- 9],大熊貓適宜生境的最低海拔將上升500m[9]。然而作為大熊貓棲息地破碎化較嚴(yán)重的邛崍山系,氣候變化對大熊貓棲息地和主食竹的研究仍是空白。

氣候變化不僅可以直接影響物種的分布,也可以通過影響植物的生長和分布,進(jìn)而影響植食性動物種群的大小與分布情況[10]。大熊貓食性單一,亞高山竹類占據(jù)其食物組成的99%以上[11]。在氣候變化背景下,秦嶺大熊貓棲息地內(nèi)的3種主食竹分布面積將大幅度減少,到21世紀(jì)末,如果竹種不能擴(kuò)散到現(xiàn)有分布以外的新氣候適宜區(qū)域,秦嶺地區(qū)幾乎所有的大熊貓棲息地都可能會消失[12]。大熊貓對主食竹的強烈依賴性,使得其分布受到主食竹的分布和生活狀態(tài)的影響[13- 15]。

從大熊貓現(xiàn)有的分布海拔以及歷史分布來看,大熊貓對氣候變化應(yīng)有一定的適應(yīng)能力。但竹子極易受氣候變化的影響,因此氣候變化對大熊貓產(chǎn)生影響主要是通過影響竹子的分布,其次才是氣候變暖本身對大熊貓產(chǎn)生的影響。因而預(yù)測氣候變化對大熊貓影響,竹種實際分布數(shù)據(jù)以及其擴(kuò)散能力是極為關(guān)鍵的。而已有的針對氣候變化對大熊貓棲息地的研究,主要是利用氣候因子和大熊貓野外痕跡分布點進(jìn)行預(yù)測[4,6- 9],也有利用氣候因子和竹種分布點進(jìn)行預(yù)測[12],少有綜合利用氣候、竹種分布(包括擴(kuò)散能力)及大熊貓痕跡點的預(yù)測研究[5],并且對主食竹分布和擴(kuò)散能力的關(guān)注還較少。本文利用最新完成的四川省第四次大熊貓野外調(diào)查的大熊貓活動痕跡點、竹類分布點數(shù)據(jù),以及2014—2015測量的各竹種生長擴(kuò)散數(shù)據(jù)；同時結(jié)合第三次大熊貓野外調(diào)查數(shù)據(jù),采用Maxent模型,綜合利用植被、氣候、地形等因子,分析在RCP8.5條件下,2050年和2070年邛崍山系大熊貓適宜棲息地和主食竹分布變化,綜合考慮大熊貓適宜棲息地以及主食竹分布變化對大熊貓棲息地的影響及變化趨勢,以期更好的了解氣候變化可能對大熊貓帶來的影響,為區(qū)域大熊貓及其生物多樣性保護(hù)提供參考。

1 研究區(qū)域概況

圖1 研究區(qū)域地理位置圖Fig.1 Location of the study area

研究區(qū)域包括邛崍山系和大相嶺山系現(xiàn)今大熊貓的主要分布區(qū)(大相嶺山系與邛崍山山系大熊貓棲息地之間無天然障礙,植被相連,故作為整體分析),行政區(qū)劃涉及四川省的成都、阿壩、雅安、眉山、樂山、甘孜6個市(州)的小金、理縣、汶川、崇州、大邑、蘆山、寶興、天全、滎經(jīng)、洪雅、峨眉山以及康定(部分地區(qū))和瀘定(部分地區(qū))13個縣,介于東經(jīng)102°1′45.9″—103°49′10.4″,北緯29°16′35.3″—31°55′6″之間(圖1)。本區(qū)域有臥龍、四姑娘山、蜂桶寨等國家級自然保護(hù)區(qū),區(qū)域最高海拔6250m,大熊貓棲息地世界遺產(chǎn)地也位于本區(qū)域。據(jù)四川省第四次大熊貓調(diào)查結(jié)果顯示：該區(qū)域有野生大熊貓558只,占四川省野生大熊貓種群數(shù)量的40.2%,占全國的29.9%[16]。

研究區(qū)內(nèi)竹類資源豐富,分布的主食竹種有6屬19種,其中分布范圍最廣的是冷箭竹(Bashaniafaberi)、短錐玉山竹(Yushaniabrevipaniculata)、拐棍竹(Fargesiarobusta)和八月竹(Chimonobambusaszechuanensis)4種,共占研究區(qū)域大熊貓棲息地所有竹種分布面積的90%以上[16]。

2 方法

2.1 數(shù)據(jù)來源

2.1.1 大熊貓痕跡點和竹種分布點數(shù)據(jù)

大熊貓痕跡點數(shù)據(jù)和主食竹分布數(shù)據(jù)來源于四川省第三、第四次大熊貓調(diào)查在野外采集的GPS位點和海拔數(shù)據(jù),共有大熊貓痕跡點數(shù)據(jù)1092個,大熊貓主食竹分布點數(shù)據(jù)8561個。

2.1.2 氣候因子數(shù)據(jù)

氣候數(shù)據(jù)來源于WorldClim網(wǎng)站(http://www.worldclim.org/)下載的當(dāng)前以及BCC-CSM1- 1模式RCP8.5下的2050、2070年的氣候數(shù)據(jù),其空間分辨率為30 arc-seconds。已有的研究顯示,由于橫斷山區(qū)域地形特殊,BCC-CSM1- 1模式對中國區(qū)域氣候的模擬效果較好[17]。在SRES的6組情景中A2情景的更為接近我國未來的可能情況[8],因此本研究選擇與A2情景當(dāng)量二氧化碳濃度更為接近的RCP8.5進(jìn)行分析(RCP8.5也是RCPs中強度最大的等級[18])。

2.1.3 其他環(huán)境因子數(shù)據(jù)

地形數(shù)據(jù)包括海拔、坡度、坡向3個環(huán)境因子,來源于中國科學(xué)院數(shù)據(jù)庫的DEM數(shù)據(jù),在ArcGIS 10.1中,通過空間分析提取獲得。河流、道路和居民點分布數(shù)據(jù)來源于四川省1∶20萬電子地圖。植被圖是以2009年的ETM衛(wèi)片(精度為1km×1km),通過非監(jiān)督分類,并以四川省第三、四次大熊貓調(diào)查的植被數(shù)據(jù)校驗獲得。

2.1.4 竹類擴(kuò)散數(shù)據(jù)

2014—2015年,分別于滎經(jīng)龍蒼溝、寶興蜂桶寨、都江堰龍池三地分別對八月竹(調(diào)查海拔1800—2200m)、冷箭竹(調(diào)查海拔2400—2800m)、短錐玉山竹(調(diào)查海拔1900—2200m)和拐棍竹(調(diào)查海拔1800—1900m)的擴(kuò)散能力進(jìn)行了調(diào)查和測量,在竹林內(nèi)隨機設(shè)置1m×1m的樣方,測量當(dāng)年生竹筍到距離該竹筍最近的成竹之間的距離(每個竹種測量不少于200個數(shù)據(jù))。從測量結(jié)果看,4個竹種每年的平均擴(kuò)散距離都很短,分別為八月竹(21.3±0.71)cm/a,短錐玉山竹(10.1±0.52)cm/a,冷箭竹(4.3±0.23)cm/a,拐棍竹(11.4±0.57)cm/a。預(yù)測年份竹種擴(kuò)散距離=竹種擴(kuò)散距離均值×(預(yù)測年份-2015)。

2.2 模型介紹

Maxent模型是基于最大熵理論開發(fā)的軟件[19]。它根據(jù)物種現(xiàn)實分布點及分布地區(qū)的環(huán)境變量得出預(yù)測模型,再利用預(yù)測模型模擬目標(biāo)物種在目標(biāo)地區(qū)及目標(biāo)情況下的可能分布情況[20]?，F(xiàn)今大量運用于物種的生境評價與預(yù)測,并表現(xiàn)出良好的預(yù)測能力。本研究采用Maxent模型對竹類和大熊貓分布進(jìn)行預(yù)測,并根據(jù)Maxent軟件生成的ROC曲線,通過AUC值對預(yù)測結(jié)果的好壞進(jìn)行判斷。其評價標(biāo)準(zhǔn)為：AUC在0.5—0.6為失??；0.6—0.7為較差；0.7—0.8為一般；0.8—0.9為好；0.9—1為非常好[21]。本研究中大熊貓訓(xùn)練集與測試集的AUC值分別為0.915和0.894,八月竹為0.963和0.944,短錐玉山竹為0.911和0.899,拐棍竹為0.950和0.936,冷箭竹為0.885和0.873。

2.3 數(shù)據(jù)分析

大熊貓適宜生境：將大熊貓痕跡點數(shù)據(jù)和氣候、植被及其他環(huán)境因子數(shù)據(jù)納入Maxent,隨機選取75%的大熊貓痕跡點用于建立模型,其余25%用于模型的驗證,其他參數(shù)均為模型的默認(rèn)值。將模型輸出結(jié)果導(dǎo)入GIS,基于模型運算出的閾值對模型預(yù)測結(jié)果進(jìn)行重分類,即獲得大熊貓適宜生境。

竹種氣候適宜區(qū)(CSAs)：將該竹種分布點數(shù)據(jù)和氣候、植被及其他環(huán)境因子數(shù)據(jù)納入Maxent,隨機選取75%的該竹種分布點用于建立模型,其余25%用于模型的驗證,其他參數(shù)均為模型的默認(rèn)值。將模型輸出結(jié)果導(dǎo)入GIS,基于模型運算出的閾值對模型預(yù)測結(jié)果進(jìn)行重分類,即獲得該竹種氣候適宜區(qū)。

主食竹氣候適宜區(qū)：將預(yù)測的4個竹種氣候適宜區(qū)相加得到。

當(dāng)前主食竹分布區(qū)：將當(dāng)前各竹種的氣候適宜區(qū)與各竹種的野外實測點疊加,沿海拔將重疊部分用平滑曲線連接,形成各竹種的當(dāng)前分布,將各個竹種的當(dāng)前分布區(qū)相加,即為當(dāng)前主食竹分布區(qū)。由于理縣和小金縣的主要竹種不是本研究的四種竹種,故分析主食竹分布及大熊貓當(dāng)前和未來棲息地時未將兩縣計算在內(nèi)。

未來竹種分布區(qū)：將當(dāng)前各個竹種分布區(qū)以擴(kuò)散距離做緩沖區(qū)后與該竹種未來氣候適宜區(qū)疊加,重疊部分為各竹種未來分布區(qū)；再將未來各竹種分布相加即為未來主食竹分布。

當(dāng)前大熊貓棲息地：將當(dāng)前大熊貓適宜生境與當(dāng)前主食竹分布區(qū)2個圖層進(jìn)行疊加,重疊的區(qū)域即為當(dāng)前大熊貓棲息地范圍。

未來大熊貓棲息地：為未來大熊貓適宜生境與未來主食竹分布區(qū)的重疊區(qū)域。

3 結(jié)果

3.1 大熊貓適宜生境和主食竹氣候適宜區(qū)變化預(yù)測

與現(xiàn)狀相比,邛崍山系中部及南部的高海拔地區(qū)將產(chǎn)生新的大熊貓適宜生境,而山系北部及低海拔地區(qū)大熊貓適宜生境面積將減少(圖2),適宜生境面積總體減少。當(dāng)前大熊貓適宜生境面積為9219 km2,2050年與當(dāng)前面積相比減少25.7%,2070年較當(dāng)前減少37.2%。

未來大熊貓主食竹氣候適宜區(qū)面積也持續(xù)減少,但減少幅度小于大熊貓適宜生境。當(dāng)前主食竹氣候適宜區(qū)面積為12764 km2,2050年較當(dāng)前氣候適宜區(qū)面積減少3.2%,到2070年主食竹氣候適宜區(qū)面積與當(dāng)前相比僅減少4.7%。減少的區(qū)域主要集中在低海拔地區(qū),而山系中部及北部的高海拔地區(qū)將產(chǎn)生新的氣候適宜區(qū)(圖2)?？傮w來看大熊貓適宜生境及主食竹氣候適宜區(qū)均有向高海拔擴(kuò)展的趨勢。

圖2 大熊貓適宜生境及主食竹氣候適宜區(qū)較當(dāng)前變化情況Fig.2 Projected changes of suitable habitat for giant panda and CSAs for bamboos, comparing with current

3.2 主食竹分布變化預(yù)測

當(dāng)前大熊貓主食竹分布面積為7353 km2,未來主食竹分布略有減少,2050年主食竹分布面積與當(dāng)前相比減少5%,到2070年將持續(xù)減少,較當(dāng)前減少8.3%。結(jié)合主食竹氣候適宜區(qū)變化看,氣候變化下主食竹分布有向高海拔擴(kuò)張的趨勢,但由于竹類擴(kuò)散能力有限,新增分布面積很少,而低海拔地區(qū)面積減少明顯(圖3),因此主食竹分布是持續(xù)減少的。而主食竹分布區(qū)域的減少可能使得大熊貓棲息地面積減少。

3.3 大熊貓棲息地變化預(yù)測

當(dāng)前11個縣中(除小金縣和理縣)大熊貓棲息地面積有6648km2, 隨著氣候變化,大熊貓適宜生境與主食竹分布均有向高海拔擴(kuò)展的趨勢,但是由于主食竹擴(kuò)散能力的限制,使得大熊貓棲息地雖然總體上有向高海拔區(qū)域擴(kuò)展的趨勢,但新增棲息地面積遠(yuǎn)小于減少的棲息地面積,因此最終預(yù)測的棲息地面積是持續(xù)減少的。2050年大熊貓棲息地分布面積減少1194km2,占當(dāng)前棲息地面積的18%,2070年較當(dāng)前棲息地減少27.2%。未來大熊貓棲息地的擴(kuò)展區(qū)域主要位于山系南部,而山系北部及低海拔區(qū)域的大熊貓棲息地一直持續(xù)減少(圖4)。

圖4 大熊貓棲息地較當(dāng)前變化情況Fig.4 Projected changes of habitat of giant panda, comparing with current giant panda habitat

4 討論

從氣候、地形條件以及植被等條件來看,氣候變化下山系中部及南部高海拔地區(qū)將產(chǎn)生新的大熊貓適宜生境,而山系北部及低海拔地區(qū)適宜生境將減少。邛崍山系的大熊貓適宜生境在未來呈現(xiàn)出持續(xù)減少并且向高海拔擴(kuò)展的變化趨勢,這與已有研究的結(jié)論相似。但大熊貓適宜生境面積減少的比例卻差異較大。一方面是由于不同山系受氣候影響程度以及不同山系大熊貓的現(xiàn)有分布數(shù)量差異較大,造成模型的基礎(chǔ)參數(shù)不同,另一方面可能是由于不同作者對模型和情景模式選取的差異。吳建國等[7]利用CART分類和回歸樹模型,分析了氣候變化對全國大熊貓適宜生境的影響,認(rèn)為氣候變化下大熊貓當(dāng)前適宜生境的東部、東北和南部將減少,新適宜生境主要向西部一些區(qū)域擴(kuò)展,并且認(rèn)為到2051—2080年目前適宜生境將減少70%。Songer等[8]運用Maxent模型對全國大熊貓適宜生境分析,認(rèn)為到2080 年大熊貓現(xiàn)有適宜生境可能喪失60%,其中邛崍山山系在CGCM3和HadCM3兩種模式下分別減少40%和67%。Fan等[9]對秦嶺地區(qū)的預(yù)測認(rèn)為,到2100年在A2和B2情景下秦嶺地區(qū)大熊貓棲息地面積將分別減少62%和37%。Jian等[22]運用物種分布模型對除秦嶺山系外的大熊貓適宜生境分析,認(rèn)為2002—2050年大熊貓適宜生境有1.5%變?yōu)椴贿m宜,同時有2.64%的區(qū)域由不適宜變?yōu)檫m宜生境,2050—2099年6.59%適宜生境變?yōu)椴贿m宜,而新增3.43%的適宜生境。本研究中,大熊貓適宜生境2050年與當(dāng)前相比面積減少25.7%,2070年較當(dāng)前減少37.2%。

邛崍山系大熊貓主食竹氣候適宜區(qū)在未來也是持續(xù)減少的趨勢,到2070年主食竹氣候適宜區(qū)面積與當(dāng)前相比減少4.7%。主食竹分布預(yù)測結(jié)果表明2050年主食竹分布面積與當(dāng)前相比減少5%,到2070年將持續(xù)減少,較當(dāng)前減少8.3%,其減少幅度大于主食竹氣候適宜區(qū)。Tuanmu等[12]對秦嶺大熊貓棲息地的竹類研究表明,氣候變化背景下未來秦嶺大熊貓棲息地內(nèi)的3個主要竹種氣候適宜區(qū)將大幅度減少,減少比例高達(dá)59%—100%,且新產(chǎn)生的氣候適宜區(qū)都與現(xiàn)有的竹分布范圍和大熊貓棲息地相隔甚遠(yuǎn)。從本研究結(jié)果看,邛崍山系竹類的氣候適宜區(qū)有一定程度的減少,但減少幅度遠(yuǎn)小于秦嶺區(qū)域。這可能與邛崍山系地理位置偏南,氣候較溫暖濕潤有關(guān)。

綜合氣候、地形、植被以及食物資源來預(yù)測邛崍山系大熊貓棲息地的變化,2050年大熊貓棲息地分布面積較當(dāng)前減少18%,2070年較當(dāng)前棲息地減少27.2%,其減少速度介于大熊貓適宜生境和主食竹分布面積之間,減少區(qū)域主要位于山系北部和低海拔地區(qū),而山系中部的核心區(qū)域減少較少。Shen等[5]采用全國第三次(1999—2003年)和第二次(1985—1988年)大熊貓調(diào)查數(shù)據(jù),運用Maxent模型預(yù)測分析岷山大熊貓棲息地變化情況,認(rèn)為到2100年大熊貓棲息地將減少(16.3±1.2)%。Shen等[5]的研究是將竹類分布與氣候及其他環(huán)境因子并列,然而本研究認(rèn)為主食竹對大熊貓棲息地分布的影響高于其他環(huán)境因子。因此如果從食物資源的角度去考慮氣候變化對大熊貓的影響,將大熊貓棲息地變化與主食竹分布變化圖層相疊加計算,2050年大熊貓棲息地加非棲息地的有竹分布區(qū)的面積,比當(dāng)前棲息地面積大4.9%,到2070年仍較當(dāng)前棲息地面積大1.5%。也即是未來大熊貓可利用棲息地面積將大于當(dāng)前大熊貓棲息地面積,表明氣候變化對邛崍山系影響較小。

總的來看,雖然氣候變化對邛崍山系影響較小,但未來邛崍山大熊貓棲息地面積會減少,減少的區(qū)域主要位于大相嶺山系的滎經(jīng)、洪雅縣以及北部的汶川縣,并且未來棲息地的破碎化程度會加劇。因此,需要進(jìn)一步加強對這些地區(qū)的保護(hù)和監(jiān)測力度,建立長期的氣候監(jiān)測站,并加強對大熊貓及主食竹變化的動態(tài)監(jiān)測,了解大熊貓及主食竹受氣候影響的趨勢；同時加強對高山區(qū)域的保護(hù)力度,因為高山區(qū)域是未來大熊貓棲息地擴(kuò)展的區(qū)域,特別是大相嶺山系的高海拔區(qū)域。邛崍山山系中部地區(qū)是受氣候變化影響較小的區(qū)域,也是現(xiàn)今大熊貓分布的核心區(qū)域,本區(qū)域的汶川、寶興、蘆山、大邑、崇州、天全6個縣(市)有大熊貓數(shù)量504只[16],占研究區(qū)域大熊貓數(shù)量的90%以上,需進(jìn)一步加強對本區(qū)域的保護(hù),特別需要加強對臥龍、蜂桶寨、鞍子河等保護(hù)區(qū)的建設(shè),加強對棲息地的保護(hù)力度,確保核心種群的長期生存安全。

致謝：四川林業(yè)廳提供四川省第三、四次大熊貓調(diào)查資料,四川省第四次大熊貓調(diào)查隊野外調(diào)查,Shane DuBay博士和吳永杰老師幫助寫作,張曼、魏玲、王彬在野外調(diào)查和數(shù)據(jù)分析中提供幫助，特此致謝。

[1] Rosenzweig C, Karoly D, Vicarelli M, Neofotis P, Wu QG, Casassa G, Menzel A, Root T L, Estrella N, Seguin B, Tryjanowski P, Liu CZ, Rawlins S, Imeson A. Attributing physical and biological impacts to anthropogenic climate change. Nature, 2008, 453(7193): 353- 357.

[2] 劉艷萍, 申國珍, 李景文. 大熊貓棲息地質(zhì)量評價研究進(jìn)展. 廣東農(nóng)業(yè)科學(xué), 2012, 39(22): 193- 198.

[3] 胡錦矗. 今昔大熊貓. 野生動物, 1981,(1): 6- 7.

[4] 王銳婷, 范雄, 劉慶, 陳文秀. 氣候變化對四川大熊貓棲息地的影響. 高原山地氣象研究, 2010, 30(4): 57- 60.

[5] Shen G Z, Stuart L, Feng C Y, Ren G F, Liu Y P, Xu W T, Li J Q, Si X F, Xie Z Q. Climate change challenges the current conservation strategy for the giant panda. Biological Conservation, 2015,190: 43- 50.

[6] 劉艷萍. 氣候變化對岷山大熊貓及棲息地的影響[D]. 北京: 北京林業(yè)大學(xué), 2012.

[7] 吳建國, 呂佳佳. 氣候變化對大熊貓分布的潛在影響. 環(huán)境科學(xué)與技術(shù), 2009, 32(12): 168- 177.

[8] Songer M, Delion M, Biggs A, Huang Q Y. Modeling impacts of climate change on giant panda habitat. International Journal of Ecology, 2012, 2012: 108752.

[9] Fan J T, Li J S, Xia R, Hu L L, WuX P, Li G. Assessing the impact of climate change on the habitat distribution of the giant panda in the Qinling Mountains of China. Ecological Modelling, 2014, 274: 12- 20.

[10] Klanderud K. Climate change effects on species interactions in an alpine plant community. Journal of Ecology, 2005, 93(1): 127- 137.

[11] 張澤鈞, 胡錦矗, 吳華. 邛崍山系大熊貓和小熊貓生境選擇的比較. 獸類學(xué)報, 2002, 22(3): 161- 168.

[12] Tuanmu M N, Vina A, Winkler J A, Li Y, Xu W H, Ouyang Z Y, Liu J G. Climate-change impacts on understorey bamboo species and giant pandas in China′s Qinling Mountains. Nature Climate Change, 2013, 3(3): 249- 253.

[13] 周世強, 吳志容, 嚴(yán)嘯, 黃金燕, 李仁貴, 黃炎, 李德生, 張和民. 自然與人為干擾對大熊貓主食竹種群生態(tài)影響的研究進(jìn)展. 竹子研究匯刊, 2015, 34(1): 1- 9.

[14] 王學(xué)志, 徐衛(wèi)華, 歐陽志云, 劉建國, 肖燚, 陳佑平, 趙聯(lián)軍, 黃俊忠. 生態(tài)位因子分析在大熊貓(Ailuropodamelanoleuca)生境評價中的應(yīng)用. 生態(tài)學(xué)報, 2008, 28(2): 821- 828.

[15] 朱靖, 龍志. 大熊貓的興衰. 動物學(xué)報, 1983, 29(1)：93- 104.

[16] 四川省林業(yè)廳. 四川的大熊貓: 四川省第四次大熊貓調(diào)查報告. 成都：四川科學(xué)技術(shù)出版社, 2015:1- 18

[17] 馮婧. 多全球模式對中國區(qū)域氣候的模擬評估和預(yù)估[D]. 南京: 南京信息工程大學(xué), 2012.

[18] 林而達(dá), 劉穎杰. 溫室氣體排放和氣候變化新情景研究的最新進(jìn)展. 中國農(nóng)業(yè)科學(xué), 2008, 41(6): 1700- 1707.

[19] 徐衛(wèi)華, 羅翀. MAXENT模型在秦嶺川金絲猴生境評價中的應(yīng)用. 森林工程, 2010, 26(2): 1- 3,26- 26.

[20] 齊增湘, 徐衛(wèi)華, 熊興耀, 歐陽志云, 鄭華, 甘德欣. 基于MAXENT模型的秦嶺山系黑熊潛在生境評價. 生物多樣性, 2011, 19(3): 343- 352.

[21] 劉振生, 高惠, 滕麗微, 蘇云, 王曉勤, 孔芳毅. 基于MAXENT模型的賀蘭山巖羊生境適宜性評價. 生態(tài)學(xué)報, 2013, 33(22): 7243- 7249.

[22] Jian J, Jiang H, Jiang Z S, Zhou G M, Yu S Q, Peng S L, Liu S Y, Wang J X. Predicting giant panda habitat with climate data and calculated habitat suitability index (HSI) map. Meteorological Applications, 2014, 21(2): 210- 217.

[23] 冉江洪, 曾宗永, 劉世昌, 王鴻加. 四川大相嶺大熊貓種群及棲息地調(diào)查. 四川大學(xué)學(xué)報自然科學(xué)版, 2006, 43(4): 889- 893.

[24] 張文廣, 唐中海, 齊敦武, 胡遠(yuǎn)滿, 胡錦矗. 大相嶺北坡大熊貓生境適宜性評價. 獸類學(xué)報, 2007, 27(2): 146- 152.

[25] 張黎明. 邛崍山系大熊貓棲息地狀況的初步研究. 四川林業(yè)科技, 1993, 14(2): 74- 76.

Climate-change impacts on bamboo distribution and giant panda habitat in

Qionglai mountains

YAN Tingting, RAN Jianghong*, ZHAO Chenhao, ZHONG Xue, LIANG Minyi

KeyLaboratoryofBio-ResourcesandEco-environmentofMinistryEducation,CollegeofLifeScience,SichuanUniversity.Chengdu610064,China

The effects of climate change on biodiversity are at the forefront of contemporary scientific research, especially regarding rare and endangered species. As an old, highly specialized and extremely rare species, the effect of global climate change on the giant panda(Ailuropodamelanoleuca) has received much attention. Due to climate change, giant panda habitats are generally becoming warmer and drier, reducing the total suitable area thereby giant pandas are migrating to higher latitudes and altitudes. Climate change not only directly affect species distribution, but also indirectly alters the size and distribution of herbivore populations by influencing the growth of staple food plants. Bamboo comprises 99% of the giant panda′s diet. This strong dependence couples the distribution of the giant panda with that of bamboo growth. Based on current locality and distributional data of giant pandas and bamboo species, respectively, from field surveys, we predicted distributional shifts of bamboo and giant panda habitats for the years 2050 and 2070 using a Maxent model with topography factors (altitude, slope and aspect), climate change data, river, road, resident community and land cover. We obtained current and future climate projections (for the years 2050 and 2070) from the WorldClim database at 30′ resolution. The model was constrained by the conditions outlined in the representative concentration pathway 8.5 of the IPCC AR5. We assessed the accuracy of predictions of the area under the ROC curve (AUC), establishing the acceptable predictions of Maxent. The results show that the climatically suitable areas of bamboo and habitat for giant pandas might decrease by 4.7% and 37.2% in 2070, respectively. The northern and lower elevation limits of current suitable habitat for giant pandas would decrease, whereas suitable habitat would primarily expand to higher elevations. Lower elevations of climatically suitable areas of bamboo would also decrease; however the distribution of bamboo would primarily expand to higher elevations. Therefore the area of giant panda habitat will decrease by 27.2% and giant pandas would move to higher altitudes. Although bamboo will spread to higher elevations, and the total proportion of bamboo distribution will decrease by 8.3% in 2070, there will be sufficient food resources for nearly all giant pandas in this region. This suggests that climate change has limited impacts on giant pandas in the Qionglai Mountains. The southern area of the Qionglai Mountains and low-altitude regions are more sensitive to climate change, which may require additional management, research, and conservation in the future. The main habitats of giant panda (Wenchuan, Baoxing, Lushan, Dayi, Chongzhou and Tianquan counties)will also required additional focus in future studies.

climate change; giant panda; qionglai mountains; bamboo; conservation

10.5846/stxb201512032424

四川省環(huán)境保護(hù)科技計劃資助項目

2015- 12- 03; 網(wǎng)絡(luò)出版日期:2016- 08- 30

晏婷婷,冉江洪,趙晨皓,鐘雪,梁敏儀.氣候變化對邛崍山系大熊貓主食竹和棲息地分布的影響.生態(tài)學(xué)報,2017,37(7):2360- 2367.

Yan T T, Ran J H, Zhao C H, Zhong X, Liang M Y.Climate-change impacts on bamboo distribution and giant panda habitat in Qionglai mountains.Acta Ecologica Sinica,2017,37(7):2360- 2367.

*通訊作者Corresponding author.E-mail: rjhong-01@163.com