京津冀地區(qū)物種多樣性保護(hù)優(yōu)先區(qū)識別研究

邢韶華,周 鑫,2,劉云強(qiáng),袁 秀

1 北京林業(yè)大學(xué)生態(tài)與自然保護(hù)學(xué)院, 北京 100083 2 拉薩市曲水縣聶當(dāng)鄉(xiāng)人民政府, 拉薩 850600 3 遼寧農(nóng)業(yè)職業(yè)技術(shù)學(xué)院園林系，營口 115009 4 中國科學(xué)院科技戰(zhàn)略咨詢研究院, 北京 100190

生物多樣性的喪失早已引起了國際社會的普遍關(guān)注,開展生物多樣性保護(hù)是一項(xiàng)長期艱巨的任務(wù)。生物多樣性保護(hù)工作者從基因、物種到群落、甚至生態(tài)系統(tǒng)的各個水平,從自然保護(hù)區(qū)、區(qū)域地帶到全球的各個尺度開展了多種形式的研究與實(shí)踐。Myers在分析熱帶雨林受威脅程度的基礎(chǔ)上,首先提出了熱點(diǎn)地區(qū)的概念[1],并根據(jù)物種特有程度和受威脅程度提出了全球范圍內(nèi)的25個生物多樣性熱點(diǎn)地區(qū)。此后生物多樣性熱點(diǎn)地區(qū)的評估在不同尺度上廣泛展開,評估指標(biāo)、評估方法也更加多樣。總體上來講,生物多樣性熱點(diǎn)地區(qū)的評估主要基于物種多樣性和生態(tài)系統(tǒng)的多樣性,如,李迪強(qiáng)等將鳥類作為生物多樣性豐富度的指示,將森林分布作為生態(tài)系統(tǒng)的指示,評價了尤溪縣生物多樣性的優(yōu)先保護(hù)地區(qū)[2];高玉平通過分析上海市植物及脊椎動物的分布熱點(diǎn),選擇出上海市物種多樣性的優(yōu)先保護(hù)地區(qū)[3]。在生態(tài)系統(tǒng)多樣性層面,除了選用最多的植被外,有研究者又選擇了一些其他的指標(biāo),如特殊的氣候地理與土壤特征、特殊生境等評價指標(biāo),進(jìn)行優(yōu)先保護(hù)地區(qū)的評價[4]。生物多樣性的優(yōu)先保護(hù)地區(qū)(熱點(diǎn)地區(qū))的形成是該地區(qū)的氣候、地形地貌以及人為活動等多種因素綜合作用的結(jié)果。在一個氣候條件相對一致的區(qū)域里,地形地貌是影響物種分布的重要因素。從全球范圍來看,生物多樣性高的熱點(diǎn)地區(qū)通常分布在地形地貌復(fù)雜、地質(zhì)條件多樣的區(qū)域[5- 7]；并且,地形復(fù)雜區(qū)域通常有更多特有種、狹域種和受威脅種,具有高的保護(hù)價值[8]。地形要素是最基本的自然地形單元,影響著地表物質(zhì)遷移與能量轉(zhuǎn)換、地表過程以及景觀格局等[9]。地形是水熱條件外的一個非常重要的、影響生物多樣性分布格局的因素。有研究表明,至少在陸地生態(tài)系統(tǒng),地形異質(zhì)性可以顯著解釋氣候的剩余變異(residual variations)[10]。地形因子的作用在于影響物質(zhì)和能量的再分配[11],地形是多種環(huán)境因子的復(fù)合,水分、溫度和土壤養(yǎng)分的變化集中體現(xiàn)在地形梯度的變化上[12]。京津冀地區(qū)作為 “我國經(jīng)濟(jì)最具活力、開放程度最高、創(chuàng)新能力最強(qiáng)、吸納外來人口最多的地區(qū)”[13]之一,同時也“面臨水土資源短缺和生態(tài)環(huán)境壓力加大”[13]等問題,因此更應(yīng)該根據(jù)《生態(tài)文明體制改革總體方案》的要求,“健全國土空間用途管制制度”[14],“將用途管制擴(kuò)大到所有自然生態(tài)空間,劃定并嚴(yán)守生態(tài)紅線,嚴(yán)禁任意改變用途,防止不合理開發(fā)建設(shè)活動對生態(tài)紅線的破壞。”[14]。本文擬通過建立典型自然保護(hù)區(qū)地形異質(zhì)性與生物多樣性關(guān)系模型,對京津冀地區(qū)物種多樣性保護(hù)優(yōu)先區(qū)進(jìn)行預(yù)測與評估,以期為京津冀地區(qū)深化生態(tài)文明體制改革、國家公園布局提供技術(shù)參考。

