張雅洲, 謝小平

曲阜師范大學(xué)地理與旅游學(xué)院，日照 276826

基于RS和GIS的南四湖生態(tài)風(fēng)險(xiǎn)評價(jià)

張雅洲, 謝小平*

曲阜師范大學(xué)地理與旅游學(xué)院，日照 276826

以南四湖為研究對象，基于遙感和地理信息系統(tǒng)技術(shù)選擇干旱、洪澇和水體污染為生態(tài)風(fēng)險(xiǎn)源，用景觀格局指數(shù)評價(jià)生態(tài)系統(tǒng)的地位和受體的易損性，并依據(jù)生態(tài)風(fēng)險(xiǎn)值的度量的原理，通過計(jì)算綜合風(fēng)險(xiǎn)概率及生態(tài)綜合損失度得到南四湖的生態(tài)風(fēng)險(xiǎn)值。研究結(jié)果顯示南四湖區(qū)域生態(tài)風(fēng)險(xiǎn)存在上下級湖分異與湖濱帶狀分異兩大空間規(guī)律，體現(xiàn)了南四湖區(qū)域生態(tài)風(fēng)險(xiǎn)的空間異質(zhì)性。

遙感影像; GIS分析; 區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià); 南四湖

區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)是將生態(tài)風(fēng)險(xiǎn)評價(jià)應(yīng)用到區(qū)域景觀上去的方法[1]，一般包括3個(gè)關(guān)鍵環(huán)節(jié)，即區(qū)域內(nèi)各生態(tài)系統(tǒng)的地位和作用評價(jià)、脆弱性評價(jià)以及可能遭受的風(fēng)險(xiǎn)概率評價(jià)[2- 3]。

經(jīng)過近20年的研究發(fā)展，國內(nèi)外學(xué)者基于風(fēng)險(xiǎn)評價(jià)的基本原理和各自的研究，建立了區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)的模型和方法步驟，或者做出了相應(yīng)完善，如Wagner[4]指出景觀生態(tài)學(xué)研究方法的優(yōu)勢在于更加注重區(qū)域內(nèi)部的異質(zhì)性和區(qū)域與區(qū)域之間的關(guān)聯(lián)性，能更好地表達(dá)空間異質(zhì)性；付在毅[5]，許學(xué)工[6]，盧宏瑋[7]分別建立了適用于遼河三角洲、黃河三角洲和洞庭湖的生態(tài)風(fēng)險(xiǎn)評價(jià)模型，并用生態(tài)指數(shù)量化了風(fēng)險(xiǎn)源，用層次分析法賦各個(gè)風(fēng)險(xiǎn)源的權(quán)重值；李謝輝[8]，李景宜[9]，陳鵬[10]在各自的研究中以GIS為基礎(chǔ)，利用插值技術(shù)描述了風(fēng)險(xiǎn)源和風(fēng)險(xiǎn)受體，同樣通過層次分析法得到了區(qū)域內(nèi)各個(gè)風(fēng)險(xiǎn)源的權(quán)重。已有的評價(jià)模型大多用層次分析法對風(fēng)險(xiǎn)源進(jìn)行加權(quán)處理，得到的針對整個(gè)研究區(qū)的某類風(fēng)險(xiǎn)的權(quán)重值，而風(fēng)險(xiǎn)源在不同區(qū)域的強(qiáng)度是有差異的，風(fēng)險(xiǎn)源權(quán)重的處理可以更加強(qiáng)調(diào)空間性，并在本研究中進(jìn)行了嘗試。

南四湖是華北最大的淡水湖泊，生態(tài)意義顯著，但是近20年來，湖區(qū)及周邊，濕地面積趨減，蓄水能力下降，生物生存環(huán)境脆弱，生態(tài)安全面臨嚴(yán)峻挑戰(zhàn)[11]。為滿足南水北調(diào)工程的水質(zhì)要求，南四湖應(yīng)制定高效的生態(tài)環(huán)境管理策略。筆者認(rèn)為生態(tài)管理更應(yīng)注重“前瞻性”，而區(qū)域生態(tài)風(fēng)險(xiǎn)的研究，恰好符合這個(gè)特點(diǎn)。

