張麗芳， 廖 雨， 楊存建， 冉丹陽

(1.四川師范大學 西南土地評價與監(jiān)測教育部重點實驗室，四川 成都 610068； 2.四川師范大學 地理與資源科學學院， 四川 成都 610068)

人類活動是導致土地利用類型發(fā)生變化的主要影響因子，對生態(tài)用地的安全性也構(gòu)成一定的威脅。自然保護區(qū)作為生物多樣性豐富、生態(tài)功能最強、水土涵養(yǎng)能力好、珍稀動植物的天然生存場所。如何協(xié)調(diào)人類活動用地與自然保護區(qū)健康可持續(xù)發(fā)展是當下面臨的可持續(xù)發(fā)展的新挑戰(zhàn)[1-2]。

近幾年來對生態(tài)安全的保護與關(guān)注度一直處于熱點，構(gòu)建生態(tài)安全格局備受中國學者的高度關(guān)注，俞孔堅[3]在Forman[4]的理論基礎(chǔ)上提出了景觀安全格局概念，創(chuàng)建識別源地、建立阻力面、構(gòu)建安全格局的模式。該模式被國內(nèi)外學者廣泛應(yīng)用于區(qū)域案例的生態(tài)安全評價構(gòu)建，并且取得了較好的效果，極大的推動了此方法的科學化、成熟化、應(yīng)用化[5-8]。針對自然保護區(qū)的景觀生態(tài)安全評價國內(nèi)也有所研究，胡艷等[9]選取景觀指數(shù)，分析了寬闊水自然保護區(qū)的景觀格局演變。付夢娣等[10]利用最小累計阻力模型對秦嶺地區(qū)的自然保護區(qū)進行了網(wǎng)絡(luò)構(gòu)建與優(yōu)化，確認出秦嶺地區(qū)的生態(tài)源地面積，對后人規(guī)劃與優(yōu)化保護區(qū)提供了參考價值。目前，針對生態(tài)安全構(gòu)建的方法多樣[11]，研究取得了一定的成果，但基于景觀生態(tài)安全理論針對自然保護區(qū)內(nèi)的人類活動變化安全與評價的研究較少。

本研究擬以壽城自然保護區(qū)2013年0.5 m空間分辨率的有人機航攝影像和2017年2 m空間分辨率的高分一號遙感影像為數(shù)據(jù)源，采用人機交互式判讀方法提取壽城自然保護區(qū)兩期土地利用分類圖層，借助GIS平臺和 Fragstat 4.2軟件，從景觀尺度、生態(tài)安全格局構(gòu)建的角度出發(fā),分析壽城自然保護區(qū)的景觀空間分布特征，并對不同安全水平下的人類活動用地進行評價與優(yōu)化。以期為壽城自然保護區(qū)的生態(tài)功能的最優(yōu)化、景觀多樣性與人類活動的區(qū)域和諧發(fā)展而服務(wù),并為壽城自然保護區(qū)合理保護生態(tài)用地、維護自然保護區(qū)安全、維系生物多樣性、合理規(guī)劃建設(shè)用地提供決策建議。

1 研究區(qū)概況與數(shù)據(jù)來源

1.1 研究區(qū)概況

壽城省級自然保護區(qū)位于廣西壯族自治區(qū)的永福縣北部和臨桂縣西北部，于1982年經(jīng)廣西壯族自治區(qū)人民政府批準建立的，經(jīng)2017年遙感解譯統(tǒng)計保護區(qū)面積范圍為70 139.59 hm2。地理坐標為東經(jīng)109°38′17″—109°53′36″, 北緯 25°02′16″—25°29′36″，最高海拔1 498 m; 最低海拔179 m，年均溫為17.2 ℃，年平均降水量約為2 500 mm[12]。

保護區(qū)內(nèi)溝壑縱橫，土地利用類型多樣，土地分類標準依據(jù)《自然保護區(qū)人類活動遙感監(jiān)測技術(shù)指南(試行)》的通知環(huán)辦〔2014〕12號文件，自然保護區(qū)土地分類一級分類為耕地、居民點、工礦用地、采石場、能源設(shè)施、旅游設(shè)施、交通設(shè)施、養(yǎng)殖場、道路、其他人工設(shè)施(為無法劃分到其余9類的設(shè)施用地)等10種一級土地分類，下分41種二級分類，本文按照實際研究需求，按照一級人類活動分類標準并疊加林地、未利用地、河流進行劃分土地利用類型，其中一級分類共計7種，并疊加林地、未利用地、河流3種土地利用類型。

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

