李海韻 王潔 燕翃翔 徐瑾

近年來(lái)，隨著人類活動(dòng)強(qiáng)度的增大和城市化的推進(jìn)，生態(tài)環(huán)境面臨的壓力與日俱增。自然保護(hù)地空間規(guī)劃中存在的一系列問(wèn)題加劇了生態(tài)保護(hù)與資源合理利用的矛盾，使得保護(hù)地外部邊界與城鎮(zhèn)開發(fā)界線沖突、保護(hù)地內(nèi)部區(qū)劃不科學(xué)等空間問(wèn)題逐漸凸顯[1]。自然保護(hù)地是由各級(jí)政府依法劃定或確認(rèn)，對(duì)重要的自然生態(tài)系統(tǒng)、自然遺跡、自然景觀及其所承載的自然資源、生態(tài)功能和文化價(jià)值實(shí)施長(zhǎng)期保護(hù)的陸域或海域[2]。雖然從2017年開始，中國(guó)國(guó)家公園體制的建立帶動(dòng)了首批重點(diǎn)自然保護(hù)地的規(guī)范化建設(shè)，但全國(guó)自然保護(hù)地?cái)?shù)量與面積數(shù)據(jù)并不能完全真實(shí)地反映有效保護(hù)現(xiàn)狀[3]。東莞市作為中國(guó)改革開放的前沿陣地和南方林區(qū)的代表區(qū)域，自20世紀(jì)開始的持續(xù)高強(qiáng)度城鎮(zhèn)建設(shè)直接導(dǎo)致了激烈的人地沖突，自然保護(hù)地面臨的空間問(wèn)題尤為突出[4]。

中國(guó)“十四五”規(guī)劃綱要明確指出應(yīng)科學(xué)劃定自然保護(hù)地保護(hù)范圍。自然保護(hù)地空間規(guī)劃的核心環(huán)節(jié)是識(shí)別出具有高保護(hù)價(jià)值的地區(qū)[5]，以該目標(biāo)為導(dǎo)向，如今的保護(hù)規(guī)劃策略傾向于兼顧生物多樣性、保護(hù)成本以及其他社會(huì)經(jīng)濟(jì)因素[6]，高效益低成本區(qū)域的識(shí)別逐漸受到相關(guān)研究者的重視[7]，如Yang等[8]利用荒野制圖研究，按照人類干擾程度將自然景觀分級(jí)，低干擾區(qū)域被劃定為高質(zhì)量荒野區(qū)，并可作為高效益低成本自然保護(hù)地。建立自然保護(hù)地的主要目的是使該區(qū)域內(nèi)的自然資源免受人為干擾[9]，而如何在以東莞市為代表的高密度建設(shè)城鎮(zhèn)開展自然保護(hù)地空間范圍優(yōu)化，是值得深入探討的問(wèn)題。相對(duì)于以實(shí)現(xiàn)全國(guó)荒野保護(hù)目標(biāo)為導(dǎo)向的大尺度荒野識(shí)別而言，東莞市生態(tài)空間的小尺度和建設(shè)用地的高密度決定了因地制宜展開自然保護(hù)地空間優(yōu)化的重要性，即建立指標(biāo)類型與權(quán)重賦值科學(xué)準(zhǔn)確的人為干擾評(píng)價(jià)指標(biāo)體系，進(jìn)行自然保護(hù)地的保護(hù)現(xiàn)狀評(píng)估和保護(hù)空缺分析，科學(xué)衡量人類活動(dòng)帶來(lái)的影響，明確自然保護(hù)地的空間優(yōu)化方向。

