鄭可君,李 琛,吳映梅,*,高彬嬪,李 嬋,武 燕

1 云南師范大學(xué)地理學(xué)部,昆明 650500 2 云南省社會(huì)科學(xué)院,昆明 650034

生態(tài)安全是經(jīng)濟(jì)社會(huì)持續(xù)健康發(fā)展的重要保障,是人類生存發(fā)展的基本條件,是國(guó)家安全的重要組成部分,在“推動(dòng)人類共享生態(tài)文明成果”的理念下,高度重視并有效防范和化解生態(tài)風(fēng)險(xiǎn),實(shí)現(xiàn)人與自然美美與共,是生態(tài)文明建設(shè)的主要內(nèi)容之一[1]。基于景觀格局的生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)從自然、人為因素影響下景觀格局與生態(tài)過程相互作用出發(fā),考量可能產(chǎn)生的不利后果從而實(shí)現(xiàn)多源風(fēng)險(xiǎn)的綜合表征及其空間可視化,對(duì)區(qū)域生態(tài)發(fā)展和生態(tài)安全格局構(gòu)建提供決策依據(jù)[2—3]。隨著經(jīng)濟(jì)社會(huì)的快速發(fā)展和人類活動(dòng)的不斷加劇,自然資源與生態(tài)系統(tǒng)的壓力日趨增大,生態(tài)環(huán)境保護(hù)面臨巨大的挑戰(zhàn)。邊境地區(qū)的景觀生態(tài)風(fēng)險(xiǎn)表征了跨境地區(qū)生態(tài)系統(tǒng)的穩(wěn)定性[4],對(duì)生態(tài)環(huán)境的全球治理意義重大。云南是中國(guó)重要的邊境省之一,受環(huán)境和經(jīng)濟(jì)基礎(chǔ)薄弱的制約云南邊境縣較為落后,正處在經(jīng)濟(jì)社會(huì)快速發(fā)展的鄉(xiāng)村振興階段,建設(shè)用地、耕地等對(duì)林地、濕地的侵占加劇了生態(tài)用地的減少,不合理的土地利用變化嚴(yán)重制約了生態(tài)系統(tǒng)功能的發(fā)揮。

基于景觀格局的生態(tài)學(xué)研究可分析生態(tài)動(dòng)態(tài)過程和空間格局的相互聯(lián)系,能夠?qū)崿F(xiàn)多元風(fēng)險(xiǎn)的綜合表征和空間可視化,對(duì)區(qū)域生態(tài)安全的維護(hù)和保障極為重要,已成為國(guó)內(nèi)外學(xué)者關(guān)注的熱點(diǎn)問題之一[1,5—6]。景觀作為一種系統(tǒng)除具有整體性外,還具有地域分異的規(guī)律性,在不同尺度上對(duì)自然景觀和人文景觀的結(jié)構(gòu)、功能和動(dòng)態(tài)發(fā)生作用[7]。20世紀(jì)90年代開始就有學(xué)者對(duì)景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)進(jìn)行研究,經(jīng)過約30年的發(fā)展,國(guó)內(nèi)學(xué)者更加關(guān)注景觀生態(tài)風(fēng)險(xiǎn),已從不同尺度探究了景觀生態(tài)風(fēng)險(xiǎn)的地域分異,并注重分析自然條件和人類活動(dòng)共同作用下對(duì)景觀生態(tài)風(fēng)險(xiǎn)地域分異的影響,評(píng)價(jià)對(duì)象以生態(tài)脆弱敏感區(qū)和人類活動(dòng)較為集中區(qū)為主,如：州[8]、市[9]、流域[10]、濕地[11]、礦區(qū)[12]等,評(píng)價(jià)方法有景觀格局指數(shù)法[13]、熵值法[14]、暴露—響應(yīng)法[15]等,影響因素有人口[4]、GDP[1]、海拔[16]、降水[5]等。國(guó)外學(xué)者對(duì)生態(tài)風(fēng)險(xiǎn)的研究多以微觀視角切入,主要圍繞污染物帶來的生態(tài)風(fēng)險(xiǎn),如達(dá)卡地區(qū)農(nóng)業(yè)土壤中的重金屬污染帶來的生態(tài)風(fēng)險(xiǎn)[17]、孟加拉國(guó)魯普薩河沉積物中重金屬的污染水平和分布及其生態(tài)和健康風(fēng)險(xiǎn)[18]等。綜上,盡管生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)已有研究取得了豐富的成果,但仍有以下兩個(gè)方面的問題需要進(jìn)一步探究：一是對(duì)邊境山區(qū)景觀生態(tài)風(fēng)險(xiǎn)的研究不足,邊境山區(qū)是生物多樣性保護(hù)的重要區(qū)域、鄉(xiāng)村振興的重點(diǎn)區(qū)域和對(duì)外開放的重要門戶,是復(fù)雜人文環(huán)境和自然環(huán)境下生態(tài)—經(jīng)濟(jì)—社會(huì)發(fā)展矛盾突出的代表,亟需針對(duì)邊境山區(qū)開展景觀生態(tài)風(fēng)險(xiǎn)的研究。二是在景觀生態(tài)風(fēng)險(xiǎn)的影響因素研究中,缺乏綜合考慮區(qū)域經(jīng)濟(jì)社會(huì)協(xié)調(diào)發(fā)展水平及經(jīng)濟(jì)社會(huì)發(fā)展模式對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響。

云南25個(gè)邊境縣是邊境山區(qū)的一個(gè)典型縮影,是中國(guó)西南生態(tài)安全的大門也是打造區(qū)域命運(yùn)共同體的通道。揭示邊境地區(qū)景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變特征與多種自然、人為干擾因素對(duì)其綜合影響,確保生態(tài)環(huán)境不斷向好發(fā)展顯得尤為重要。本研究基于2000、2010及2020年土地利用數(shù)據(jù),計(jì)算區(qū)域景觀格局指數(shù),構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)模型,在系統(tǒng)分析景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變特征的基礎(chǔ)上,進(jìn)一步探究自然、可達(dá)性和社會(huì)經(jīng)濟(jì)因素對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響,以期分析區(qū)域潛在景觀生態(tài)風(fēng)險(xiǎn)及風(fēng)險(xiǎn)區(qū)域分異為風(fēng)險(xiǎn)防范和評(píng)價(jià)提供借鑒,對(duì)豐富邊境地區(qū)生態(tài)系統(tǒng)穩(wěn)定性的研究具有一定的理論與實(shí)踐意義。