本研究采用目前國內(nèi)外較認(rèn)可的生態(tài)風(fēng)險(xiǎn)評價(jià)過程，即受體分析、風(fēng)險(xiǎn)源分析、暴露與危害分析和綜合風(fēng)險(xiǎn)評價(jià)[4- 10,12- 13]，充分利用GIS技術(shù)，對南四湖完成了生態(tài)風(fēng)險(xiǎn)評價(jià)，研究結(jié)果能為風(fēng)險(xiǎn)管理者和防災(zāi)減災(zāi)部門提供風(fēng)險(xiǎn)管理和環(huán)境決策依據(jù)。

1 研究區(qū)概況

圖1 研究區(qū)范圍示意圖Fig.1 Location of study area

南四湖位于山東省西南部，北緯34°27′—35°20′，東經(jīng)116°34′—117°21′之間，是南陽湖、獨(dú)山湖、昭陽湖和微山湖4個(gè)相互連貫湖泊的總稱。南四湖位于地質(zhì)構(gòu)造凹陷區(qū)，由古泗水河道被黃河泥沙淤塞而成，屬河跡洼地湖，地處暖溫帶、半濕潤季風(fēng)氣候區(qū)[11]，土壤絕大部分為水稻土和砂姜黑土，水稻土面積約2800 hm2，主要分布在湖西濱湖洼地上，砂姜黑土主要分布在湖東濱湖洼地附近，面積約1300 hm2[14]。

本文區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)的范圍為南四湖湖區(qū)，即南四湖東大堤到西大堤之間，總面積約1197.47 km2。研究中采用了表征相對量的景觀格局指數(shù)，格局指數(shù)對空間尺度的選擇十分敏感[15- 17]，不同的空間尺度會對應(yīng)不同的計(jì)算結(jié)果。在本研究中因?yàn)楫?dāng)景觀格局指數(shù)的計(jì)算范圍僅為湖區(qū)時(shí)，其結(jié)果并不能很好的表達(dá)該區(qū)域景觀斑塊的特點(diǎn)，故計(jì)算了除湖區(qū)外還有湖區(qū)周邊的濟(jì)寧市任城區(qū)、微山縣和魚臺縣的景觀格局指數(shù)，涉及區(qū)域總面積約3350 km2，但在實(shí)際評價(jià)時(shí)僅選用計(jì)算結(jié)果的東大堤到西大堤之間的部分以突出湖區(qū)范圍內(nèi)的景觀特點(diǎn)(圖1藍(lán)色部分)。

2 研究方法

風(fēng)險(xiǎn)(risk)往往被定義為一個(gè)不理想事件發(fā)生的概率及其導(dǎo)致的后果[18]，進(jìn)行區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)，還需強(qiáng)調(diào)空間性[18- 19]。本研究對南四湖湖區(qū)進(jìn)行生態(tài)風(fēng)險(xiǎn)評價(jià)，為體現(xiàn)區(qū)域差異，首先劃分風(fēng)險(xiǎn)小區(qū)，每一個(gè)風(fēng)險(xiǎn)小區(qū)均為面積1 km2的方格，對應(yīng)著不同的風(fēng)險(xiǎn)源和風(fēng)險(xiǎn)受體，針對風(fēng)險(xiǎn)源，評估各類風(fēng)險(xiǎn)源的概率和強(qiáng)度，針對風(fēng)險(xiǎn)受體，引入景觀格局指數(shù)評估其地位和易損性，最終用綜合生態(tài)風(fēng)險(xiǎn)值表征區(qū)域生態(tài)風(fēng)險(xiǎn)的相對大小，每一個(gè)風(fēng)險(xiǎn)小區(qū)對應(yīng)一個(gè)綜合生態(tài)風(fēng)險(xiǎn)值，用以下模型計(jì)算：

Rk=PkDk

式中，Rk是第k個(gè)風(fēng)險(xiǎn)小區(qū)的綜合生態(tài)風(fēng)險(xiǎn)值；Pk是風(fēng)險(xiǎn)小區(qū)的綜合風(fēng)險(xiǎn)概率；Dk是風(fēng)險(xiǎn)小區(qū)的生態(tài)綜合損失度。其中：

式中，Pkj是第k個(gè)風(fēng)險(xiǎn)小區(qū)內(nèi)j類風(fēng)險(xiǎn)源不同級別生態(tài)風(fēng)險(xiǎn)源的概率；?kj是第k個(gè)小區(qū)內(nèi)j類風(fēng)險(xiǎn)源的權(quán)重；j是風(fēng)險(xiǎn)源；n為風(fēng)險(xiǎn)源數(shù)；N為景觀類型的數(shù)量；Ai為小區(qū)內(nèi)第i類景觀組分的面積；A為小區(qū)景觀總面積；Ri為第i種景觀所對應(yīng)的生態(tài)損失度。