1 研究區(qū)域概況

京津冀地區(qū)地處華北平原北部,包括北京市、天津市以及河北省(北緯 36°03′—42°40′,東經(jīng) 113°27′— 119°50′),總面積約21.72萬km2,北靠燕山山脈,南面華北平原,西倚太行山脈,東臨渤海,東部和南部地形較為平坦,屬于華北平原,是典型的農(nóng)業(yè)區(qū),北部和西北部地勢相對較高[15]。地形由西北向的燕山-太行山山系構(gòu)造向東南逐步過渡為平原,呈現(xiàn)出西北高東南低的地形特點(diǎn)[16]。最高峰為位于河北省的小五臺山,海拔2882 m。除西北壩上高原、圍場山地、張家口、宣化河谷盆地、蔚縣盆地、豐寧山地屬溫帶大陸性季風(fēng)氣候區(qū)外,其余皆為暖溫帶大陸性季風(fēng)氣候區(qū),年平均氣溫0—13℃,年平均降水量為300—800 mm[17]。

京津冀地區(qū)的地帶性植被是以櫟類、楊樺林等為主的暖溫帶落葉闊葉林,間或有溫性松林分布；山地高海拔地區(qū)和河北省的北部地區(qū)分布有山地草甸,亦或是草甸草原。東南部平原區(qū)是主要農(nóng)作區(qū),以小麥、玉米等農(nóng)作物為主,自然植被很少。地帶性土壤主要有棕壤和褐土,此外有少量的鹽堿土和沼澤土。自山地、山麓至濱海大致是：棕壤(棕色森林土)-淋溶褐土-潮褐土-潮土和褐土化潮土-濱海鹽堿土。

2 材料與方法

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

(1)自然保護(hù)區(qū)野生動植物信息數(shù)據(jù)資料

在研究區(qū)域內(nèi)選擇了野生動植物信息數(shù)據(jù)齊全的35個自然保護(hù)區(qū)(表1),其中屬于森林生態(tài)系統(tǒng)類型的有24個,內(nèi)陸濕地類型的有7個,草原草甸類型的有4個。野生動植物信息數(shù)據(jù)包括野生動植物種類及被IUCN物種紅色名錄、生物多樣性紅色名錄、中國物種紅色名錄、國家重點(diǎn)保護(hù)野生動植物名錄的收錄情況。數(shù)據(jù)來源于公開出版的自然保護(hù)區(qū)綜合科學(xué)考察報告[18- 26]、學(xué)位論文[27- 29]、專著[30- 33]等,或者是自然保護(hù)區(qū)的內(nèi)部資料。

(2)植被分布數(shù)據(jù)

來源于1∶1000000 的矢量化中國植被圖[34]。

(3)地形高程數(shù)據(jù)

來源于91衛(wèi)圖助手下載的分辨率為60 m的DEM數(shù)據(jù)。

(4)行政區(qū)區(qū)劃、居民點(diǎn)等其他數(shù)據(jù)

京津冀地區(qū)行政區(qū)邊界、居民點(diǎn)位置信息來源于我國行政區(qū)劃圖[35]。

2.2 研究方法

2.2.1自然保護(hù)區(qū)物種多樣性保護(hù)價值量化