圖1 研究區(qū)位置Fig.1 Study area location該圖基于國(guó)家測(cè)繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載審圖號(hào)為GS(2020)4619 號(hào)的標(biāo)準(zhǔn)地圖制作,底圖無(wú)修改。DEM: 數(shù)字高程模型 Digital elevation model; 1隴川縣；2芒市； 3騰沖市； 4貢山縣；5龍陵縣；6鎮(zhèn)康縣；7瀘水市；8耿馬縣；9滄源縣；10孟連縣；11西盟縣；12瀾滄縣；13勐海縣；14景洪市；15勐臘縣；16江城縣；17綠春縣；18金平縣；19河口縣；20馬關(guān)縣；21麻栗坡縣；22富寧縣；23瑞麗市；24盈江縣；25福貢縣

1 研究區(qū)概況

云南25個(gè)邊境縣(市)[19]位于中國(guó)西南部(21°8′ —28°23′N,97°31′ —106°11′E),與緬甸、老撾、越南3個(gè)國(guó)家毗鄰,邊境線長(zhǎng)約4060km,如圖1所示,其分布在滇西北高山峽谷地區(qū)、滇西邊境山區(qū)、滇南水資源豐富地區(qū)、滇東南巖溶石漠化地區(qū),總面積約90838km2。2020年,云南邊境縣國(guó)內(nèi)生產(chǎn)總值約2484.8億元,戶籍人口約674.6萬(wàn)人。25 個(gè)縣中有 21 個(gè)為民族地區(qū),少數(shù)民族人口近60%,分布有云南特有9個(gè)直接由原始社會(huì)形態(tài)過渡到社會(huì)主義社會(huì)的民族和 8 個(gè)人口較少民族,有16個(gè)與境外同一民族毗鄰而居的少數(shù)民族[20]。云南邊境縣在自然條件、資源分布、人文環(huán)境、社會(huì)經(jīng)濟(jì)發(fā)展等的空間分布上差異較大。作為鄉(xiāng)村振興和保護(hù)生態(tài)系統(tǒng)良性循環(huán)的重要區(qū)域,云南邊境縣的發(fā)展既要全面推進(jìn)沿邊開放城鎮(zhèn)帶建設(shè),提高沿邊地區(qū)發(fā)展能力和開放水平,同時(shí)也肩負(fù)著構(gòu)建生物多樣性保護(hù)網(wǎng)絡(luò)的重要使命。進(jìn)一步協(xié)調(diào)好“發(fā)展”與“保護(hù)”的關(guān)系是云南邊境縣發(fā)展的主要目標(biāo)之一。

2 數(shù)據(jù)來源與研究方法

2.1 數(shù)據(jù)來源及處理

土地利用數(shù)據(jù)來源于 Globeland30 全球地表利用數(shù)據(jù)庫(kù)(http://www.globallandcover.com),包括2000、2010 和 2020 年 3 期數(shù)據(jù),空間分辨率為 30 m×30 m,該數(shù)據(jù)集數(shù)據(jù)較新、精度較高。全區(qū)土地利用類型共 9種,包括耕地、林地、草地、灌木地、濕地、水體、人造地表、裸地、冰川永久積雪。高程數(shù)據(jù)來自地理空間數(shù)據(jù)云(http://www.gscloud.cn)；年降水量和年平均氣溫?cái)?shù)據(jù)來自中國(guó)科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心(http://www.resdc.cn)；道路數(shù)據(jù)來源于全國(guó)地理信息資源目錄服務(wù)系統(tǒng) (https://www.webmap.cn/main.do?method=index)；縣域人口和GDP數(shù)據(jù)來源于《中國(guó)縣域統(tǒng)計(jì)年鑒(縣市卷)》；縣域經(jīng)濟(jì)發(fā)展聚類數(shù)據(jù)來源于“云南邊境縣經(jīng)濟(jì)社會(huì)發(fā)展探析”[20]文獻(xiàn)研究。

2.2 評(píng)價(jià)單元?jiǎng)澐?/h3>

風(fēng)險(xiǎn)評(píng)價(jià)單元?jiǎng)澐质巧鷳B(tài)風(fēng)險(xiǎn)評(píng)價(jià)和生態(tài)風(fēng)險(xiǎn)空間可視化表達(dá)的重要步驟。結(jié)合研究區(qū)實(shí)際情況和斑塊大小,采用10km×10km的格網(wǎng),將云南邊境縣劃分成1212個(gè)評(píng)價(jià)單元,利用 Fragstats 4.2軟件,計(jì)算每個(gè)單元的生態(tài)風(fēng)險(xiǎn)值,將其作為屬性值賦予各評(píng)價(jià)單元的中心點(diǎn)[21]。

2.3 景觀生態(tài)風(fēng)險(xiǎn)指數(shù)構(gòu)建

景觀格局指數(shù)高度濃縮了景觀格局信息,從景觀格局的角度出發(fā),構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)指數(shù)的計(jì)算公式為[22]：

(1)

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

(2)

式中Si為景觀結(jié)構(gòu)指數(shù),Fi為景觀脆弱度指數(shù),景觀脆弱度的大小與其在景觀自然演替過程中所處的階段有關(guān),根據(jù)前人研究經(jīng)驗(yàn)[23—24]結(jié)合研究區(qū)特點(diǎn),本文將9類景觀類型按其脆弱度賦值并進(jìn)行打分,歸一化得到各自的脆弱度指數(shù)：耕地0.11、林地0.04、草地0.09、灌木地0.07、濕地0.16、水體0.13、人造地表0.02、裸地0.18、冰川永久積雪0.20。

Si=aCi+bNi+cDi

(3)