在評價(jià)過程中，結(jié)合遙感數(shù)據(jù)、實(shí)測數(shù)據(jù)和歷史統(tǒng)計(jì)數(shù)據(jù)，以ARCGIS9.3為平臺，采用反距離權(quán)重和普通克里格法進(jìn)行空間插值并進(jìn)行對比，發(fā)現(xiàn)反距離權(quán)重法插值易受數(shù)據(jù)點(diǎn)集群的影響，而克里格法插值由于不僅考慮距離，而且通過變異函數(shù)和結(jié)構(gòu)分析，考慮已知樣本點(diǎn)的空間分布與未知樣本點(diǎn)的空間方位關(guān)系，其插值結(jié)果能較好反應(yīng)真實(shí)情況，故在本研究中確定采用克里格法進(jìn)行空間插值去模擬風(fēng)險(xiǎn)源的強(qiáng)度與空間分布以及風(fēng)險(xiǎn)受體的易損性等，并充分運(yùn)用了字段計(jì)算、柵格計(jì)算等其他的GIS功能。除此之外，在數(shù)據(jù)處理時(shí)還用到了數(shù)據(jù)量綱和級別統(tǒng)一的標(biāo)準(zhǔn)化法和統(tǒng)計(jì)風(fēng)險(xiǎn)頻率的概率法等。

3 受體分析

受體就是風(fēng)險(xiǎn)的承受者，是風(fēng)險(xiǎn)評價(jià)中的生態(tài)系統(tǒng)里可能受到來自風(fēng)險(xiǎn)源干擾的不利作用的組成部分[20]。在區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)中，受體則往往是多個(gè)生態(tài)系統(tǒng)類型，不同的生態(tài)系統(tǒng)類型在區(qū)域整體的生態(tài)功能方面所發(fā)揮的作用是不相同的[13]。本研究選定水域生態(tài)系統(tǒng)和濕地生態(tài)系統(tǒng)為主要的生態(tài)受體。

在具有不確定性的風(fēng)險(xiǎn)源作用下，風(fēng)險(xiǎn)受體可能受到的損害，以及由此而發(fā)生的區(qū)域生態(tài)系統(tǒng)結(jié)構(gòu)和功能的損傷就是生態(tài)終點(diǎn)[13]。南四湖湖區(qū)，生態(tài)終點(diǎn)可能是珍稀物種的滅絕、生物種群數(shù)量的減少、濕地的退化和消失、植被的退化或演替中斷、生態(tài)系統(tǒng)功能的損傷和人類利用價(jià)值的喪失等。

4 風(fēng)險(xiǎn)源分析

4.1 風(fēng)險(xiǎn)源的識別

風(fēng)險(xiǎn)源就是可能受到的脅迫。區(qū)域生態(tài)風(fēng)險(xiǎn)源往往包括自然因素和人為因素[9,21]?？紤]發(fā)生的概率、強(qiáng)度和范圍，本研究選擇干旱和洪澇為主要的自然生態(tài)風(fēng)險(xiǎn)源；人為生態(tài)風(fēng)險(xiǎn)源指可能導(dǎo)致危害或嚴(yán)重干擾生態(tài)系統(tǒng)的人為活動[8]，在本研究中，主要選取污廢水排放和水體化學(xué)污染為主要人為生態(tài)風(fēng)險(xiǎn)源，統(tǒng)一稱水體污染。

表1 濟(jì)寧市水旱災(zāi)害發(fā)生概率表

Table 1 The probability of ecological risk sources in the Nansi Lake

位置Location洪澇Flooddisaster級別Grade概率Probability干旱Droughtdisaster級別Grade概率Probability上級湖中等0.4680中等0.3404Highlevellake嚴(yán)重0.0425嚴(yán)重0.1063下級湖中等0.3523中等0.2127Lowlevellake嚴(yán)重0.1063嚴(yán)重0.0180

來源于濟(jì)寧市1945至2005年間后營水文站統(tǒng)計(jì)資料和《濟(jì)寧市地方志》，本文中級別以50年一遇及以上為嚴(yán)重，以下為中等

4.2 風(fēng)險(xiǎn)源描述

