基于景觀格局的干旱區(qū)綠洲生態(tài)風險評價與管理
——以青銅峽市為例

韓曉佳, 劉小鵬,2, 王亞娟, 溫勝強, 葉均艷

(1.寧夏大學 資源環(huán)境學院, 銀川 750021; 2.寧夏(中阿)旱區(qū)資源評價與環(huán)境調(diào)控重點實驗室,銀川 750021)

基于景觀格局的干旱區(qū)綠洲生態(tài)風險評價與管理
——以青銅峽市為例

韓曉佳1, 劉小鵬1,2, 王亞娟1, 溫勝強1, 葉均艷1

(1.寧夏大學 資源環(huán)境學院, 銀川 750021; 2.寧夏(中阿)旱區(qū)資源評價與環(huán)境調(diào)控重點實驗室,銀川 750021)

干旱區(qū)綠洲抗干擾能力弱，開展其生態(tài)風險評價具有的重要意義。以景觀生態(tài)學理論為指導(dǎo)，選取青銅峽市作為研究對象，采用三期遙感數(shù)據(jù)，運用景觀格局指數(shù)評估生態(tài)風險，并對生態(tài)風險指數(shù)進行空間化，繪制研究區(qū)2006年、2010年、2016年的生態(tài)風險空間分布圖。在此基礎(chǔ)上，結(jié)合空間統(tǒng)計法分析了研究區(qū)生態(tài)風險的聚集模式。結(jié)果表明：(1) 研究期間，青銅峽市景觀格局變化明顯。耕地、水域和未利用土地的景觀破碎度增加，大斑塊破碎化趨勢明顯，區(qū)域的景觀穩(wěn)定性降低。景觀格局指數(shù)的空間分布顯示景觀破碎度、損失度較高區(qū)分布在人口和經(jīng)濟活動聚集的綠洲核心區(qū)。(2) 研究區(qū)的生態(tài)風險時空變化差異較大。其中高風險區(qū)面積減少最為顯著，并隨著風險等級之間的轉(zhuǎn)換，中等風險區(qū)成為研究區(qū)的主要風險類型。生態(tài)風險呈現(xiàn)緩和趨勢。(3) 研究區(qū)的生態(tài)風險存在顯著空間正相關(guān)，“低低”和“高高”聚集是主要的生態(tài)風險聚集模式，低值為高值包圍或高值為低值包圍的空間孤立點表現(xiàn)不顯著，高低風險之間嚴重分化是研究區(qū)亟待解決的問題。

景觀格局; 生態(tài)風險; 干旱區(qū)綠洲; 風險管理; 青銅峽市

在一定區(qū)域內(nèi)，不確定性的事故和災(zāi)害對生態(tài)系統(tǒng)造成的不利影響稱為生態(tài)風險，具有潛在性、不確定性和危害性等特點[1-2]。隨著土地退化、物種滅絕以及環(huán)境污染等生態(tài)問題相繼出現(xiàn)，嚴重影響和制約了社會經(jīng)濟發(fā)展。因此對環(huán)境中潛在的不利效應(yīng)加以評估和識別，并作為風險管理依據(jù)應(yīng)用于風險預(yù)警和防范之中[3]，是降低生態(tài)風險對生態(tài)系統(tǒng)損傷的重要途徑。生態(tài)風險評價始于20世紀80年代。歷經(jīng)30多年的發(fā)展，風險評價的風險源識別由單一的化學污染物擴展到自然災(zāi)害、人類活動等宏觀因子[4-6]，風險受體也從人體健康擴展到對種群、景觀類型以及整個生態(tài)系統(tǒng)的研究[7-9]；評價方法由簡單的定性描述發(fā)展為利用模型法[10]、熵值法[11]并結(jié)合3S技術(shù)[12]完成定量的計算和空間模擬分析；評價對象的尺度，由流域研究擴展至城市等行政區(qū)域[13-14]。因此，區(qū)域生態(tài)風險評價便具有了動態(tài)性、空間異質(zhì)性以及綜合性等特征。景觀格局，綜合反映了自然與人為相互作用下產(chǎn)生的區(qū)域生態(tài)環(huán)境，充分體現(xiàn)了生態(tài)過程作用在不同研究尺度上的最終結(jié)果?；诰坝^格局的生態(tài)風險評價，從景觀要素的多樣性和空間關(guān)系中揭示了區(qū)域景觀的安全格局，直觀地反映出風險壓力因子的空間分布規(guī)律和風險值梯度演變特征[15]，是目前生態(tài)風險研究與評估熱點。

