重慶武隆區(qū)景觀類型轉(zhuǎn)移的生態(tài)風(fēng)險(xiǎn)時(shí)空演化研究

黃九松　劉婷　王鐘書　張傳華　鄧煒

摘要［目的］分析武隆區(qū)景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變特征。［方法］以武隆區(qū)2010、2015、2020年土地利用數(shù)據(jù)為基礎(chǔ)，運(yùn)用ArcGIS 10.2和Fragstats 4.3軟件，通過構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)模型測(cè)算和分析景觀生態(tài)風(fēng)險(xiǎn)大小及時(shí)空演變特征。［結(jié)果］2010—2020年武隆區(qū)景觀格局變化明顯，耕地面積先增后減，景觀優(yōu)勢(shì)度降低；林地面積先減后增，破碎度和分離度呈先升高后降低趨勢(shì)；建設(shè)用地面積先減后增，景觀優(yōu)勢(shì)度升高；濕地和裸地景觀分離度和損失度先降低后升高。2010、2015和2020年都以中高生態(tài)風(fēng)險(xiǎn)等級(jí)為主，低風(fēng)險(xiǎn)生態(tài)區(qū)占比較少；空間分布以烏江為軸線，風(fēng)險(xiǎn)等級(jí)距離城區(qū)呈環(huán)狀逐級(jí)遞減狀態(tài)。2010—2015、2015—2020年低等級(jí)向高等級(jí)風(fēng)險(xiǎn)轉(zhuǎn)換的面積分別占總面積的3.62%、4.46%，高等級(jí)向低等級(jí)風(fēng)險(xiǎn)轉(zhuǎn)換的面積占比分別為0.52%、17.13%，即研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)等級(jí)呈下降趨勢(shì)。［結(jié)論］武陵山區(qū)應(yīng)在生態(tài)環(huán)境保護(hù)中尋求社會(huì)經(jīng)濟(jì)發(fā)展，建立和完善生態(tài)環(huán)境保護(hù)責(zé)任制，提升生態(tài)系統(tǒng)服務(wù)功能和抗風(fēng)險(xiǎn)能力。

關(guān)鍵詞景觀格局；景觀生態(tài)風(fēng)險(xiǎn)；時(shí)空演變；武隆區(qū)

中圖分類號(hào)X718.55;X826文獻(xiàn)標(biāo)識(shí)碼A

文章編號(hào)0517-6611（2023）11-0056-06

doi：10.3969/j.issn.0517-6611.2023.11.013開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Study on Temporal and Spatial Evolution of Ecological Risk of Landscape Type Transfer in Wulong District of Chongqing

HUANG Jiu-song LIU Ting WANG Zhong-shu et al（1. Chongqing Planning and Natural Resources Information Center，Chongqing 401147;2. College of Management， Chongqing University of Technology， Chongqing 400054;3.Chongqing Institute of Geology and Mineral Resources， Chongqing 400042）

Abstract［Objective］ To analyze the temporal and spatial evolution characteristics of landscape ecological risk in Wulong District.［Method］Based on the land use data of Wulong District in 2010， 2015， and 2020， using ArcGIS 10.2 and Fragstats 4.3 software， the landscape ecological risk model was constructed to calculate and analyze the magnitude and spatiotemporal evolution characteristics of landscape ecological risk.［Result］From 2010 to 2020， the landscape pattern changed obviously， the area of arable land increased first and then decreased， and the landscape dominance decreased; the area of forest land decreased first and then increased， and the fragmentation and separation increased first and then decreased; the area of construction land decreased first and then increased， and the landscape dominance increased;the landscape separation and loss of wetland and bare land decreased at first and then increased.In 2010， 2015 and 2020， the ecological risk levels were dominated by medium and high ecological risk levels， and the low-risk ecological areas accounted for a small proportion; the spatial distribution took the Wujiang River as the axis， and the risk level from the urban area was in a state of gradual decrease in a ring shape.During 2010-2015 and 2015-2020， the area converted from low-grade to high-grade risk accounted for 3.62% and 4.46% of the total area， and the area converted from high-grade risk to low-grade risk accounted for 0.52% and 17.13% of the total area， that was， the landscape ecological risk level in the study area showed a downward trend. ［Conclusion］Seek social and economic development in ecological environmental protection in Wuling District，establish and improve the ecological environmental protection responsibility system， and enhance the ecosystem service function and anti-risk ability.

Key wordsLandscape pattern;Landscape ecological risk;Spatiotemporal evolution;Wulong District