風(fēng)險(xiǎn)源表現(xiàn)為風(fēng)險(xiǎn)發(fā)生的概率和強(qiáng)度兩大要素，在區(qū)域尺度上要體現(xiàn)空間的差異性[6]。南四湖生態(tài)風(fēng)險(xiǎn)源的發(fā)生概率和其強(qiáng)度的空間分布情況分別如表1和圖2所示。由于水體污染風(fēng)險(xiǎn)的發(fā)生概率難以獲得，考慮到水體污染給生態(tài)系統(tǒng)帶來的巨大影響，賦其概率為0.5。

4.2.1 干旱

南四湖幾乎每一年都有不同程度的旱情發(fā)生，特別是20世紀(jì)80年代至今，發(fā)生過4次非常嚴(yán)重的旱災(zāi)。最近兩次分別是：1997年，濟(jì)寧市40多條河道斷流，南四湖低于死水位，農(nóng)作物受災(zāi)達(dá)43.29萬 hm2，其中23.4萬 hm2重災(zāi)；2002年，南四湖入湖河道斷流53條，南四湖7月中旬干涸，生態(tài)環(huán)境遭受嚴(yán)重破壞[11]。

參考《濟(jì)寧市地方志》和歷年降水量資料，收集統(tǒng)計(jì)湖區(qū)旱情的分布范圍，并通過GIS插值得到湖區(qū)旱情分布柵格圖，以其像元值度量干旱風(fēng)險(xiǎn)源的強(qiáng)度(圖2)。

4.2.2 洪澇

隨著建國后南四湖的一系列控制性工程的建成，南四湖的洪澇問題得到了有效控制。濟(jì)寧市1949—2005年氣象統(tǒng)計(jì)資料顯示，年降雨量在650 mm以上的年份為33a，其中1964、1990和2004年降雨量分別為1150.1、1071、1723.7 mm，2003年更是高達(dá)2068.1 mm，均是大澇之年。

解譯湖區(qū)高水位2004年6月和低水位2010年9月兩時(shí)陸地衛(wèi)星5號搭載的主題成像傳感器(TM，Thematic Mapper)影像中的水體信息，以面積差異度量湖區(qū)洪澇的空間分布，并認(rèn)為強(qiáng)度是一致的(圖2)。

4.2.3 水體污染

南四湖湖水中污染物較高的項(xiàng)目是有機(jī)物、氨氮(NH3-N、NO3-N、NO2-N)、酚、汞和砷，有機(jī)物的普遍超標(biāo)，顯示湖域局部已經(jīng)遭到了較嚴(yán)重的有機(jī)物污染[22]。南四湖入湖河流的水質(zhì)情況很大程度上決定著南四湖的水質(zhì)優(yōu)劣情況，南陽湖和獨(dú)山湖是4個(gè)湖區(qū)中有機(jī)物污染比較嚴(yán)重的區(qū)域[20]。

結(jié)合濟(jì)寧市環(huán)保局在南四湖湖中11個(gè)長期監(jiān)測點(diǎn)多年監(jiān)測資料和南四湖主要入湖河流的廢污排放量[11]，插值得到湖區(qū)水體污染空間分布柵格圖，并以像元值度量水體污染的強(qiáng)度(圖2)。

圖2 風(fēng)險(xiǎn)源強(qiáng)度空間分布圖Fig.2 Distribution of the risk sources

5 暴露與危害分析

暴露分析就是研究區(qū)各風(fēng)險(xiǎn)源與生態(tài)受體的接觸關(guān)系，而與之相聯(lián)系的危害分析則是要確定風(fēng)險(xiǎn)源對風(fēng)險(xiǎn)受體的損害程度[2- 3,12- 13]。

5.1 暴露在風(fēng)險(xiǎn)中的受體

本研究中主要的生態(tài)風(fēng)險(xiǎn)受體是水域生態(tài)系統(tǒng)、濕地生態(tài)系統(tǒng)。選用2010年9月18日，條帶號122，行編號為35、36的兩景空間分辨率30 m的TM影像，以ENVI48為平臺，進(jìn)行拼接裁剪得到研究區(qū)，用人工神經(jīng)網(wǎng)絡(luò)法提取生態(tài)受體信息，并用2011年12月，2012年9月實(shí)地考察所獲得的實(shí)際地表覆被資料進(jìn)行精度檢驗(yàn)，其中訓(xùn)練樣本分離度均在1.85以上，檢驗(yàn)精度為97.5%，Kappa系數(shù)為0.9。