綠洲是干旱區(qū)內(nèi)人類生產(chǎn)生活發(fā)展的主要空間，是自然與人文要素共同作用下的復(fù)合生態(tài)系統(tǒng)[16]，屬干旱區(qū)內(nèi)特有的生態(tài)景觀現(xiàn)象。干旱區(qū)綠洲的生態(tài)安全，直接關(guān)系到干旱區(qū)的生存環(huán)境安全與社會經(jīng)濟發(fā)展。因此，對綠洲生態(tài)系統(tǒng)的穩(wěn)定性研究成為近年來干旱區(qū)的研究熱點，大多集中在綠洲的形成機制與演化規(guī)律[17]、土地覆被變化[18]以及景觀格局[19]研究上。其中，基于綠洲景觀的高度異質(zhì)性和強烈人為干擾下形成的景觀獨特性，從景觀格局視角研究綠洲生態(tài)安全，能充分揭示綠洲景觀特征，是評價綠洲生態(tài)安全的有效方法。

青銅峽市位于寧夏平原中部，屬西北內(nèi)陸干旱區(qū)，境內(nèi)黃河過境水量約400億m3，有得天獨厚的引黃灌溉條件，是典型干旱區(qū)綠洲地帶，脆弱性與環(huán)境敏感性是其生態(tài)環(huán)境典型特點。而近年來，隨著人類活動的擴張，加之實施的寧夏“十二五”中南部生態(tài)移民工程帶來的新生態(tài)壓力，預(yù)防和解決青銅峽市面臨的生態(tài)風險刻不容緩。本研究立足于這一典型區(qū)域，選取2006年、2010年和2016年三期TM遙感影像為數(shù)據(jù)基礎(chǔ)，依據(jù)景觀生態(tài)學原理，結(jié)合研究區(qū)當?shù)貙嶋H情況構(gòu)建景觀生態(tài)風險指標，運用景觀格局分析軟件Fragstats 3.4計算景觀指數(shù)。在ArcGIS平臺上，揭示生態(tài)風險的空間分異，為研究區(qū)的生態(tài)風險管理提供科學依據(jù)，對于保障整個干旱區(qū)生態(tài)系統(tǒng)功能與穩(wěn)定性具有重要意義。

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

1.1 研究區(qū)概況

青銅峽市位于黃河上游，寧夏平原的中部，隸屬吳忠市，地處105°22′—109°09′E，37°41′—38°33′N，南北長60 km，東西寬30 km，市轄8個鎮(zhèn)和2個場。地勢由西南向東北自高而低呈現(xiàn)階梯狀分布，境內(nèi)地貌可分為山地、丘陵、洪積扇、黃河沖積平原以及庫區(qū)5大類地貌類型，屬中溫帶大陸性氣候，年降水量約260.7 mm，四季分明，晝夜溫差較大，水資源較為豐富，是寧夏回族自治區(qū)重要的工業(yè)基地之一?！笆濉逼陂g，寧夏實施了大規(guī)模的生態(tài)移民工程，移民人口約7.65萬戶(32.9萬人)，涉及南部山區(qū)原州、西吉、隆德、涇源、彭陽、同心、海原等縣區(qū)。其中青銅峽市，共承接2.75萬來自南部山區(qū)的移民，建成有4個生態(tài)移民安置區(qū)。

1.2 研究方法