生態(tài)文明建設(shè)已納入中國(guó)特色社會(huì)主義事業(yè)的總體布局，景觀格局風(fēng)險(xiǎn)預(yù)控是實(shí)現(xiàn)區(qū)域生態(tài)風(fēng)險(xiǎn)綜合防范、景觀生態(tài)管理和生態(tài)保護(hù)的重要手段之一。景觀生態(tài)風(fēng)險(xiǎn)是指景觀環(huán)境與生態(tài)環(huán)境在自然因素或人為活動(dòng)的干擾下相互作用而產(chǎn)生的負(fù)面影響［1］，是基于經(jīng)濟(jì)、社會(huì)和自然生態(tài)層面對(duì)景觀格局演變?cè)斐傻牟焕蠊M(jìn)行區(qū)域性和綜合性的預(yù)判［2-3］。

國(guó)內(nèi)外學(xué)者對(duì)景觀生態(tài)風(fēng)險(xiǎn)的研究成果豐碩，學(xué)者們通過景觀指數(shù)法［4］、風(fēng)險(xiǎn)“源-匯”法［5-6］、空間自相關(guān)分析法［7-8］等方法，從景觀風(fēng)險(xiǎn)評(píng)價(jià)［9-10］、格局優(yōu)化［11-12］、時(shí)空分異［13-14］等視角，對(duì)河湖流域［15-16］、城市地域［17-18］、自然保護(hù)區(qū)［19］、水利樞紐區(qū)［20］等不同尺度進(jìn)行研究?？v觀已有研究，鮮有學(xué)者關(guān)注到山地景觀生態(tài)風(fēng)險(xiǎn)防控。我國(guó)山地丘陵地貌分布廣，且生態(tài)環(huán)境較平原和河湖地區(qū)更為脆弱，山區(qū)受到自然因素和人為活動(dòng)的干擾不斷加劇，伴隨產(chǎn)生的一系列生態(tài)安全問題嚴(yán)重威脅區(qū)域生態(tài)環(huán)境的持續(xù)發(fā)展和生態(tài)安全格局的合理構(gòu)建［21］。因此，以山地為視角研究景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變有重要的現(xiàn)實(shí)意義。

武陵山脈貫穿黔東、湘西、鄂西、渝東南地區(qū)，是烏江、沅江、澧水的分水嶺，是長(zhǎng)江流域的重要生態(tài)安全屏障。武隆區(qū)位于武陵山脈西部，地處烏江下游流域，是重慶的生態(tài)功能區(qū)。武隆區(qū)崇山峻嶺，溝壑縱橫，自然資源豐富，對(duì)當(dāng)?shù)氐淖匀簧鷳B(tài)進(jìn)行保護(hù)式利用有利于當(dāng)?shù)刈匀簧鷳B(tài)旅游的發(fā)展，對(duì)武隆區(qū)的社會(huì)經(jīng)濟(jì)發(fā)展、當(dāng)?shù)孛撠毘晒柟獭?shí)現(xiàn)經(jīng)濟(jì)持續(xù)增長(zhǎng)具有重要意義。因此，該研究選取武隆區(qū)為研究對(duì)象，通過構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)模型，對(duì)武隆區(qū)的景觀格局指數(shù)時(shí)序變化和景觀生態(tài)風(fēng)險(xiǎn)時(shí)空分異情況進(jìn)行研究，以期為武隆區(qū)生態(tài)景觀格局優(yōu)化及風(fēng)險(xiǎn)防范、山區(qū)自然生態(tài)環(huán)境保護(hù)與管理提供理論基礎(chǔ)和依據(jù)。

1資料與方法

1.1研究區(qū)概況武隆區(qū)位于重慶東南部（107°14′～108°05′E、29°02′～29°40′N），處于武陵山與大婁山的結(jié)合部，東鄰彭水，西接南川和涪陵，南靠貴州道真，北抵豐都，位于重慶“一圈兩翼”的渝東南翼，有“渝黔門屏”之稱（圖1）。武隆區(qū)土地面積為288 947.49 hm²，最低海拔160 m，最高海拔2 033 m，全區(qū)崇山峻嶺，溝壑縱橫，自然資源豐富，土壤以紫色土、黃壤土、黃棕壤土和水稻土為主，主要植被為人工林和天然次生林。該區(qū)屬于亞熱帶濕潤(rùn)季風(fēng)氣候，全年平均氣溫17.4 ℃，年降水量1 199.7 mm，年日照時(shí)數(shù)939.0 h。截至2021年底，全年地區(qū)生產(chǎn)總值突破262.14億元，同比增長(zhǎng)7.8%；工業(yè)增加值59.9億元，增長(zhǎng)7.4%；一般公共預(yù)算收入15.8億元，增長(zhǎng)18.1%；固定資產(chǎn)投資86.8億元，增長(zhǎng)3.3%；社會(huì)零售總額136.8億元，增長(zhǎng)21.7%；城鄉(xiāng)居民人均可支配收入分別達(dá)44 194、17 175元，分別增長(zhǎng)9.1%、11.0%。武隆有豐富的旅游資源，包括芙蓉江國(guó)家重點(diǎn)風(fēng)景名勝區(qū)、仙女山國(guó)家森林公園景區(qū)、天生三橋、千里烏江畫廊等旅游勝地，被譽(yù)為中國(guó)西部地質(zhì)之鄉(xiāng)。