5.2 生態(tài)系統(tǒng)的地位與易損性

南四湖區(qū)域的各類生境所處的位置不同，組合方式更是有所差異。不同的生境類型在維護(hù)生物多樣性、保護(hù)物種、完善生態(tài)系統(tǒng)的結(jié)構(gòu)和功能、促進(jìn)景觀結(jié)構(gòu)自然演替等方面的作用是不一樣的[11]。同強(qiáng)度、同一風(fēng)險(xiǎn)源作用于不同的生境類型時(shí)，可能對整個(gè)區(qū)域的生態(tài)結(jié)構(gòu)與功能產(chǎn)生不同強(qiáng)度的危害[13]。景觀分布在遙感中的表現(xiàn)具有直觀性和準(zhǔn)確性，它可以準(zhǔn)確的顯示出各種生態(tài)影響的空間分布和梯度變化特征，使各種空間分析的手段成為可能[15]。本文引入景觀格局指數(shù)和生態(tài)損失度來表達(dá)生態(tài)系統(tǒng)的地位與易損性[8- 10,23]。景觀格局指數(shù)是一個(gè)相對量，本研究計(jì)算該指數(shù)時(shí)將空間尺度范圍擴(kuò)展到湖區(qū)周邊的縣域(表2)。

表2 生態(tài)系統(tǒng)類型景觀格局指數(shù)Table 2 Calculation methods of landscape pattern indices

表中公式中,ni為生態(tài)系統(tǒng)類型i的斑塊數(shù)量；Ai為生態(tài)系統(tǒng)i的總面積，A是區(qū)域總面積；Ii是景觀類型的距離指數(shù)；Hmax是生態(tài)系統(tǒng)多樣性指數(shù)的最大值；Pi為i類生態(tài)系統(tǒng)所占面積的比例；N為生態(tài)系統(tǒng)類型的總數(shù)；a、b、c為相應(yīng)指數(shù)的權(quán)重。生態(tài)損失度Ri的大小由結(jié)構(gòu)指數(shù)Si和脆弱度Fi決定。Si反映了相對于整個(gè)研究區(qū)域而言的各個(gè)生態(tài)系統(tǒng)受到人為活動干擾的程度，Si由破碎度Ci、分離度Ni和優(yōu)勢度Di疊加得到，三者權(quán)重a+b+c=1，通過分析權(quán)衡并結(jié)合已有研究經(jīng)驗(yàn)[6,8- 10,24- 26]，賦a為0.6，b為0.3，c為0.1。Fi是生態(tài)系統(tǒng)脆弱度，反映的是各個(gè)生態(tài)系統(tǒng)本身的抗干擾能力，亦可以反映其對環(huán)境改變的敏感性[8]，景觀的自然演替過程與生態(tài)系統(tǒng)的脆弱度關(guān)系緊密，一般來說處于初級演替階段、食物鏈單一并且生物多樣性差的生態(tài)系統(tǒng)脆弱性較高[2,8]，結(jié)合已有研究經(jīng)驗(yàn)[5- 9,18- 20,24- 25]，并根據(jù)專家打分并歸一化結(jié)果得到。

根據(jù)表2中公式，計(jì)算得表3結(jié)果。

表3 湖區(qū)2010年景觀格局指數(shù)Table 3 Landscape pattern indices of Nansi lake in 2010

Ci，Ni，Di，Si，F(xiàn)i，Ri的意義及算法參照表(2)

6 綜合生態(tài)風(fēng)險(xiǎn)評價(jià)

綜合生態(tài)風(fēng)險(xiǎn)評價(jià)就是要結(jié)合風(fēng)險(xiǎn)源分析、受體分析和暴露與危害分析對區(qū)域的生態(tài)風(fēng)險(xiǎn)大小給出綜合的評價(jià)，本文即是計(jì)算在研究方法部分中提到的Pk、Dk、Rk值。