研究數(shù)據(jù)來源于研究區(qū)2006年、2010年和2016年三期分辨率為30 m的LandsatTM遙感影像(L129R34)。以ENVI 5.1為圖像處理平臺，采用監(jiān)督分類、目視解譯和實地驗證相結(jié)合的方法對研究區(qū)的遙感影像進行解譯。結(jié)合研究區(qū)的土地利用現(xiàn)狀，將景觀類型劃分為耕地、林地、草地、水域、建設(shè)用地和未利用土地6種作為風險受體。選取樣本分離度均高于1.87，三期影像解譯Kappa指數(shù)分別是0.84，0.85，0.83，達到解譯精度需求。最后分別得到三期景觀類型圖。借鑒鞏杰等[20]學者的相關(guān)研究，構(gòu)建包括景觀結(jié)構(gòu)指數(shù)、景觀脆弱度指數(shù)和生態(tài)風險綜合指數(shù)在內(nèi)的景觀格局風險指標體系。

1.2.1 景觀結(jié)構(gòu)分析 人類活動通過土地利用影響區(qū)域的生態(tài)系統(tǒng)。在區(qū)域生態(tài)系統(tǒng)受到干擾后，風險受體的變化顯著地表現(xiàn)在景觀結(jié)構(gòu)的空間組成上。而景觀結(jié)構(gòu)指數(shù)用以反映景觀生態(tài)系統(tǒng)遭受到外界壓力干擾和脅迫程度?？紤]青銅峽市的實際情況，本研究選取景觀破碎度(Ci)、景觀分離度(Di)和面積周長分維度(Fi)衡量景觀的空間結(jié)構(gòu)變化，參考相關(guān)研究文獻，將三項指數(shù)通過附加權(quán)重后可得到景觀結(jié)構(gòu)指數(shù)(Si)。

Si=aCi+bDi+cFi

式中：a，b，c為各個指標的權(quán)重，權(quán)重的大小，依據(jù)已有學者的相關(guān)研究[21]，并運用專家打分法，分別賦值為0.5，0.3，0.2。

1.2.2 脆弱性分析 不同景觀類型其本身的脆弱程度不同，在受到風險源干擾后呈現(xiàn)出不同程度上的損害。景觀脆弱度指數(shù)(Vi)表明該類型景觀對外部干擾的抵抗能力。對于易損數(shù)值的確定，常采用專家賦值歸一化處理和層次分析法。本研究在借鑒已有的研究成果和經(jīng)驗[22]，結(jié)合研究區(qū)景觀的實際情況，采用專家賦值法，對6類景觀脆弱性分級打分，最后進行歸一化處理得到各類景觀的脆弱度指數(shù)Vi，對其由高到低依次排序為未利用土地(0.285 7)、水域(0.238 1)、耕地(0.190 5)、草地(0.142 9)、林地(0.095 2)、建設(shè)用地(0.047 6)。而景觀損失度(Ri)指數(shù)，用以描述不同景觀類型在遭受到干擾后所受到生態(tài)損失的多少，是由一種土地類型景觀的景觀結(jié)構(gòu)指數(shù)(Si)和景觀脆弱度指數(shù)(Vi)相乘而得到：

Ri=Si×Vi

1.2.3 綜合生態(tài)風險分析 基于景觀結(jié)構(gòu)的變化分析，僅能從表面體現(xiàn)區(qū)域生態(tài)風險的特征變化，為了更好描述研究區(qū)綜合生態(tài)風險，構(gòu)建起景觀結(jié)構(gòu)與生態(tài)風險的聯(lián)系，需利用每一類型景觀的損失度指數(shù)與該類型景觀的面積比重來定量化區(qū)域綜合生態(tài)風險(ERI)。

式中：Aji是風險小區(qū)內(nèi)第i類景觀的面積；Aj是風險小區(qū)的面積；Ri是第i類型景觀的生態(tài)損失度；n是景觀類型；ERI是第i個風險的生態(tài)風險指數(shù)。

依據(jù)景觀生態(tài)學中對景觀樣本的基本要求，即樣本面積應(yīng)介于平均斑塊面積的2～5倍[23]，將青銅峽市以3km×3km網(wǎng)格劃分單元進行采樣，共243個風險小區(qū)。分別計算各個風險小區(qū)的綜合生態(tài)風險值，利用ArcGIS中的Kriging空間插值法，將量化后的生態(tài)風險指數(shù)轉(zhuǎn)化成為空間分布形式。

