梁野山國(guó)家級(jí)自然保護(hù)區(qū)景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)

李 崢

（福建省林業(yè)調(diào)查規(guī)劃院 自然保護(hù)地監(jiān)測(cè)中心，福建 福州 350003）

自然保護(hù)區(qū)建設(shè)作為自然保護(hù)地體系建設(shè)的重要組成部分，面臨著自然資源保護(hù)和開發(fā)和諧統(tǒng)一的問(wèn)題。人類活動(dòng)在一定范圍和程度上影響著景觀格局和生態(tài)過(guò)程的變化，在不同尺度下，通過(guò)景觀生態(tài)風(fēng)險(xiǎn)的評(píng)價(jià)和空間表征，直觀地表達(dá)出人類活動(dòng)因素所產(chǎn)生的風(fēng)險(xiǎn)影響[1]，張?zhí)餥2]等研究不同尺度下景觀生態(tài)風(fēng)險(xiǎn)關(guān)聯(lián)的顯著性變化，提出多尺度綜合風(fēng)險(xiǎn)防范的城市管理需求；李謝輝等[3]基于不同尺度下研究渭河下游河流沿線區(qū)域的景觀生態(tài)特征，揭示生態(tài)風(fēng)險(xiǎn)的空間分布特征,建立在人類干擾下景觀結(jié)構(gòu)與區(qū)域生態(tài)風(fēng)險(xiǎn)的有機(jī)聯(lián)系,為區(qū)域環(huán)境管理提供數(shù)量化的決策依據(jù)和理論支持；張瑩等[4]研究不同空間尺度下扎龍自然保護(hù)區(qū)景觀格局時(shí)空變化特征和演變規(guī)律，構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)指數(shù)，對(duì)研究區(qū)不同時(shí)期空間尺度變化過(guò)程中景觀生態(tài)風(fēng)險(xiǎn)時(shí)空變化特征進(jìn)行分析，為自然保護(hù)區(qū)的管理提供有效的科學(xué)支撐，為制定自然資源管理和降低生態(tài)風(fēng)險(xiǎn)措施提供有力的決策支持，對(duì)自然保護(hù)區(qū)自然資源保護(hù)、利用和可持續(xù)發(fā)展具有重要的意義[5-8]。

以梁野山國(guó)家級(jí)自然保護(hù)區(qū)為例，通過(guò)分析研究區(qū)不同空間粒度下景觀指數(shù)的變化情況，確定適宜的空間粒度；以景觀指數(shù)為基礎(chǔ)，構(gòu)建研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)評(píng)指數(shù)；通過(guò)對(duì)比不同空間尺度下景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果和空間分布情況，確定適宜的空間尺度，為制定自然保護(hù)區(qū)保護(hù)管理工作奠定堅(jiān)實(shí)的理論基礎(chǔ)。

1 研究區(qū)與研究方法

1.1 研究區(qū)概況

梁野山國(guó)家級(jí)自然保護(hù)區(qū)位于福建武平縣境內(nèi)，山系為武夷山脈最南端，南嶺山脈東頭，屬中亞熱帶、南亞熱帶過(guò)渡區(qū)域，地理坐標(biāo)為：東經(jīng)116°07'～116°19'，北緯25°04'～25°20'，研究區(qū)總面積16 640 hm2。最高海拔1 538.4 m，最低海拔273 m，相對(duì)高差1 265.4 m。

1.2 數(shù)據(jù)來(lái)源

依據(jù)研究區(qū)2018年12月的資源三號(hào)衛(wèi)星（ZY-3）高分影像，按照國(guó)家土地利用現(xiàn)狀分類標(biāo)準(zhǔn)，對(duì)2018年研究區(qū)森林資源二類數(shù)據(jù)進(jìn)行土地利用分類，生成2018年研究區(qū)土地利用現(xiàn)狀本地?cái)?shù)據(jù)庫(kù)；30 m分辨率的數(shù)字高程數(shù)據(jù)，下載自地理空間數(shù)據(jù)云網(wǎng)站。