1.2數(shù)據(jù)來源與處理該研究所采用的數(shù)據(jù)主要包括：①土地利用類型來源于2010、2015、2020年3個(gè)時(shí)期土地利用遙感監(jiān)測(cè)數(shù)據(jù)，運(yùn)用Erdas imagine 9.0軟件對(duì)各期影像進(jìn)行幾何校正和配準(zhǔn)，采用人工目視解譯、判讀和識(shí)別提取3期土地利用信息，參照GB/T 21010—2017全國(guó)土地利用分類體系將武隆區(qū)景觀劃分為6種基本類型：耕地、林地、草地、濕地、建設(shè)用地、裸地（圖2）。②數(shù)字高程模型（digital elevation model）來源于ASTER GDEM數(shù)據(jù)，運(yùn)用ArcGIS 10.2處理得到武隆區(qū)坡度、高程和地形數(shù)據(jù)。以上數(shù)據(jù)獲取后，通過ArcGIS 10.2建立空間數(shù)據(jù)庫(kù)進(jìn)行整合，并統(tǒng)一各專題圖件的空間投影坐標(biāo)系（Gauss_Kruger，CGCS2000_3_Degree_GK_Zone_35），統(tǒng)一柵格數(shù)據(jù)空間分辨率為30 m×30 m。③景觀格局指數(shù)運(yùn)用 Fragstats 4.3計(jì)算得出。④人口、社會(huì)、經(jīng)濟(jì)數(shù)據(jù)來源于歷年《武隆區(qū)統(tǒng)計(jì)年鑒》以及武隆區(qū)國(guó)民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)。

1.3研究方法

1.3.1風(fēng)險(xiǎn)小區(qū)劃分。參考已有研究［22］，該研究利用ArcGIS 10.2等間距系統(tǒng)采樣法對(duì)武隆區(qū)2010、2015和2020年的土地利用數(shù)據(jù)進(jìn)行5 km×5 km的正方形網(wǎng)格化處理，采集劃分風(fēng)險(xiǎn)小區(qū)155個(gè)（圖1），風(fēng)險(xiǎn)小區(qū)面積為研究區(qū)景觀斑塊平均面積的2～5倍［23］。以風(fēng)險(xiǎn)小區(qū)為評(píng)價(jià)單元，計(jì)算風(fēng)險(xiǎn)小區(qū)內(nèi)各景觀類型的景觀格局指數(shù)，利用半方差分析法進(jìn)行生態(tài)風(fēng)險(xiǎn)等級(jí)空間特征分析。

1.3.2景觀生態(tài)風(fēng)險(xiǎn)模型。

1.3.3半方差分析法。

2結(jié)果與分析

2.1景觀格局指數(shù)時(shí)序變化利用Excel和Fragstats 4.2軟件的統(tǒng)計(jì)分析功能，根據(jù)表1的景觀格局指數(shù)計(jì)算公式計(jì)算得到2010、2015、2020年6種基本類型的景觀格局指數(shù)（表2）。

研究期內(nèi)，耕地和林地的面積有變化，但其在研究區(qū)景觀類型中的主要地位沒有發(fā)生改變。2007年武隆喀斯特地貌成功申遺，武隆區(qū)進(jìn)入快速發(fā)展階段，建設(shè)用地?cái)U(kuò)張速度急劇加快，但得益于2008年開始實(shí)施的重慶地票制度，武隆區(qū)大量廢棄閑置的農(nóng)村建設(shè)用地復(fù)墾為耕地?？傮w來看，2010—2020年建設(shè)用地增加了3 104 hm²，耕地減少了16 299 hm²；其中前5年建設(shè)用地面積略有減少，耕地面積略有增加；而后5年建設(shè)用地面積快速增加，耕地面積快速減少。隨著武隆區(qū)旅游業(yè)的快速發(fā)展和城鎮(zhèn)化進(jìn)程的不斷加快，人為的開發(fā)利用活動(dòng)對(duì)景觀格局的影響不斷加強(qiáng)。