1.2.4 空間自相關(guān)分析 探究區(qū)域內(nèi)生態(tài)風險的空間異質(zhì)性規(guī)律一般采用空間統(tǒng)計學方法進行分析[24-25]。其中，運用空間自相關(guān)分析法可以分析風險小區(qū)的風險取值與相鄰風險區(qū)值的相似程度，即空間正相關(guān)為某一風險小區(qū)的取值與其相鄰區(qū)取值有相同的變化趨勢，相反則為空間負相關(guān)?？臻g自相關(guān)的相似度一般采用Moran指數(shù)進行刻畫，一般分為全局Moran指標和局部Moran指標。其中局部Moran指數(shù)可揭示研究區(qū)局部空間變異情況以及與相鄰間風險小區(qū)的空間關(guān)系，其計算公式為：

2 結(jié)果與分析

2.1 青銅峽市景觀格局特征

由表1可知，在2006—2016年青銅峽市的景觀格局有明顯的變化。耕地、水域和未利用土地呈現(xiàn)景觀破碎化趨勢，區(qū)域景觀的穩(wěn)定性有所降低。由于人類活動的擴張，耕地和水域的破碎化較為明顯，破碎度呈上升趨勢，耕地的分維度變小，說明耕地和水域這兩類景觀受人類活動干擾明顯，較大斑塊被分割，數(shù)量隨之增多，景觀的破碎度加大。其中建設(shè)用地的破碎度在2006—2010年增加明顯，是由于此期間處于研究區(qū)城鎮(zhèn)迅速擴張階段，土地利用管理較為松散，城鎮(zhèn)建成區(qū)尚未形成，破碎化較為嚴重。之后的2010—2016年，城鎮(zhèn)擴張速度下降，加強對城鎮(zhèn)土地的集約利用，破碎度隨之明顯下降，景觀類型的分維度也隨之變小。未利用土地的破碎度有所增加，且變化較為明顯，是因為除城鎮(zhèn)擴張開墾了部分未利用土地外，研究區(qū)還承接了寧夏“十二五”期間生態(tài)移民工程搬遷的27 577人，建成4個規(guī)模較大的生態(tài)移民安置區(qū)，且均在未利用土地的荒灘上新建而成，一定程度上加大了未利用地的破碎度。但較耕地、建設(shè)用地等利用率高的景觀類型，破碎化程度和受人為干擾程度還較低。草地、林地的破碎化趨勢和分維度均有所下降。究其原因是研究區(qū)較好地響應(yīng)政府退耕還林還草政策，大面積地種植人工林地和草地，景觀形狀較為規(guī)則，受人類活動影響明顯。景觀的生態(tài)損失度顯示：除耕地和草地有較為明顯的上升趨勢外，其余景觀類型整體上變化不明顯。

從景觀格局指數(shù)的空間分布來看，景觀破碎度、損失度較高區(qū)分布在研究區(qū)人口和經(jīng)濟活動聚集的東北部，且分維度較低，該景觀類型的斑塊形狀趨于規(guī)則，人類活動的干擾表現(xiàn)明顯。分布在研究區(qū)西部和東南部的未利用土地景觀破碎度、損失度較低，景觀的分維度較大，說明未利用土地較好地保存了原有的自然狀態(tài)，受到人為干擾較小。從景觀格局指數(shù)的時間變化看，研究區(qū)東北部和中南部的破碎度和損失度范圍在2006—2016年呈擴大趨勢，分維度范圍也隨之增大，說明景觀受人為干擾影響在不斷加大。這些變化多分布在東北部(綠洲核心區(qū))的邊緣，屬生態(tài)敏感區(qū)，需加強生態(tài)風險防范。

表1 2006-2016年青銅峽市景觀格局指數(shù)

2.2 青銅峽市生態(tài)風險時空分異