壽城自然保護區(qū)2013，2017年土地利用類型研究數(shù)據(jù)來源于廣西省保護區(qū)遙感監(jiān)測項目。數(shù)據(jù)基于GIS軟件平臺，采用了2013年0.5 m空間分辨率的有人機航攝影像、2017年2 m空間分辨率的高分一號影像，經(jīng)過圖像預(yù)處理以及人工目視解譯后獲得。全球數(shù)字高程數(shù)據(jù)下載于地理空間數(shù)據(jù)云網(wǎng)站(http:∥www.gscloud.cn/)，數(shù)據(jù)經(jīng)過拼接、裁剪得到壽城自然保護區(qū)30 m分辨率的DEM數(shù)據(jù)，運用GIS軟件計算出坡度并進行重分類處理；此外，壽城自然保護區(qū)行政邊界、自然保護區(qū)邊界等數(shù)據(jù)來源于廣西省自然保護區(qū)遙感監(jiān)測項目。

2 研究方法

研究根據(jù)壽城自然保護區(qū)的地理特征，首先選取景觀指數(shù)對保護區(qū)內(nèi)的土地利用類型進行景觀格局演變分析。然后以景觀安全格局為理論背景，以MCR模型為依托，識別源地，選取距居民點距離、距河流距離、距道路距離、海拔、坡度等5個阻力因子創(chuàng)建阻力面，利用最小阻力模型得出保護區(qū)的不同生態(tài)安全水平下的等級最小累積阻力面，與人類活動用地疊加，以此來對保護區(qū)內(nèi)的人類活動用地進行評價與優(yōu)化。

2.1 景觀格局指數(shù)選取

景觀指數(shù)是景觀生態(tài)學研究的最重要內(nèi)容，選取高度濃縮景觀格局信息、反映其結(jié)構(gòu)組成和空間配置某些方面特征的景觀指標[13-15]，來研究景觀結(jié)構(gòu)組成特征和空間配置關(guān)系，是景觀生態(tài)學研究中使用廣泛且較為成熟的定量研究方法。研究結(jié)合前人經(jīng)驗與自然保護區(qū)的區(qū)域特點，堅持常用、全面、系統(tǒng)等原則，選取了斑塊數(shù)量(NP)、斑塊密度(PD)、最大斑塊指數(shù)(LPI)、景觀形狀指數(shù)(LSI)、聚集度(AI)、邊緣密度(ED)、香農(nóng)多樣性指數(shù)(SHDI)、香農(nóng)均勻度指數(shù)(SHEI)這8個景觀指標,從土地利用景觀斑塊個數(shù)及密度、斑塊形狀、斑塊空間結(jié)構(gòu)等方面來描述壽城自然保護區(qū)土地利用景觀格局變化特征[16-17]。景觀指數(shù)的計算運用GIS軟件平臺，將壽城自然保護區(qū)土地利用類型數(shù)據(jù)轉(zhuǎn)換為TIFF格式，并將數(shù)據(jù)加載至Fragstats4.2軟件，得出壽城自然保護區(qū)各年各景觀格局指數(shù)。

2.2 生態(tài)安全格局構(gòu)建

2.2.1 保護區(qū)源地確定 生態(tài)源地是指作為區(qū)域生態(tài)功能最佳、能量與物質(zhì)流通性好、水源涵養(yǎng)功能強的區(qū)域[18-19]，是景觀生態(tài)安全質(zhì)量最佳的核心地帶。本研究是基于生態(tài)安全格局下對保護區(qū)人類活動進行的評價優(yōu)化。因此，保護區(qū)源地的確定優(yōu)先選取生態(tài)服務(wù)功能最佳、生態(tài)質(zhì)量較好、受人類活動影響較小的林地與水域。