式中a、b、c為相應(yīng)景觀指數(shù)的權(quán)重,a+b+c=1,根據(jù)前人研究結(jié)果[22,25]結(jié)合研究區(qū)實(shí)際情況,a、b、c分別賦值為0.5、0.3、0.2。Ci為景觀破碎度指數(shù),景觀破碎化是由于自然或人為干擾所導(dǎo)致的景觀由單一、均質(zhì)和連續(xù)的整體趨向于復(fù)雜、異質(zhì)和不連續(xù)的斑塊鑲嵌體的過程,景觀破碎化是生物多樣性喪失的重要原因之一,其值越大必然導(dǎo)致生態(tài)風(fēng)險(xiǎn)越大；Ni為景觀分離度指數(shù),表示景觀類型中不同斑塊個(gè)體分布的分離程度,其值越大景觀分布越復(fù)雜,景觀內(nèi)不穩(wěn)定性越低,對(duì)應(yīng)的景觀生態(tài)穩(wěn)定性越低,生態(tài)風(fēng)險(xiǎn)越高；Di為景觀優(yōu)勢(shì)度指數(shù),是衡量斑塊在景觀中重要地位的一種指標(biāo),其大小直接反映了斑塊對(duì)景觀格局形成和變化影響的大小,其值越大,斑塊對(duì)景觀格局形成、變化影響越大對(duì)應(yīng)的生態(tài)風(fēng)險(xiǎn)就越大[25]。

(4)

式中Qi為斑塊i出現(xiàn)的格網(wǎng)數(shù)/總格網(wǎng)數(shù)；Mi為斑塊i的數(shù)目/斑塊的總數(shù)目；Li為斑塊i的面積/樣方總面積；ni為景觀類型i的斑塊數(shù)；A為景觀總面積；Ai為景觀類型i的面積。

2.4 探索性空間數(shù)據(jù)分析

運(yùn)用 ArcGIS結(jié)合空間自相關(guān)分析對(duì)研究區(qū)進(jìn)行空間化表達(dá),通過地統(tǒng)計(jì)分析模塊對(duì)點(diǎn)數(shù)據(jù)進(jìn)行半變異函數(shù)擬合,得到最佳擬合模型[21]并選取指數(shù)克里金插值法分別對(duì)2000、2010、2020年數(shù)據(jù)進(jìn)行空間插值,得到景觀生態(tài)風(fēng)險(xiǎn)空間分布圖,采用自然間斷點(diǎn)法劃分景觀生態(tài)風(fēng)險(xiǎn)等級(jí),進(jìn)一步利用轉(zhuǎn)移矩陣模型分析不同時(shí)期景觀生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)移情況[1]。

通過空間自相關(guān)Moran′sI指數(shù)確認(rèn)相鄰或近鄰區(qū)域風(fēng)險(xiǎn)評(píng)價(jià)單元屬性值是否具有關(guān)聯(lián)性[21]。Moran′sI值在-1—1,大于0表示正相關(guān),小于0表示負(fù)相關(guān),等于0表示不相關(guān)[26]。

2.5 地理探測(cè)器歸因分析

綜合前人研究[1,4—5,16]結(jié)果和研究區(qū)自然環(huán)境、人文環(huán)境、資源空間分布及社會(huì)經(jīng)濟(jì)發(fā)展差距較大的特點(diǎn),利用地理探測(cè)器既可以探測(cè)數(shù)值型數(shù)據(jù),也可以探測(cè)定性數(shù)據(jù)的優(yōu)勢(shì)[27]。選取社會(huì)經(jīng)濟(jì)[28]、自然環(huán)境、可達(dá)性[29]三個(gè)維度10個(gè)影響因子,包括人口地理集中度、經(jīng)濟(jì)地理集中度、人口經(jīng)濟(jì)地理集中度不一致指數(shù)[30]、發(fā)展模式[20]、年降水量、年均氣溫、高程、地形起伏度、距道路距離、距城鎮(zhèn)距離。在 ArcGIS 中將所有影響因子?xùn)鸥窕?、統(tǒng)一投影坐標(biāo)系并運(yùn)用自然斷點(diǎn)法將其分類,利用地理探測(cè)器軟件中的因子探測(cè),測(cè)度景觀生態(tài)風(fēng)險(xiǎn)與各影響因子之間的關(guān)系及影響程度[5],同時(shí)兩個(gè)因子共同作用下對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響。兩兩因子的影響可分為5類,如表1所示。運(yùn)用Origin軟件對(duì)交互探測(cè)結(jié)果作圖。

表1 雙因子交互作用結(jié)果類型[27]

2.6 地理集中度

地理集中度是綜合考慮區(qū)域人口、經(jīng)濟(jì)、面積等因素,是衡量人口、經(jīng)濟(jì)空間分布的有效指標(biāo),一般應(yīng)用地理集中度指數(shù)來表征。綜合考慮云南25個(gè)邊境縣在人口和經(jīng)濟(jì)總量與區(qū)域面積上發(fā)展不平衡的問題,用人口地理集中度和經(jīng)濟(jì)地理集中度衡量云南邊境縣人口與經(jīng)濟(jì)的空間分布及集聚狀況,人口地理集中度和經(jīng)濟(jì)地理集中度計(jì)算公式如下：

(5)

式中,Rpopit為人口地理集中度,RGDPit為經(jīng)濟(jì)地理集中度；popit、POPt分別表示i地區(qū)t時(shí)刻人口數(shù)量和區(qū)域總?cè)丝跀?shù)量；gdpit、GDPt分別表示i地區(qū)t時(shí)刻GDP數(shù)值和區(qū)域的總 GDP 數(shù)值；landit、LANDt分別表示i地區(qū)t時(shí)刻國(guó)土面積和區(qū)域的總國(guó)土面積[31]。

在此基礎(chǔ)上,人口與經(jīng)濟(jì)空間分布的不一致性可以用人口地理集中度Rpopit和經(jīng)濟(jì)地理集中度RGDPit的相對(duì)比例來衡量,則地區(qū)i在t時(shí)間人口與經(jīng)濟(jì)地理集中度的不一致指數(shù)Iit的計(jì)算公式如下[32]：

(6)

2.7 經(jīng)濟(jì)社會(huì)發(fā)展模式