1.3 研究方法

1.3.1 空間粒度的選取

空間粒度的選擇要以研究區(qū)景觀類型的特征為依據(jù)，空間粒度選擇過(guò)小，不利用對(duì)景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的總體把握；空間粒度選擇過(guò)大，容易丟失相關(guān)景觀類型信息，景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)準(zhǔn)確性降低。因此，需要在不同的空間粒度大小下，分析景觀指數(shù)的變化趨勢(shì)，從而確定適宜的空間粒度，是保證景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)準(zhǔn)確性的重要前提[9-10]。

1.3.2 多尺度的選擇

選取樣方尺度、生態(tài)單元尺度和行政村尺度3種空間尺度。樣方尺度的評(píng)價(jià)單元，利用GIS軟件將研究區(qū)生成正方形的漁網(wǎng)格作為評(píng)價(jià)單元，由于樣方單元?jiǎng)澐值拇笮∨c其景觀生態(tài)風(fēng)險(xiǎn)值存在一定的相關(guān)性，因此，需要對(duì)不同大小的樣方單元進(jìn)行適宜度分析，分析其景觀生態(tài)風(fēng)險(xiǎn)區(qū)的分布情況，進(jìn)而確定適宜的樣方單元大小；利用GIS軟件的水文分析工具，以研究區(qū)內(nèi)子流域的分水嶺的邊界范圍作為生態(tài)單元尺度的評(píng)價(jià)單元；以研究區(qū)內(nèi)涉及的26個(gè)行政村的范圍作為行政村尺度的評(píng)價(jià)單元。依據(jù)不同空間尺度下評(píng)價(jià)單元代表的景觀信息特征，計(jì)算評(píng)價(jià)單元的景觀生態(tài)風(fēng)險(xiǎn)值，利用地統(tǒng)計(jì)方法進(jìn)行空間表征，得出3種空間尺度下的景觀生態(tài)風(fēng)險(xiǎn)區(qū)的分布狀況。

1.3.3 景觀生態(tài)風(fēng)險(xiǎn)指數(shù)

根據(jù)已有研究成果[11-14]，選取景觀破碎度（Ci）、景觀分離度（Ni）、景觀優(yōu)勢(shì)度（Di）3個(gè)景觀指標(biāo)來(lái)構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)指數(shù)，用于反映研究區(qū)各景觀類型受干擾的風(fēng)險(xiǎn)程度，其計(jì)算公式為：

式中:ERIk為第k個(gè)評(píng)價(jià)單元的景觀生態(tài)風(fēng)險(xiǎn)指數(shù)；Aki為第k個(gè)風(fēng)險(xiǎn)小區(qū)內(nèi)景觀類型i的面積；Ak為第k個(gè)風(fēng)險(xiǎn)小區(qū)的總面積；n為評(píng)價(jià)單元內(nèi)景觀類型種數(shù)；Ri為i類景觀的景觀損失度指數(shù)。

Ri的計(jì)算公式為：式中：Si為景觀干擾度指數(shù)；Vi為景觀脆弱度指數(shù)。

Si的計(jì)算公式為：

式中：a、b、c為景觀指數(shù)的權(quán)重，分別為0.5、0.3、0.2。

Vi表示不同景觀類型對(duì)外界干擾的敏感性，值越大，生態(tài)風(fēng)險(xiǎn)越大。采用專家打分歸一化方法獲取[4，11-14]。結(jié)合研究區(qū)特點(diǎn)，將研究區(qū)景觀類型分別賦值：建設(shè)用地1、林地2、草地3、耕地4、園地5、水域6、其他用地7，并進(jìn)行歸一化處理。

1.3.4 地統(tǒng)計(jì)分析