2010—2020年，耕地的斑塊數(shù)量先減后增，但面積先增后減，耕地的破碎度指數(shù)（C_i）和分離度指數(shù)（N_i）呈先減后增趨勢(shì)，耕地被農(nóng)村居民點(diǎn)、城鎮(zhèn)用地等景觀類型用地分割，最終形成小塊零散的分布狀態(tài)，耕地優(yōu)勢(shì)度指數(shù)（D_i）總體呈降低趨勢(shì)，損失度指數(shù)（Ri）總體呈升高趨勢(shì)，由2010年的0.016 7升高至2020年的0.020 0。與耕地相反，林地斑塊數(shù)量先增后減，但面積先減后增，C_i和N_i呈先升高后降低趨勢(shì)，Ri沒有變化。建設(shè)用地斑塊數(shù)量增加，面積呈先減后增趨勢(shì)，C_i和N_i呈先升高后降低趨勢(shì)，Di不斷升高，由2010年的0.313 2升高至2020年的0.342 4，R_i先升高后降低。草地斑塊數(shù)量和面積不斷減少，C_i和D_i降低，N_i和R_i升高。濕地D_i先升高后降低，C_i先降低后升高，N_i和R_i先降低后升高。裸地斑塊數(shù)量和面積均呈先增后減趨勢(shì)，C_i升高，D_i先升高后降低，N_i和R_i先降低后升高。草地和裸地轉(zhuǎn)換為其他景觀類型是造成上述兩類景觀指數(shù)變化的主要原因。

2.2景觀生態(tài)風(fēng)險(xiǎn)時(shí)空分異借鑒高賓等［34］和黃木易等［36］的研究，結(jié)合武隆區(qū)實(shí)際，運(yùn)用手動(dòng)斷點(diǎn)法將各風(fēng)險(xiǎn)小區(qū)的景觀生態(tài)風(fēng)險(xiǎn)指數(shù)劃分為5個(gè)等級(jí)：低風(fēng)險(xiǎn)區(qū)（ERI≤0.006 4）、較低風(fēng)險(xiǎn)區(qū)（0.006 4＜ERI≤0.008 0）、中風(fēng)險(xiǎn)區(qū)（0.008 0＜ERI≤0.009 4）、較高風(fēng)險(xiǎn)區(qū)（0.009 4＜ERI≤0.010 8）、高風(fēng)險(xiǎn)區(qū)（ERI＞0.010 8），統(tǒng)計(jì)不同等級(jí)生態(tài)風(fēng)險(xiǎn)面積比例（圖3），進(jìn)行生態(tài)風(fēng)險(xiǎn)面積動(dòng)態(tài)變化分析（表3），繪制得出2010、2015、2020年武隆區(qū)景觀生態(tài)風(fēng)險(xiǎn)等級(jí)分布圖（圖4）。

從時(shí)間序列分析（圖3和表3），2010、2015和2020年研究區(qū)均以中高風(fēng)險(xiǎn)等級(jí)為主，分別占研究區(qū)面積的69.39%、70.01%、68.63%，且3個(gè)年份較高生態(tài)風(fēng)險(xiǎn)所占比重都是最高的，分別是28.08%、27.61%和30.28%；低生態(tài)風(fēng)險(xiǎn)所占比重都是最低的，分別是11.31%、11.11%和10.98%；較低風(fēng)險(xiǎn)區(qū)和中風(fēng)險(xiǎn)區(qū)所占比重維持在20%左右；高風(fēng)險(xiǎn)區(qū)在2020年的比重明顯降低，為13.90%。2010—2015年，只有高風(fēng)險(xiǎn)區(qū)面積增加，動(dòng)態(tài)度為7.06%，其他4種風(fēng)險(xiǎn)區(qū)面積均減少，動(dòng)態(tài)度為-2.15%～-1.33%；2015—2020年，高風(fēng)險(xiǎn)區(qū)面積大幅減少，動(dòng)態(tài)度為-33.67%，中風(fēng)險(xiǎn)區(qū)面積大幅增加，動(dòng)態(tài)度為14.00%。總的來看，后5年較前5年各個(gè)風(fēng)險(xiǎn)區(qū)面積和動(dòng)態(tài)度變化加劇。