根據(jù)游琰[20]對(duì)云南 25 個(gè)邊境縣經(jīng)濟(jì)社會(huì)發(fā)展的聚類分析結(jié)果,將4種不同經(jīng)濟(jì)社會(huì)發(fā)展模式作為影響研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)的影響因素進(jìn)行探測(cè)。該聚類分析綜合考慮了經(jīng)濟(jì)發(fā)展和社會(huì)事業(yè)中的36個(gè)變量。結(jié)果如下：聚類1,經(jīng)濟(jì)社會(huì)發(fā)展領(lǐng)軍型,包括勐?？h、盈江縣、騰沖市、芒市、景洪市、龍陵縣、耿馬縣；聚類2,經(jīng)濟(jì)社會(huì)發(fā)展均衡型,包括河口縣、瑞麗市、勐臘縣；聚類3,經(jīng)濟(jì)社會(huì)發(fā)展緩慢型,包括江城縣、瀘水市、臨滄縣、馬關(guān)縣、金平縣、麻栗坡縣、隴川縣、孟連縣、綠春縣、富寧縣、滄源縣、鎮(zhèn)康縣；聚類4,經(jīng)濟(jì)社會(huì)發(fā)展落后型,包括福貢縣、西蒙縣、貢山縣。

3 結(jié)果與分析

3.1 景觀類型的時(shí)序變化

研究區(qū)景觀類型以林地和耕地為主,林地面積占研究區(qū)總面積的62.77%以上,2000年、2010年、2020年林地面積分別為58511.63km2、57318.19km2、56985.43km2,面積共減少了1526.2km2；耕地作為研究區(qū)第二大景觀類型,面積占比22.71%以上且持續(xù)增加,20年間共增加2396.78km2；草地的變化表現(xiàn)為面積持續(xù)減少,草地面積占比5.9%以上,面積從2000年的7238.83km2減少至2020年的5397.69km2,減少幅度較大,共減少 1841.14km2；灌木地面積占比3.73%以上,面積呈現(xiàn)先上升后下降的變化趨勢(shì),總體增加201.42km2；裸地面積占比在0.04%-0.35%之間,總體面積翻兩番共增加126.16km2；人造地表面積占比在0.35%—0.84%之間,20年間面積持續(xù)增加,共增加了443.63km2；水體面積20年間面積翻一番,共增加237.19km2；其他占比小的濕地面積呈現(xiàn)先增加后減少的變化趨勢(shì),冰川和永久積雪面積呈現(xiàn)先減少后增加的變化趨勢(shì),總體均減少。

從景觀類型的時(shí)序變化可以看出,20年間經(jīng)濟(jì)欠發(fā)達(dá)的邊境縣得益于國(guó)家一系列興邊富民的政策保障,耕地和人造地表面積增加一定程度上反映了邊境地區(qū)人口增加和生活條件得到改善,同時(shí)伴隨著林地和草地大面積減少,景觀格局發(fā)生了較大改變。國(guó)家對(duì)生態(tài)環(huán)境保護(hù)的重視和規(guī)范化建設(shè),灌木地、水體面積總體增加明顯?？傮w來看,研究區(qū)景觀類型的時(shí)序變化一定程度上反映了政策干預(yù)下社會(huì)經(jīng)濟(jì)發(fā)展和生態(tài)環(huán)境保護(hù)對(duì)景觀格局變化的影響[33]。

3.2 景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變特征

3.2.1景觀生態(tài)風(fēng)險(xiǎn)的空間化

云南邊境縣2000年、2010年、2020年景觀生態(tài)風(fēng)險(xiǎn)的全局Moran′sI指數(shù)均大于0且P值均小于等于0.001,分別為0.559、 0.585、0.531(圖2)。反映出研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)在空間上表現(xiàn)為聚集分布,正相關(guān)性較強(qiáng),相鄰空間單元的景觀生態(tài)風(fēng)險(xiǎn)程度相似。全局Moran′sI指數(shù)先上升后下降,總體下降,表明空間相鄰單元相互影響在減弱,空間趨同性逐漸降低。局部空間自相關(guān)分析結(jié)果顯示(圖3),研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)值高高聚集主要分布在研究區(qū)北部貢山縣、騰沖市、龍陵縣、芒市和東部的馬關(guān)縣、麻栗坡縣、富寧縣；低低聚集主要分布在研究區(qū)北部福貢縣、瀘水市,中部的耿馬縣、滄源縣、綠春縣、景洪市、江城縣、勐臘縣；研究區(qū)還存在較少的低高聚集和高低聚集的現(xiàn)象,零星分布具有較強(qiáng)的空間異質(zhì)性,20年間變化較為穩(wěn)定。

圖2 2000—2020 年云南邊境縣景觀生態(tài)風(fēng)險(xiǎn)Moran′s I散點(diǎn)圖Fig.2 Landscape ecological risk Moran′s I scatter diagram of Ecological Protection Area in Border counties in Yunnan, 2000—2020

圖3 2000—2020 年云南邊境縣景觀生態(tài)風(fēng)險(xiǎn)局部空間自相關(guān)LISA結(jié)果Fig.3 Local spatial autocorrelation LISA results for landscape ecological risk Area in Border counties in Yunnan, 2000—2020

通過普通克里金插值法對(duì)景觀生態(tài)風(fēng)險(xiǎn)值進(jìn)行空間插值,得到2000—2020年研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)空間分布情況,再通過自然斷點(diǎn)法將2000年景觀生態(tài)風(fēng)險(xiǎn)分為5個(gè)等級(jí),低風(fēng)險(xiǎn)(ERI<0.1002)、較低風(fēng)險(xiǎn)(0.1002≤ERI<0.1102)、中等風(fēng)險(xiǎn)(0.1102≤ERI<0.1198)、較高風(fēng)險(xiǎn)(0.1198≤ERI<0.1316)、高生態(tài)風(fēng)險(xiǎn)(ERI≥0.1316),為便于數(shù)據(jù)比較分析,2010年、2020年數(shù)據(jù)均采用2000年的分級(jí)區(qū)間(圖4)。從三期數(shù)據(jù)來看,研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)值空間分布較為相似,總體呈東西高、南北低的規(guī)律。

3.2.2景觀生態(tài)風(fēng)險(xiǎn)等級(jí)分布時(shí)空特征