根據(jù)林業(yè)行業(yè)相關(guān)技術(shù)標(biāo)準(zhǔn)[36],對35個自然保護(hù)區(qū)綜合科學(xué)考察報告中的維管束植物、陸地脊椎動物與鳥類名錄進(jìn)行整理,將動植物的瀕危性和保護(hù)等級作為主要參考指標(biāo),從中篩選出IUCN物種紅色名錄、生物多樣性紅色名錄、中國物種紅色名錄中近危種及以上的動植物物種,國家二級及以上的重點(diǎn)保護(hù)動植物物種(除國家重點(diǎn)保護(hù)野生動植物名錄以外的名錄,以下統(tǒng)稱為其他名錄),共計(jì)218種。按保護(hù)等級為這些動植物物種賦分：國家一級重點(diǎn)保護(hù)物種、其他名錄中極危種賦值為8分；國家二級保護(hù)物種、其他名錄中瀕危種賦值為4分；其他名錄中易危種賦值為2分；其他名錄中近危種和無危種賦值為1分[37- 38]。一個物種在不同名錄中保護(hù)級別不同時,按照最高級別賦分的方式。在此基礎(chǔ)上算出每個自然保護(hù)區(qū)的物種多樣性保護(hù)價值得分,記作保護(hù)價值,并以此作為判斷生物多樣性熱點(diǎn)分布的指引依據(jù)。

2.2.2物種多樣性保護(hù)優(yōu)先區(qū)預(yù)測

1)以這24個森林生態(tài)系統(tǒng)類型自然保護(hù)區(qū)的平均面積為參考,利用ArcGIS的創(chuàng)建漁網(wǎng)(fishnet)功能,將京津冀地區(qū)劃分為12 km×12 km大小的格網(wǎng),共計(jì)1651個網(wǎng)格單元(含不完整網(wǎng)格單元)。

2)利用GIS軟件將京津冀地區(qū)的DEM數(shù)據(jù)、植被分布數(shù)據(jù)與1651個網(wǎng)格單元(含不完整網(wǎng)格單元)進(jìn)行疊加,計(jì)算出每個網(wǎng)格單元內(nèi)不同海拔段(0—200 m,201—400 m,401—600 m,…, >2601 m)、不同植被類型(參考《中國植被圖集》的分類)、坡度(平坡(0°—5°)、緩坡(5°—15°)、斜坡(15°—25°)、陡坡(25°—35°)、急坡(35°—45°)和險坡(大于45°))、坡向(平地、陰坡(N、NW)、陽坡(S、SW)、半陰坡(W、NE)、半陽坡(SE、E))等因子的面積比例以及最高海拔、最低海拔和海拔差等,共計(jì)32個因子。

4)利用SPSS軟件對2)中得到的各地形因子、植被分布因子等32個因子進(jìn)行主成分分析(PCA)降維,獲得影響自然保護(hù)區(qū)保護(hù)價值的主成分變量。

5)提取24個森林生態(tài)系統(tǒng)類型自然保護(hù)區(qū)所對應(yīng)的單元網(wǎng)格的主要成分變量,利用SPSS軟件將24個自然保護(hù)區(qū)的保護(hù)價值與對應(yīng)的單元網(wǎng)格的主成分變量進(jìn)行多元回歸分析,確定回歸系數(shù)及回歸方程,判斷其顯著性。將24個自然保護(hù)區(qū)涵蓋范圍外的網(wǎng)格單元的主成分變量代入回歸方程,從而獲得每個網(wǎng)格單元的預(yù)測值,預(yù)測值即為每一個網(wǎng)格單元的保護(hù)價值。

多元線性回歸模型的建立：

假設(shè)某一因變量y受k個自變量x1,x2,…,xk的影響,其觀測值為(y,x1,x2,…,xk),多元線性回歸模型的結(jié)構(gòu)為：

y=β0+β1x1+β2x2+,…,+βkxk+ε

(1)

式中,β0,β1,…,βk為待定參數(shù),ε為常量

如果b0,b1,…,bk分別是β0,β1,β2,…,βk的擬合值,則回歸方程為：

γ′=b0+b1x1+b2x2+,…,+bkxk

(2)