計算已劃分的生態(tài)風險小區(qū)風險值，在ArcGIS平臺利用克里金空間插值法將量化后的生態(tài)風險轉(zhuǎn)化至空間上，能夠較為直觀地表現(xiàn)出研究區(qū)內(nèi)生態(tài)風險的空間分布狀況。統(tǒng)計結(jié)果顯示，在2006—2016年各風險小區(qū)的風險值為0.001 3～0.651 8，高低風險值之間差值較大，但3期數(shù)據(jù)風險平均值變化不明顯。在ArcGIS中，按等間隔重分類將生態(tài)風險值由小到大依次劃分5個等級，即低風險區(qū)(0.001 3～0.131 4)、較低風險區(qū)(0.131 4～0.261 5)、中等風險區(qū)(0.261 5～0.391 6)、較高風險區(qū)(0.391 6～0.521 7)和高風險區(qū)(0.5217～0.651 8)。

由表2可知，青銅峽在研究期間的10 年內(nèi)生態(tài)風險有較大的變化。2006年研究區(qū)的生態(tài)風險以較低風險區(qū)和高風險區(qū)為主，生態(tài)風險等級兩極分化嚴重。2010年高風險區(qū)面積開始減少，中等風險區(qū)和較高風險區(qū)面積加大，生態(tài)風險表現(xiàn)出一定的緩和趨勢。2016年高風險區(qū)大面積減少，變化劇烈，較低風險區(qū)和中等風險區(qū)占研究區(qū)面積的73.75%，成為青銅峽市的主要風險等級類型。其中，城鎮(zhèn)面積的迅速擴張和積極響應(yīng)退耕還林還草政策是生態(tài)風險降低的重要原因。

生態(tài)風險的空間分布顯示，高風險和較高風險區(qū)分布在未利用土地面積較大的西部和東南部，這部分是自然景觀保存完好卻易于受到人為干擾的區(qū)域。隨著城鎮(zhèn)擴張，大部分未利用土地被開墾利用，一定程度上增加了高風險區(qū)景觀的空間異質(zhì)性，原有未利用土地的生態(tài)損失度下降，高生態(tài)風險也隨之降低。而綠洲核心區(qū)，盡管人為干擾度較大，但由于林地和草地的面積增加，生態(tài)風險類型以較低和低風險為主。從整體來看，中等生態(tài)風險逐漸成為研究區(qū)主要的生態(tài)風險類型。

2.3 青銅峽市生態(tài)風險聚集模式

計算研究區(qū)2006年、2010年、2016年生態(tài)風險的全局Moran指數(shù)均為正值，說明在研究區(qū)內(nèi)，生態(tài)風險空間分布上存在顯著的空間正相關(guān)，即生態(tài)風險表現(xiàn)出相似值之間的空間聚集趨勢。運用空間自相關(guān)分析軟件OpenGeoDa生成Moran散點圖(圖1)和LISA聚類圖(圖2)。由圖1可知，大部分風險小區(qū)位于第一、第三象限，空間正相關(guān)聯(lián)系顯著，屬于“低低”聚集(低值包圍低值區(qū))和“高高”聚集(高值包圍高值區(qū))風險聚集模式。第二象限的“低高”聚集模式和第三象限的“高低”聚集模式也有零星分布，但數(shù)量上明顯要小于“低低”和“高高”聚集模式。從LISA聚類圖中(圖2)，可以更為清晰地看到在2006—2016年生態(tài)風險聚集模式的空間分布狀況。“低低”風險聚集區(qū)分布在綠洲核心區(qū)和中南部低風險區(qū)上，且面積有所擴展，位置向南偏移。“高高”風險聚集區(qū)隨時間變化有縮小趨勢，位置向西北、西南偏移。說明研究區(qū)的生態(tài)風險有一定的緩和趨勢，但高風險和低風險之間分化嚴重。

表2 2006-2016年青銅峽市生態(tài)風險級別面積統(tǒng)計

圖12006-2016年研究區(qū)生態(tài)風險小區(qū)Moran散點圖

圖22006-2016年研究區(qū)生態(tài)風險小區(qū)LISA聚類圖

3 結(jié)論與討論

3.1 結(jié) 論