2.2.2 阻力面構(gòu)建 生態(tài)源地間生物水平方向的空間運動與生態(tài)競爭演變過程，主要是克服人類活動與自然環(huán)境的各個阻力實現(xiàn)的，阻力面反映了生態(tài)用地克服阻力的空間連通性[19]。研究利用地理學中的表面模型——最小累積阻力(MCR)模型建立最小累計阻力面，它反映了從源地出發(fā)到達空間的某一個點所克服的最小阻力、成本耗費最小的距離[7]。其公式如下：

(1)

式中：f——區(qū)域中任意一點的最小阻力與其到所有源地的距離和景觀特征的正相關(guān)關(guān)系[8,20]；Dij——從空間中的某一景觀i到源地j的距離；Ri——景觀i的阻力值。最小累積阻力模型是根據(jù) Knaapen等人提出的模型和GIS中的成本距離發(fā)展得來的[21-23]。該模型主要利用GIS軟件中成本距離工具分析完成。

2.2.3 安全格局構(gòu)建下的人類活動識別 本研究根據(jù)各阻力面與源地進行成本距離分析，劃分出不同等級的生態(tài)安全水平用地。然后，將生態(tài)安全格局水平圖層與人類活動用地(耕地、道路、居民點、養(yǎng)殖場、采石場、能源用地、其他人工設(shè)施用地)斑塊圖層疊加，得出不同安全水平的人類活動景觀安全格局，并對人類活動的發(fā)展建設(shè)提出相應(yīng)的優(yōu)化建議。

3 結(jié)果與分析

3.1 景觀指數(shù)動態(tài)變化

3.1.1 土地利用景觀格局總體變化 通過計算壽城自然保護區(qū)2013年和2017年包括斑塊個數(shù)、斑塊邊緣密度、聚集度指數(shù)、香農(nóng)多樣性指數(shù)、香農(nóng)均勻性指數(shù)等5個景觀指數(shù)，來把握壽城自然保護區(qū)土地利用景觀格局總體變化特征，計算結(jié)果如表1所示。

表1 壽城自然保護區(qū)景觀多樣性指數(shù)

從表1中可以看出壽城自然保護區(qū)土地利用斑塊個數(shù)和斑塊邊緣密度均呈上升的趨勢，與2013年相比，2017年壽城自然保護區(qū)土地利用斑塊個數(shù)增加了83個，斑塊邊緣密度增加了0.304 5 m/hm2，表明壽城自然保護區(qū)土地利用斑塊形狀趨于復雜，斑塊破碎化程度越來越高。此外，斑塊聚集度指數(shù)呈下降的趨勢，表明壽城自然保護區(qū)土地利用斑塊在空間分布上趨于分散。壽城自然保護區(qū)土地利用景觀香農(nóng)多樣性指數(shù)和香農(nóng)均勻性指數(shù)均呈上升趨勢，說明壽城自然保護區(qū)土地利用景觀復雜程度和均勻程度均有所增加。

3.1.2 不同土地利用類型景觀格局變化 為進一步探討壽城自然保護區(qū)土地利用景觀格局動態(tài)變化趨勢，計算了2013年—2017年壽城自然保護區(qū)各土地利用景觀類型斑塊數(shù)量及密度、最大斑塊指數(shù)、斑塊邊緣密度、斑塊形狀指數(shù)、斑塊聚集度等6個指標用于研究保護區(qū)的景觀動態(tài)變化特征(表2)。

從表2中可以看出，壽城自然保護區(qū)道路、居民點、水域等土地利用類型斑塊個數(shù)最多，耕地和林地次之，其余地類斑塊個數(shù)較少，與2013年相比，2017年壽城自然保護區(qū)居民點斑塊大幅度增加，其次為道路、其他人工設(shè)施。從斑塊密度變化幅度來看，居民點的斑塊密度增長最為明顯(0.08)，道路斑塊次之(0.02)。從最大斑塊指數(shù)來看，林地的最大斑塊指數(shù)最大，而其余地類最大斑塊指數(shù)均不足3，表明林地在壽城自然保護區(qū)土地利用類型中占據(jù)絕對優(yōu)勢。從斑塊形狀變化來看，各地類斑塊邊緣密度和除道路外的各地類景觀形狀指數(shù)均有所增加，其中居民點斑塊邊緣密度和景觀形狀指數(shù)變化最大，分別增加了0.25和1.12，道路的景觀形狀指數(shù)下降了0.16，表明道路的通達性有所增強。壽城自然保護區(qū)各土地利用類型聚集度指數(shù)均呈下降趨勢，表明各土地利用類型斑塊空間分布上趨于分散，斑塊破碎化程度增加。

