氣候變暖情景下東北丹頂鶴繁殖地植被生境變化分析

李明陽 ，張稱意 ，吳 軍 ，徐 婷

氣候變暖情景下東北丹頂鶴繁殖地植被生境變化分析

李明陽1，張稱意2，吳 軍3，徐 婷1

（1. 南京林業大學 森林資源與環境學院, 江蘇 南京 210037；2. 國家氣候中心，北京 100081；3. 環保部南京環境科學研究所，江蘇 南京 210042）

明晰氣候變化情景下瀕危物種生境變化的趨勢及其驅動因素，是進行瀕危物種種群數量預測和制定野生動物保護規劃的前提。以東北地區丹頂鶴繁殖地主要植被類型蘆葦為主要研究對象，以東北25個丹頂鶴自然保護區邊界數據、1999～2003年保護區蘆葦空間分布數據、國家氣候中心2005年與2030年及2050年模擬數據為主要信息源，采用氣候包絡法Bioclim、Domain模型，在模型精度評價的基礎上，進行2005～2050年間丹頂鶴繁殖地蘆葦的動態變化分析、空間格局分析。研究表明：(1) 與2005年相比，2030年、2050年研究區域中部、東部的蘆葦生境將大面積消失；(2) 與2005年相比， 2030年質量好、質量中等的蘆葦生境面積比例分別下降了3.812‰、6.617‰，2050年則分別下降了7.121‰、8.057‰，高質量蘆葦生境空間聚集性減弱，破碎化趨勢增強；(3) 2030年、2050年蘆葦地理分布中心將由東部分別向西北、東北方向偏移，與2005年相比，2050年地理分布中心向西北方向平均每年偏移了19.756 km。

丹頂鶴；氣候變暖；生境分析；東北地區；氣候包絡法

野生動物保護專家認為，保護瀕危物種的最好方法是保護其棲息環境[1-2]。確定瀕危物種的地理空間分布、明晰影響瀕危物種生境變化的主要驅動因素，是進行瀕危物種種群數量預測和制定野生動物保護規劃的前提[3]。丹頂鶴Grus japonensis是我國一級保護鳥類，2001 年被列入《亞洲鳥類紅皮書》，屬于世界級瀕危物種[4]。丹頂鶴在濕地生態系統中處于食物鏈頂層，是濕地環境動態變化最敏感和最明顯的生物指示物種[5]。東北地區是丹頂鶴最重要的集中繁殖棲息地，也是氣候變暖最顯著的區域之一[6]。蘆葦Phragmites australis、小葉樟Deyeuxia langsdorffii和多種苔草屬Carex沼澤植被，是丹頂鶴的隱蔽所和筑巢的主要原料。研究氣候變化情景下丹頂鶴繁殖地沼澤植被生境動態，對于丹頂鶴種群保護和生境管理具有重要的理論意義和實際參考價值。