人為干擾是在人為行動(dòng)影響下對(duì)自然造成的干擾，即人類由于生產(chǎn)、生活和其他社會(huì)活動(dòng)，在有目的的行為指導(dǎo)下對(duì)自然環(huán)境和生態(tài)系統(tǒng)施加的各種影響[10]。與自然干擾相比，人為干擾的方式更為復(fù)雜[11-12]，因此，合理辨識(shí)人為干擾的作用方式和強(qiáng)度大小，分析人為干擾加劇所導(dǎo)致的生態(tài)系統(tǒng)異化過(guò)程，對(duì)于維持生態(tài)系統(tǒng)結(jié)構(gòu)與功能的持續(xù)性和完整性具有重要意義[13]。人為干擾度（hemeroby）的概念于1955年由植物學(xué)家Jalas首次引入，用于量化人類對(duì)植被的影響[14]。20世紀(jì)70年代，這個(gè)概念被應(yīng)用于整個(gè)生態(tài)系統(tǒng)的人類影響評(píng)估，人為干擾度成為衡量所有人類干擾對(duì)生態(tài)系統(tǒng)影響的綜合指標(biāo)。

隨著衛(wèi)星圖像數(shù)據(jù)集的開放使用[15]和景觀指數(shù)軟件包的大量開發(fā)[16]，景觀生態(tài)學(xué)視角下的人為干擾研究得到了廣泛開展，通常以地理空間技術(shù)定量并圖示人為干擾的影響范圍及程度。當(dāng)前，基于土地覆被類型的生態(tài)環(huán)境人為干擾度分析，被認(rèn)為可以綜合評(píng)估各種人為影響的累積性結(jié)果，并準(zhǔn)確顯示出人為影響程度的空間分布和梯度變化特征[17]。其中最為常用的人為干擾量化方法來(lái)源于Steinhardt等[18]提出的人為干擾度計(jì)算公式，其原理是土地覆被類型的空間分異可直觀反映人為干擾度的空間分異。目前已有研究將其應(yīng)用于遼寧省營(yíng)口市南部海岸[19]、山東省煙臺(tái)市廟島群島北五島[20]、北京城區(qū)[21]、湖南省長(zhǎng)株潭主城區(qū)[22]等區(qū)域尺度生態(tài)環(huán)境的人為干擾度評(píng)價(jià)研究，并表明人為干擾度考慮了不同人類活動(dòng)的影響程度差異，同時(shí)具有計(jì)算簡(jiǎn)便、數(shù)據(jù)可獲取性強(qiáng)以及普適性等優(yōu)點(diǎn)。然而現(xiàn)有人為干擾評(píng)價(jià)方法主要是從對(duì)海岸、灘涂等海拔平緩地區(qū)的研究演變而來(lái)，且諸多研究表明地形因子的分異對(duì)人為活動(dòng)具有約束性影響[23-25]，如低海拔地區(qū)土地利用程度相對(duì)較高，人為影響因素較強(qiáng)；坡度越陡、地形起伏度越大的區(qū)域人類活動(dòng)越少等[26-27]，因此現(xiàn)有方法在地形地貌復(fù)雜且海拔坡度跨度較大的地區(qū)具有一定局限性，有必要嘗試借助三維空間來(lái)量化地形因子對(duì)人為干擾的地理影響。東莞市生態(tài)空間多分布在山地區(qū)域，林地、草地、水體所處的地形位置也在一定程度上決定了其是否易受人為干擾，所以地形因子對(duì)人為干擾的地理影響不容忽視[28]。

本研究以東莞市現(xiàn)有以及潛在的自然保護(hù)地為研究對(duì)象，綜合研判研究區(qū)特征，從人為干擾評(píng)價(jià)的視角出發(fā)，以地形因子修正基于土地覆被類型的單一評(píng)價(jià)體系，通過(guò)三維視角評(píng)估東莞市人為干擾的空間分布強(qiáng)度，科學(xué)評(píng)估自然保護(hù)地保護(hù)現(xiàn)狀，并結(jié)合空間自相關(guān)分析方法，識(shí)別可能存在的保護(hù)空缺區(qū)域，為有效保護(hù)生態(tài)環(huán)境、合理開發(fā)自然資源、推動(dòng)自然經(jīng)濟(jì)社會(huì)良性發(fā)展、調(diào)控人類行為提供重要的決策依據(jù)。

1 研究區(qū)概況