(1) 研究期間，青銅峽市的景觀格局變化明顯。耕地、水域和未利用土地三種景觀類型的破碎度增加，大斑塊破碎化趨勢明顯，研究區(qū)的景觀穩(wěn)定性降低。景觀格局指數(shù)的空間分布顯示景觀破碎度、損失度較大的區(qū)域分布在研究區(qū)人口和經(jīng)濟活動聚集的綠洲核心區(qū)。

(2) 青銅峽市生態(tài)風險的時空變化差異較大。其中高風險區(qū)面積減少顯著，并隨著風險等級之間的轉(zhuǎn)換，中等風險區(qū)成為青銅峽市的主要風險類型。其中，積極響應(yīng)退耕還林還草政策是近10 年研究區(qū)生態(tài)風險緩和的主要原因。

(3) 研究區(qū)的生態(tài)風險存在顯著空間正相關(guān)，LISA聚類圖顯示，“低低”和“高高”聚集是主要生態(tài)風險聚集模式，低值為高值包圍或高值為低值包圍的空間孤立點不顯著。高低風險之間嚴重分化是研究區(qū)亟待解決的問題。

3.2 討 論

近些年，生態(tài)風險評價的側(cè)重點由污染因子對生態(tài)系統(tǒng)的破壞向人類活動對生態(tài)系統(tǒng)干擾轉(zhuǎn)變，如開始探究土地利用變化、氣候變化等對生態(tài)系統(tǒng)的脅迫。研究對象尺度從種群擴展至整個區(qū)域景觀研究。這種變化有利于將生態(tài)風險評價與區(qū)域的環(huán)境與發(fā)展之間相結(jié)合，能為區(qū)域發(fā)展監(jiān)測和預(yù)測生態(tài)環(huán)境變化[26]。依據(jù)青銅峽生態(tài)風險評價結(jié)果，結(jié)合不同等級風險區(qū)的特點和實際情況，提出相應(yīng)的風險管理對策，對降低和減少研究區(qū)的生態(tài)風險損失具有重要意義。針對研究區(qū)各級生態(tài)風險區(qū)的空間分布情況并結(jié)合研究區(qū)的實際情況提出相應(yīng)的風險管理對策建議如下：

(1) 應(yīng)對高生態(tài)風險和較高生態(tài)風險的管理建議。青銅峽市的高生態(tài)風險區(qū)和較高生態(tài)風險區(qū)主要分布在未利用土地面積廣泛、地勢較高的西部和東南地區(qū)。這些地區(qū)本身的脆弱性較高，且為水土流失、滑坡等地質(zhì)災(zāi)害高發(fā)區(qū)。雖受到的人類活動和干擾較少，但卻是生態(tài)系統(tǒng)恢復(fù)和重建的困難地區(qū)，一旦受到人類活動的破壞，將會加劇地質(zhì)災(zāi)害的發(fā)生。因此，該地區(qū)的風險管理重點應(yīng)放在加強生態(tài)環(huán)境的建設(shè)上，嚴格控制對該地區(qū)未利用土地的隨意利用，加強地表植被種植，優(yōu)化配置山區(qū)間的土地利用，遵循自然生態(tài)適應(yīng)性規(guī)律，以提高生態(tài)系統(tǒng)的服務(wù)能力。

(2) 應(yīng)對中等生態(tài)風險的管理建議。研究區(qū)的中等風險區(qū)主要分布在綠洲核心區(qū)邊緣與未利用地廣布的交界處。建設(shè)用地、草地、水域是該風險區(qū)的主要景觀類型，植被覆蓋度中等，并隨著城鎮(zhèn)綠洲區(qū)擴張逐漸成為青銅峽市的主要風險等級。該區(qū)的風險管理對策應(yīng)從維持現(xiàn)狀和降低生態(tài)風險入手，加強因地制宜建設(shè)，合理地布局工業(yè)和農(nóng)業(yè)生產(chǎn)活動，增加草地和林地等植被的覆蓋度，積極推進研究區(qū)內(nèi)的產(chǎn)業(yè)向環(huán)保產(chǎn)業(yè)轉(zhuǎn)型。