模型中Ai、A由文中4.1處遙感解譯后統(tǒng)計(jì)得到，Ri由表3得到，Pkj由表1得到，關(guān)于?kj即風(fēng)險(xiǎn)源權(quán)重，以往的許多生態(tài)風(fēng)險(xiǎn)研究，基本上都是采用層次分析法得到針對整個(gè)研究區(qū)域的權(quán)重值[5- 10]，本研究將論文第4部分獲得的風(fēng)險(xiǎn)源插值結(jié)果置于風(fēng)險(xiǎn)小區(qū)內(nèi)，并用其在小區(qū)內(nèi)分布面積的比例與該風(fēng)險(xiǎn)的強(qiáng)度等級的乘積作為該小區(qū)的風(fēng)險(xiǎn)源權(quán)重，這樣，同一風(fēng)險(xiǎn)源在整個(gè)區(qū)域內(nèi)沒有統(tǒng)一的權(quán)重，只存在針對風(fēng)險(xiǎn)小區(qū)的風(fēng)險(xiǎn)源的權(quán)重。將Pk插值得圖3，可以看到這種權(quán)重處理表現(xiàn)出了良好的空間性。

利用GIS字段計(jì)算功能最終得到1178個(gè)生態(tài)風(fēng)險(xiǎn)值Rk，統(tǒng)計(jì)得表4，以Rk為風(fēng)險(xiǎn)小區(qū)中心點(diǎn)的屬性值插值得到生態(tài)風(fēng)險(xiǎn)空間分布圖(圖3)。

南四湖湖區(qū)1178個(gè)風(fēng)險(xiǎn)小區(qū)的生態(tài)風(fēng)險(xiǎn)值差異明顯，標(biāo)準(zhǔn)差1.662，最大值與最小值相差90倍以上(表4)。統(tǒng)計(jì)發(fā)現(xiàn)風(fēng)險(xiǎn)值在0—3之間出現(xiàn)了733次，3—6之間出現(xiàn)了409次，6—10之間出現(xiàn)了46次，由于生態(tài)風(fēng)險(xiǎn)值本身表達(dá)的是相對整個(gè)研究區(qū)域的風(fēng)險(xiǎn)大小，是一個(gè)相對量，故不再進(jìn)行歸一化處理。

表4 南四湖生態(tài)風(fēng)險(xiǎn)值統(tǒng)計(jì)表Table 4 The statistics of Nansi lake ecological risk value

圖3 綜合風(fēng)險(xiǎn)概率(PK)、生態(tài)綜合損失度(DK)、綜合生態(tài)風(fēng)險(xiǎn)(RK)空間分布圖Fig.3 The spatial distribution of integrated risk probability (PK), integrated eco-loss (DK) and integrated eco-risk (RK)

從空間分布上看，生態(tài)風(fēng)險(xiǎn)表現(xiàn)出上級湖大于下級湖，濱湖地帶大于湖心敞水區(qū)的分異性(圖3)，具體而言，在南四湖湖區(qū)，相對生態(tài)風(fēng)險(xiǎn)較高的區(qū)域主要集中在：洸府河、泗河、白馬河和東魚河河口附近，南陽鎮(zhèn)附近水域以及二級壩周邊地區(qū)。其主要原因可能是：1)與注入上下級湖的河流有關(guān)：目前注入南四湖的河流的主要功能之一便是納污[11]，而注入上級湖的河流的數(shù)量與密度都要大于注入下級湖的河流，大量河流排污排廢增加了上級湖的生態(tài)風(fēng)險(xiǎn)。2)由于二級壩的存在，下級湖受到洪澇和干旱災(zāi)害機(jī)率都要小于上級湖，水位變動對下級湖的生態(tài)擾動相對比較小。3)與研究區(qū)降水有關(guān)：南四湖多年平均降水由南向北遞減。4)湖濱地區(qū)受水位影響大：南四湖水旱災(zāi)害頻繁，湖濱地區(qū)在豐水期為澇，枯水則旱，生境的頻繁變化增加了該區(qū)域的生態(tài)風(fēng)險(xiǎn)。5)南四湖流速很慢，換水周期為503 d，由河流帶來的生活污水，工業(yè)廢水以及農(nóng)業(yè)灌溉回水中的有機(jī)污染物、有毒物質(zhì)等累積在河流入湖口附近。

7 結(jié)論與討論

研究結(jié)果顯示了南四湖生態(tài)風(fēng)險(xiǎn)的空間分異特征，具體找到了湖區(qū)生態(tài)風(fēng)險(xiǎn)較高的地區(qū)，同時(shí)也量化了風(fēng)險(xiǎn)受體以及湖區(qū)的干旱、洪澇和水體污染3個(gè)風(fēng)險(xiǎn)源，并且找到了他們的空間分布規(guī)律。

