王金鳳,徐基良,李建強,*,周春發(fā),鄧文洪

1 北京林業(yè)大學(xué)生態(tài)與自然保護學(xué)院,北京 100083 2 北京師范大學(xué)生命科學(xué)學(xué)院,北京 100875

2018年修訂的《中華人民共和國野生動物保護法》明確規(guī)定要加強野生動物棲息地保護,并提出建立野生動物重要棲息地名錄制度[1]。已有實踐表明,建立自然保護地是就地保護野生動物及其棲息地的有效手段[2]。截至2020年底,我國已建成各級各類自然保護地約1.18萬余處,覆蓋了我國陸域面積的18.0%和管轄海域的4.1%[3]。這些自然保護地在保護我國重要生態(tài)系統(tǒng)、珍稀瀕危動植物及其棲息環(huán)境上發(fā)揮了重要作用,但也存在著空間布局不合理、與物種重要分布區(qū)域在空間上吻合不足等問題,進而影響其保護效果的發(fā)揮[4],逐漸引起了社會關(guān)注[5—7]。

衡量自然保護地保護成效的重要指標(biāo)之一是評估對野生動物重要棲息地的保護效果[8—9],目前已經(jīng)有許多學(xué)者開展了這方面的研究[10—13]。然而,相關(guān)研究多數(shù)集中于對單一物種的分析[14—15]。這雖然可以反映自然保護地對特定物種或特定類群及其棲息地的保護效果[16],但是可能并不足以反映其對不同類群的整體保護效果[17—18]。因此,有研究建議以某一區(qū)域的多個物種為對象進行評估,可以更加客觀地反映區(qū)域內(nèi)自然保護地的保護成效[9]。

北京市(115.7°—117.4°E,39.4°—41.6°N)地處燕山山脈和太行山脈的交叉地帶,野生動植物資源豐富[19],擁有陸生脊椎野生動物約581種,包括褐馬雞(Crossoptionmantchuricum)、黑鸛(Ciconianigra)、麋鹿(Elaphurusdavidianus)和狼(Canislupus)等國家重點保護野生動物81種[20]。北京市雖然是一個高度城市化的地區(qū),但隨著生態(tài)環(huán)境的不斷改善,北京市野生動物種群數(shù)量和棲息地分布范圍近年來呈現(xiàn)明顯的擴張趨勢[21—22]。如今,北京市已經(jīng)建立了大量的自然保護地以加強對野生動物及其棲息地的保護。一些研究評估了這些自然保護地對部分野生動物及其棲息地的保護成效[22—23],然而相關(guān)研究多關(guān)注于單一物種或單一類群的保護成效。本文選取了北京地區(qū)六種較為典型的野生動物作為研究對象,利用MaxEnt模型分析了其適宜棲息地的分布,并結(jié)合北京市自然保護地的類型及其空間布局,深入分析其適宜棲息地在北京市的分布狀況和保護空缺,評估北京市自然保護地保護有效性,旨在為加強北京市野生動物資源保護、劃定野生動物重要棲息地[24]、推動自然保護地整合優(yōu)化提供重要參考。

1 材料與方法

1.1 研究區(qū)域

截至2019年底,北京市共有5類79處自然保護地,其中自然保護區(qū)21處(國家級2處、市級12處、區(qū)級7處)、風(fēng)景名勝區(qū)11處、森林公園31處、地質(zhì)公園6處、濕地公園10處[19](圖1)。自然保護地主要分布在西北山區(qū)、淺山地帶以及平原的河流水系地區(qū),總面積4648.05km2(含交叉重疊面積),空間上實際面積共3470.61km2,占市域行政面積的21.15%。

圖1 北京地區(qū)各類自然保護地及其分布Fig.1 The distribution of various types of the natural protected areas in Beijing

1.2 研究對象