東莞市位于廣東省中南部，地處珠江口東岸，全市陸地面積2 460.1 km2，海域面積82.6 km2，屬于典型的南亞熱帶季風(fēng)氣候。東莞市自然山水格局可概括為“東山、西水、南海、北江”，地勢(shì)東南高、西北低，地貌以丘陵臺(tái)地、沖積平原為主，最高峰銀瓶山海拔898.2 m，位于東莞市東部山區(qū)（圖1）。全市森林生態(tài)系統(tǒng)分布廣泛，森林覆蓋率41.83%，集中分布在東部和南部的低山丘陵區(qū)域。現(xiàn)有區(qū)、縣級(jí)以上自然保護(hù)地36處（圖2），批復(fù)總面積46 958.91 hm2，占東莞市國(guó)土面積的18.46%，其中各個(gè)自然保護(hù)地之間的空間重疊面積高達(dá)10 339.12 hm2，自然保護(hù)地空間規(guī)劃體系亟待更新完善。

1 東莞市海拔分布Altitude distribution in Dongguan City

2 東莞市現(xiàn)有自然保護(hù)地分布Distribution of existing protected areas in Dongguan City

2 研究方法

2.1 數(shù)據(jù)來(lái)源與處理

2.1.1 土地覆被分類與整合

本研究使用中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心（www.resdc.cn）提供的2018年?yáng)|莞市土地覆被數(shù)據(jù)（空間分辨率為30 m），該數(shù)據(jù)依據(jù)國(guó)家土地覆被分類方法，結(jié)合已建立的土地利用和土地覆被變化（land use and land cover change, LUCC）分類系統(tǒng)，通過(guò)衛(wèi)星圖像遙感解譯將中國(guó)土地覆被類型劃分為林地、草地、耕地、建設(shè)用地、水域、未利用地6種一級(jí)類型和25種二級(jí)類型，東莞市土地覆被類型包括5種一級(jí)類型和14種二級(jí)類型。

由于東莞市位于亞熱帶常綠闊葉林區(qū)，有林地（郁閉度＞30%的天然林和人工林）分布面積廣泛，因此為得到盡可能精確的研究結(jié)果，本研究保留已獲取的東莞市土地覆被數(shù)據(jù)中草地、耕地、建設(shè)用地、水域的分類，并參照前人對(duì)南方林區(qū)的研究[29]，結(jié)合東莞市土地覆被特點(diǎn)，將該數(shù)據(jù)中的林地分為有林地和其他林地（即非有林地），有林地按照最大似然法分為常綠闊葉林、針葉林、針闊混交林，其他林地直接保留現(xiàn)有數(shù)據(jù)分類，即分為灌木林、疏林和園地、苗圃。有林地分類使用的遙感影像圖獲取網(wǎng)站為地理空間數(shù)據(jù)云（www.gscloud.cn），首先根據(jù)研究時(shí)間、尺度需要，選擇地面站接收時(shí)間為2020年2月且云量低于10%的Landsat 8遙感數(shù)據(jù)，其次利用ENVI 5.3軟件將有林地進(jìn)行二次分類，分類步驟包括圖像預(yù)處理（精確配準(zhǔn)、鑲嵌裁剪、正射校正、輻射定標(biāo)以及大氣校正等），并結(jié)合目視解譯進(jìn)行監(jiān)督分類，最后根據(jù)現(xiàn)場(chǎng)調(diào)研進(jìn)行修正，得到有林地分類結(jié)果，經(jīng)整合得到最終的東莞市土地覆被分類結(jié)果。

2.1.2 地形因子選取與分級(jí)