表2 壽城自然保護區(qū)不同用地類型的景觀指數(shù)

注：NP為斑塊數(shù)量；PD為斑塊密度；LPI為最大斑塊指數(shù)；LSI為景觀形狀指數(shù)；AI為聚集度；ED為邊緣密度；SHDI為香農(nóng)多樣性指數(shù)；SHEI為香農(nóng)均勻度指數(shù)。

總之，壽城自然保護區(qū)受人類活動的影響較為強烈，主要表現(xiàn)為道路、居民點、其他人工設(shè)施等人工景觀的破碎化程度增加，空間規(guī)模在擴大，而林地、水域等自然景觀趨于萎縮，斑塊個數(shù)減少。各土地利用類型空間分布趨于分散，斑塊形狀更加復雜，景觀多樣性和景觀均勻性均有所提高。

3.2 生態(tài)安全格局構(gòu)建

3.2.1 保護區(qū)源地識別 本研究根據(jù)2017年壽城自然保護區(qū)的各個地類的實際分布情況，將坡度分為0°～5°,5～10°,10～15°,15～25°,>25°等5個級別，篩選出保護區(qū)內(nèi)分布在坡度>25°較為連續(xù)的林地與水域作為保護區(qū)源地。用GIS軟件提取相應(yīng)圖層，源地面積為29 046.53 hm2，占保護區(qū)總面積的41.41%，因其獨特的自然保護區(qū)地理特性，生態(tài)源地所占比例較大。

3.2.2 阻力面構(gòu)建研究 基于2017年數(shù)據(jù)，選取到居民點距離、到道路距離、海拔、坡度、到河流距離5項阻力因子構(gòu)建綜合阻力面。在參照前人相關(guān)研究的基礎(chǔ)上[19-20，24]，結(jié)合專家經(jīng)驗將各阻力因子的相對阻力值劃分為5個等級，并賦予阻力系數(shù)。各個因子的相對阻力值系數(shù)與權(quán)重通過專家打分并結(jié)合實際情況獲得(見表3)，其中阻力值越大、對生態(tài)用地產(chǎn)生的阻力越大，反之，阻力值越小，對生態(tài)用地產(chǎn)生的阻力越小。最后，依據(jù)表3所建立的各阻力因子權(quán)重與指標體系，運用空間分析工具對不同要素的阻力值進行加權(quán)疊加分析，最終得到壽城自然保護區(qū)綜合阻力分布圖(附圖1)。

表3 壽城自然保護區(qū)各阻力面體系、權(quán)重及阻力系數(shù)

3.2.3 不同安全水平的人類活動用地評價 根據(jù)阻力面構(gòu)建模型所述的阻力模型構(gòu)建方法，利用GIS軟件將源地、綜合阻力面添加至成本距離工具進行分析，得到綜合阻力面中每個像元到源地的最小累計成本距離，即壽城自然保護區(qū)的最小累積阻力面。將最小累積阻力面根據(jù)自然斷點法重分類為4個不同安全水平，從而得到壽城自然保護區(qū)生態(tài)安全優(yōu)化布局圖(附圖2)。其中，生態(tài)安全水平為1級的表示生態(tài)安全水平質(zhì)量最高，為保護區(qū)生態(tài)功能價值最高的地帶，應(yīng)該作為保護區(qū)的生態(tài)核心地帶，優(yōu)先加強管理與保護。

提取2003年與2017年兩期的耕地、居民點、養(yǎng)殖場、采石場、能源設(shè)施、道路、其他人工設(shè)施等人類活動矢量斑塊與重分類后的最小累計阻力面進行疊加分析，劃分出4個基于生態(tài)安全的不同安全水平的人類活動等級：發(fā)展核心區(qū)、發(fā)展緩沖區(qū)、發(fā)展優(yōu)化區(qū)、發(fā)展限制區(qū)(如附圖3所示)。其中生態(tài)安全等級為1級的為生態(tài)安全用地較好的區(qū)域，不適宜人類活動景觀的分布與建設(shè)，為人類活動用地的發(fā)展限制區(qū)；生態(tài)安全級別為4的表示生態(tài)用地影響較弱，可以優(yōu)先進行人類活動建設(shè)，為發(fā)展核心區(qū)。