依據(jù)北京市野生動物種類和分布現(xiàn)狀,通過咨詢專家,結(jié)合野外分布資料的可獲得性,選取在北京地區(qū)具有代表性的6個國家重點保護野生動物,包括國家I級重點保護野生動物黑鸛、褐馬雞、大鴇(Otistarda)、金雕(Aquilachrysaetos),國家II級重點保護野生動物鴛鴦(Aixgalericulata)和斑羚 (Naemorhedusgriseus)[25]。其中黑鸛作為近些年來在北京繁殖、越冬的涉禽,可體現(xiàn)山谷濕地生態(tài)狀況[26]；褐馬雞作為中國特有種[27]、典型的地棲性山地森林鳥類[28],近些年來在北京的種群及其分布范圍穩(wěn)步上升[22]；大鴇是典型的草地生態(tài)指示物種[29]；金雕作為北京地區(qū)的大型猛禽,可以指示區(qū)域生態(tài)系統(tǒng)健康[30]；鴛鴦則是重要的評價濕地生態(tài)系統(tǒng)的代表物種之一[31]；斑羚為典型的林棲獸類[15]。這6個物種的棲息地包括森林、濕地、草地、濕地水域、高山峭壁等類型,涉及到北京市濕地公園、森林公園、自然保護區(qū)、風(fēng)景名勝區(qū)、地質(zhì)公園等5類自然保護地,具有較強的典型性和代表性。

1.3 研究方法

1.3.1物種數(shù)據(jù)來源及處理