地統(tǒng)計(jì)學(xué)，是指以具有空間分布特點(diǎn)的區(qū)域化變量理論為基礎(chǔ)，研究自然現(xiàn)象的空間變異與空間結(jié)構(gòu)的一門學(xué)科。地統(tǒng)計(jì)是一系列檢測(cè)、模擬和估計(jì)變量在空間上的相關(guān)關(guān)系和格局的統(tǒng)計(jì)方法[15-16]。景觀生態(tài)風(fēng)險(xiǎn)指數(shù)作為一種典型的區(qū)域化變量，依據(jù)評(píng)價(jià)單元及其風(fēng)險(xiǎn)值，利用GIS軟件的克里金插值功能模塊對(duì)景觀生態(tài)風(fēng)險(xiǎn)指數(shù)進(jìn)行空間插值運(yùn)算，得出研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)的空間分布。

1.4 數(shù)據(jù)處理

景觀類型劃分：按照國(guó)家土地利用類型分類標(biāo)準(zhǔn)，將研究區(qū)的景觀類型進(jìn)行劃分為建設(shè)用地、林地、草地、耕地、園地、水域、其他用地等7種景觀類型。

行政村界線：依據(jù)研究區(qū)的本地?cái)?shù)據(jù)庫(kù)中的村代碼，提取行政村界線。

地形地貌分布圖：按照福建省地形地貌分類，依據(jù)30 m分辨率數(shù)字高程數(shù)據(jù)，將研究區(qū)的地貌形態(tài)進(jìn)行劃分，得到地形地貌分布圖。

分水嶺分布圖：依據(jù)30 m分辨率數(shù)字高程數(shù)據(jù)，在GIS軟件中進(jìn)行水文分析，得到子流域的分水嶺分布圖。

2 結(jié)果與分析

2.1 空間粒度的適宜度分析

選擇反應(yīng)景觀水平的8個(gè)景觀指數(shù)：斑塊數(shù)量(NP)、斑塊密度(PD)、最大斑塊數(shù)(LPI)、景觀形狀指數(shù)(LSI)、聚集度(AI)、香農(nóng)多樣性指數(shù)(SHDI)、周長(zhǎng)面積分維數(shù)(PAFRAC)、蔓延度指數(shù)(CONTAG)作為空間粒度適宜度分析的指標(biāo)。在12種空間粒度(5、10、15、20、25、30、35、40、45、50、55、60 m)下計(jì)算8個(gè)景觀指數(shù)，并分析景觀指數(shù)的變化趨勢(shì)，結(jié)果如表1和圖1所示。最大斑塊數(shù)(LPI)在空間粒度10 m上波動(dòng)很大；斑塊數(shù)量(NP)和斑塊密度(PD)在空間粒度10 m后值變得平穩(wěn)變化，隨后在20 m后出現(xiàn)規(guī)律性遞減；香農(nóng)多樣性指數(shù)(SHDI)基本沒(méi)有變化；其他景觀指數(shù)出現(xiàn)線性規(guī)律性變化，沒(méi)有明顯的轉(zhuǎn)折點(diǎn)。因此，選擇10 m作為本次研究適宜的空間粒度，同時(shí)，研究區(qū)屬于森林生態(tài)系統(tǒng)類型的自然保護(hù)區(qū)，人工活動(dòng)少，人工設(shè)施建設(shè)受到一定限制，選擇此空間粒度，可較好的保留研究區(qū)的人為干擾景觀信息。利用GIS軟件，把研究區(qū)的地類本地?cái)?shù)據(jù)庫(kù)按劃分的景觀類型轉(zhuǎn)成柵格數(shù)據(jù)，像元大小為10 m，生成研究區(qū)景觀類型分布圖。

表1 12種空間粒度下的景觀指數(shù)Tab.1 Landscape index under 12 spatial granularity

圖1 12種空間粒度下的景觀指數(shù)變化趨勢(shì)Fig.1 Trend of landscape index changeunder 12 spatial granularity

2.2 多尺度下的景觀風(fēng)險(xiǎn)評(píng)價(jià)

2.2.1 樣方尺度