從空間分布分析（圖4），高風(fēng)險(xiǎn)區(qū)集中分布在武隆區(qū)中部和東南部，且以烏江為軸線的南北兩岸，風(fēng)險(xiǎn)等級(jí)距離城區(qū)呈環(huán)狀逐級(jí)遞減狀態(tài)。該風(fēng)險(xiǎn)區(qū)地勢(shì)平坦，經(jīng)濟(jì)發(fā)達(dá)，人口密度大，建設(shè)用地分布廣，人類活動(dòng)干擾程度大，建設(shè)用地的破碎度和分離度顯著高于耕地和林地，因此處于高生態(tài)風(fēng)險(xiǎn)等級(jí)；距烏江越遠(yuǎn)，海拔越高，地形起伏越大，社會(huì)經(jīng)濟(jì)條件越差，人類活動(dòng)干擾程度越低，生態(tài)風(fēng)險(xiǎn)等級(jí)越低。

2.3景觀生態(tài)風(fēng)險(xiǎn)轉(zhuǎn)換分析運(yùn)用ArcGIS 10.2的疊加分析功能，將武隆區(qū)2010、2015、2020年的生態(tài)風(fēng)險(xiǎn)等級(jí)分布圖進(jìn)行疊加，得到2010—2015和2015—2020年2個(gè)時(shí)間段不同等級(jí)生態(tài)風(fēng)險(xiǎn)轉(zhuǎn)換矩陣（表4～5）。

武隆區(qū)正向（高風(fēng)險(xiǎn)—低風(fēng)險(xiǎn)）景觀生態(tài)風(fēng)險(xiǎn)的轉(zhuǎn)換類型包括較低—低風(fēng)險(xiǎn)、中—較低風(fēng)險(xiǎn)、較高—中風(fēng)險(xiǎn)、高—較高風(fēng)險(xiǎn)4種類型；武隆區(qū)逆向（低風(fēng)險(xiǎn)—高風(fēng)險(xiǎn)）景觀生態(tài)風(fēng)險(xiǎn)的轉(zhuǎn)換類型包括低—較低風(fēng)險(xiǎn)、較低—中風(fēng)險(xiǎn)、中—較高風(fēng)險(xiǎn)、較高—高風(fēng)險(xiǎn)4種類型。從表4可以看出，2010—2015年正向轉(zhuǎn)換的面積為1 490.75 hm²，占研究區(qū)總面積的0.52%；逆向轉(zhuǎn)換的面積為10 454.50 hm²，占研究區(qū)總面積的3.62%。正向轉(zhuǎn)換的面積和比例低于逆向轉(zhuǎn)換，說明2010—2015年研究區(qū)整體景觀生態(tài)風(fēng)險(xiǎn)程度升高，生態(tài)環(huán)境有惡化趨勢(shì)。由于2010—2015年武隆區(qū)處于城鎮(zhèn)化快速發(fā)展階段，人類活動(dòng)對(duì)景觀干擾強(qiáng)度加大，建設(shè)用地迅速擴(kuò)張，耕地、林地被無序占用，政府管理力度小，沒有迅速采取措施應(yīng)對(duì)環(huán)境變化，生態(tài)環(huán)境遭到破壞卻無法及時(shí)恢復(fù)，生態(tài)風(fēng)險(xiǎn)等級(jí)升高。

由表5可知，2015—2020年生態(tài)風(fēng)險(xiǎn)等級(jí)的轉(zhuǎn)換類型與2010—2015年的一致，正向轉(zhuǎn)換的面積為49 495.75 hm²，占研究區(qū)總面積的17.13%；逆向轉(zhuǎn)換的面積為12 877.75 hm²，占研究區(qū)總面積的4.46%。正向轉(zhuǎn)換的面積和比例遠(yuǎn)高于逆向轉(zhuǎn)換，說明2015—2020年研究區(qū)整體景觀生態(tài)風(fēng)險(xiǎn)程度大幅下降，生態(tài)環(huán)境明顯改善。一是因?yàn)榻ㄔO(shè)用地集聚利用，景觀破碎度和分離度有所下降；二是政府為了促進(jìn)旅游業(yè)的發(fā)展而加大了生態(tài)環(huán)境保護(hù)力度，促進(jìn)了高等級(jí)生態(tài)風(fēng)險(xiǎn)范圍面積縮小。