根據美國大氣研究國家中心(NCAR)(http://www.ncar.ucar.edu/)的研究成果，隨著溫室氣體排放量的增加和臭氧層的破壞，地球表面溫度呈上升趨勢。從1800年開始，地球表面溫度上升了0.7℃。從20世紀70年代至今,這種上升趨勢呈現加速狀態。據NCAR預測,到2010和2050年，地球表面溫度將分別上升0.6℃～2.5℃、1.4℃～3.8℃。中國國家氣候中心的研究表明，過去100年，我國的平均氣溫上升了0.5℃～0.8℃，比全球略高。研究表明，近百年來東北年平均溫度增加了1.43℃，是全球增溫率的2倍，全國增溫率的3倍[7]。隨著全球氣候變化，東北丹頂鶴繁殖地沼澤植被生境也會發生動態變化。

近年來，關于氣候變化對野生動物種群空間分布格局影響，以哺乳動物的研究較多，而鳥類方面的研究報道相對較少，而且以文獻綜述為主[8]。氣候變化對丹頂鶴繁殖地沼澤植被的影響研究，重點集中在過去數十年氣候變化對植被物候特征的影響方面[9-10]，而關于未來氣候變化情景下丹頂鶴繁殖地植被生境空間分布格局動態方面相對較少。由于大區域及長時間尺度上土地利用、氣候變化、物種分布綜合監測數據的缺乏，目前國際上尚缺乏系統的、操作性強的評估工具，氣候變暖情境下丹頂鶴繁殖地沼澤植被生境變化的量化分析國內研究成果，尚未見到公開報道。

本文以東北地區丹頂鶴繁殖地主要沼澤植被類型蘆葦為主要研究對象，以東北25個丹頂鶴頂鶴自然保護區邊界數據、1999～2003年保護區蘆葦空間分布數據、國家氣候中心2005年、2030年、2050年模擬數據為主要信息源，采用生物氣候包絡模型，在模型精度評價的基礎上，進行2005～2050年間丹頂鶴繁殖地植被的動態變化分析、空間格局分析，以期總結出氣候變化對大區域尺度上珍稀物種棲息生境影響評估的一般方法和基本途徑，從而對我國丹頂鶴繁殖地生境保護提供科學依據。

1 研究區域基本概況

東北地區（包括遼寧、吉林、黑龍江三省，以及內蒙古自治區東部呼倫貝爾盟、興安盟、赤峰市、通遼市、錫林郭勒盟東經115°以東部分），位于中國的東北部，東與北朝鮮接壤、北與俄羅斯一衣帶水，地理范圍為東經118°0′0″～134°46′29″E，北緯 38°42′57″～ 53°33′38″N，土地總面積135.61萬km2，占全國的14.12%，2010年人口11 022萬人，占全國的8.2%。受緯度、海陸位置、地勢等因素的影響，東北地區跨寒溫帶、溫帶、暖溫帶的濕潤、半濕潤地區，溫度差異明顯，主要氣候特征是冬季嚴寒，夏季冷濕。遼寧南部＞10℃積溫可達3 400℃以上，中部大部分地區＞10℃積溫在2 400 ℃～3 400 ℃之間，降水量自東向西遞減，長白山東南側鴨綠江流域年降水量可達1 000 mm以上,長白山西側為600～700 mm，松嫩平原400～600 mm，大興安嶺東側500 mm，西側為350 mm。東北地形總的特點是東、北、西三面環山，中部為丘陵和遼闊的平原。

東北水系比較發達, 是我國濕地最多的地區，主要分布在三江平原、大小興安嶺、長白山地區、松嫩平原、遼河下游平原沿海河口和東部山地，包括了中國濕地的濱海、江河、湖泊、沼澤、人工濕地5大類，具有蘚類沼澤、草本沼澤、灌叢沼澤、森林沼澤、內陸鹽沼、水庫等14 種濕地類型。自20 世紀中期中國東北地區出現了持續而顯著的增溫現象。20 世紀80～90 年代的平均氣溫與20世紀60～70 年代相比己經上升了1.0℃～2.5℃, 而且在未來幾十年或更長的時間將繼續升溫。隨著氣溫的升高，濕地蒸騰蒸發量和降水變率增大，極端降水事件頻率增加，導致大片的蘆葦、苔草濕地退化甚至消失，大幅增加了棲息在東北沼澤濕地的丹頂鶴、灰鶴Grus grus等珍稀水禽種群的滅絕風險。

2 數據來源及研究方法

2.1 數據來源

論文所用的主要數據有：① 中國氣象局國家氣候中心提供的溫室氣體中等排放情景下（SRESA1B）生成的研究區域2005年、2030年、2050年跨學科研究全球模式驅動下區域氣候模型（MIROC-RegCM）模式模擬數據（http://www.climatechange-data.cn），空間分辨率為1 000 m，包括月平均氣溫（℃）、月平均最高溫度（℃）、月平均最低溫度（℃）、月平均降水（mm）、月平均凈輻射（W/m2）5個氣候因子；②由全球土地覆蓋數據中心（GLCF,http://glcf.umiacs.umd.edu/）提供的研究區域1999～2003年Landsat ETM遙感數據，多光譜波段空間分辨率為30m，全色波段空間分辨率為15m；③ 國家環保部南京環境科學研究所提供的東北地區25個丹頂鶴自然保護區邊界矢量文件；④ 由國家基礎地理信息系統（http://nfgis.nsdi.gov.cn/）提供的研究區域行政區劃基礎數據，包括省、市、縣三級行政區劃界線。

2.2 蘆葦空間分布數據的生成

在25個丹頂鶴保護區中，選取具有1∶50000地形圖的20個自然保護區，分別以1999-2003年Landsat ETM遙感數據為主要信息源，以保護區地形圖為地理參照，運用遙感影像處理軟件ERDAS IMAGINE 9.2，在對 Landsat ETM遙感圖像進行幾何精校正、輻射增強、空間子集運算等預處理的基礎上，通過監督分類，將保護區的土地利用類型分為森林、農田、灌叢、建筑用地、水域、沼澤草地、裸地7類，每個自然保護區的總體分類精度均在85%以上。在東北丹頂鶴繁殖地沼澤濕地中，沼澤草地的主要植被類型是蘆葦。除集中分布的葦叢外，散生的蘆葦往往與多種苔草屬植被混生在一起。將沼澤草地從監督分類柵格文件中提取出來，通過ArcGIS平臺上Spatial Analyst模塊的“柵格轉矢量”工具，轉化為Point Shape 文件。在此基礎上，將沼澤草地植被矢量點與研究區域1999～2003年Landsat ETM彩色合成遙感影像疊加。為消除混生的苔草屬草本植物對蘆葦空間分布點的影響，借助于各個自然保護區的地形圖、土地利用類型圖，通過目視解譯，剔除散生的沼澤草地矢量點，在位于沼澤中間、沼澤邊緣且集中分布的沼澤草地中選取典型的草本植被空間分布點477個，構成1999-2003年蘆葦空間分布數據庫，以此作為東北丹頂鶴繁殖地沼澤植被蘆葦生態位建模的數據基礎。

2.3 氣候包絡法生態位建模

生物氣候模型Bioclim、領域模型Domain都只需要物種分布的存在數據就可以對物種的潛在分布區進行估測。由于在生物多樣性保護實踐中物種不存在空間數據很難獲得，因此這2種模型在瀕危物種潛在生境空間建模中得到較為廣泛應用[11]。 Bioclim和Domain模型都只考慮氣候因素對物種分布的影響, 這種方法被稱作氣候包絡法(Climate envelope)[12]。Bioclim模型原理主要是利用已有的物種分布資料和環境數據產生以生態位為基礎的物種生態需求，探索物種已知分布區的環境特征與研究區域的非隨機關系。它采用分布百分比算法(Percentile distribution)對分布區域中每個網格的每個環境變量進行多元單邊分析, 計算每個網格的環境適生度[13]。如果網格中所有的環境變量都在物種分布的氣候包絡之內，就認為這個地方是適合該物種生存的。Domain模型是在建立Gower矩陣基礎上, 將分布點的環境變量與它周圍地區的環境變量進行相似性比較，來對物種的分布區進行預測，相似性的最大值為100，相似性值越高表示與已知分布區的氣候相似性越大。

2.4 蘆葦生境氣候因子的提取

研究表明，在國家和大區域尺度上，影響植被空間分布的生態環境因子主要是溫度和降水等氣候因子[12]。預測物種分布的Bioclim模型和Domain模型是基于國際馬鈴薯中心開發的生物地理信息系統軟件Diva-GIS。在該軟件安裝目錄的Environ子文件夾下，已經內嵌了1950～2000年來自世界各地氣象站的氣候信息采用插值法生成的全球氣候數據，該數據由WORLDCLIM ( http: / /worldclim. org)提供，空間分辨率為2.5′（赤道處1′= 1852 m）。直接采用WORLDCLIM 的全球氣候數據進行東北地區丹頂鶴繁殖地蘆葦生境空間建模，存在著數據時相性差、空間分辨率低、區域針對性弱的弊端。在論文中，利用Diva-GIS軟件的氣候數據庫開放功能，將研究區域2005年模式模擬數據中的月平均最高溫度（℃）、月平均最低溫度（℃）、月平均降水（mm）3個因子轉換成CLM氣候文件格式，導入Environ子文件夾下，即可由軟件自動生成包括年均溫(Annual mean temperature)、平均周溫度變化范圍(Mean diurnal temperature range)、平均年溫度變化范圍(Mean annual temperature range)、極端最高溫(Maximum temperature of the warmest month)在內的19個氣候因子，以此作為蘆葦生境空間建模的環境因子。

2.5 氣候變暖情景分析

區域氣候模型模式模擬數據對比分析表明，2005年、2030年、2050年3個時段的年平均溫度分別為2.752℃、4.418℃、5.794℃。與2005年相比，2030年、2050年的平均溫度上升了1.666℃、3.042℃。在Bioclim、Dommain 2個氣候包絡模型性能評價的基礎上，選取預測精度較高的模型，在Diva-GIS軟件生境建模菜單下的“氣候變化”分析模塊支持下，利用與2005年相同的蘆葦空間分布數據庫，分別年平均氣溫升高1.666℃、3.042℃ 2種氣候變暖情景，進行氣候變化情景下蘆葦潛在生境變化分析。

3 結果與分析

3.1 2005年蘆葦生境空間分布分析

將1999-2003年丹頂鶴生境植被蘆葦的空間分布數據分為2部分：70%（334個）樣本點用于建模，30%（143）個用于模型驗證。采用Diva-GIS軟件提供的受試者工作特征曲線（ROC）下面積（AUC）計算功能來評價2個氣候包絡預測模型的性能。AUC是一個與參考閾值無關的統計量，通過百分比的方法計算提升置信區間(bootstrap confidence interval)，來評價模型區分存在點和背景點的診斷性能[14]。一般認為AUC值為0.5～0.7時模型診斷價值較低；為0.7～0. 9時診斷價值中等；大于0.9時診斷價值優秀。通過ArcGIS中ArcToolbox數據轉換功能， 分別將Diva-GIS軟件輸出的Bioclim、Domain模型預測結果從ASCII文件轉換為Raster文件格式，生成基于Bioclim模型的2005年丹頂鶴繁殖地植被適應性概率空間分布圖（圖1）、基于Domain模型的2005年丹頂鶴繁殖地植被適應性概率空間分布圖（圖2）。

圖1 基于Bioclim模型的2005年蘆葦適應性概率Fig. 1 Bioclim-based suitable probability of reed habitat in 2005

圖2 基于Domain模型的2005年蘆葦生境適應性概率Fig. 2 Domain-based suitable probability of reed habitat in 2005

Bioclim、Domain模型預測精度AUC分別為0.904、0.489，表明Bioclim模型的預測精度超過Domain，達到優良性能，因此采用Bioclim模型進行蘆葦生境空間建模、氣候變化情景分析。從圖1可以看出，丹頂鶴繁殖地沼澤植被蘆葦生境主要集中在5大區域，其中2個大片區域分別位于黑龍江省東部的佳木斯市、雙鴨山市、雞西市（簡稱東部片區），黑龍江、吉林、內蒙古3省（區）交界處的齊齊哈爾、大慶市、白城市、興安盟東南部、通遼市北部（簡稱中部片區），3個小片區域分布在錫林郭勒盟東北部（簡稱錫盟片區）、呼倫貝爾盟西南部（簡稱呼盟片區）、黑龍江省牡丹江市（簡稱牡丹江片區）。按照蘆葦生境適應性概率從高到低的順序排列，排在前10位的縣（區、縣級市）全部位于黑龍江省境內，它們依次是：寶清縣（0.618）、友誼縣（0.603）、集賢縣（0.516）、齊齊哈爾市轄區(0.455)、杜爾伯特縣（0.450）、龍江縣(0.438)、虎林縣（0.412）、泰來縣（0.395）、七臺河市區（0.372）、樺川縣（0.366）。

3.2 氣候變化情景下植被生境變化分析

按照2.5中的方法，采用Bioclim模型，進行2030年、2050年氣候變暖情境下蘆葦生境動態分析，生成2030年蘆葦生境適應性概率空間分布圖（圖3）、2050年蘆葦生境適應性概率空間分布圖（圖4）。對照圖2、圖3可以看出，與2005年相比，2030年平均溫度從2.752℃上升到4.418℃，包括扎龍濕地在內的中部片區生境大面積消失。對照圖3、圖4可以看出，隨著2050年平均溫度從2030年的4.418℃上升到5.794 ℃，降水量較為豐富的東部片區、牡丹江片區生境也大面積消失，只剩下了位于研究區西部的錫盟片區、呼盟片區2小塊生境。對比圖2、圖3、圖4可以看出，2005～2050年，東北丹頂鶴繁殖地沼澤植被蘆葦的生境在空間分布上呈現出從中部到東部大面積消失的過程。從變化的程度來看，太陽輻射適中、降水量較少的中部片區的蘆葦生境變化最為劇烈，降水量較為豐富的東部片區次之，降水稀少、太陽輻射強烈的西部最為穩定。這一變化趨勢可能與研究對象蘆葦喜光、耐濕的生物學特征有關。

圖3 2030年蘆葦生境適應性概率空間分布Fig. 3 Suitable probability of reed habitat in 2030

圖4 2050年蘆葦生境適應性概率空間分布Fig.4 Suitable probability of reed habitat in 2050

為進一步分析2005～2050年東北沼澤植被蘆葦生境變化，在生成2030年、2050年適應性概率空間分布圖基礎上，進行不同質量等級蘆葦生境面積統計分析、蘆葦生境空間分布格局分析。為便于進行不同等級蘆葦生境面積比例分析，按照適應性概率的大小將東北地區蘆葦生境分為3級：好（＞0.7）、中（0.5～0.7）、差（＜0.5）。將適生概率＞0.7的像素點作為蘆葦高質量生境點，采用地統計學的Moron I 系數、Z 兩個指標進行蘆葦空間聚集性分析。Moran I 系數是一種最常用的自相關系數,系數的取值在-1 和1之間:小于0 表示負相關，等于0 表示不相關，大于0 表示正相關[15]。在統計學中,Z 是測量標準偏差的一個統計量，等于偏離適應性概率平均值的標準偏差的倍數。當可信度P=0.95、Z 位于區間范圍[-1.96,1.96]時，表征了一種統計變量隨機分布的空間格局。當Z 值落在區間范圍之外，則表示統計變量呈現出離散或聚集的分布格局[16]。

表1 2005-2050年蘆葦生境變化分析Table 1 Habitat change analysis of reed from 2005 to 2050

從表1可以看出，從生境質量來看，2005～2050年，蘆葦生境質量呈迅速下降趨勢，質量好的生境面積比例、質量中等的生境面積比例呈顯著減少趨勢，質量差的生境面積比例呈緩慢增加趨勢。與2005年比，2030年、2050年質量好的生境面積比例分別下降了3.812‰、7.121‰，質量中等的生境面積比例分別減少了6.617‰、8.057‰，質量差的生境面積比例則分別增加了10.429‰、15.178‰。從空間分布來看，2005年、2030年、2050年3個時段，Moran I均大于0，Z均大于1.96，表明2005～2050年，沼澤植被呈空間聚集性分布。從表1可以看出，與2005年相比，2030年、2050年，Moran I、Z均呈大幅度下降趨勢，表明隨著氣溫升高、降水量減少，高質量沼澤植被生境面積減少、空間聚集性減弱、破碎化趨勢增強。

3.3 高質量植被生境空間分布中心變化分析

提取適應性概率＞0.7 的像素點作為高質量蘆葦植被生境，利用ArcGIS的空間統計模塊，生成3個時段的高質量蘆葦生境地理分布中心，并與研究地區邊界矢量文件疊加，生成2005～2050年高質量蘆葦生境地理分布中心變化圖（圖5）。

圖5 2005～2050年蘆葦地理分布中心Fig.5 Geographical distribution center of marsh vegetation from 2005 to 2050

從圖5可以看出，2005年沼澤植被的地理分布中心位于黑龍江省雙鴨山市集賢縣東南，2030年分布中心則向西偏北方向偏移到內蒙古自治區呼倫貝爾盟鄂溫克自治縣北部。與2005年相比，2030年地理分布中心向西偏北方向偏移了900.3 km，平均每年偏移36.012 km。與2030年相比，2050年沼澤植被的地理分布中心從內蒙古呼倫貝爾盟鄂溫克自治縣境內向東偏北方向移到了呼倫貝爾盟海拉爾市境內，東偏北方向偏移了35.3km，平均每年偏移1.765 km。與2005年相比，2050年的沼澤植被地理分布中心向西偏北方向偏移了889.031 km，平均每年向西偏北方向偏移了19.756 km。

2005-2050年，東北丹頂鶴繁殖地蘆葦地理分布中心呈現出由東部向西北方向、由西北方向向東北方向偏移的復雜過程。從圖5可以看出，研究地區的西部為光照較為強烈的內蒙古高原，陽光充足，東部靠近日本海，水熱條件優越。作為一種喜光、耐濕的沼澤植被，在氣候變化影響情境下，蘆葦的地理分布中心趨向于向環境條件改善的區域偏移。

4 結論與討論

在溫室氣體中等排放情景下，全球模式驅動下區域氣候模型模式模擬數據分析表明，2005～2050年，隨著氣溫的升高、降水量的減少，蘆葦的高質量生境比例下降、空間聚集性減弱、分布中心向西北方向偏移。一些野生動物保護專家的研究表明，草本沼澤蘆葦的面積大小及其空間分布是影響丹頂鶴巢址選擇的主要生態因子 。因此，我們可以推定，在氣候變暖情景下，位于東北地區中部、東部的丹頂鶴繁殖地潛在生境將大面積消失，高質量生境面積縮小、空間聚集性減弱，丹頂鶴繁殖地分布中心向生態環境改善的西北方向偏移。

在東北區域尺度上，缺乏帶有地理坐標的丹頂鶴繁殖地沼澤植被蘆葦空間分布數據。在論文中，研究者借助于Landsat ETM遙感數據、丹頂鶴自然保護區邊界矢量文件，通過一系列遙感圖像處理、GIS空間分析、目視解譯等數據預處理過程，篩選出1999～2003年蘆葦空間分布數據。在提取的蘆葦空間分布數據中，難免存在著少量的小葉樟、苔草等植被類型。不同沼澤植被類型對氣候因子的反應機制有所不同，在一定程度上會影響研究結論科學性。其次，研究區域幅員遼闊，面積高達135.61萬km2，而一景Landsat ETM遙感影像的空間覆蓋范圍只有170 km×183 km=31 110 km2，很難獲取20個自然保護區同一時段遙感數據，遙感影像時相的不一致性也會影響分析結果的準確性。第三，蘆葦空間分布數據與生態位建模氣候數據時段的不一致性，也會影響分析結果的精度。

論文的研究結果是建立在溫室氣體中等排放情景下，跨學科研究全球模式驅動下區域氣候模型模式模擬數據基礎之上的[18]。溫室氣體中等排放情景假定：世界人口趨于穩定，高新技術廣泛應用，全球合作，經濟快速發展，各種能源平衡發展，在此前提下CO2排放濃度加倍。溫室氣體排放的情景不同、采用的全球驅動模式不同，研究區域2005～2050年的氣候變化情景有所不同。由于目前對于未來氣候變化的預估還存在很大的不確定性，土地管理部門、野生動物保護者在引用論文的研究成果時，尚需謹慎。

[1] Ferrier S. Mapping spatial pattern in biodiversity for regional conservation planning: where to from here[J]. Systems Biology,2002, 51: 331-363.

[2] Rushton S P, Ormerod S J, Kerby G. New paradigms for modeling species distributions[J]. Journal of Applied Ecology, 2004, 41:193-200.

[3] Ricklefs R. A comprehensive framework for global patterns in biodiversity [J]. Ecology Letters, 2004, 7: 1-15.

[4] 鄒紅菲，吳慶明.扎龍自然保護區丹頂鶴 (Grus japonensis)巢的內分布型及巢域[J]. 生態學報，2009，29(4): 1710-1718.

[5] 馬逸清，李曉民.丹頂鶴研究[M].上海：上海科技教育出版社，2002：17-18.

[6] 沃曉棠，黃智超，孫彥坤，等. 扎龍濕地近 50 年溫度和降水年際變化分析[J].東北林業大學學報，2010，41(5)：61-65.

[7] 潘響亮，鄧 偉，張道勇，等. 東北地區濕地的水文景觀分類及其對氣候變化的脆弱性[J].環境科學研究，2003，16（1）：14-18.

[8] 孫全輝，張正旺.氣候變暖對我國鳥類分布的影響[J].動物學雜志, 2000, 35(6): 45-48.

[9] 趙串串，楊曉陽，張鳳臣，等.氣候變化對濕地植被生物量影響分析——以三江源區為例[J].干旱區資源與環境，2008,(9):88—91.

[10] 李榮平，周 莉，謝艷兵，等. 1957～2004年盤錦蘆葦濕地的氣候變化特征分析[J]. 氣象與環境學報，2006,22(4):13-17.

[11] 徐曉婷，楊 永，王利松. 白豆杉的地理分布及潛在分布區估計[J].植物生態學報，2008，32 (5)：1134-1145.

[12] Hijmans R J, Graham C H. The ability of climate envelope models to predict the effect of climate change on species distributions [J].Global Change Biology, 2006, 12, 2272-2281.

[12] 胥 勇，周平陽，王 揚，等.應用BIOCLIM生態位模型對松材線蟲傳入云南的風險評估[J].云南農業大學學報，2008，23(6):746-753.

[14] Hanley J A, McNeil B J. The meaning and use of the area under a Receiver Operating Characteristic (ROC) curve [J]. Radiology,1982, 143: 29-36

[15] 王勁峰，廖一蘭，劉 鑫.空間數據分析教程[M].北京：科學出版社，2010.

[16] Arthur G, Ord J K. The analysis of spatial association by use of distance statistics [J]. Geographical Analysis,1992, 24(3):189-206.

[17] 劉學昌.三江平原地區丹頂鶴（Grus japonensis）巢址選擇[D]. 哈爾濱：東北林業大學，2006：13-45.

[18] 石 英. RegCM3對21世紀中國區域氣候變化的高分辨率數值模擬[D].北京: 中國科學院大氣物理研究所，2010:118-120.

Vegetation dynamics analysis in northeastern breeding habitat of Grus japonensis under scenaries of climate warming

LI Ming-yang1, ZHANG Cheng-yi2, WU Jun3, XU Ting1
(1.College of Forest Resources and Environment, Nanjing Forestry University, Nanjing 210037; Jiangsu, China; 2.National Climate Center of China, Beijing 100081,China; 3.Nanjing Institute of Environmental Sciences, Ministry of Environmental Protection, Nanjing 210042, Jiangsu, China)

The detailed knowledge of habitat change trends and drivers of endangered species under scenaries of climate change is the premise to predict endangered species population and the plan making for wildlife conservation. Phragmites australis which is the major vegetation type in northeastern breeding marshes of Grus japonensis was chosen as the research object, while the simulated data in 2005,2030,2050 from the National Climate Center, together with the boundary vector data of 25 red-crowned nature reserves were collected as the main information source to build climate envelop models of Bioclim and Domain to analyze dynamics of reed spatial distribution and spatial pattern from 2005 to 2050. Research results show that: (1) compared with in 2005, the large reed habitat in central part and eastern part of research region will disappear in 2030 and 2050 respectively; (2)compared with 2005, the area proportion of good and medium quality habitat in 2030 and 2050 will decrease by 3.812‰ and 6.617‰ , 7.121‰ and 8.057‰ respectively, with less aggregated and more fragmented spatial distribution of high quality reed potential habitat; (3) under impact of regional climate warming, compared with in 2005, the geographical distribution center of reed habitat in 2030 and 2050 will move from eastern part to northwest and northeast part of research region respectively, with annual average moving speed of 19.756 km to word to northwest direction.

Grus japonensis (red-crowned crane); climate warming; habitat analysis; northeast of China; climate envelop model

2012-01-24

環保部重大公益性行業科研專項“氣候變化對東北野生動植物影響的評估技術研究”（200909070）

李明陽(1967—)，男，河南三門峽人，教授，博士，從事森林資源與生態環境監測、風景林調查規劃、3S應用研究；

E-mail：lmy196727@126.com

S718.5

A

1673-923X(2012)08-0058-06

[本文編校：吳 彬]