式中,b0為常數(shù),b1,b2,…,bk為偏回歸系數(shù)。

該回歸模型中,y即為觀測值,x1,x2,…,xk即為生成主成分變量,β0,β1,…,βk,即為系數(shù)矩陣中的常數(shù),ε為殘差統(tǒng)計(jì)矩陣中的殘差,γ′即為預(yù)測值。

6)將網(wǎng)格單元預(yù)測值導(dǎo)入ArcGIS中對應(yīng)網(wǎng)格單元,即該網(wǎng)格的保護(hù)價值,形成京津冀地區(qū)的物種多樣性保護(hù)價值分布圖,最后根據(jù)保護(hù)價值的分布情況確定保護(hù)優(yōu)先區(qū)。

3 結(jié)果與分析

3.1 自然保護(hù)區(qū)物種多樣性保護(hù)價值情況

通過對京津冀地區(qū)35個自然保護(hù)區(qū)的物種多樣性保護(hù)價值進(jìn)行量化評分,結(jié)果見表1。

由表1可見,國家級自然保護(hù)區(qū)普遍得分較高,平均得分為204,比參與評分的35個自然保護(hù)區(qū)平均分高出40分。在國家級自然保護(hù)區(qū)中河北衡水湖自然保護(hù)區(qū)的保護(hù)價值得分最高,為306,也是唯一一個得分超過300分的自然保護(hù)區(qū),其次是河北駝梁自然保護(hù)區(qū)、塞罕壩自然保護(hù)區(qū),得分分別為288、266。得分較低的國家級自然保護(hù)區(qū)分別是北京松山自然保護(hù)區(qū)、北京百花山和天津八仙山自然保護(hù)區(qū)、得分分別為152、143和139。在省級自然保護(hù)區(qū)中也存在一些得分較高的自然保護(hù)區(qū),如唐海濕地自然保護(hù)區(qū)和河北南大港自然保護(hù)區(qū),這兩個濕地類型的自然保護(hù)區(qū)的保護(hù)價值得分都在200以上,且超過了參評國家級自然保護(hù)區(qū)保護(hù)價值得分的平均值,因此應(yīng)引起注意。

表1 京津冀地區(qū)自然保護(hù)區(qū)動植物保護(hù)等級、數(shù)量及物種多樣性保護(hù)價值統(tǒng)計(jì)表

從不同類型自然保護(hù)區(qū)的保護(hù)價值的得分情況看,森林類型的自然保護(hù)區(qū)得分最高的為河北駝梁自然保護(hù)區(qū),得分為288；得分最低的為河北興隆六里坪獼猴省級自然保護(hù)區(qū),得分僅為40。濕地類型的自然保護(hù)區(qū)得分最高的為河北衡水湖國家級自然保護(hù)區(qū),得分為306,得分最低的為河北白洋淀濕地省級自然保護(hù)區(qū),得分為45。草原草甸類型的自然保護(hù)區(qū)得分最高的為河北圍場紅松洼國家級自然保護(hù)區(qū),得分為243；得分最低的是河北豐寧滿族灤河源省級自然保護(hù)區(qū),得分為112。由此可以看出不同類型自然保護(hù)區(qū)的保護(hù)價值之間不存在顯著差異,同一類型自然保護(hù)區(qū)的保護(hù)價值差別卻很大,尤其是森林生態(tài)系統(tǒng)類型自然保護(hù)區(qū)和內(nèi)陸濕地類型的自然保護(hù)區(qū)。

3.2 物種多樣性保護(hù)優(yōu)先區(qū)預(yù)測及保護(hù)優(yōu)先等級劃分

3.2.1回歸模型的構(gòu)建

對32個環(huán)境變量因子進(jìn)行主成分分析,選取特征值大于1的根作為主成分?jǐn)?shù)量,共提取了8個主成分,即λ1=11.998、λ2=3.593、λ3=3.045……λ8=1.053,分別對應(yīng)的第1個主成分、第2個主成分、第3個主成分……第8個主成分,其方差貢獻(xiàn)率分別為37.495%、11.229%、9.516%……3.291%,累計(jì)方差貢獻(xiàn)率為80.196%,見表2。