本研究需綜合考量并選擇合理的地形因子修正基于土地覆被類型的人為干擾評(píng)價(jià)。基于前人對(duì)人為干擾與地形因子的關(guān)系研究[23-25]，并通過(guò)對(duì)東莞市地形地貌的研判可知，海拔、坡度等單一因子雖然能夠在一定程度上反映人為干擾的難易程度，但單一因子缺乏評(píng)價(jià)的針對(duì)性，無(wú)法直觀地反映地貌類型，因而無(wú)法直觀地體現(xiàn)人為干擾影響。綜合以上判斷，結(jié)合有關(guān)專家建議，本研究選擇坡位和地形起伏度2種與人為干擾直接相關(guān)的地形因子。坡位是指研究點(diǎn)所處地形縱剖面的上下位置[30]，該因子在地形位置的基礎(chǔ)上綜合了坡度信息，通過(guò)地形、地貌劃分，消除了絕對(duì)高程特征在區(qū)域尺度范圍的影響[31]，分為山谷、下坡、平坡、中坡、上坡、山脊6種類型；以地貌受人為干擾的難易程度來(lái)說(shuō)，平坡、下坡、山谷極易受干擾，中坡較易受干擾，上坡不易受干擾，山脊極不易受干擾。地形起伏度是指在所指定的分析區(qū)域內(nèi)所有柵格中最大高程與最小高程的差[26]，其數(shù)值取決于柵格大小，并在宏觀層面反映了研究區(qū)地貌被人為開發(fā)建設(shè)的難易程度，在等面積的柵格中，數(shù)值越大，代表該柵格所在區(qū)域越不易被人為開發(fā)建設(shè)。

本研究擬構(gòu)建三維人為干擾評(píng)價(jià)體系，數(shù)據(jù)來(lái)源于國(guó)家測(cè)繪局的數(shù)字高程模型（digital elevation model, DEM），分辨率為30 m×30 m。具體利用ArcGIS 10.2軟件從DEM數(shù)據(jù)中分別提取研究所需的地形因子，根據(jù)劃分尺度和研究區(qū)地形特征，結(jié)合自然間斷點(diǎn)分級(jí)法確定數(shù)量分級(jí)梯度。

2.2 人為干擾度計(jì)算

2.2.1 人為干擾評(píng)價(jià)指標(biāo)體系構(gòu)建

人為干擾指數(shù)的確定依據(jù)是人類活動(dòng)頻率和人類活動(dòng)程度對(duì)生態(tài)系統(tǒng)或景觀的綜合影響程度[9]。本研究基于對(duì)各項(xiàng)土地覆被指標(biāo)天然性的整體判斷并參考前人研究結(jié)果[9,21,32]構(gòu)建人為干擾評(píng)價(jià)指標(biāo)體系，該指標(biāo)體系的權(quán)重賦值主要采用層次分析法結(jié)合專家打分法[33]。首先構(gòu)造判斷矩陣，分為目標(biāo)層、準(zhǔn)則層、因素層、指標(biāo)層，利用1～9標(biāo)度選擇方法建立指標(biāo)體系權(quán)重賦值調(diào)查問(wèn)卷，并發(fā)放給12位具有自然保護(hù)地研究背景的專家；其次將調(diào)查結(jié)果導(dǎo)入Yaahp 10.2軟件并進(jìn)行層次總排序及一致性檢驗(yàn)；最后由此確定各項(xiàng)評(píng)價(jià)指標(biāo)權(quán)重。

2.2.2 研究網(wǎng)格構(gòu)建

本研究主要應(yīng)用網(wǎng)格分析法進(jìn)行人為干擾度的計(jì)算。網(wǎng)格應(yīng)用中需要強(qiáng)調(diào)指出的是，景觀格局中尺度效應(yīng)的存在已被諸多研究證實(shí)[34-35]，即網(wǎng)格選擇的科學(xué)性直接關(guān)系到人為干擾度的計(jì)算結(jié)果，因此由地圖空間模式分析軟件Fragstats 4.2來(lái)確定網(wǎng)格尺度[36]。具體以東莞市土地覆被分類圖進(jìn)行景觀斑塊面積的計(jì)算，并根據(jù)景觀指數(shù)計(jì)算中獲得的斑塊平均面積，以其2～5倍設(shè)置自然保護(hù)地的人為干擾分析網(wǎng)格面積。

2.2.3 人為干擾度計(jì)算和分級(jí)

人為干擾度計(jì)算方法來(lái)源于Steinhardt等[18]提出的人為干擾度M計(jì)算公式。根據(jù)式（1）計(jì)算自然保護(hù)地每個(gè)網(wǎng)格單元的人為干擾度[21]，并繪制空間分異圖：

式中：M為人為干擾度，Hi為第i類土地覆被類型或地形因子類型的干擾指數(shù)，Si為第i類土地覆被類型或地形因子類型的面積，S為網(wǎng)格總面積。