物種分布位點數(shù)據(jù)是以北京市野生動物資源調(diào)查中獲取的坐標(biāo)信息為基礎(chǔ),并將在中國觀鳥記錄中心(http://www.birdreport.cn/)、全球生物多樣性信息網(wǎng)絡(luò)數(shù)據(jù)庫(https://www.gbif.org/)和有明確坐標(biāo)信息的相關(guān)文獻中所獲取的準(zhǔn)確的物種分布點坐標(biāo)信息作為補充。由于研究物種活動范圍較大,因此采用緩沖區(qū)分析對收集的物種分布點進行篩選和校對,排除人為因素導(dǎo)致的分布點重合以及不符合實際情況的分布點[32]。根據(jù)6個目標(biāo)物種的平均活動距離[31,33—37],在ArcGIS 10.4中以物種分布點為圓心,建立以500m為半徑的緩沖區(qū)；當(dāng)兩個分布點之間的距離小于1km時,只保留其中一點。最終確定進入分析的各目標(biāo)物種分布點,其中黑鸛105個、褐馬雞24個、大鴇21個、鴛鴦97個、金雕33個、斑羚11個。

1.3.2環(huán)境變量的收集、處理

根據(jù)現(xiàn)有文獻資料中對各目標(biāo)物種的棲息地需求[15,37],選擇海拔、坡度、坡向、土地利用類型、歸一化植被指數(shù)、到居民點(鄉(xiāng)鎮(zhèn)及鄉(xiāng)鎮(zhèn)級別以上的)的距離、到水源(河流湖泊、水庫)的距離、到不同等級的道路(包括鐵路；高速公路、國道、省道；市縣級道路；人行道四個等級)的距離等共11個環(huán)境因子作為初步變量進行分析。

地形數(shù)據(jù)來源于中國科學(xué)院計算機網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云平臺(http://www.gscloud.cn/)下載的數(shù)字高程數(shù)據(jù)(DEM),空間分辨率為30m?；贏rcGIS 10.4軟件提供的空間分析功能,運用表面分析工具從DEM數(shù)據(jù)中提取出坡度和坡向兩種地形因子的柵格數(shù)據(jù)。土地利用分類數(shù)據(jù)來源于國家基礎(chǔ)地理信息中心(http://www.globeland30.org/),空間分辨率為30m。在ArcGIS 10.4中,使用重分類工具從土地利用類型中得出研究區(qū)范圍內(nèi)的關(guān)鍵地類,然后通過焦點統(tǒng)計工具,以平均活動距離作為半徑得出各關(guān)鍵地類占比的柵格數(shù)據(jù),以各關(guān)鍵地類在活動區(qū)域內(nèi)的占比來限制物種的出現(xiàn),得到影響物種分布的環(huán)境變量。歸一化植被指數(shù)(NDVI)數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心(https://www.resdc.cn/),空間分辨率為30m。在ArcGIS 10.4中,使用腌膜提取工具得到研究區(qū)域的植被覆蓋數(shù)據(jù)。通過查閱文獻,針對研究對象選擇到居民點的距離、到不同等級道路(包括鐵路；高速公路、國道、省道；市縣級道路；人行道4個等級)的距離來反應(yīng)人為干程度。其中,鄉(xiāng)鎮(zhèn)及鄉(xiāng)鎮(zhèn)級別以上的居民點、河流湖泊和水庫、不同等級的道路數(shù)據(jù)來源于全國地理信息資源目錄服務(wù)系統(tǒng)(http://www.webmap.cn)提供的1∶2500000的全國基礎(chǔ)地理數(shù)據(jù)。運用ArcGIS 10.4中的歐氏距離工具來獲取到居民點的距離、到水源的距離和到不同等級道路的距離的柵格數(shù)據(jù)。北京市行政邊界數(shù)據(jù)來自于國家基礎(chǔ)地理信息系統(tǒng)(http://nfgis.nsdi.gov.cn/)1∶4000000的中國行政區(qū)劃矢量地圖,自然保護地邊界來源于北京市園林綠化局。

1.3.3模型構(gòu)建、精度檢驗及相關(guān)性檢驗

使用MaxEnt v3.4.1軟件[38]評估各物種的適宜棲息地。首先檢驗環(huán)境變量兩兩之間的Pearson′s相關(guān)系數(shù),對于|r|≥0.8的一對變量,保留對初步模型結(jié)果貢獻大的一個[39]；然后針對不同物種,根據(jù)刀切法(Jackknife test)檢驗系數(shù)的結(jié)果,剔除貢獻率低于0.4的環(huán)境變量,最終確定用于預(yù)測研究對象棲息地的環(huán)境因子[40]。

將6個物種的分布點數(shù)據(jù)和環(huán)境因子數(shù)據(jù)分別導(dǎo)入MaxEnt模型中,針對不同研究對象隨機選取該物種75%的樣本數(shù)據(jù)作為訓(xùn)練數(shù)據(jù)構(gòu)建最大熵模型,剩余25%的樣本用于模型驗證[32]。使用AUC評價指標(biāo)對模型的效果和各個環(huán)境特征變量的重要性進行檢測,來驗證最大熵模型的準(zhǔn)確度[41]。AUC的取值在0—1之間,數(shù)值越接近于1表示預(yù)測精確度越高[32]。對于每個物種,使用模型中的自舉重復(fù)功能運算10次,取AUC值較大者作為模型模擬的最終結(jié)果[39]。

1.4 適宜棲息地分析

MaxEnt模型輸出的預(yù)測結(jié)果是物種在研究區(qū)域的存在概率,適宜棲息地的值越高,代表物種存在的條件越有利[42]。模型輸出的預(yù)測結(jié)果通過閾值來劃分適宜生境和不適宜生境,閾值的選擇對于確定最終的物種適宜棲息地分布范圍十分重要[32]。在模型運行過程中,得出的10%存在點訓(xùn)練集邏輯斯蒂閾值(10%TPLT)可以確保在全部已知樣點中至少包含90%的樣點是適宜的[43]。所以選擇10%存在點訓(xùn)練集邏輯斯蒂閾值將預(yù)測的概率圖轉(zhuǎn)化為Presence/Absence二值圖[32]。利用ArcGIS 10.4對其進行賦值,Presence賦值為1,即適宜棲息地；Absence賦值為0,即不適宜。最終得到6個代表性物種的二值圖[39]。通過疊加6個物種的二值圖得到單個像元值最高為6,取像元值≥3者為分布熱點區(qū),得到物種適宜棲息地的分布熱點[40]。

1.5 保護成效與空缺分析

將各目標(biāo)物種的適宜棲息地分布范圍與自然保護地進行疊加,將每個物種在自然保護地內(nèi)的分布面積占其分布總面積的比例視為該物種的受保護率[9]。在ArcGIS 10.4中,分別計算物種受保護率以及北京市自然保護地對各目標(biāo)物種適宜棲息地的覆蓋比例,用于從物種和自然保護地類型兩個方面反映北京市自然保護地的保護成效。將6個物種適宜棲息地的分布熱點圖與北京市自然保護地分布圖進行空間疊加,識別需要保護但目前尚不屬于自然保護地保護范圍的區(qū)域,即為保護空缺[40]。

2 研究結(jié)果

2.1 預(yù)測模型精度檢驗

6個野生動物模型模擬結(jié)果的相應(yīng)ROC值均大于0.90(圖2),表明MaxEnt模型預(yù)測能力優(yōu)秀,影響6個物種分布的環(huán)境因子和模型預(yù)測的對應(yīng)物種重要棲息地分布之間存在比較好的相關(guān)性。

圖2 6個物種的MaxEnt模型預(yù)測結(jié)果的受試者工作特征曲線(ROC曲線)Fig.2 Receiver Operating Characteristic Curve (ROC) of MaxEnt model prediction results of the 6 species

2.2 適宜棲息地分布

黑鸛在北京地區(qū)的適宜棲息地面積為3559.67km2,占研究區(qū)總面積的21.22%。主要分布在密云區(qū)密云水庫及周邊地區(qū)、房山區(qū)的中部地區(qū)以及西南部、門頭溝區(qū)的東部和中部、海淀區(qū)的中北部、昌平區(qū)北部、懷柔區(qū)的中部和南部等地。大鴇在北京地區(qū)的適宜棲息地面積為95.70km2,占北京市行政總面積的0.57%。主要分布在密云區(qū)的密云水庫周邊水系附近、平谷區(qū)-順義區(qū)-通州區(qū)北部的河流水系沿線、昌平區(qū)東北部、延慶區(qū)的官廳水庫及周邊地區(qū)等地。鴛鴦在北京地區(qū)的適宜棲息地面積為2993.90km2,占北京市行政總面積的17.85%。主要分布在中心城區(qū)的城市公園、密云區(qū)密云水庫周邊區(qū)域、房山區(qū)西南部、門頭溝區(qū)東部和中部地區(qū)、延慶區(qū)官廳水庫周邊等地,多靠近河流水系分布(圖3)。

褐馬雞在北京地區(qū)的適宜棲息地面積為19.70km2,共占北京市行政總面積的0.12%。主要集中在門頭溝區(qū)西南部的山林地區(qū)。金雕在北京地區(qū)的適宜棲息地面積為642.75km2,占北京市行政總面積的3.83%。主要集中在門頭溝區(qū)和房山區(qū)西北部的山區(qū),另外在延慶區(qū)和密云區(qū)也均有分布。斑羚在北京地區(qū)的適宜棲息地面積為711.02km2,占北京市行政總面積的4.24%。主要集中在北京西部的山地區(qū)域,包括懷柔區(qū)、延慶區(qū)、昌平區(qū)、門頭溝區(qū)以及房山區(qū)等地,另外在密云區(qū)也有少量分布(圖3)。

圖3 6個物種在北京市的適宜棲息地分布Fig.3 Distribution of suitable habitats of the 6 species in Beijing

2.3 適宜棲息地的熱點區(qū)域

6個物種適宜棲息地分布的熱點區(qū)域面積為280.01km2(圖4),占北京市行政面積的1.71%。熱點區(qū)域主要包括密云區(qū)密云水庫周邊地區(qū)、延慶區(qū)官廳水庫及周邊區(qū)域,門頭溝區(qū)靈山自然風(fēng)景區(qū)-小龍門森林公園、永定河-清水河沿線和房山區(qū)西北部百花山國家級自然保護區(qū)、拒馬河沿線、大石河沿線等地。

圖4 6個物種適宜棲息地在北京地區(qū)的熱點分布(物種種數(shù)≥3者即為熱點分布區(qū)域)Fig.4 The suitable habitat hotspots of the 6 species in Beijing (the area with ≥ 3 species is classified as the hotspot area)

2.4 自然保護地對物種適宜棲息地的保護成效

6個研究物種在現(xiàn)有的自然保護地體系中均得到了不同程度的保護,平均受保護率為41.00%。其中褐馬雞的受保護率最高(92.82%),其次為斑羚(52.43%)、金雕(43.79%)、大鴇(23.04%)、黑鸛(19.94%),鴛鴦的受保護率最低(13.66%)。

北京市自然保護地覆蓋了141.54km2的熱點區(qū)域,有49.45%的熱點區(qū)域未被覆蓋,未被覆蓋的熱點區(qū)域即為保護空缺。各類自然保護地對典型物種熱點地區(qū)的覆蓋率存在差異,其中風(fēng)景名勝區(qū)相對覆蓋比例較高(14.74%)；其余類型的覆蓋率由高到低依次為自然保護區(qū)(14.17%)、地質(zhì)公園(10.49%)、森林公園(4.10%)和濕地公園(0.83%)(表1)。

表1 北京市不同類型自然保護地對6個典型物種適宜棲息地?zé)狳c分布的覆蓋情況

2.5 保護空缺

物種適宜棲息地?zé)狳c地區(qū)的保護空缺主要分布在密云水庫各支流的沿線及周邊區(qū)域；延慶區(qū)南部、懷柔區(qū)南部與昌平區(qū)北部相接壤的地區(qū)；門頭溝區(qū)中部永定河-清水河沿線及周邊地區(qū)；房山區(qū)的北部、西部及中部地區(qū)；順義區(qū)的潮白河沿線等地。其中,門頭溝區(qū)和房山區(qū)的保護空缺面積最大(圖5)。