表2 主成分分析的解釋總方差

利用24個自然保護(hù)區(qū)的保護(hù)價值和相對應(yīng)的8個(主成分)環(huán)境變量進(jìn)行回歸分析,確定回歸系數(shù),建立回歸方程,見公式(3)。并通過方差分析,對整個回歸方程做顯著性檢驗(yàn),結(jié)果顯示Sig=0.048,見表3。說明主成分變量與保護(hù)價值存在顯著相關(guān)性,可以用這個模型來預(yù)測。

表3 顯著性分析表

回歸方程如下：

γ′=204.81-23.00x1+23.42x2+39.46x3-16.46x4+29.06x5-14.27x6-14.41x7-46.97x8

(3)

3.2.2保護(hù)價值的計(jì)算

將京津冀地區(qū)所有網(wǎng)格單元的主成分環(huán)境因子代入以上回歸方程,計(jì)算出所有網(wǎng)格的保護(hù)價值,將網(wǎng)格單元保護(hù)價值導(dǎo)入ArcGIS中對應(yīng)網(wǎng)格單元,結(jié)果如圖1所示,圖中顏色越深,說明該網(wǎng)格的保護(hù)價值得分越高,即保護(hù)價值越高。

圖1 網(wǎng)格化的京津冀地區(qū)物種多樣性保護(hù)價值分布圖 Fig.1 Biodiversity conservation value distribution in latticed BTHr

對24個森林類型自然保護(hù)區(qū)的保護(hù)價值得分進(jìn)行“有序樣品聚類”分析(圖2),將24個自然保護(hù)區(qū)的保護(hù)價值分為了三組,其中編號為23,24的自然保護(hù)區(qū)組成一組,該組的保護(hù)價值最高；編號為1,2,3,4的自然保護(hù)區(qū)組成一組,該組的保護(hù)價值最低。剩余的劃為一組,其保護(hù)價值居中,取這一組保護(hù)價值中的低值(133)作為保護(hù)優(yōu)先區(qū)和非保護(hù)優(yōu)先區(qū)的分界線；取這一組保護(hù)價值中的高值(202)作為一級保護(hù)優(yōu)先區(qū)和二級保護(hù)優(yōu)先區(qū)的分界線；這一組又可分為兩小組,取保護(hù)價值較低一組的高值(156)作為二級保護(hù)優(yōu)先區(qū)和三級保護(hù)優(yōu)先區(qū)的分界線。最終確定以133分、156分、202分3個值作為劃分優(yōu)先保護(hù)級別的臨界值。133分以下的區(qū)域不作為保護(hù)區(qū)域,133分—156分的區(qū)域劃分為三級保護(hù)優(yōu)先區(qū),157分—202分的區(qū)域劃分為二級保護(hù)優(yōu)先區(qū),203分以上的區(qū)域劃分為一級保護(hù)優(yōu)先區(qū)。

圖2 自然保護(hù)區(qū)保護(hù)價值得分聚類圖Fig.2 Distribution trend of conservation value of nature reserve

為了增加保護(hù)優(yōu)先區(qū)的精度,將預(yù)測到的物種多樣性保護(hù)優(yōu)先區(qū)等級分布圖與京津冀地區(qū)自然植被圖進(jìn)行疊加分析,去除其中農(nóng)田、水域、鹽場等非自然植被,疊加分析后得到圖3。統(tǒng)計(jì)表明,劃定的三個等級保護(hù)優(yōu)先區(qū)總面積為36791.35 km2,占京津冀總面積的16.94%。其中,一級保護(hù)優(yōu)先區(qū)總面積為4611.57 km2,占總優(yōu)先區(qū)域面積的12.53%；二級保護(hù)優(yōu)先區(qū)總面積為16045.79 km2,占總優(yōu)先區(qū)域面積的43.61%；三級保護(hù)優(yōu)先區(qū)總面積為16133.98 km2,占總優(yōu)先區(qū)域面積的43.85%。