(3) 應(yīng)對低生態(tài)風險和較低生態(tài)風險的管理建議。低風險區(qū)和較低風險區(qū)分布在研究區(qū)東北部，即綠洲區(qū)核心。雖然人口和人類活動在這里高度聚集，但該區(qū)草地、林地等植被覆蓋度較高，河流、水庫以及沼澤等水域濕地廣布，生物量相應(yīng)較大，景觀類型的異質(zhì)性較強，生態(tài)風險等級分布呈現(xiàn)出以較低風險和低風險為主。針對該區(qū)的風險管理應(yīng)注重對現(xiàn)階段低生態(tài)風險的維持，合理地規(guī)劃和利用土地，有效調(diào)控人類的生產(chǎn)生活活動，提高人們的環(huán)保和可持續(xù)發(fā)展理念，加強公眾參與生態(tài)風險防范意識和途徑，以減少區(qū)域環(huán)境污染，降低不合理的人類活動對生態(tài)系統(tǒng)的破壞，積極地維護并提高土地的生態(tài)系統(tǒng)服務(wù)價值。

[1] 彭建,黨威雄,劉焱序,等.景觀生態(tài)風險評價研究進展與展望[J].地理學報,2015,70(4):664-677.

[2] 鄧飛,于云江,全占軍.區(qū)域生態(tài)風險評價研究進展[J].環(huán)境科學與技術(shù),2011,34(S1):141-147.

[3] 周平,蒙吉軍.區(qū)域生態(tài)風險管理研究進展[J].生態(tài)學報,2009,29(4):2097-2106.

[4] Mohamed M, Mohamed S, Mehdi M, et al. Environmental and ecological risk assessment of heavy metals in sediments of Nador lagoon, Morocco[J]. Ecological Indicators, 2015,48(1):616-626.

[5] 許學工,顏磊,徐麗芬,等.中國自然災(zāi)害生態(tài)風險評價[J].北京大學學報:自然科學版,2011,47(5):901-908.

[6] 石朋浩,于開芹,馮永軍.基于景觀格局的城鄉(xiāng)結(jié)合部生態(tài)風險分析:以泰安市岱岳區(qū)為例[J].應(yīng)用生態(tài)學報,2013,24(3):705-712.

[7] 張思鋒,劉晗夢.生態(tài)風險評價方法述評[J].生態(tài)學報,2010,30(10):2735-2744.

[8] 王亞娟,米文寶,李建華.寧夏中部干旱帶農(nóng)業(yè)景觀格局變化研究[J].水土保持研究,2013,20(4):203-206.

[9] 康鵬,陳衛(wèi)平,王美娥.基于生態(tài)系統(tǒng)服務(wù)的生態(tài)風險評價研究進展[J].生態(tài)學報,2016,36(5):1192-1203.

[10] Petra T B, Douglas P. An ecological risk management and capacity building model[J]. Human Ecology, 2014,42(4):591-603.

[11] 周婷,蒙吉軍.區(qū)域生態(tài)風險評價方法研究進展[J].生態(tài)學雜志,2009,28(4):762-747.

[12] 王文杰,張永福,王慧杰.基于GIS干旱區(qū)綠洲縣域土地利用變化生態(tài)風險分析:以新疆澤普縣為例[J].水土保持研究,2016,23(6):216-220.

[13] 李博,石培基,金淑婷.石羊河流域植被生態(tài)系統(tǒng)生態(tài)風險評價研究[J].水土保持通報,2013,33(1):201-205.

[14] 王美娥,陳衛(wèi)平,彭馳.城市生態(tài)風險評價研究進展[J].應(yīng)用生態(tài)學報,2014,25(3):911-918.

[15] Echeverría C, Newton A, Nahuelhual L, et al. How landscapes change: Integration of spatial patterns and human processes in temperate landscapes of southern Chile[J]. Applied Geography, 2012,32(2):822-831.

[16] 趙文智,莊艷麗.中國干旱區(qū)綠洲穩(wěn)定性研究[J].干旱區(qū)研究,2008,25(2):155-162.

[17] 王永興,張小雷,闞耀平.綠洲地域系統(tǒng)及其演變規(guī)律的初步研究[J].干旱區(qū)地理,1999,22(1):62-68.