根據(jù)圖1可知，2013年—2017年保護區(qū)內(nèi)的人類活動總面積主要呈現(xiàn)增加趨勢，2013年發(fā)展限制區(qū)面積為232.84 hm2，發(fā)展優(yōu)化區(qū)面積為758.10 hm2，發(fā)展過渡區(qū)為1 203.06 hm2，發(fā)展核心區(qū)為3 652.66 hm2；2017年發(fā)展限制區(qū)面積為232.75 hm2，發(fā)展優(yōu)化區(qū)面積為763.35 hm2，發(fā)展過渡區(qū)為1 212.21 hm2，發(fā)展核心區(qū)為3 670.14 hm2?？傮w來看，2013—2017年，除發(fā)展限制區(qū)人類活動面積減少外，其余等級區(qū)域的人類活動面積均有增加。為進一步分析2013—2017年的不同安全等級的人類活動用地的具體情況，現(xiàn)對4個級別的人類活動用地進行了統(tǒng)計分析(如表4所示)。

圖1 壽城自然保護區(qū)各等級下的人類活動用地面積特征

表4 壽城自然保護區(qū)各等級下的不同人類活動用地面積分布情況 hm2

(1) 發(fā)展核心區(qū)。2013年發(fā)展核心區(qū)內(nèi)人類活動總面積為3 652.66 hm2，占當年壽城自然保護區(qū)人類活動總面積的62.43%。其中人類活動用地主要以耕地、居民點用地和道路用地為主，其他人類活動用地面積較小。與2013年相比，2017年發(fā)展核心區(qū)內(nèi)人類活動總面積增加了17.48 hm2，居民點用地和其他人工設(shè)施用地分別增加了12.53，11.64 hm2，耕地減少了7.98 hm2，表明核心區(qū)內(nèi)人類活動用地變化主要以新建房屋和增加人工設(shè)施為主，居民點和其他人工設(shè)施的擴張占用耕地的現(xiàn)象明顯。發(fā)展核心區(qū)主要分布在壽城自然保護區(qū)南部以及北部和中部沿河流兩岸，該地區(qū)地勢相對平坦、人類活動較為成熟，主要表現(xiàn)為耕地面積廣，人口分布集中，對生態(tài)用地影響較小，可在此區(qū)域優(yōu)先進行居民點與工礦用地等人類活動用地建設(shè)，作為人類活動重點發(fā)展的區(qū)域，不建議進行生態(tài)功能服務(wù)的規(guī)劃投入。

(2) 發(fā)展過渡區(qū)。2013年發(fā)展過渡區(qū)內(nèi)人類活動總面積為1 203.06 hm2，占當年壽城自然保護區(qū)人類活動總面積的20.58%。其中人類活動主要為耕地、道路、居民點、其他能源設(shè)施、養(yǎng)殖場。與2013年相比，2017年發(fā)展過渡區(qū)總面積增加了9.15 hm2，耕地和居民用地面積增長明顯，分別為4.12，4.71 hm2，可以看出，發(fā)展過渡區(qū)人類活動用地變化主要以耕地和居民點的擴張為主，人口向過渡區(qū)轉(zhuǎn)移明顯。發(fā)展過渡區(qū)主要分布在發(fā)展核心區(qū)邊緣，是人類活動的過渡地帶，同時也是生態(tài)源地恢復、擴展的緩沖地帶，應(yīng)做好生態(tài)修復與環(huán)境保護等工作，對于離生態(tài)用地較近的采石場、養(yǎng)殖場進行監(jiān)測與管理，避免人類活動對生態(tài)用地產(chǎn)生污染，該區(qū)域作為發(fā)展優(yōu)化區(qū)域與發(fā)展限制區(qū)域的搬遷安置點，應(yīng)合理開發(fā)、構(gòu)建生態(tài)屏障、保障人類活動的合理與生態(tài)用地的安全。