云南25個(gè)邊境縣2000—2020年景觀生態(tài)風(fēng)險(xiǎn)平均值持續(xù)上升,2000年為0.111、2010年為0.113、2020年為0.115。景觀生態(tài)風(fēng)險(xiǎn)等級(jí)主要以較低生態(tài)風(fēng)險(xiǎn)、中等生態(tài)風(fēng)險(xiǎn)為主(圖4)。20年間低生態(tài)風(fēng)險(xiǎn)區(qū)面積占比16.92%以上,面積呈現(xiàn)出先增加后減少的變化規(guī)律,總體減少了4896.47km2,主要分布在研究區(qū)北部的福貢縣、瀘水市,南部的孟連縣、勐?？h、景洪市、勐臘縣、江城縣、綠春縣,這些地方有較好的自然條件加上政策支持和人為保護(hù),對(duì)該類地區(qū)的生態(tài)系統(tǒng)安全發(fā)展有較好的影響,比如高黎貢山國(guó)家級(jí)自然保護(hù)區(qū)和南部邊境生態(tài)屏障區(qū)均分布于此；較低生態(tài)風(fēng)險(xiǎn)區(qū)面積占比在24.12%以上,面積先減少后增加總體減少2746.32km2,主要分布在研究區(qū)低生態(tài)風(fēng)險(xiǎn)區(qū)周圍,多集中于研究區(qū)中段；中等生態(tài)風(fēng)險(xiǎn)區(qū)占比19.22%以上,面積先減少后增加總體增加1959.05km2,其分布范圍較廣,在研究區(qū)的大部分縣均有分布,主要分布在較低生態(tài)風(fēng)險(xiǎn)區(qū)和較高生態(tài)風(fēng)險(xiǎn)區(qū)周圍；較高生態(tài)風(fēng)險(xiǎn)區(qū)面積占比17.14%以上,呈現(xiàn)先減少后增加的變化規(guī)律,總體增加5745.28km2,主要分布在研究區(qū)中部騰沖縣、盈江縣、芒市、龍陵縣和南部的勐?？h、瀾滄縣及東部的馬關(guān)縣、富寧縣；高生態(tài)風(fēng)險(xiǎn)區(qū)面積占比7.11%以上,其變化較大,2000—2010年面積增加4720.40km2,2010—2020年面積減少4781.94km2,總體面積減少了61.54km2,主要分布在貢山縣、騰沖市、麻栗坡縣、富寧縣,其中麻栗坡縣和富寧縣是石漠化地區(qū),該類地區(qū)自然條件不適宜植被生長(zhǎng),除此之外傳統(tǒng)的農(nóng)耕方式對(duì)土地的利用率較低,使得土地利用類型改變較大,破壞了生態(tài)系統(tǒng)的整體性和穩(wěn)定性,使得以耕地和林地為主的景觀破碎度及分離度增大,最終導(dǎo)致景觀生態(tài)風(fēng)險(xiǎn)增大。

圖4 2000—2020 年云南邊境縣景觀生態(tài)風(fēng)險(xiǎn)等級(jí)變化空間分布Fig.4 Spatial distribution of landscape ecological risk levels in Border counties of Yunnan province from 2000 to 2020

3.2.3景觀生態(tài)風(fēng)險(xiǎn)的等級(jí)轉(zhuǎn)移時(shí)空特征

2000—2020年云南邊境縣景觀生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)移情況如表2所示,除等級(jí)沒有變化的區(qū)域外,低生態(tài)風(fēng)險(xiǎn)等級(jí)以轉(zhuǎn)為較低生態(tài)風(fēng)險(xiǎn)等級(jí)為主,轉(zhuǎn)移面積7104.13km2；較低生態(tài)風(fēng)險(xiǎn)等級(jí)以轉(zhuǎn)為中等生態(tài)風(fēng)險(xiǎn)等級(jí)為主,轉(zhuǎn)移面積10021.21km2；中等生態(tài)風(fēng)險(xiǎn)等級(jí)以轉(zhuǎn)為較高生態(tài)風(fēng)險(xiǎn)為主,轉(zhuǎn)移面積7036.40km2,以上三個(gè)等級(jí)的轉(zhuǎn)移顯示研究區(qū)局部景觀生態(tài)風(fēng)險(xiǎn)有不斷增加趨勢(shì),應(yīng)重點(diǎn)關(guān)注該區(qū)域的景觀結(jié)構(gòu)變化,維護(hù)區(qū)域生態(tài)的穩(wěn)定性,避免因景觀破碎度、分離度、優(yōu)勢(shì)度指數(shù)增大從而導(dǎo)致景觀生態(tài)風(fēng)險(xiǎn)進(jìn)一步增加。除此之外,較高生態(tài)風(fēng)險(xiǎn)等級(jí)以轉(zhuǎn)為中等生態(tài)風(fēng)險(xiǎn)等級(jí)為主,高生態(tài)風(fēng)險(xiǎn)等級(jí)以轉(zhuǎn)為較高生態(tài)風(fēng)險(xiǎn)等級(jí)為主,高生態(tài)風(fēng)險(xiǎn)等級(jí)的面積得到了控制,可以看出在科學(xué)發(fā)展觀和新發(fā)展理念的指導(dǎo)下,區(qū)域發(fā)展更加注重綠色高效,現(xiàn)階段的總體發(fā)展情況對(duì)生態(tài)系統(tǒng)良性循環(huán)較為有利,但同時(shí)應(yīng)加強(qiáng)對(duì)較高生態(tài)風(fēng)險(xiǎn)轉(zhuǎn)為高生態(tài)風(fēng)險(xiǎn)區(qū)和中等生態(tài)風(fēng)險(xiǎn)轉(zhuǎn)為較高生態(tài)風(fēng)險(xiǎn)區(qū)的監(jiān)控。

表2 2000—2020年云南邊境縣景觀生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)移矩陣/km2

3.3 景觀生態(tài)風(fēng)險(xiǎn)時(shí)空分異的影響因素

3.3.1單因子探測(cè)