圖5 北京市6個典型物種的保護空缺分布Fig.5 Spatial distribution of the GAPs of the 6 species in Beijing

3 討論

保護以典型物種為代表的生物多樣性是設(shè)立自然保護地的首要目標(biāo)之一,因此對物種保護成效的評估可以反映自然保護地對物種及其棲息地的保護效果[16]。以國家重點保護的、珍稀瀕危的野生動植物物種作為評價生物多樣性熱點區(qū)域、物種多樣性保護優(yōu)先區(qū)的指標(biāo)則是當(dāng)前物種多樣性保護和分布格局的重要方法[13,44]。熱點區(qū)域的識別則可以幫助我們以最小的成本、最大限度地保護區(qū)域生物多樣性[45]。

褐馬雞、金雕等6個物種適宜棲息地的熱點區(qū)域主要分布于密云區(qū)、懷柔區(qū)與昌平區(qū)的交界處、延慶區(qū)西部、門頭溝區(qū)和房山區(qū)等地。其中,黑鸛、大鴇和鴛鴦主要分布在淺山地帶以及平原的濕地區(qū)域,與人類活動空間重疊范圍較大[26,31,34]。褐馬雞作為一種典型的森林雉類[35],在北京主要分布于西北山地森林區(qū)域,受到的人為干擾較小[22]；金雕的適宜棲息地主要集中在海拔較高的山林區(qū)域,人為干擾也較弱[30]；斑羚對人為干擾極為敏感,其適宜棲息地主要分布在北京市北部和西南部中高海拔、植被豐富、靠近水源的區(qū)域[15,46]。