在計(jì)算網(wǎng)格單元人為干擾度的基礎(chǔ)上，本研究使用ArcGIS 10.2軟件的自然間斷點(diǎn)分級(jí)法進(jìn)行人為干擾度的數(shù)量分級(jí)。

2.3 空間自相關(guān)分析

確定保護(hù)優(yōu)先區(qū)的方法因研究區(qū)尺度而異，利用空間自相關(guān)分析確定小尺度范圍內(nèi)保護(hù)優(yōu)先區(qū)是目前較為公認(rèn)的方法[7,37]?？臻g自相關(guān)是指同一變量在不同空間位置上的相關(guān)性，可使用局部空間自相關(guān)分析計(jì)算每一單元與其鄰近單元就某一屬性的相似性，用以確定和識(shí)別空間聚集的分布[7]。本研究使用局部空間自相關(guān)分析中的Getis-Ord Gi*系數(shù)，該系數(shù)是基于距離權(quán)矩陣識(shí)別高值或低值要素在空間上發(fā)生聚類的位置，即具有統(tǒng)計(jì)顯著性的熱點(diǎn)與冷點(diǎn)空間聚類的指標(biāo)[7]。熱點(diǎn)應(yīng)具有高值，且被其他同樣具有高值的要素所包圍，冷點(diǎn)反之。該系數(shù)可用來(lái)檢測(cè)人為干擾在空間上表現(xiàn)的聚集程度[38]。本研究利用ArcGIS 10.2軟件對(duì)人為干擾熱點(diǎn)與冷點(diǎn)的區(qū)域分布特征進(jìn)行識(shí)別與分析，并與東莞市現(xiàn)有自然保護(hù)地空間分布進(jìn)行對(duì)比，從而提出東莞市自然保護(hù)地范圍調(diào)整的相關(guān)建議（圖3）。

3 研究技術(shù)路線Technical route for this research

3 結(jié)果與分析

3.1 數(shù)據(jù)處理結(jié)果

3.1.1 土地覆被分類與整合結(jié)果

通過(guò)數(shù)據(jù)的分類與整合，得到東莞市土地覆被類型有6個(gè)大類和16個(gè)小類（圖4）。由分類結(jié)果可知，東莞市北部、中部、西南、東南遍布城鎮(zhèn)用地，而林地集中分布在南部及東部，西部則多河渠、水田、園地等類型，并與工交建設(shè)用地交叉分布。

3.1.2 地形因子提取與分級(jí)結(jié)果通過(guò)ArcGIS 10.2軟件對(duì)東莞市坡位和地形起伏度分布進(jìn)行提取與分級(jí)（圖5、6）。對(duì)照東莞市海拔分布圖（圖1）可知，坡位成片分布地區(qū)與地形起伏度高值區(qū)域集中在中西部和東部2座山脈，即大嶺山和銀瓶山；此外，2種地形因子的分布與土地覆被類型分布模式基本相似，尤其是天然程度高的類型，如有林地一般分布在中坡、上坡、山脊集中的區(qū)域和地形起伏度高的區(qū)域。

5 東莞市坡位分布Distribution of slope position in Dongguan City

6 東莞市地形起伏度分布Distribution of topographic relief in Dongguan City

3.2 東莞市人為干擾評(píng)價(jià)指標(biāo)體系及權(quán)重賦值

本研究在前人研究基礎(chǔ)上，針對(duì)東莞市自然特征，建立2個(gè)維度、4個(gè)層次（目標(biāo)層、準(zhǔn)則層、因素層、指標(biāo)層）的人為干擾評(píng)價(jià)指標(biāo)體系（表1）。