由圖5可知，2010—2015和2015—2020年2個(gè)時(shí)段的不同等級(jí)生態(tài)風(fēng)險(xiǎn)轉(zhuǎn)換類型均發(fā)生在相鄰2個(gè)等級(jí)的邊緣地帶，后5年（圖5b）比前5年（圖5a）生態(tài)風(fēng)險(xiǎn)等級(jí)轉(zhuǎn)換范圍更廣。前5年的轉(zhuǎn)換范圍是11 945.25 hm²，后5年的轉(zhuǎn)換范圍是62 373.50 hm²，是前5年的5.22倍。說明隨著時(shí)間的推移，人類活動(dòng)的干擾強(qiáng)度加大。2010—2020年（圖5c）高等級(jí)向低等級(jí)風(fēng)險(xiǎn)轉(zhuǎn)換的面積為50 986.50 hm²，低等級(jí)向高等級(jí)風(fēng)險(xiǎn)轉(zhuǎn)換的面積為23 332.25 hm²，前者是后者的2.19倍，說明生態(tài)風(fēng)險(xiǎn)等級(jí)呈下降趨勢(shì)。

3結(jié)論與討論

3.1結(jié)論該研究以風(fēng)險(xiǎn)小區(qū)為評(píng)價(jià)單元進(jìn)行網(wǎng)格化定量研究，運(yùn)用ArcGIS 10.2和Fragstats 4.3軟件，通過構(gòu)建景觀生態(tài)風(fēng)險(xiǎn)模型測(cè)算和分析武隆區(qū)景觀生態(tài)風(fēng)險(xiǎn)大小及時(shí)空演變特征，主要結(jié)論如下：

（1）景觀格局指數(shù)時(shí)序特征：2010—2020年景觀格局變化明顯，耕地面積先增后減，景觀優(yōu)勢(shì)度降低；林地面積先減后增，破碎度和分離度呈先升高后降低趨勢(shì)；建設(shè)用地面積先減后增，景觀優(yōu)勢(shì)度升高；濕地和裸地景觀分離度和損失度先降低后升高。

（2）景觀生態(tài)風(fēng)險(xiǎn)時(shí)空特征：2010、2015和2020年研究區(qū)都以中高風(fēng)險(xiǎn)等級(jí)為主，分別占研究區(qū)總面積的69.39%、70.01%和68.63%；低風(fēng)險(xiǎn)等級(jí)占比很小，3個(gè)年份都維持在10%左右；空間分布上，大致以烏江為軸線，風(fēng)險(xiǎn)等級(jí)距離城區(qū)呈環(huán)狀向外圍逐級(jí)遞減狀態(tài)。

（3）景觀生態(tài)風(fēng)險(xiǎn)轉(zhuǎn)換特征：正向轉(zhuǎn)換和逆向轉(zhuǎn)換各有4種類型；2010—2015、2015—2020年低等級(jí)向高等級(jí)風(fēng)險(xiǎn)（逆向）轉(zhuǎn)換的面積分別占總面積的3.62%、4.46%，高等級(jí)向低等級(jí)風(fēng)險(xiǎn)（正向）轉(zhuǎn)換的面積占比分別為0.52%、17.13%，即研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)等級(jí)呈下降趨勢(shì)。說明社會(huì)經(jīng)濟(jì)發(fā)展中人類活動(dòng)的干擾以及有關(guān)生態(tài)文明建設(shè)的政策對(duì)于景觀生態(tài)風(fēng)險(xiǎn)轉(zhuǎn)移有重大貢獻(xiàn)。

3.2討論景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)可為武陵山區(qū)生態(tài)景觀格局優(yōu)化及風(fēng)險(xiǎn)防范、自然生態(tài)環(huán)境保護(hù)與管理提供有力的依據(jù)。武陵山區(qū)應(yīng)加強(qiáng)生態(tài)環(huán)境保護(hù)，依托當(dāng)?shù)刎S富的自然資源發(fā)展生態(tài)旅游，打造武隆喀斯特文旅品牌，在現(xiàn)有的生態(tài)保護(hù)旅游區(qū)的基礎(chǔ)上，形成特有的自然環(huán)境保護(hù)體系和區(qū)域連接鏈條，建立和完善環(huán)境保護(hù)責(zé)任制，全面整改生態(tài)環(huán)境保護(hù)督察反饋問題，形成系統(tǒng)的生態(tài)保護(hù)體系，提升生態(tài)系統(tǒng)服務(wù)功能和抗風(fēng)險(xiǎn)能力。

該研究構(gòu)建的景觀生態(tài)風(fēng)險(xiǎn)模型沒有綜合考慮武隆區(qū)當(dāng)?shù)氐纳鐣?huì)、經(jīng)濟(jì)和生態(tài)環(huán)境等其他因素，只簡(jiǎn)單表征生態(tài)問題可能發(fā)生的綜合性概率；該研究缺乏生態(tài)風(fēng)險(xiǎn)轉(zhuǎn)移的驅(qū)動(dòng)力分析，沒有探討不同等級(jí)生態(tài)風(fēng)險(xiǎn)的應(yīng)對(duì)和管理措施。