然而,北京市自然保護地對這6個物種適宜棲息地?zé)狳c區(qū)域的覆蓋率整體較低,其中對褐馬雞、金雕和斑羚的保護效率較高,而對黑鸛、大鴇和鴛鴦的保護效率則相對較低。這與北京市自然保護地空間布局和發(fā)展過程密切相關(guān)。京郊山區(qū)是北京市最早開展自然保護地建設(shè)工作的區(qū)域,擁有著北京市70%以上的自然保護地[47]。與此相應(yīng)的,北京市平原及城市中心地區(qū)則是城市發(fā)展的重點區(qū)域[48],相關(guān)政策難以直接惠及該區(qū)域的生物多樣性保護[23]。并且,過去北京市的自然保護地建設(shè)管理具有我國自然保護地發(fā)展的共同特點,即分部門管理[47]。各部門在自然資源相對豐富、有生態(tài)旅游潛力的優(yōu)勢地區(qū),集中建立了很多森林公園、濕地公園、風(fēng)景名勝區(qū)和地質(zhì)公園,這也促使北京市自然保護地呈現(xiàn)聚集分布、保護范圍交叉重疊嚴(yán)重的現(xiàn)象[49]。

在當(dāng)前建立以國家公園為主體的自然保護地體系的背景下,整合優(yōu)化現(xiàn)有自然保護地體系、提升整體保護成效是迫切需要解決的問題,也是重構(gòu)北京市自然保護地體系的重要契機[47]。根據(jù)本研究結(jié)果,為了提升北京市自然保護地的保護成效,應(yīng)因地制宜,針對不同的保護對象采取不同的策略來完善自然保護地體系：對于門頭溝區(qū)和房山區(qū)山林地帶、人為干擾較少的野生動物適宜棲息地保護空缺,可以根據(jù)實際受威脅程度、利用強度等情況,增設(shè)新的自然保護地；對于懷柔區(qū)與昌平區(qū)交界處的未能劃入自然保護地的野生動物適宜棲息地保護空缺區(qū)域,可通過擴大自然保護地范圍將其劃入就近的自然保護地內(nèi)；對于平原區(qū)面積小且分布零散的野生動物適宜棲息地保護空缺,可以通過建設(shè)生態(tài)走廊、設(shè)定保護小區(qū)或納入生態(tài)紅線等形式予以保護；對于分布在城市建成區(qū)內(nèi)的野生動物,也可以在濕地、城市公園等區(qū)域生態(tài)關(guān)鍵節(jié)點處,根據(jù)保護對象的棲息地需求營造微生境,為野生保護動物提供棲息場所。