指標(biāo)層的權(quán)重賦值均有相關(guān)依據(jù)，如常綠闊葉林作為東莞市最主要的植被類型，大多為天然植被；而針葉林多為人工種植的馬尾松林等次生林，抗干擾能力低于常綠闊葉林；其他林地整體人為干擾指數(shù)高于有林地，而灌木林由于多為天然形成的矮林地和灌叢林地，且郁閉度大于40%，因此人為干擾指數(shù)較低；東莞市高山丘陵地區(qū)的草地主要為天然形成，而平原區(qū)域大面積草地多為高爾夫球場(chǎng)內(nèi)人工種植的草地，存在一定人為干擾；東莞市水系發(fā)達(dá)，西北區(qū)域天然河網(wǎng)密布，中部人工開挖水庫(kù)、坑塘較多，總體來(lái)說(shuō)河渠的人為干擾指數(shù)小于水庫(kù)、坑塘。結(jié)合專家打分結(jié)果，利用Yaahp 10.2軟件計(jì)算各層次判斷矩陣，目標(biāo)層判斷矩陣的最大特征值λ=2，一致性比率CR=0（＜0.1），一致性檢驗(yàn)通過(guò)。最終建立人為干擾評(píng)價(jià)指標(biāo)體系（表1）。

3.3 東莞市人為干擾評(píng)價(jià)結(jié)果

根據(jù)Fragstats 4.2軟件的計(jì)算，研究區(qū)景觀斑塊平均面積為0.527 km2，取其2倍的整數(shù)值1 km2，利用ArcGIS 10.2軟件將東莞市劃分為1 km×1 km的等面積柵格后對(duì)各評(píng)價(jià)指標(biāo)進(jìn)行賦值，再疊加權(quán)重得到最終人為干擾度分布結(jié)果，按照自然間斷點(diǎn)分級(jí)法進(jìn)行人為干擾度7級(jí)劃分（圖7）。

7 東莞市人為干擾度分布Distribution of human disturbance in Dongguan City

整體來(lái)看，東莞市人為干擾度呈現(xiàn)北高南低的分布模式。根據(jù)東莞市行政區(qū)域（圖2）可知，最為明顯的差異是分布在北部大片城區(qū)的城鎮(zhèn)用地和分布于南部大嶺山鎮(zhèn)、樟木頭鎮(zhèn)、清溪鎮(zhèn)等山地丘陵區(qū)的林地，其人為干擾度分別為最高和最低水平；西部麻涌鎮(zhèn)、沙田鎮(zhèn)天然水域較多，但由于其主要處于下坡和平坡等易受人為干擾的坡位，且地形起伏度較低，因此人為干擾度基本處于中間水平；南部低人為干擾度區(qū)域主要呈現(xiàn)出東西向的連續(xù)成片分布，其中每一片區(qū)域人為干擾度由中心向外圍呈現(xiàn)由低到高的趨勢(shì)?？傮w上，地形起伏度高、坡位復(fù)雜的山區(qū)人為干擾度低于平原地區(qū)，城市南部人為干擾度明顯高于北部，城市中部人為干擾度高于東西兩側(cè)的城市邊緣。

3.4 基于人為干擾評(píng)價(jià)的東莞市自然保護(hù)地空間優(yōu)化策略

就自然保護(hù)地的保護(hù)目的而言，可以認(rèn)為其人為干擾越少，保護(hù)狀態(tài)越好。在得到東莞市人為干擾度分布結(jié)果的基礎(chǔ)上，進(jìn)一步通過(guò)局部空間自相關(guān)分析得出人為干擾冷熱點(diǎn)區(qū)域分布，并將其與東莞市現(xiàn)有自然保護(hù)地進(jìn)行空間疊加（圖8），進(jìn)行自然保護(hù)地的保護(hù)狀態(tài)分析與空間優(yōu)化策略探討。

8 東莞市現(xiàn)有自然保護(hù)地邊界疊加人為干擾冷、熱點(diǎn)分布Distribution of the boundaries of existing protected areas superimposed with areas featuring little and much human disturbance in Dongguan City