由表3 可以看出,研究期內(nèi)自然、區(qū)域可達(dá)性、經(jīng)濟(jì)社會(huì)因素對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響均較大。在自然影響因素中年降水量解釋力最大,分別為2000年32.17、2010年24.28、2020年16.25,表明研究區(qū)受降水影響植被空間分異明顯,從而影響景觀結(jié)構(gòu)最終作用于景觀生態(tài)風(fēng)險(xiǎn),然而體現(xiàn)云南邊境縣復(fù)雜地形地貌的地形起伏度對(duì)景觀生態(tài)風(fēng)險(xiǎn)的解釋力較小。在區(qū)域可達(dá)性因素中,距城鎮(zhèn)距離的解釋力較大,20年間解釋力增大了2.32,表明城鎮(zhèn)化的過程中人造地表在快速擴(kuò)張改變了區(qū)域景觀結(jié)構(gòu),對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響逐漸增大。在經(jīng)濟(jì)社會(huì)因素中,人口與經(jīng)濟(jì)地理集中度不一致指數(shù)、人口地理集中度指數(shù)解釋力最大,反映出人口與經(jīng)濟(jì)發(fā)展的不協(xié)調(diào)對(duì)區(qū)域景觀結(jié)構(gòu)的影響較大,人口多、經(jīng)濟(jì)不發(fā)達(dá)導(dǎo)致土地的利用效率降低,資源消耗較大使得景觀的破碎度增大,景觀生態(tài)風(fēng)險(xiǎn)也隨之增大；人口集中度低于經(jīng)濟(jì)集中度則資源和土地利用率較高,經(jīng)濟(jì)發(fā)展向外擴(kuò)散,綠色高效的發(fā)展對(duì)景觀結(jié)構(gòu)的保護(hù)較為有益從而使得景觀生態(tài)風(fēng)險(xiǎn)隨之降低。人口地理集中度對(duì)景觀生態(tài)風(fēng)險(xiǎn)解釋力較大,但在研究期內(nèi)逐漸減小,表明人類活動(dòng)對(duì)生態(tài)系統(tǒng)的擾動(dòng)較為明顯,隨著城鎮(zhèn)化的有序推進(jìn),建設(shè)用地相對(duì)穩(wěn)定,加之生態(tài)文明建設(shè)的不斷推進(jìn),人類活動(dòng)對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響得到有效控制。總體來看,景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變是在多種因素共同影響下形成的復(fù)雜過程,只有充分發(fā)揮好云南邊境縣自然環(huán)境和資源優(yōu)勢(shì),高效利用能源,嚴(yán)守生態(tài)保護(hù)紅線,才能有效防控景觀生態(tài)風(fēng)險(xiǎn)進(jìn)一步惡化。

表3 云南邊境縣景觀生態(tài)風(fēng)險(xiǎn)影響因素地理探測(cè)器單因子探測(cè)結(jié)果

3.3.2交互探測(cè)

交互探測(cè)體現(xiàn)因子之間的共同作用相對(duì)于單因子作用時(shí)對(duì)景觀生態(tài)風(fēng)險(xiǎn)影響的差異。探測(cè)結(jié)果顯示研究期內(nèi)驅(qū)動(dòng)因子間交互探測(cè)結(jié)果均表現(xiàn)出非線性增強(qiáng)和雙因子增強(qiáng)作用,不存在相互獨(dú)立或減弱的情況,說明因子間交互作用對(duì)景觀生態(tài)風(fēng)險(xiǎn)的解釋力相對(duì)于單因子作用均有不同程度的增強(qiáng),同時(shí)也印證了景觀生態(tài)風(fēng)險(xiǎn)變化是一種復(fù)雜的因子交互作用過程。如圖5所示,2000年以人口地理集中度指數(shù)∩人口經(jīng)濟(jì)地理集中度不一致指數(shù)(0.5793)、年降水量∩人口地理集中度指數(shù)(0.5691)、年降水量∩年均氣溫(0.5438)的解釋力最大,反映出云南邊境縣人口與經(jīng)濟(jì)發(fā)展不協(xié)調(diào)對(duì)景觀結(jié)構(gòu)的改變較大,較為粗放的經(jīng)濟(jì)發(fā)展方式對(duì)生態(tài)系統(tǒng)良性循環(huán)造成一定威脅,尤其是人口在地理空間的規(guī)?；奂瘜?dǎo)致土地利用類型的快速改變,使得景觀分離度和優(yōu)勢(shì)度迅速提高,對(duì)景觀生態(tài)風(fēng)險(xiǎn)產(chǎn)生不利影響。降水是影響不同土地覆被空間分異的重要影響因素,尤其是不同生態(tài)用地受其影響劇烈,如林地、草地、灌木地等,對(duì)景觀生態(tài)風(fēng)險(xiǎn)影響十分明顯。人口地理集中度指數(shù)與降水的交互作用對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響更為顯著。同時(shí)受不同氣候類型的影響研究區(qū)氣溫差異較大,如熱帶雨林區(qū)和寒溫性針葉林區(qū)是在不同的氣候條件下形成的不同景觀類型,水量是否適宜植被生長(zhǎng)一定程度上決定了林地的景觀優(yōu)勢(shì)度及破碎度,因此當(dāng)降水和氣溫共同作用下對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響更顯著。與2000年相比,2010年人口地理集中度指數(shù)、人口經(jīng)濟(jì)地理集中度不一致指數(shù)、年均氣溫、年降水量仍然是影響區(qū)域景觀生態(tài)風(fēng)險(xiǎn)的重要交互因子但交互的解釋力均有所下降,所有交互中僅經(jīng)濟(jì)地理集中度指數(shù)∩發(fā)展模式(0.2322)解釋力略有上升,表明發(fā)展模式對(duì)資源和土地的利用效率不同,作用于景觀分離度和破碎度就會(huì)影響景觀生態(tài)風(fēng)險(xiǎn)。與2010年相比,2020年距城鎮(zhèn)距離∩人口與經(jīng)濟(jì)地理集中度不一致指數(shù)(0.3872)、人口與經(jīng)濟(jì)地理集中度不一致指數(shù)∩發(fā)展模式(0.3679)兩組因子交互探測(cè)解釋力有明顯提升,人口經(jīng)濟(jì)發(fā)展的不協(xié)調(diào)與距城鎮(zhèn)距離是人類活動(dòng)對(duì)生態(tài)系統(tǒng)擾動(dòng)的量化因子,二者解釋力的增強(qiáng),表明隨著社會(huì)經(jīng)濟(jì)發(fā)展需求的日益增長(zhǎng),人類活動(dòng)不斷加劇,自然因素影響在不斷減弱,經(jīng)濟(jì)社會(huì)及區(qū)域可達(dá)性的影響在不斷增強(qiáng)。