圖3 京津冀地區(qū)物種多樣性保護(hù)優(yōu)先區(qū)分布圖Fig.3 Priority area for biodiversity conservation in BTHr

保護(hù)優(yōu)先區(qū)內(nèi)分布有森林植被8204.47 km2,灌叢植被14457.72 km2,灌草叢植被3427.88 km2,草甸、草原、草叢等草本植被10701.25 km2,分別占總面積的29.09%,9.32%,39.30%,22.30%。森林植被主要有白樺林、蒙古櫟林、華北落葉松林、油松林、山楊林等,灌叢植被主要有荊條酸棗灌叢、繡線菊灌叢、虎榛子灌叢、黃櫨灌叢等,草本植被類型較復(fù)雜,可以分為草原、草甸和一般草地,草原主要有羊草、叢生禾草草原、西北針茅草原、貝加爾針茅、雜類草草原等,草甸主要有苔草、雜類草草甸、野古草、大油芒、雜類草草甸、芨芨草草甸等。

物種多樣性保護(hù)優(yōu)先區(qū)共涉及京津冀地區(qū)56個區(qū)縣,主要有河北省北部的圍場滿族蒙古族自治縣、豐寧滿族自治縣、隆化縣、赤城縣、沽源縣、康保縣、灤平縣、平泉縣、承德縣、涿鹿縣、張北縣；北京市的延慶區(qū)、密云區(qū)、房山區(qū)、昌平區(qū)等。所有區(qū)、縣(市)中圍場滿族蒙古族自治縣內(nèi)保護(hù)優(yōu)先區(qū)面積最大,約為5153.48 hm2,其次是豐寧滿族自治縣,約為4044.21 hm2。

4 討論

(1)以國家重點(diǎn)保護(hù)的、珍稀瀕危的的野生動植物物種作為評價生物多樣性熱點(diǎn)區(qū)域的指標(biāo)是當(dāng)前研究的重要方法。如欒曉峰等人在進(jìn)行東北地區(qū)生物多樣性熱點(diǎn)區(qū)域評價和保護(hù)空缺研究時,將國家一級和二級重點(diǎn)保護(hù)物種、東北地區(qū)特有物種、東北地區(qū)罕見或數(shù)量極為稀少的珍稀種類以及國際上較為關(guān)注的物種,作為生物多樣性熱點(diǎn)地區(qū)評價的指示物種,并最終確定了64種瀕危和特有動植物][39]。盧怡萌在進(jìn)行中國野生蘭科植物優(yōu)先保護(hù)研究時,僅用了《中國生物多樣性紅色名錄-高等植物卷》中的受脅迫的蘭科植物作為中國蘭科植物優(yōu)先保護(hù)物種[40]。馬琳在進(jìn)行長白山森林生物多樣性體系研究時,選取了國家級、省級重點(diǎn)保護(hù)物種,國際上較關(guān)注的物種(IUCN紅色名錄物種、CITES附錄物種)、對生態(tài)系統(tǒng)或生態(tài)過程具有重要意義的物種等,共確定141種優(yōu)先保護(hù)對象作為優(yōu)先保護(hù)物種選擇的標(biāo)準(zhǔn)[41]。可見將重點(diǎn)保護(hù)動植物名錄作為評價指標(biāo)物種,并以此作為研判生物多樣性熱點(diǎn)區(qū)域、物種多樣性保護(hù)優(yōu)先區(qū)已成為了研究生物多樣性保護(hù)與分布格局的重要方法之一,只是在對于重點(diǎn)保護(hù)物種選擇的名錄范圍上可能不同。本研究也采用當(dāng)前這種評價指標(biāo)選擇方法,以京津冀地區(qū)自然保護(hù)區(qū)內(nèi)的國家重點(diǎn)保護(hù)野生動植物為評價指標(biāo),共選擇物種218種,可以看出在選擇物種的數(shù)量上相對于以上幾個研究是比較多的,也應(yīng)更能代表京津冀地區(qū)的生物多樣性。