區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)是一個(gè)多學(xué)科交叉的研究領(lǐng)域，涉及內(nèi)容較多，本文研究存在一些不足。關(guān)于風(fēng)險(xiǎn)源，在有條件的情況下應(yīng)盡量全面地反映研究區(qū)的生態(tài)脅迫，南四湖湖區(qū)的圍墾現(xiàn)象和湖泊的泥沼化亦是南四湖濕地退化的重要原因，在以后的研究中因受重視。關(guān)于生態(tài)受體，單從景觀學(xué)上說明受體的地位和脆弱性顯得不夠全面，景觀破碎化與生物多樣性的關(guān)系尚待進(jìn)一步的深入研究。同許多區(qū)域生態(tài)風(fēng)險(xiǎn)的研究一樣，本文最終獲得的綜合生態(tài)風(fēng)險(xiǎn)值在指示生態(tài)事件或生態(tài)終點(diǎn)方面存在不足。

研究中基于GIS技術(shù)對風(fēng)險(xiǎn)源權(quán)重的確定進(jìn)行了新的嘗試，很好的反映出了風(fēng)險(xiǎn)源的空間性，與層次分析法得到的風(fēng)險(xiǎn)源權(quán)重相比，本研究考慮到了風(fēng)險(xiǎn)源在空間上有不同強(qiáng)度表現(xiàn)的可能性，是對區(qū)域生態(tài)風(fēng)險(xiǎn)研究技術(shù)的探討和區(qū)域風(fēng)險(xiǎn)源權(quán)重計(jì)算方面的提高。

[1] Hunsaker C T, Grahm R L, Suter G W, O′Neill R V, Barnthouse L W, Gardner R H. Assessing ecological risk on a regional scale. Environmental Management, 1990, 14(3): 325- 332.

[2] 肖楊, 毛顯強(qiáng). 區(qū)域景觀生態(tài)風(fēng)險(xiǎn)空間分析. 中國環(huán)境科學(xué), 2006, 26(5): 623- 626.

[3] 周婷, 蒙古軍. 區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)方法研究進(jìn)展. 生態(tài)學(xué)雜志, 2009, 28(4): 762- 767.

[4] Wagner H H, Fortin M J. Spatial analysis of landscapes: concepts and statistics. Ecology, 2005, 86(8): 1975- 1987.

[5] 付在毅, 許學(xué)工, 林輝平, 王憲禮. 遼河三角洲濕地區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià). 生態(tài)學(xué)報(bào), 2001, 21(3): 365- 373.

[6] 許學(xué)工, 林輝平, 付在毅. 黃河三角洲濕地區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià). 北京大學(xué)學(xué)報(bào): 自然科學(xué)版, 2001, 37(1): 111- 120.

[7] 盧宏瑋, 曾光明, 謝更新, 張碩輔, 黃國和, 金相燦, 劉鴻亮. 洞庭湖流域區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià). 生態(tài)學(xué)報(bào), 2003, 23(12): 2521- 2530.

[8] 李謝輝, 王磊, 李景宜. 基于GIS的渭河下游河流沿線區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià). 生態(tài)學(xué)報(bào), 2009, 29(10): 5524- 5534.

[9] 李謝輝, 李景宜. 基于GIS的區(qū)域景觀生態(tài)風(fēng)險(xiǎn)分析——以渭河下游河流沿線區(qū)域?yàn)槔? 干旱區(qū)研究, 2008, 25(6): 899- 903.

[10] 陳鵬, 潘曉玲. 干旱區(qū)內(nèi)陸流域區(qū)域景觀生態(tài)風(fēng)險(xiǎn)分析——以阜康三工河流域?yàn)槔? 生態(tài)學(xué)雜志, 2003, 22(4): 116- 120.

[11] 沈吉, 張祖陸, 楊麗原, 孫慶義. 南四湖——環(huán)境與資源研究. 北京: 地震出版社, 2009: 4- 5.

[12] 李自珍, 何俊紅. 生態(tài)風(fēng)險(xiǎn)評價(jià)與風(fēng)險(xiǎn)決策模型及應(yīng)用——以河西走廊荒漠綠洲開發(fā)為例. 蘭州大學(xué)學(xué)報(bào): 自然科學(xué)版, 1999, 35(3): 150- 155.

[13] 蒙吉軍, 趙春紅. 區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)指標(biāo)體系. 應(yīng)用生態(tài)學(xué)報(bào), 2009, 20(4): 983- 990.