圖5 景觀生態(tài)風(fēng)險(xiǎn)影響因子交互作用探測(cè)結(jié)果Fig.5 Detection results of interaction between factors affecting landscape ecological risk

4 討論

4.1 景觀生態(tài)風(fēng)險(xiǎn)對(duì)土地利用變化的響應(yīng)

研究期間,云南邊境山區(qū)景觀生態(tài)風(fēng)險(xiǎn)的平均值持續(xù)上升,從景觀生態(tài)風(fēng)險(xiǎn)的演變過程來看相較2000—2010年2010—2020年景觀生態(tài)高風(fēng)險(xiǎn)區(qū)和低風(fēng)險(xiǎn)區(qū)面積減少明顯,較高、中等、較低生態(tài)風(fēng)險(xiǎn)區(qū)面積增加明顯。景觀生態(tài)風(fēng)險(xiǎn)的變化主要反映了景觀格局的改變,不同景觀類型構(gòu)成了景觀格局,從區(qū)域景觀類型的變化來看,林地、草地、灌木地在前十年的減少幅度較大,被大量占為耕地,生態(tài)系統(tǒng)的平衡狀態(tài)被打破,后十年通過天然林保護(hù)等生態(tài)環(huán)境保護(hù)工程的實(shí)施,林地面積減少的速度放緩,灌木地面積增加,高生態(tài)風(fēng)險(xiǎn)減少,生態(tài)系統(tǒng)逐漸趨于穩(wěn)定。但由于天然林的減少破壞了天然生態(tài)系統(tǒng),難以被修復(fù),導(dǎo)致低生態(tài)風(fēng)險(xiǎn)區(qū)面積減少且很難增加,較高生態(tài)風(fēng)險(xiǎn)區(qū)面積增加較快。除此之外,由于云南邊境縣經(jīng)濟(jì)發(fā)展較為緩慢,在后十年的脫貧攻堅(jiān)中,易地搬遷,產(chǎn)業(yè)布局等人類活動(dòng)的加劇對(duì)區(qū)域景觀類型的改變較大,資源的開發(fā)利用對(duì)景觀的結(jié)構(gòu)也造成了較大影響,使得較高生態(tài)風(fēng)險(xiǎn)區(qū)面積增加明顯。

4.2 區(qū)域協(xié)調(diào)水平及發(fā)展模式對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響

景觀生態(tài)風(fēng)險(xiǎn)的演變是受復(fù)雜的地理環(huán)境和人類活動(dòng)綜合作用的動(dòng)態(tài)變化過程。因此選取能夠反映研究區(qū)自然環(huán)境和人類活動(dòng)特征的因子,進(jìn)一步分析景觀生態(tài)風(fēng)險(xiǎn)動(dòng)態(tài)變化的影響因素尤為重要。綜合前人研究[20,31—32,34],在探討影響因素時(shí),經(jīng)濟(jì)因素作為主要的指標(biāo)常被考慮在內(nèi),但現(xiàn)有研究較少結(jié)合研究區(qū)特點(diǎn)考慮社會(huì)因素,而社會(huì)因素和經(jīng)濟(jì)因素是人類活動(dòng)的兩個(gè)重要因素,缺一不可。為了更全面地反映研究區(qū)人類活動(dòng)的社會(huì)、經(jīng)濟(jì)因素,結(jié)合邊境山區(qū)經(jīng)濟(jì)發(fā)展較緩慢,社會(huì)發(fā)展較滯后且空間分異較大的現(xiàn)狀,引入人口地理集中度、經(jīng)濟(jì)地理集中度、人口經(jīng)濟(jì)地理集中度不一致指數(shù)和發(fā)展模式,在一定程度上反映了區(qū)域發(fā)展協(xié)調(diào)水平和經(jīng)濟(jì)社會(huì)不同的發(fā)展模式,相較于僅以公里網(wǎng)格或行政區(qū)域的人口數(shù)和GDP等指標(biāo)表征區(qū)域發(fā)展水平來看,更加注重區(qū)域內(nèi)部經(jīng)濟(jì)與社會(huì)發(fā)展的協(xié)調(diào)關(guān)系。從研究結(jié)果來看,相較區(qū)域發(fā)展模式,協(xié)調(diào)水平對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響更大,其中人口地理集中度和人口經(jīng)濟(jì)地理集中度不一致指數(shù)的影響尤為顯著,反映了人口在空間上的集中分布和落后的經(jīng)濟(jì)增長(zhǎng)方式對(duì)區(qū)域景觀格局的改變有較大影響。伴隨著社會(huì)經(jīng)濟(jì)發(fā)展和技術(shù)進(jìn)步,人類對(duì)土地利用的效率在不斷提高,對(duì)生態(tài)環(huán)境保護(hù)的意識(shí)在增強(qiáng),人口地理集中度和發(fā)展模式對(duì)區(qū)域景觀生態(tài)風(fēng)險(xiǎn)的影響在減弱。然而區(qū)域發(fā)展不協(xié)調(diào)的問題更加突出,人口經(jīng)濟(jì)地理集中度不一致指數(shù)的影響增大,反映出協(xié)調(diào)發(fā)展對(duì)生態(tài)系統(tǒng)保護(hù)的重要性。因此,在研究過程中,不能只注重區(qū)域經(jīng)濟(jì)發(fā)展水平,更應(yīng)關(guān)注經(jīng)濟(jì)社會(huì)的發(fā)展模式,并且綜合考慮區(qū)域協(xié)調(diào)發(fā)展水平,從而找到景觀生態(tài)風(fēng)險(xiǎn)空間分異的根源,最終實(shí)現(xiàn)以發(fā)展模式為抓手,具有可操作性和目標(biāo)明確的區(qū)域調(diào)控方案。

4.3 邊境山區(qū)景觀生態(tài)風(fēng)險(xiǎn)調(diào)控對(duì)策