參考文獻(xiàn)

［1］ 賈艷艷，唐曉嵐，任宇杰.長(zhǎng)江流域安徽段生態(tài)系統(tǒng)服務(wù)價(jià)值與景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變及其關(guān)聯(lián)分析［J］.南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2022，46（3）：31-40.

［2］ LIU S L，CUI B S，DONG S K，et al. Evaluating the influence of road networks on landscape and regional ecological risk： A case study in Lancang River Valley of Southwest China［J］.Ecological engineering，2008，34（2）： 91-99.

［3］ 曹祺文，張曦文，馬洪坤，等.景觀生態(tài)風(fēng)險(xiǎn)研究進(jìn)展及基于生態(tài)系統(tǒng)服務(wù)的評(píng)價(jià)框架：ESRISK［J］.地理學(xué)報(bào)，2018，73（5）：843-855.

［4］ 傅伯杰.地理學(xué)綜合研究的途徑與方法：格局與過程耦合［J］.地理學(xué)報(bào)，2014，69（8）：1052-1059.

［5］ WU J S，ZHU Q L，QIAO N，et al.Ecological risk assessment of coal mine area based on “source-sink” landscape theory-A case study of Pingshuo mining area［J］. Journal of cleaner production，2021，295：1-18.

［6］ 程晰鈺.基于“源-匯”景觀格局理論的南四湖流域非點(diǎn)源污染控制研究［D］.曲阜：曲阜師范大學(xué)，2020.

［7］ 杜春艷，張丁軒，付梅臣，等.基于景觀格局的土地生態(tài)系統(tǒng)完整性空間分析研究［J］.安徽農(nóng)業(yè)科學(xué)，2012，40（27）：13466-13469，13483.

［8］ 王權(quán)，唐芳，李陽兵，等.巖溶地區(qū)景觀格局演變及其生態(tài)安全的時(shí)空分異：以貴州省東北部槽谷為例［J］.生態(tài)學(xué)報(bào)，2021，41（18）：7273-7291.

［9］ 劉珍環(huán)，張國(guó)杰，付鳳杰.基于景觀格局-服務(wù)的景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)：以廣州市為例［J］.生態(tài)學(xué)報(bào)，2020，40（10）：3295-3302.

［10］ 傅微，呂一河，傅伯杰，等.陜北黃土高原典型人類活動(dòng)影響下景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)［J］.生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)，2019，35（3）：290-299.

［11］ 李青圃，張正棟，萬露文，等.基于景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的寧江流域景觀格局優(yōu)化［J］.地理學(xué)報(bào)，2019，74（7）：1420-1437.

［12］ 潘竟虎，劉曉.疏勒河流域景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)與生態(tài)安全格局優(yōu)化構(gòu)建［J］.生態(tài)學(xué)雜志，2016，35（3）：791-799.

［13］ 陳丹，蔣貴國(guó)，張妍，等.基于生態(tài)服務(wù)價(jià)值的沱江流域土地景觀生態(tài)風(fēng)險(xiǎn)時(shí)空分異研究［J］.水土保持通報(bào)，2019，39（3）：223-230.

［14］ 李少玲，謝苗苗，李漢廷，王回茴，許萌，周偉.資源型城市景觀生態(tài)風(fēng)險(xiǎn)的時(shí)空分異：以烏海市為例［J］.地學(xué)前緣，2021，28（4）：100-109.

［15］ XIE H L，WEN J M，CHEN Q R，et al.Evaluating the landscape ecological risk based on GIS： A case-study in the Poyang Lake region of China［J］. Land degradation & development，2021，32（9）：2762-2774.

［16］ ZHANG X M，DU H M，WANG Y，et al.Watershed landscape ecological risk assessment and landscape pattern optimization： Take Fujiang River Basin as an example［J］. Human and ecological risk assessment： An international journal，2021，27（9/10）：2254-2276.

［17］ 王瑩，張淑梅，張煜，等.鄭州市綠地景觀破碎化對(duì)人為活動(dòng)的時(shí)空響應(yīng)［J］.西北林學(xué)院學(xué)報(bào)，2021，36（2）：231-239.

［18］ 王飛，葉長(zhǎng)盛，華吉慶，等.南昌市城鎮(zhèn)空間擴(kuò)展與景觀生態(tài)風(fēng)險(xiǎn)的耦合關(guān)系［J］.生態(tài)學(xué)報(bào)，2019，39（4）：1248-1262.