以景觀類型分布圖為基礎(chǔ)，利用GIS軟件創(chuàng)建4種漁網(wǎng)類型 （0.5 km×0.5 km，1.0 km×1.0 km，1.5 km×1.5 km，2 km×2 km）作為樣方尺度，以生成的漁網(wǎng)格作為樣方單元，分別計(jì)算其景觀生態(tài)風(fēng)險(xiǎn)值，以漁網(wǎng)格的幾何中心點(diǎn)作為地學(xué)統(tǒng)計(jì)分析的插值點(diǎn)。為了便于各樣方尺度下景觀生態(tài)風(fēng)險(xiǎn)程度的對(duì)比，以0.5 km×0.5 km的樣方尺度為基礎(chǔ)，采用自然間斷點(diǎn)法，進(jìn)行景觀生態(tài)風(fēng)險(xiǎn)程度劃分，并確定劃分區(qū)間：低生態(tài)風(fēng)險(xiǎn)區(qū)[0.135 4,0.156 0）、較低生態(tài)風(fēng)險(xiǎn)區(qū)[0.156 0,0.179 0）、中生態(tài)風(fēng)險(xiǎn)區(qū)[0.179 0,0.209 9）、較高生態(tài)風(fēng)險(xiǎn)區(qū)[0.209 9,0.252 0）、高生態(tài)風(fēng)險(xiǎn)區(qū)[0.252 0,0.337 6]。對(duì)其他樣方尺度的景觀生態(tài)風(fēng)險(xiǎn)程度進(jìn)行分類，結(jié)果如圖2所示。

由圖2可得，隨著漁網(wǎng)格面積的增大，其表征的景觀生態(tài)風(fēng)險(xiǎn)的精度隨之降低。在1.0 km×1.0 km的樣方尺度下，高生態(tài)風(fēng)險(xiǎn)區(qū)的信息丟失，北部的中生態(tài)風(fēng)險(xiǎn)區(qū)和較低風(fēng)險(xiǎn)區(qū)信息丟失嚴(yán)重；在1.5 km×1.5 km的樣方尺度下，北部的中生態(tài)風(fēng)險(xiǎn)區(qū)和較低風(fēng)險(xiǎn)區(qū)信息丟失，南部的較低風(fēng)險(xiǎn)區(qū)和中生態(tài)風(fēng)險(xiǎn)區(qū)信息丟失嚴(yán)重；在2.0 km×2.0 km的樣方尺度下，高生態(tài)風(fēng)險(xiǎn)區(qū)和較高生態(tài)風(fēng)險(xiǎn)區(qū)的信息丟失。所以采用0.5 km×0.5 km的漁網(wǎng)格作為樣方尺度，能較好的保證景觀生態(tài)風(fēng)險(xiǎn)區(qū)空間表征的信息量和精度。

在0.5 km×0.5 km的樣方尺度下，研究區(qū)高生態(tài)風(fēng)險(xiǎn)區(qū)面積768.18 hm2，占總面積的4.62%；較高生態(tài)風(fēng)險(xiǎn)區(qū)面積1 430.82 hm2，占總面積的8.60%；中生態(tài)風(fēng)險(xiǎn)區(qū)面積3 392.43 hm2，占總面積的20.39%；較低生態(tài)風(fēng)險(xiǎn)區(qū)面積5 212.73 hm2，占總面積的31.33%；低生態(tài)風(fēng)險(xiǎn)區(qū)面積5 835.60 hm2，占總面積的35.06%。研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)平均值為0.180 6，處于中生態(tài)風(fēng)險(xiǎn)區(qū)間。高生態(tài)風(fēng)險(xiǎn)區(qū)和較高生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在居民聚集區(qū)及周邊耕地集中區(qū)，低生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在連片的森林分布區(qū)，中生態(tài)風(fēng)險(xiǎn)區(qū)和較低生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在以上3個(gè)區(qū)的中間過(guò)渡帶。

圖2 4種樣方尺度下景觀生態(tài)風(fēng)險(xiǎn)區(qū)域分布Fig.2 Regional distribution of landscape ecological risk at four sample scales

2.2.2 生態(tài)單元尺度

研究區(qū)地勢(shì)從北至南中間高，分別向東西兩方向逐漸降低。在GIS軟件中，對(duì)30 m分辨率的數(shù)字高程數(shù)據(jù)進(jìn)行水文分析，得到分級(jí)的河網(wǎng)和河流鏈接數(shù)據(jù)，依據(jù)河網(wǎng)生成子流域節(jié)點(diǎn)，依據(jù)河流鏈接生成分水嶺。以分水嶺為生態(tài)單元計(jì)算景觀生態(tài)風(fēng)險(xiǎn)值，并將生態(tài)單元的值賦值給子流域節(jié)點(diǎn)，并進(jìn)行研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)的空間表征（圖3）。

圖3 生態(tài)單元尺度下景觀生態(tài)風(fēng)險(xiǎn)區(qū)域分布Fig.3 Regional distribution of landscape ecological risk areas at ecological scale

生態(tài)單元尺度下，研究區(qū)包含5種景觀生態(tài)風(fēng)險(xiǎn)類型：高生態(tài)風(fēng)險(xiǎn)區(qū)面積624.73 hm2，占總面積的3.75%；較高生態(tài)風(fēng)險(xiǎn)區(qū)面積1 602.53 hm2，占總面積的9.63%；中生態(tài)風(fēng)險(xiǎn)區(qū)面積4 090.86 hm2，占總面積的24.58%；較低生態(tài)風(fēng)險(xiǎn)區(qū)面積5 330.78 hm2，占總面積的32.04%；低生態(tài)風(fēng)險(xiǎn)區(qū)面積4 990.86 hm2，占總面積的30.00%。全區(qū)景觀生態(tài)風(fēng)險(xiǎn)平均值為0.180 8，處于中生態(tài)風(fēng)險(xiǎn)區(qū)間。高生態(tài)風(fēng)險(xiǎn)區(qū)和較高生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在研究區(qū)南部的居民聚集區(qū)及周邊耕地集中區(qū)，低生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在北部和中部連片的森林分布區(qū)，中生態(tài)風(fēng)險(xiǎn)區(qū)和較低生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在以上3個(gè)區(qū)的中間過(guò)渡帶。

2.2.3 行政村尺度

研究區(qū)涉及武平縣26個(gè)行政村，以研究區(qū)內(nèi)的村界為評(píng)價(jià)單元邊界，對(duì)研究區(qū)進(jìn)行評(píng)價(jià)單元的劃分，得到26個(gè)評(píng)價(jià)單元，并利用GIS軟件生成各評(píng)價(jià)單元的幾何中心點(diǎn)，將各評(píng)價(jià)單元的景觀生態(tài)風(fēng)險(xiǎn)值賦值給幾何中心點(diǎn)，進(jìn)行地統(tǒng)計(jì)空間表征。如圖4所示。

圖4 行政村尺度下景觀生態(tài)風(fēng)險(xiǎn)區(qū)域分布Fig.4 Regional distribution of landscape ecological risk areas at administrative village scale

在行政村尺度下，研究區(qū)包含較低生態(tài)風(fēng)險(xiǎn)區(qū)和低生態(tài)風(fēng)險(xiǎn)區(qū)2種景觀生態(tài)風(fēng)險(xiǎn)類型。較低生態(tài)風(fēng)險(xiǎn)區(qū)面積313.56 hm2，占總面積的1.88%；低生態(tài)風(fēng)險(xiǎn)區(qū)面積16 326.20 hm2，占總面積的98.12%。全區(qū)景觀生態(tài)風(fēng)險(xiǎn)平均值為0.143 7，處于低生態(tài)風(fēng)險(xiǎn)區(qū)間。較低生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在靠近研究區(qū)東南邊界附近。

2.3 多尺度下景觀生態(tài)風(fēng)險(xiǎn)結(jié)果對(duì)比分析

研究區(qū)3種尺度下的景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果如表2所示，可以看出，在不同類型尺度下，景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的結(jié)果存在一定的差異。樣方尺度和生態(tài)單元尺度的風(fēng)險(xiǎn)區(qū)劃分種類和面積基本一致，從綜合風(fēng)險(xiǎn)等級(jí)劃分來(lái)看，樣方尺度和生態(tài)單元尺度下一致，屬于中生態(tài)風(fēng)險(xiǎn)區(qū)，而行政村尺度屬于低生態(tài)風(fēng)險(xiǎn)區(qū)。所以行政村尺度的結(jié)果可靠性低、科學(xué)性差，不適用于研究區(qū)。

表2 3種空間尺度下景觀生態(tài)風(fēng)險(xiǎn)區(qū)面積及其占比Tab.2 Area and proportion of landscape ecological risk area in three spatial scales

對(duì)樣方種尺和生態(tài)單元尺度的景觀生態(tài)風(fēng)險(xiǎn)區(qū)的分布情況進(jìn)行空間疊加，結(jié)果如表3所示?？梢钥闯觯瑯臃匠叨扰c生態(tài)單元尺度在5種景觀生態(tài)風(fēng)險(xiǎn)區(qū)類型種都有重疊區(qū)域。提取景觀生態(tài)風(fēng)險(xiǎn)區(qū)一致的部分，其重疊面積7 944.84 hm2，占研究區(qū)總面積的47.75%。

表3 樣方尺度與生態(tài)單元尺度下景觀生態(tài)風(fēng)險(xiǎn)區(qū)分布對(duì)比情況Tab.3 Comparison of the distribution of landscape ecological risk areas between quadrat scale and eco unit scale hm2

3 討論

通過(guò)3種尺度下結(jié)果的對(duì)比分析，樣方尺度（0.5 km×0.5 km）和生態(tài)單元尺度比較符合研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)，而行政村尺度以村行政邊界作為評(píng)價(jià)單元進(jìn)行評(píng)價(jià)，采樣點(diǎn)密度不夠，其信息丟失嚴(yán)重，不適宜用于研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)。如圖2所示，在樣方尺度下，隨著評(píng)價(jià)單元逐級(jí)增大（采樣點(diǎn)密度降低），其景觀生態(tài)風(fēng)險(xiǎn)信息丟失情況逐漸加劇，評(píng)價(jià)單元的大?。ú蓸狱c(diǎn)的密度）與景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)精度存在相關(guān)性，需要對(duì)評(píng)價(jià)單元的大小進(jìn)行適宜度分析和選擇。所以對(duì)于研究區(qū)，采用生態(tài)單元尺度比樣方尺度和行政村尺度更加適宜：①避免了樣方尺度下人為選取樣方單元大小是否適宜的不確定因素；②避免了行政村尺度下評(píng)價(jià)單元?jiǎng)澐值娜藶橐蛩?，取而代之的是以自然地形地貌為特點(diǎn)劃分的評(píng)價(jià)單元。