由人為干擾冷、熱點(diǎn)的區(qū)域分布可知，人為干擾熱點(diǎn)區(qū)域集中成片分布在東莞市主城區(qū)以及周邊寮步鎮(zhèn)、橫瀝鎮(zhèn)，加上西南部厚街鎮(zhèn)、虎門鎮(zhèn)、長(zhǎng)安鎮(zhèn)和東南部塘廈鎮(zhèn)等城鎮(zhèn)建設(shè)區(qū)域；冷點(diǎn)區(qū)域集中于南部大嶺山鎮(zhèn)、大朗鎮(zhèn)、黃江鎮(zhèn)、樟木頭鎮(zhèn)、清溪鎮(zhèn)等山林地區(qū)，總體上熱點(diǎn)區(qū)域在冷點(diǎn)區(qū)域四周環(huán)繞分布。結(jié)合東莞市及其周邊的土地覆被格局分析可知，東莞市西部為平原地區(qū)，雖然水域分布較廣，可帶來(lái)生態(tài)效益，但其地形特點(diǎn)及地理位置決定了周邊建設(shè)用地極易擴(kuò)張，人為干擾可能進(jìn)一步向東蔓延；此外，雖然東莞市南部、東部山脈向外延伸，可形成一大生態(tài)屏障，但其周圍遍布廣州、深圳、惠州等一系列粵港澳大灣區(qū)的經(jīng)濟(jì)發(fā)展重鎮(zhèn)，城鎮(zhèn)擴(kuò)張無(wú)疑會(huì)加大生態(tài)屏障承載的環(huán)境壓力，使得自然保護(hù)地邊界進(jìn)一步受到蠶食。

通過(guò)疊加?xùn)|莞市現(xiàn)有自然保護(hù)地邊界可知，目前的自然保護(hù)地空間分布基本位于人為干擾冷點(diǎn)區(qū)域，說(shuō)明目前自然保護(hù)地的保護(hù)現(xiàn)狀基本良好。但仍存在部分未重疊的空間，主要集中于大朗鎮(zhèn)東部、樟木頭鎮(zhèn)南部的山林區(qū)域中間位置，總面積約60 km2，占東莞市現(xiàn)有自然保護(hù)地總面積的12%，是東莞市單個(gè)自然保護(hù)地平均面積的4.38倍，若可納入自然保護(hù)地范圍，將為筑牢區(qū)域生態(tài)屏障提供重要價(jià)值。進(jìn)一步分析可知，未處于自然保護(hù)地范圍的人為干擾冷點(diǎn)區(qū)域現(xiàn)狀用地為廣東省樟木頭林場(chǎng)等國(guó)有林場(chǎng)，且基本已納入生態(tài)公益林，管理與保護(hù)條件較成熟；現(xiàn)狀林地類型主要為南亞熱帶季風(fēng)常綠闊葉林，具有較高保護(hù)價(jià)值；另外，該區(qū)域周邊緊鄰巍峨山森林公園和黃牛埔森林公園，二者均位于黃江鎮(zhèn)，且部分邊界范圍與黃江鎮(zhèn)行政邊界重合，可見早期自然保護(hù)地的范圍劃定在一定程度上受行政邊界影響。

此外，本研究還疊加得出一些未處于人為干擾冷點(diǎn)區(qū)域的現(xiàn)有自然保護(hù)地，一方面與研究尺度有關(guān)，由于本研究采用的網(wǎng)格為1 km×1 km，因此零星散布的小面積自然保護(hù)地所處網(wǎng)格無(wú)法成為人為干擾冷點(diǎn)區(qū)域；另一方面可能表示現(xiàn)有自然保護(hù)地存在一定問(wèn)題，如中西部和中南部少量區(qū)域以及東南部帶狀區(qū)域都位于山地邊緣，且緊鄰大面積城鎮(zhèn)建設(shè)用地，其處于人為干擾冷點(diǎn)區(qū)域邊緣以外，說(shuō)明目前可能已存在建設(shè)用地邊界擴(kuò)張倒逼自然保護(hù)地邊界后退的風(fēng)險(xiǎn)。