[14] 濟(jì)寧市土壤. 濟(jì)寧市土壤肥料工作站. 濟(jì)南: 山東省地圖出版社, 1990: 39- 190.

[15] 陳利頂,劉洋,呂一河,馮曉明,傅伯杰. 景觀生態(tài)學(xué)中的格局分析: 現(xiàn)狀、困境與未來. 生態(tài)學(xué)報(bào), 2008, 28(11): 5522- 5531.

[16] 布仁倉, 李秀珍, 胡遠(yuǎn)滿, 常禹, 賀紅士. 尺度分析對景觀格局指標(biāo)的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2003, 14(12): 2181- 2186.

[17] 申衛(wèi)軍, 鄔建國, 任海, 林永標(biāo), 李明輝. 空間幅度變化對景觀格局分析的影響. 生態(tài)學(xué)報(bào), 2003, 23(11): 2219- 2231.

[18] 付在毅, 許學(xué)工. 區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià). 地球科學(xué)進(jìn)展, 2001, 16(2): 267- 271.

[19] 曾勇. 區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)——以呼和浩特市區(qū)為例. 生態(tài)學(xué)報(bào), 2010, 30(3): 668- 673.

[20] 高學(xué)平,李文猛,張晨,趙世新.入湖河流對南四湖水質(zhì)的影響.水資源保護(hù), 2013, 29(2): 1- 5.

[21] 陳立項(xiàng), 傅伯杰. 黃河三角洲地區(qū)人類活動對景觀結(jié)構(gòu)的影響分析——以山東省東營市為例. 生態(tài)學(xué)報(bào), 1996, 22(4): 116- 120.

[22] 武周虎, 張可, 金玲仁, 楊連寬, 張建. 南四湖水質(zhì)空間分布特征分析與改善效果評估. 水資源保護(hù), 2012, 28(6): 1- 7.

[23] 鞠強(qiáng). 博斯騰湖區(qū)域生態(tài)風(fēng)險(xiǎn)評價(jià)研究 [D]. 烏魯木齊: 新疆大學(xué), 2005.

[24] 許妍, 高俊峰, 高永年. 基于土地利用動態(tài)變化的太湖地區(qū)景觀生態(tài)風(fēng)險(xiǎn)評價(jià). 湖泊科學(xué), 2011, 23(4): 642- 648.

[25] 林揚(yáng)碧, 陳秋波, 鄒偉, 鐘同暢. 區(qū)域景觀結(jié)構(gòu)的生態(tài)風(fēng)險(xiǎn)評價(jià)方法與案例研究. 環(huán)境科學(xué)與管理, 2009, 34(4): 186- 189.

[26] 盧遠(yuǎn), 蘇文靜, 華璀, 湯傳勇. 左江上游流域景觀生態(tài)風(fēng)險(xiǎn)評價(jià). 熱帶地理, 2010, 30(5): 496- 502.

Regional ecological risk assessment in Nansi lake based on RS and GIS

ZHANG Yazhou, XIE Xiaoping*

SchoolofGeography&Tourism,QufuNormalUniversity,Rizhao276826,China

Nansi Lake was selected as our study site for this study. Droughts, floods and water pollution for this site were used as a source of the ecological risk (ER) with information by the technologies of remote sensing images and geographic information system. Landscape indices were developed to evaluate ecosystem status and fragile ER receptors. Based on the principles of the measurement of ER, the ER values for Nansi Lake were obtained by integrated comprehensive loss probability and the ER degree. The results showed that the ER for Nansi Lake area had significant differences between upper and lower lakes with spatial distributions. It also showed that ER for Nansi Lake had a spatial and regional heterogeneity.

remote-sensing image; GIS; regional ecological risk; Nansi Lake

國家自然科學(xué)基金(41072164); 曲阜師范大學(xué)南四湖濕地生態(tài)與環(huán)境保護(hù)省高校重點(diǎn)實(shí)驗(yàn)室資助

2013- 05- 10;

日期:2014- 04- 17

10.5846/stxb201305101012

*通訊作者Corresponding author.E-mail: xp.xie@263.net

張雅洲, 謝小平.基于RS和GIS的南四湖生態(tài)風(fēng)險(xiǎn)評價(jià).生態(tài)學(xué)報(bào),2015,35(5):1371- 1377.

Zhang Y Z, Xie X P.Regional ecological risk assessment in Nansi lake based on RS and GIS.Acta Ecologica Sinica,2015,35(5):1371- 1377.