結(jié)合研究區(qū)景觀類型分布圖，生態(tài)單元尺度下景觀生態(tài)風(fēng)險(xiǎn)區(qū)的分布與人類活動(dòng)頻繁程度存在較高的相關(guān)性。高生態(tài)風(fēng)險(xiǎn)區(qū)和較高生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在景觀格局穩(wěn)定性低的區(qū)域，其景觀破碎化程度高，耕地、住宅、園地和水域景觀類型分布比較密集，人類活動(dòng)頻繁程度最高；低生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在景觀格局穩(wěn)定性高的區(qū)域，其景觀破碎化程度低，主要以林地景觀類型為主，研究區(qū)內(nèi)森林保護(hù)管理嚴(yán)格，人類活動(dòng)頻繁程度最低；較低生態(tài)風(fēng)險(xiǎn)區(qū)和中生態(tài)風(fēng)險(xiǎn)區(qū)主要分布在其他風(fēng)險(xiǎn)區(qū)的過(guò)渡地帶，其景觀格局穩(wěn)定性一般，景觀破碎化程度較低，主要以林地、耕地景觀類型為主，人類活動(dòng)頻繁程度中等。同時(shí)，人類活動(dòng)頻繁程度與人類居住地聚集程度有關(guān)，居住地的選址又與生產(chǎn)條件（水源、土地、氣候等）是否滿足人類生活需求相關(guān)，生產(chǎn)條件又與研究區(qū)的地形地貌（匯水區(qū)、坡度、坡向、海拔）有直接關(guān)系，本文采用的生態(tài)單元尺度是基于研究區(qū)地形地貌中地表自然界線（河網(wǎng)子流域和分水嶺）進(jìn)行景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)，既能充分反映自然成因?yàn)橹鞯木坝^類型（林地、草地等）特征，也能間接反映出以人為活動(dòng)成因?yàn)橹鞯木坝^類型（住宅、耕地等）特征。綜上所述，生態(tài)單元尺度下的景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)更加符合研究區(qū)景觀類型特點(diǎn)和實(shí)際情況。