[18] 吳文婕,石培基,張學斌,等.干旱區(qū)內(nèi)陸河流域綠洲水土資源生態(tài)風險動態(tài)研究[J].水土保持通報,2012,33(3):249-255.

[19] 宋冬梅,肖篤寧,張志城,等.甘肅民勤綠洲的景觀格局變化及驅(qū)動力分析[J].應(yīng)用生態(tài)學報,2003,14(4):535-539.

[20] 鞏杰,趙彩霞,謝余初,等.基于景觀格局的甘肅白龍江流域生態(tài)風險評價與管理[J].應(yīng)用生態(tài)學報,2014,25(7):2041-2048.

[21] 謝花林.基于景觀結(jié)構(gòu)的土地利用生態(tài)風險空間特征分析:以江西興國縣為例[J].中國環(huán)境科學,2011,31(4):688-695.

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

[23] 鄔建國.景觀生態(tài)學[M].北京:高等教育出版社,2012.

[24] 魏偉,石培基,雷莉,等.基于景觀結(jié)構(gòu)和空間統(tǒng)計方法的綠洲區(qū)生態(tài)風險分析[J].自然資源學報,2014,29(12):2023-2035.

[25] 謝花林.基于景觀結(jié)構(gòu)和空間統(tǒng)計學的區(qū)域生態(tài)風險分析[J].生態(tài)學報,2008,8(10):5020-5026.

[26] 巫麗蕓.區(qū)域景觀生態(tài)風險評價及生態(tài)風險管理研究[D].福州:福建師范大學，2004.

EcologicalRiskAssessmentandItsManagementofOasisinAridAreaBasedonLandscapePattern—ACaseStudyinQingtongxia

HAN Xiaojia1, LIU Xiaopeng1,2, WANG Yajuan1, WEN Shengqiang1, YE Junyan1

(1.SchoolofResourcesandEnvironment,NingxiaUniversity,Yinchuan750021,China; 2.KeyLaboratory(China-Arab)ofResourceEvaluationandEnvironmentalRegulationofAridRegioninNingxia,Yinchuan750021,China)

It is of great significance to carry out ecological risk assessment because of the weak anti-interference ability of oasis in arid area. Based on the theory of landscape ecology, we selected Qingtongxia City as the research sample, and used remote sensing data in three periods, and the landscape pattern index to assess the ecological risk, and achieved the spatialization of index, as well as made the spatial distribution maps of ecological risk in 2006, 2010 and 2016. On the basis above, the accumulation mode of ecological risk in the study area was analyzed by spatial statistics. The results showed that: (1) during the study period, the change of landscape pattern was obvious in Qingtongxia City; the fragmentation of cultivated land, water area and unused land increased, the fragmentation trend of large patches was obvious, and the landscape stability decreased; the spatial distribution of landscape pattern index show that the landscape fragmentation and the loss degree of the higher areas distributed in the oasis core area which was gathered by population and economic activity; (2) the spatial and temporal changes of ecological risk in the study area were different; among them, the reduction of high risk area was the most remarkable, and with the conversion of the risk grade, the medium risk area became the main risk type; ecological risk presented the mitigative trend; (3) there was prominent spatial positive correlation between the ecological risks in teh study area, the ‘low—low’ and ‘high—high’ accumulations were the main ecological risk accumulation mode, the spatial outliers that low values were surrounded by high values or high values were surrounded by low values were not significant, the serious differentiation between high and low risks is the urgent problem to be solved in the study area.

landscape pattern; ecological risk; arid area oasis; risk management; Qingtongxia City

2017-03-05

：2017-03-31

國家自然科學基金資助項目“寧夏生態(tài)移民安置區(qū)生態(tài)風險及其空間分異研究”(41461039)

韓曉佳(1992—),女,山西長治人,碩士,主要研究生態(tài)經(jīng)濟與生態(tài)評估。E-mail:nxdxhxj@163.com

劉小鵬(1973—),男(滿族),寧夏海原人,博士,教授,主要研究生態(tài)經(jīng)濟發(fā)展。E-mail:2641996568@qq.com

X171.1

：A

：1005-3409(2017)05-0285-06