(2)區(qū)域尺度上綜合考慮地形地貌對物種多樣性分布格局的影響更有意義。地貌單元和地形的局部變化嚴(yán)重制約著光照、溫度、水分甚至土壤資源的再分配[42]。在大多數(shù)有關(guān)地形地貌與生物多樣性分布格局的研究中,都是將地形地貌分解成幾個單因素分別進(jìn)行研究,其中研究較多的是海拔梯度與生物多樣性分布格局的關(guān)系,其次是坡向、坡度、凹凸度等地形地貌因子與生物多樣性分布格局的關(guān)系,只有很少的研究是將地表褶皺度,或者是地表粗糙度等綜合性的地形地貌因子作為影響因素,開展其與生物多樣性分布格局研究的[43]。然而地形地貌對生物多樣性分布的影響,尤其是植物分布的影響是綜合的,因此研究綜合的地形地貌因子對生物多樣性分布的影響會更有效。

另外,與地形有關(guān)的生物多樣性的研究多是在群落或樣地這樣小尺度上開展的,以群落樣方數(shù)據(jù)為依據(jù),開展不同坡度、坡向,海拔等不同地形條件下物種多樣性的變化,如張謐、王家鳴等人開展的相關(guān)研究[40,44],都是基于大樣地開展的。物種多樣性的空間分布格局與取樣區(qū)域、取樣面積、生境異質(zhì)性、干擾以及生物區(qū)系等有關(guān)[45- 46]。區(qū)域尺度上單位面積物種多樣性的分布格局與群落尺度明顯不同,“尺度效應(yīng) ”在塑造物種多樣性地理分布格局中具有重要重要作用[47]。因此,有研究認(rèn)為僅采取小尺度樣方調(diào)查,難以反映區(qū)域的物種多樣性[48],Shannon多樣性均值及多樣性指數(shù)也不能反映區(qū)域物種多樣性全貌[49- 50]。因此以京津冀地區(qū)自然保護(hù)區(qū)的平均面積這樣一個中等尺度開展地形異質(zhì)性與生物多樣性的分布格局研究,更有利于生物多樣性保護(hù)工作的開展。

5 結(jié)論

(1)以國家重點(diǎn)保護(hù)野生動植物物種、珍稀瀕危動植物種為評價指標(biāo),評價了京津冀地區(qū)自然保護(hù)區(qū)的物種多樣性保護(hù)價值。國家級自然保護(hù)區(qū)普遍得分較高,平均得分為204分,比參與評分的35個自然保護(hù)區(qū)平均分高出40分；在省級自然保護(hù)區(qū)中也存在一些得分較高的自然保護(hù)區(qū),如,唐海濕地自然保護(hù)區(qū)和河北南大港自然保護(hù)區(qū),且得分超過了參評國家級自然保護(hù)區(qū)保護(hù)價值得分的平均值,應(yīng)引起注意。不同類型自然保護(hù)區(qū)之間的保護(hù)價值得分不存在顯著差異,同一類型自然保護(hù)區(qū)的保護(hù)價值差別卻很大,尤其是森林生態(tài)系統(tǒng)類型自然保護(hù)區(qū)和內(nèi)陸濕地類型的自然保護(hù)區(qū)。

(2)以京津冀地區(qū)的24個森林類型自然保護(hù)區(qū)為參照,結(jié)合地形地貌、自然植被等因素,利用多元回歸分析方法預(yù)測了京津冀陸域范圍內(nèi)的生物多樣性熱點(diǎn)分布格局,劃定了三級保護(hù)優(yōu)先區(qū),總面積為36791.35 km2,占京津冀總面積的16.94%。其中,一級保護(hù)優(yōu)先區(qū)總面積為4611.57 km2,占總優(yōu)先區(qū)域面積的12.53%；二級保護(hù)優(yōu)先區(qū)總面積為16045.79 km2,占總優(yōu)先區(qū)域面積的43.61%；三級保護(hù)優(yōu)先區(qū)總面積為16133.98 km2,占總優(yōu)先區(qū)域面積的43.85%。