［19］ GAINES K F， PORTER D E， DYER S A， et al. Using wildlife as receptor species：A landscape approach to ecological risk assessment［J］. Environmental management， 2004，34（4）：528-545.

［20］ 任金銅，莫世江，陳群利，等.貴州夾巖水利樞紐區(qū)域景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究［J］.環(huán)境科學(xué)與技術(shù)，2018，41（4）：182-189.

［21］ 張傳華，張鳳太，王鐘書，等.生態(tài)敏感區(qū)村域尺度土地生態(tài)經(jīng)濟(jì)系統(tǒng)的耦合協(xié)調(diào)發(fā)展研究：以重慶市江津區(qū)為例［J］.長(zhǎng)江流域資源與環(huán)境，2021，30（7）：1547-1558.

［22］ 李程程，南忠仁，王若凡，等.基于景觀結(jié)構(gòu)和3S技術(shù)的干旱區(qū)綠洲生態(tài)風(fēng)險(xiǎn)分析：以高臺(tái)縣為例［J］.干旱區(qū)資源與環(huán)境，2012，26（11）：31-35.

［23］ 徐維藝，闞瑤川，王國(guó)文，等.基于價(jià)值—風(fēng)險(xiǎn)的阜平縣生態(tài)分區(qū)研究［J］.林業(yè)與生態(tài)科學(xué)，2022，37（1）：46-56.

［24］ 許鳳嬌. 基于土地利用變化的江蘇沿海地區(qū)生態(tài)風(fēng)險(xiǎn)格局［D］.曲阜：曲阜師范大學(xué)，2018.

［25］ 肖篤寧，李秀珍.當(dāng)代景觀生態(tài)學(xué)的進(jìn)展和展望［J］.地理科學(xué)，1997，17（4）：356-364.

［26］ 岳啟發(fā)，趙筱青，李思楠，等.“一帶一路”背景下博多河流域景觀格局變化及生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究［J］.世界地理研究，2021，30（4）：839-850.

［27］ 謝花林.基于景觀結(jié)構(gòu)的土地利用生態(tài)風(fēng)險(xiǎn)空間特征分析：以江西興國(guó)縣為例［J］.中國(guó)環(huán)境科學(xué)，2011，31（4）：688-695.

［28］ 賈艷艷，唐曉嵐，劉振威，等.長(zhǎng)江沿岸蕪湖區(qū)段景觀生態(tài)風(fēng)險(xiǎn)時(shí)空演變分析［J］.中南林業(yè)科技大學(xué)學(xué)報(bào)，2019，39（11）：78-87.

［29］ 張文靜，孫小銀，單瑞峰，等.1975—2018年南四湖流域景觀生態(tài)風(fēng)險(xiǎn)時(shí)空變化及其驅(qū)動(dòng)因素研究［J］.生態(tài)科學(xué)，2020，39（3）：172-181.

［30］ 徐麗麗. 鹽城市道路網(wǎng)絡(luò)對(duì)景觀格局的影響和累積效應(yīng)分析［D］.南京：南京信息工程大學(xué)，2015.

［31］ 謝花林.基于景觀結(jié)構(gòu)和空間統(tǒng)計(jì)學(xué)的區(qū)域生態(tài)風(fēng)險(xiǎn)分析［J］.生態(tài)學(xué)報(bào)，2008，28（10）：5020-5026.

［32］ 李謝輝，李景宜.基于GIS的區(qū)域景觀生態(tài)風(fēng)險(xiǎn)分析：以渭河下游河流沿線區(qū)域?yàn)槔跩］.干旱區(qū)研究，2008，25（6）：899-903.

［33］ 許學(xué)工，林輝平，付在毅，等.黃河三角洲濕地區(qū)域生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)［J］.北京大學(xué)學(xué)報(bào)（自然科學(xué)版），2001，37（1）：111-120.

［34］ 高賓，李小玉，李志剛，等.基于景觀格局的錦州灣沿海經(jīng)濟(jì)開發(fā)區(qū)生態(tài)風(fēng)險(xiǎn)分析［J］.生態(tài)學(xué)報(bào)，2011，31（12）：3441-3450.

［35］ 劉春艷，張科，劉吉平.1976—2013年三江平原景觀生態(tài)風(fēng)險(xiǎn)變化及驅(qū)動(dòng)力［J］.生態(tài)學(xué)報(bào)，2018，38（11）：3729-3740.

［36］ 黃木易，何翔.近20年來巢湖流域景觀生態(tài)風(fēng)險(xiǎn)評(píng)估與時(shí)空演化機(jī)制［J］.湖泊科學(xué)，2016，28（4）：785-793.