4 結(jié)論

以景觀指數(shù)在不同空間粒度下的變化趨勢(shì)為依據(jù)進(jìn)行空間粒度的確定；利用景觀指標(biāo)構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)指數(shù)模型，在3種空間尺度下進(jìn)行景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果的分析與對(duì)比，得到如下結(jié)論:

(1)空間粒度適宜度分析：通過(guò)計(jì)算不同空間粒度下景觀格局指數(shù)，分析其變化情況，結(jié)合研究區(qū)景觀類型的特點(diǎn)，確定10 m為研究區(qū)適宜的空間粒度。

(2)樣方尺度下樣方單元適宜度分析：計(jì)算和分析4種評(píng)價(jià)單元下景觀生態(tài)風(fēng)險(xiǎn)結(jié)果，根據(jù)所表達(dá)信息的完整性和準(zhǔn)確度來(lái)看，確定0.5 km×0.5 km樣方尺度為適宜尺度。

(3)通過(guò)3種尺度下的景觀生態(tài)風(fēng)險(xiǎn)結(jié)果對(duì)比分析，風(fēng)險(xiǎn)類型、面積和空間分布有所不同。行政村尺度，由于采樣密度不夠，信息流失嚴(yán)重，此方法不適合研究區(qū)的景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)；樣方尺度 （0.5 km×0.5 km）與生態(tài)單元尺度空間重疊區(qū)域最大，風(fēng)險(xiǎn)類型較全，風(fēng)險(xiǎn)區(qū)空間分布更加合理，研究區(qū)整體景觀生態(tài)風(fēng)險(xiǎn)等級(jí)為中生態(tài)風(fēng)險(xiǎn)。但是，除0.5 km×0.5 km樣方尺度與生態(tài)單元尺度下風(fēng)險(xiǎn)分布結(jié)果相似度最高。其他樣方尺度信息量丟失嚴(yán)重和空間分布差異較大，所以，采用生態(tài)單元尺度進(jìn)行研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)可避免人為選擇評(píng)價(jià)單元大小所帶來(lái)的影響，同時(shí)，以地表特征（河網(wǎng)子流域和分水嶺）為評(píng)價(jià)單元的生態(tài)單元尺度，即直接反映出了自然因素成因的景觀類型特征，又間接反映出了人類因素成因的景觀類型特征，具有一定的全面性。

通過(guò)空間粒度適宜度分析，確定了研究區(qū)適宜的空間粒度；通過(guò)多尺度下景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果的分析，確定了適宜的空間尺度。研究方法和成果對(duì)同類型的研究區(qū)開展景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)具有指導(dǎo)意義，對(duì)制定保護(hù)管理措施提供堅(jiān)實(shí)的理論依據(jù)。本研究以景觀為風(fēng)險(xiǎn)受體，人類干擾為風(fēng)險(xiǎn)源，在景觀脆弱度指數(shù)的量化上，采用以往學(xué)者的研究成果，可能不夠確切。同時(shí)，本文只考慮了空間尺度，未對(duì)時(shí)間尺度展開研究，這將在后續(xù)的研究中完善，并進(jìn)行綜合的分析與討論。