(3) 發(fā)展優(yōu)化區(qū)。2013年發(fā)展優(yōu)化區(qū)人類活動總面積為758 hm2，占當年壽城自然保護區(qū)人類活動總面積的12.97%。主要有耕地、居民點用地、道路用地、其他能源設(shè)施用地。與2013年相比，2017年發(fā)展優(yōu)化區(qū)人類活動總面積增加了5.24 hm2，主要是由居民點用地面積增加造成，表明發(fā)展優(yōu)化區(qū)內(nèi)人類活動有所增強，人口規(guī)模呈逐漸擴大的趨勢。發(fā)展優(yōu)化區(qū)主要分布在發(fā)展核心區(qū)和發(fā)展過渡區(qū)的外圍，作為人類活動的發(fā)展優(yōu)化地帶，是生態(tài)源地恢復、擴展的緩沖地帶，該區(qū)域建議適當調(diào)整優(yōu)化人類活動用地，禁止耕地的過度開墾，避免因人類的干擾與破壞造成的水土流失、石漠化等問題，在不影響生態(tài)用地穩(wěn)定性的情況下，適當優(yōu)化土地資源配置。

(4) 發(fā)展限制區(qū)。2013年發(fā)展限制區(qū)人類活動總面積為232.84 hm2，占當年壽城自然保護區(qū)人類活動總面積的3.98%，主要為耕地、道路、居民點用地。2017年發(fā)展限制區(qū)人類活動總面積為232.75 hm2，較2013年減少0.08 hm2，可以看出發(fā)展限制區(qū)人類活動占地面積變化幅度較小，總體上呈略微下降的趨勢。發(fā)展限制區(qū)零散分布在自然保護區(qū)內(nèi)，遠離其他發(fā)展區(qū)，人口較為稀少，該區(qū)受人類活動的影響相對較弱，經(jīng)濟發(fā)展方式單一，主要以農(nóng)業(yè)耕作為主。發(fā)展限制區(qū)作為人類活動的限制地帶，適宜作為動植物的棲息地，是生態(tài)源地重點保護的地帶，該區(qū)域建議禁止人類活動的開發(fā)與建設(shè)，對于該區(qū)域的耕地等人類活動用地，應(yīng)做好土地規(guī)劃，向發(fā)展優(yōu)化區(qū)進行土地轉(zhuǎn)移，營造該區(qū)域優(yōu)質(zhì)的生態(tài)環(huán)境。

4 結(jié) 論

(1) 壽城自然保護區(qū)2013—2017年受人類活動的影響景觀的破碎化程度增加，景觀類型空間分布趨于分散，景觀多樣性和景觀均勻性均有所提高。人類活動用地面積在發(fā)展核心區(qū)、過渡區(qū)、優(yōu)化區(qū)均有增加趨勢，僅在生態(tài)功能最強的限制區(qū)內(nèi)有所下降?？傮w來看，人類活動限制區(qū)內(nèi)控制的較好，但仍有不足之處，應(yīng)當繼續(xù)加強監(jiān)管與規(guī)劃，對優(yōu)化區(qū)、過渡區(qū)的人類活動進行適當?shù)膬?yōu)化與調(diào)整，避免人類活動用地對生態(tài)用地與自然保護區(qū)的可持續(xù)發(fā)展帶來干擾與不利影響。

(2) 由于壽城自然保護區(qū)自然環(huán)境面積范圍較小，從而選取的各個阻力因子受到了一定的限制。但本文利用了兩期較高分辨率的遙感影像，解譯的土地利用數(shù)據(jù)精度較為精準，研究對保護區(qū)進行了景觀指數(shù)演變分析，同時借助MCR模型與景觀安全格局理論對自然保護區(qū)的人類活動進行了較為合理的布局評價優(yōu)化，可為其他自然保護區(qū)的人類活動評價提供科學的參考。

(3) 本文指標創(chuàng)建的過程中，阻力系數(shù)與生態(tài)系統(tǒng)服務(wù)價值指標的確定基本是參考前人的研究經(jīng)驗與專家打分法[25-27]。該方法雖然在某種程度上可以得到較為科學、準確的結(jié)果，但主觀性強。因此，在今后的研究中要綜合考慮更加科學的賦值模型與數(shù)學方法來完善指標的創(chuàng)建。