因此，對(duì)東莞市自然保護(hù)地提出3個(gè)方面的空間優(yōu)化策略。1）可遵循“應(yīng)保盡?！钡脑瓌t，將具備自然保護(hù)地設(shè)立條件的南部大面積人為干擾冷點(diǎn)區(qū)域納入自然保護(hù)地體系，使其連通東西兩側(cè)的大嶺山地區(qū)和銀瓶山地區(qū)，為區(qū)域生態(tài)屏障的功能加強(qiáng)和自然保護(hù)地的整合優(yōu)化提供基礎(chǔ)；2）應(yīng)進(jìn)一步綜合評(píng)估現(xiàn)有自然保護(hù)地邊界的保護(hù)情況，并對(duì)接國(guó)土空間規(guī)劃進(jìn)行進(jìn)一步研究，疊加生態(tài)紅線、基本農(nóng)田等劃定，優(yōu)化調(diào)整自然保護(hù)地邊界，嚴(yán)控城鎮(zhèn)開發(fā)邊界；3）對(duì)于散布的小面積自然保護(hù)地，可進(jìn)一步評(píng)估其人為干擾程度，并確定其自然保護(hù)等級(jí)。

4 結(jié)論與展望

本研究構(gòu)建了三維層面的東莞市人為干擾評(píng)價(jià)體系，并將人為干擾冷、熱點(diǎn)區(qū)域與現(xiàn)有、自然保護(hù)地邊界進(jìn)行空間疊加，為東莞市自然保護(hù)地空間規(guī)劃的更新、完善提供了新的思路。本研究的創(chuàng)新點(diǎn)在于將普遍應(yīng)用的人為干擾研究方法精細(xì)化、三維化，并基于此進(jìn)行自然保護(hù)地空間優(yōu)化。1）基于對(duì)研究區(qū)域特性的研判，進(jìn)行了林地二級(jí)分類與圖像解譯；2）基于對(duì)地形因子的研究，確定了坡位和地形起伏度2種影響人為干擾的指標(biāo)，從而構(gòu)建了更為完善的人為干擾評(píng)價(jià)指標(biāo)體系，科學(xué)權(quán)衡并進(jìn)行指標(biāo)賦值，得到研究結(jié)果；3）基于人為干擾冷點(diǎn)區(qū)域分布，對(duì)自然保護(hù)地保護(hù)現(xiàn)狀進(jìn)行科學(xué)評(píng)估，提出自然保護(hù)地空間優(yōu)化策略。

本研究的不足之處在于以土地覆被類型為主要依據(jù)的人為干擾評(píng)價(jià)并未將建設(shè)用地類型進(jìn)行具體細(xì)分，如未將道路交通及其輻射范圍等因素納入考慮，人為干擾評(píng)價(jià)指標(biāo)的綜合性、全面性有待提升；同時(shí)該指標(biāo)體系在應(yīng)用于大面積林地時(shí)存在林下人為干擾無(wú)法識(shí)別的情況（例如為林蔭所覆蓋的道路、居民點(diǎn)、景點(diǎn)等），需進(jìn)一步結(jié)合實(shí)地調(diào)研，細(xì)化空間網(wǎng)格單元，區(qū)分人為干擾的類型、位點(diǎn)和影響范圍，進(jìn)一步提高人為干擾評(píng)價(jià)的精度和準(zhǔn)確度。此外，盡管本研究證明了人為干擾評(píng)價(jià)可以作為自然保護(hù)地范圍劃定的重要輔助工具，但仍然存在使用單一評(píng)價(jià)體系無(wú)法識(shí)別空間生態(tài)價(jià)值高低的弊端。未來(lái)研究需綜合考慮生態(tài)系統(tǒng)完整性、生態(tài)系統(tǒng)服務(wù)價(jià)值、生態(tài)系統(tǒng)脆弱性等與自然保護(hù)地生態(tài)保護(hù)與資源利用密切相關(guān)的評(píng)價(jià)指標(biāo)，進(jìn)一步完善自然保護(hù)地空間優(yōu)化的科學(xué)依據(jù)。

致謝(Acknowledgments)：

感謝東莞市林業(yè)局對(duì)本研究的大力支持。

圖表來(lái)源(Sources of Figures and Table)：

文中圖表均由作者繪制，其中圖1、2、4～8的地圖底圖來(lái)源于天地圖平臺(tái)的東莞市標(biāo)準(zhǔn)地圖，審圖號(hào)為粵S（2018）016，圖1的高程數(shù)據(jù)由國(guó)家測(cè)繪局提供，圖2、8的自然保護(hù)地?cái)?shù)據(jù)由東莞市林業(yè)局提供。