與利用單一物種或單一類群評估區(qū)域自然保護地保護成效的研究相比,本研究表明從多個類群物種的角度評估自然保護地保護成效更具有全面性和科學(xué)性,同時也可以為野生動物重要棲息地的識別與劃定提供可行的方案。然而,本研究也存在一些不足。比如受到數(shù)據(jù)和文獻資料的限制,本研究僅選取了在北京地區(qū)比較典型的國家級重點保護動物作為北京市自然保護地保護成效評估的研究對象,也沒有考慮到北京市自然保護地對植物棲息地的保護成效。因此,今后的工作中,應(yīng)盡量考慮增加更多的物種,從物種的種類和數(shù)量變化、種群的空間分布格局與范圍、種群生存力分析等方面深入評估北京市自然保護地的保護成效,進而更好地提升北京市自然保護地對區(qū)域生物多樣性的保護成效。

4 結(jié)論

(1)研究所選取的6個目標(biāo)物種中,黑鸛、大鴇和鴛鴦主要分布于北京市淺山地帶和平原上的濕地區(qū)域；褐馬雞、金雕和斑羚的棲息地則集中分布于市域西南部和北部中高海拔的山林區(qū)域。由于北京市自然保護地分布不夠均衡,所以總體上黑鸛、大鴇和鴛鴦受自然保護地保護率低,而褐馬雞、金雕和斑羚受自然保護地保護率高。

(2)6個物種適宜棲息地分布的熱點區(qū)域主要分布于密云區(qū)、懷柔區(qū)與昌平區(qū)的交界處、延慶區(qū)西部、門頭溝區(qū)和房山區(qū)等地。北京市現(xiàn)有的自然保護地對物種適宜棲息地分布熱點區(qū)域的覆蓋率相對較低,這與北京市自然保護地空間布局和發(fā)展過程密切相關(guān)。

(3)針對于現(xiàn)存的保護空缺,建議政府相關(guān)部門因地制宜,針對不同的保護對象采取不同的保護策略來完善現(xiàn)有自然保護地體系。如擴大或合并現(xiàn)有自然保護地、增設(shè)新的自然保護地、納入生態(tài)保護紅線或營造微生境等方式來提升北京市自然保護地的保護成效。