邊境山區(qū)由于特殊的地理區(qū)位和復(fù)雜的地形地貌及人文社會(huì)環(huán)境對(duì)生態(tài)系統(tǒng)保護(hù)有著更為迫切的需求,是國(guó)家生態(tài)安全的一道重要防線,也是展示國(guó)家生態(tài)形象的重要窗口,其景觀生態(tài)風(fēng)險(xiǎn)提高直接威脅到國(guó)家邊境地區(qū)的生態(tài)安全,因此對(duì)邊境山區(qū)景觀生態(tài)風(fēng)險(xiǎn)的防控意義重大。為防止景觀生態(tài)風(fēng)險(xiǎn)進(jìn)一步惡化,應(yīng)打破行政區(qū)劃找到研究區(qū)不同經(jīng)濟(jì)社會(huì)增長(zhǎng)模式適合的發(fā)展路徑,從綠色產(chǎn)業(yè)布局和生態(tài)環(huán)境保護(hù)兩方面入手,構(gòu)筑邊境一體化生態(tài)安全網(wǎng)絡(luò)。綜合考慮不同等級(jí)景觀生態(tài)風(fēng)險(xiǎn)區(qū)域的特點(diǎn),出臺(tái)有針對(duì)性的生態(tài)環(huán)境保護(hù)措施,加強(qiáng)對(duì)當(dāng)?shù)鼐用竦慕逃c培訓(xùn),避免為種植經(jīng)濟(jì)林破壞天然林的行為發(fā)生。對(duì)高生態(tài)風(fēng)險(xiǎn)區(qū)加強(qiáng)生態(tài)保護(hù)、修復(fù)和監(jiān)測(cè),制定生態(tài)紅線合理分配“三生空間”,修復(fù)破碎景觀,降低其景觀生態(tài)風(fēng)險(xiǎn)等級(jí)；對(duì)較高生態(tài)風(fēng)險(xiǎn)區(qū)加強(qiáng)調(diào)控,防止其景觀生態(tài)風(fēng)險(xiǎn)等級(jí)增加,制定以生態(tài)為導(dǎo)向的政府績(jī)效考核指標(biāo)；對(duì)中等生態(tài)風(fēng)險(xiǎn)區(qū),提高土地利用效率,綜合考量城鎮(zhèn)和生態(tài)用地范圍,提前規(guī)劃構(gòu)建區(qū)域生態(tài)安全網(wǎng)絡(luò),提高區(qū)域生態(tài)系統(tǒng)韌性；對(duì)較低生態(tài)風(fēng)險(xiǎn)區(qū)和低生態(tài)風(fēng)險(xiǎn)區(qū)繼續(xù)開展好生態(tài)環(huán)境保護(hù)和監(jiān)測(cè)工作,總結(jié)并推廣低生態(tài)風(fēng)險(xiǎn)區(qū)的發(fā)展模式。

5 結(jié)論

本研究以云南25個(gè)邊境縣為例,利用2000、2010、2020年三期土地利用數(shù)據(jù),通過構(gòu)建基于景觀格局指數(shù)的景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)模型,分析邊境山區(qū)景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變特征,采用地理探測(cè)器對(duì)研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)的解釋力進(jìn)行因子探測(cè)和因子交互探測(cè),結(jié)論如下：(1)2000—2020年研究區(qū)景觀類型主要以林地和耕地為主,20年間景觀類型變化較為明顯,林地、草地面積減少最多,耕地面積增加最多,灌木地、裸地、人造地表、水體面積增加,濕地冰川和永久積雪面積減少。林地、草地與濕地是維系生態(tài)系統(tǒng)循環(huán)的重要景觀類型應(yīng)在發(fā)展中重點(diǎn)關(guān)注其變化。(2)研究期間,云南邊境縣景觀生態(tài)風(fēng)險(xiǎn)平均值持續(xù)上升。研究區(qū)整體以較低生態(tài)風(fēng)險(xiǎn)、中等生態(tài)風(fēng)險(xiǎn)區(qū)為主,低生態(tài)風(fēng)險(xiǎn)、較低生態(tài)風(fēng)險(xiǎn)、高生態(tài)風(fēng)險(xiǎn)面積減少,中等生態(tài)風(fēng)險(xiǎn)、較高生態(tài)風(fēng)險(xiǎn)面積增加。景觀生態(tài)風(fēng)險(xiǎn)指數(shù)在空間上表現(xiàn)出顯著的正相關(guān)性,大多以高高、低低集聚為主,分布較為穩(wěn)定。(3)2000—2020年云南邊境縣景觀生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)移以低生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)為較低生態(tài)風(fēng)險(xiǎn)等級(jí)、較低生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)為中等生態(tài)風(fēng)險(xiǎn)等級(jí)、中等生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)為較高生態(tài)風(fēng)險(xiǎn)等級(jí)、較高生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)為中等生態(tài)風(fēng)險(xiǎn)等級(jí)、高生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)為較高生態(tài)風(fēng)險(xiǎn)等級(jí)為主。(4)景觀生態(tài)風(fēng)險(xiǎn)受自然環(huán)境和人類活動(dòng)共同作用產(chǎn)生空間分異,從因子探測(cè)結(jié)果來看,年降水量、年均氣溫、距城鎮(zhèn)距離、人口地理集中度指數(shù)、人口與經(jīng)濟(jì)地理集中度不一致指數(shù)是影響區(qū)域景觀生態(tài)風(fēng)險(xiǎn)的主要因素；從因子交互探測(cè)結(jié)果來看,兩兩交互均呈現(xiàn)非線性增強(qiáng)和雙因子增強(qiáng),研究期內(nèi)年降水量與年均氣溫交互、人口地理集中度指數(shù)與人口經(jīng)濟(jì)地理集中度不一致指數(shù)交互、年降水量與人口地理集中度指數(shù)交互對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響最大。邊境地區(qū)的景觀生態(tài)風(fēng)險(xiǎn)表征了跨境地區(qū)生態(tài)系統(tǒng)的穩(wěn)定性,對(duì)生態(tài)環(huán)境的全球治理意義重大,探究人類活動(dòng)對(duì)景觀生態(tài)風(fēng)險(xiǎn)的影響能夠進(jìn)一步豐富“社會(huì)-經(jīng)濟(jì)-自然”復(fù)合生態(tài)系統(tǒng)的研究,為區(qū)域協(xié)調(diào)發(fā)展理論的實(shí)踐提供參考。