2005～2014年平潭縣土地利用變化及生態(tài)風(fēng)險評價

施婷婷，許章華,2

(1. 福州大學(xué) 環(huán)境與資源學(xué)院，福建 福州 350116；2. 福建省資源環(huán)境監(jiān)測與可持續(xù)經(jīng)營利用重點實驗室，福建 福州 350002)

?

2005～2014年平潭縣土地利用變化及生態(tài)風(fēng)險評價

施婷婷1，許章華1,2

(1. 福州大學(xué) 環(huán)境與資源學(xué)院，福建 福州 350116；2. 福建省資源環(huán)境監(jiān)測與可持續(xù)經(jīng)營利用重點實驗室，福建 福州 350002)

為了研究區(qū)域土地利用變化對生態(tài)環(huán)境造成的影響，以2005年Landsat 5 TM與2014年Landsat 8 OLI多光譜遙感影像數(shù)據(jù)為主要數(shù)據(jù)源，采用決策樹分類方法輔以目視解譯進(jìn)行平潭縣土地利用分類；從多風(fēng)險源、多風(fēng)險受體、多風(fēng)險效應(yīng)3個方面構(gòu)建平潭縣生態(tài)風(fēng)險評價指標(biāo)體系和模型，分析2005～2014年平潭縣生態(tài)風(fēng)險變化.結(jié)果表明：2005～2014年來，填海造陸、海島開發(fā)建設(shè)等導(dǎo)致平潭縣土地利用類型發(fā)生著較大變化，建設(shè)用地增加3 405.6 hm2，變化幅度最大；較之于2005年，2014年總體平均生態(tài)風(fēng)險值上升0.03，大部分鄉(xiāng)鎮(zhèn)平均生態(tài)風(fēng)險值上升；極低生態(tài)風(fēng)險、低生態(tài)風(fēng)險和中等生態(tài)風(fēng)險等級的面積減少，高生態(tài)風(fēng)險面積增加.平潭縣生態(tài)風(fēng)險的定量化綜合評價可為區(qū)域生態(tài)環(huán)境與社會經(jīng)濟(jì)可持續(xù)發(fā)展提供科學(xué)依據(jù)，同時為土地利用生態(tài)風(fēng)險指標(biāo)體系與模型的科學(xué)合理構(gòu)建提供參考.

區(qū)域生態(tài)風(fēng)險評價； 土地利用變化； 風(fēng)險管理； 平潭縣

土地利用/覆蓋變化(LUCC)反映自然條件改變與人類活動對土地生態(tài)系統(tǒng)的綜合影響, 是當(dāng)前全球變化研究的重要內(nèi)容之一.目前，關(guān)于土地利用的研究日趨完善并向變化監(jiān)測[1-3]、驅(qū)動力機(jī)制[4]、數(shù)據(jù)庫構(gòu)建[5]和資源環(huán)境效應(yīng)[6]等綜合化的方向發(fā)展.近年來人類社會的快速發(fā)展引起的土地利用變化導(dǎo)致區(qū)域生態(tài)系統(tǒng)的結(jié)構(gòu)和功能造成較大影響，故基于土地利用變化的區(qū)域生態(tài)風(fēng)險評價研究開始受到關(guān)注[6-7].區(qū)域生態(tài)風(fēng)險評價是在區(qū)域尺度上對不同類型的生態(tài)系統(tǒng)進(jìn)行環(huán)境污染、人為活動或自然災(zāi)害等多風(fēng)險源對評價終點造成的負(fù)面效應(yīng)及危害程度的風(fēng)險評估[8].眾多研究表明[9-10]，土地利用變化與區(qū)域生態(tài)風(fēng)險有著密不可分的關(guān)系，土地利用引起的生態(tài)風(fēng)險具有區(qū)域性和積累性，能夠直接反映區(qū)域中各種潛在的生態(tài)影響及其空間異質(zhì)性.因此，基于土地利用的區(qū)域生態(tài)風(fēng)險評價能夠預(yù)測一定區(qū)域范圍內(nèi)生態(tài)風(fēng)險發(fā)生的概率及可能產(chǎn)生的負(fù)面效應(yīng)，提出預(yù)防風(fēng)險發(fā)生的措施，對于改善生態(tài)環(huán)境質(zhì)量、提高生態(tài)系統(tǒng)價值、促進(jìn)區(qū)域環(huán)境管理與風(fēng)險決策等具有重要意義.目前國內(nèi)外學(xué)者針對不同的研究區(qū)在評價指標(biāo)體系、方法與模型等方面開展了大量研究[11-13]，但仍缺乏科學(xué)統(tǒng)一合適的區(qū)域生態(tài)風(fēng)險評價指標(biāo)體系.如Angela[11]等構(gòu)建基于壓力因子-受體-生境三者相互關(guān)系的模型來研究Codorus Creek流域的生態(tài)風(fēng)險，其選擇的壓力因子包括土地利用、水土流失、地表覆蓋和廢水排放等；評價端點有水量、水質(zhì)、本地魚類數(shù)量減少、水生生物減少、恢復(fù)損害和暴雨危害等；傅麗華[12]等選取景觀及土地利用變化類型、土地利用分級程度、土地利用變化率和不同景觀的生態(tài)服務(wù)價值作為主要評價指標(biāo)，建立生態(tài)風(fēng)險評價模型，對長株潭城市群核心區(qū)土地利用進(jìn)行了生態(tài)風(fēng)險評價；孫洪波[13]等采用相對風(fēng)險模型，從風(fēng)險源、風(fēng)險受體和風(fēng)險效應(yīng)3個方面選擇指標(biāo)構(gòu)建生態(tài)風(fēng)險指標(biāo)體系，分析南京市地區(qū)土地利用生態(tài)風(fēng)險空間分異規(guī)律及影響因素.由此可見，不同的研究者對于生態(tài)風(fēng)險評價指標(biāo)體系與模型的構(gòu)建各不相同，目前仍未形成科學(xué)規(guī)范的評價標(biāo)準(zhǔn)；同時由于不同研究區(qū)的各類詳細(xì)監(jiān)測數(shù)據(jù)的獲取比較困難，在對各種生態(tài)系統(tǒng)的生態(tài)過程和作用機(jī)理尚未清楚的情況下構(gòu)建評價指標(biāo)體系與模型是較為困難的.因此，進(jìn)一步完善基于土地利用變化的區(qū)域生態(tài)風(fēng)險評價定量研究方法，構(gòu)建合適的生態(tài)風(fēng)險評價指標(biāo)體系與模型具有重要的意義.筆者從多風(fēng)險源、多風(fēng)險受體、多風(fēng)險效應(yīng)3個方面來構(gòu)建土地利用變化下的生態(tài)風(fēng)險評價指標(biāo)體系與模型，綜合定量地評價平潭縣生態(tài)風(fēng)險，旨在為平潭縣生態(tài)環(huán)境的綜合管理和社會經(jīng)濟(jì)的可持續(xù)發(fā)展提供科學(xué)依據(jù)，同時為土地利用生態(tài)風(fēng)險指標(biāo)體系與模型的科學(xué)合理構(gòu)建提供參考.

1 研究區(qū)概況

圖1 平潭縣地理區(qū)位示意圖

平潭縣位于福建省東部，與臺灣隔海相望，地理位置介于25°15′ ～25°45′ N，119°32′ ～120°06′ E之間(如圖1所示)，地勢低平，中部略高，主要地貌類型為丘陵和海積平原.平潭縣屬于亞熱帶海洋性季風(fēng)氣候，溫?zé)釢駶?，冬暖夏?2009年設(shè)立平潭綜合實驗區(qū)，根據(jù)平潭綜合實驗區(qū)2013年國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計公報數(shù)據(jù)顯示，平潭縣2013年總?cè)丝诩s41.8萬，全年實現(xiàn)生產(chǎn)總值155.39億元，比上年增長13.0%.近十年來，平潭縣正處于大力發(fā)展的時期，隨著區(qū)域人口增長和社會經(jīng)濟(jì)的快速發(fā)展，憑借著國家和地方政策的大力支持，基礎(chǔ)設(shè)施建設(shè)、填海造陸工程等人類活動引起全區(qū)土地利用類型與景觀格局發(fā)生劇烈變化，造成環(huán)境污染、森林銳減、水土流失等區(qū)域性環(huán)境問題日益突出，制約區(qū)域的可持續(xù)發(fā)展.

2 研究方法

2.1 數(shù)據(jù)收集與影像預(yù)處理 收集的數(shù)據(jù)包括：平潭縣2005年Landsat 5 TM與2014年Landsat 8 OLI多光譜遙感影像，并對影像進(jìn)行輻射校正、幾何校正、剪裁等處理；平潭縣1 ∶10 000地形圖，用于影像的配準(zhǔn)；行政區(qū)劃邊界圖，用于影像的裁剪；2005年與2014年土地利用現(xiàn)狀圖，主要從Google Earth中獲取，用于決策樹分類.投影類型為Transverse Mercator，橢球體為Krasovsky，誤差控制在0.5個像元內(nèi).2.2 決策樹分類 決策樹分類是按照一定規(guī)則把遙感數(shù)據(jù)集逐級細(xì)分以得到具有不同屬性的各個子類別[14].RuleGen是ENVI軟件下關(guān)于決策樹分類的插件，運行原理基于CART算法，基于ENVI自帶Decision Tree進(jìn)行制作與執(zhí)行決策樹規(guī)則.首先將歸一化差值植被指數(shù)(NDVI)與預(yù)處理后的原始影像波段進(jìn)行組合，結(jié)合Google Earth進(jìn)行一定數(shù)量的訓(xùn)練區(qū)樣本選取，然后利用訓(xùn)練區(qū)樣本獲取專家知識規(guī)則，進(jìn)行RuleGen決策樹分類,最后對分類后的遙感影像輔以目視解譯，結(jié)合Google Earth對錯誤的分類進(jìn)行糾正，提高分類精度.

2.3 土地利用變化分析 采用定性與定量相結(jié)合的方式，選取土地利用變化面積、變化幅度和變化速率等3個方面進(jìn)行土地利用變化分析.將2005年和2014年2期土地利用分類圖，在ERDAS中進(jìn)行空間建模，生成各類土地利用變化圖，并對此進(jìn)行統(tǒng)計計算生成2005～2014年平潭縣土地利用轉(zhuǎn)移矩陣，來定量描述各土地利用類型間的轉(zhuǎn)化.

2.4 平潭縣生態(tài)風(fēng)險評價指標(biāo)體系的構(gòu)建 結(jié)合國內(nèi)外關(guān)于區(qū)域生態(tài)風(fēng)險評價研究成果，廣泛征求各方面專家意見，借鑒相對風(fēng)險模型(RRM)，從多風(fēng)險源、多風(fēng)險受體、多風(fēng)險效應(yīng)3個方面構(gòu)建土地利用變化下的平潭縣生態(tài)風(fēng)險評價指標(biāo)體系.多風(fēng)險源強(qiáng)反映區(qū)域可能面臨的多種不利影響，如自然災(zāi)害、人類活動等，主要考慮與土地利用活動有關(guān)的壓力因子以及土地利用強(qiáng)度.多風(fēng)險受體指在風(fēng)險評價中生態(tài)系統(tǒng)可能受到來自風(fēng)險源不利作用的組成部分.從景觀生態(tài)學(xué)的角度出發(fā)，認(rèn)為人類活動引發(fā)的土地利用變化在很大程度上反映生態(tài)系統(tǒng)景觀破碎化以及景觀格局變化，故指標(biāo)多選取與風(fēng)險受體空間承載力相關(guān)的景觀格局指數(shù).多風(fēng)險效應(yīng)指標(biāo)反映風(fēng)險受體對風(fēng)險源的響應(yīng)，指標(biāo)多選用自然生態(tài)系統(tǒng)及復(fù)合生態(tài)系統(tǒng)的功能性指標(biāo).

2.4.1 多風(fēng)險源指標(biāo)選取及計算

1) 建設(shè)用地比例是城市擴(kuò)張和人口增長的重要標(biāo)志，比例越高，意味人類活動對原有的生態(tài)系統(tǒng)的破壞越大，帶來一些如熱島效應(yīng)、生物多樣性銳減、有毒有害物質(zhì)劇增等負(fù)面效應(yīng).因此，建設(shè)用地比例越高，對生態(tài)環(huán)境的影響越大.建設(shè)用地成為區(qū)域生態(tài)環(huán)境惡化的主要風(fēng)險源.計算公式

P=D/L×100% ，

(1)

其中，P是建設(shè)用地比例，D為建設(shè)用地面積，L為土地總面積.

2)耕地墾殖指數(shù)是保障區(qū)域人類生存和發(fā)展的基礎(chǔ)條件.耕地對生態(tài)環(huán)境的負(fù)面效應(yīng)主要表現(xiàn)在水土流失、水體污染、土地退化和土地污染等.耕地的開墾導(dǎo)致地表植被對土壤的保護(hù)性降低，一旦暴雨沖刷，易造成水土流失.同時，耕地意味著大量的農(nóng)藥和化肥流入土地，一些難降解的持久性污染物將留在土壤里，形成土地污染.因此，耕地墾殖指數(shù)越高，意味著其對生態(tài)風(fēng)險的負(fù)面影響程度越高.計算公式[15]

K=F/L×100% ，

(2)

其中，K是耕地墾殖指數(shù)，F(xiàn)為耕地面積，L為土地總面積.

3)土地利用綜合指數(shù)反映土地利用的廣度和深度[16]，在一定程度上體現(xiàn)自然因素和人文因素等綜合影響對土地利用方式和結(jié)構(gòu)的作用，其值越大，代表土地利用強(qiáng)度越大，對生態(tài)風(fēng)險的影響程度越高.參考劉紀(jì)遠(yuǎn)[17]等提出的土地利用綜合指數(shù)，計算公式

(3)其中，La為土地利用綜合指數(shù)，Ai為第i類土地利用程度分級指數(shù)，Ci為第i類土地面積占總面積的比重.

2.4.2 多風(fēng)險受體指標(biāo)選取及計算

1)破碎度指數(shù)反映景觀空間結(jié)構(gòu)的復(fù)雜性，在一定程度上反映人為因素對景觀的干擾程度，其值越大說明景觀受到較為強(qiáng)烈的干擾，造成抵抗干擾的能力越弱，景觀易遭受破壞，生態(tài)風(fēng)險越高.計算公式[18]

Ci=Ni/Ai，

(4)

其中，Ci為景觀i的破碎度，Ni為景觀i的斑塊數(shù)，Ai為景觀i的總面積.

2)景觀形狀指數(shù)體現(xiàn)生態(tài)系統(tǒng)整體結(jié)構(gòu)的復(fù)雜性，其值越大，景觀生態(tài)系統(tǒng)的復(fù)雜度越高，受干擾程度越高，生態(tài)系統(tǒng)穩(wěn)定性越低，生態(tài)風(fēng)險越高.計算公式[18]

(5)

其中，LSI為景觀形狀指數(shù)，E為景觀中所有斑塊邊界的總長度，A為景觀面積.

3)香農(nóng)多樣性指數(shù)是指景觀元素在結(jié)構(gòu)、功能以及隨時間變化方面的多樣性，反映景觀類型的豐富度和復(fù)雜度.土地利用越豐富，破碎度越高，其香農(nóng)多樣性指數(shù)越高，生態(tài)風(fēng)險越高.計算公式[18]

(6)

其中，SHDI為香農(nóng)多樣性指數(shù)，Pi是景觀類型i所占面積的比例，n為景觀類型數(shù)目.

4)聚集度指數(shù)是指景觀中不同斑塊類型的非隨機(jī)性或聚集程度，反映景觀要素的空間配置特征，其值越小，景觀越離散，生態(tài)風(fēng)險越高.計算公式[19]

(7)

其中，C為聚集度指數(shù)，Cmax為聚集度指數(shù)的最大值(2ln n)，n為景觀中斑塊類型總數(shù)，Pij為斑塊類型i與j相鄰的概率.

2.4.3 多風(fēng)險效應(yīng)指標(biāo)選取及計算

1)森林覆蓋率體現(xiàn)區(qū)域內(nèi)森林資源豐富程度以及實現(xiàn)綠化的程度，是反映生態(tài)平衡狀況的重要指標(biāo).森林有著調(diào)節(jié)氣候、水土保持、維持生物多樣性等作用，區(qū)域內(nèi)森林覆蓋率的變化決定著生態(tài)環(huán)境質(zhì)量.因此其值越大，生態(tài)風(fēng)險越低.計算公式

C=W/L×100% ，

(8)

其中，C是森林覆蓋率，W為林地面積，L為土地總面積.

2)生物豐度指數(shù)反映區(qū)域內(nèi)生物豐度的豐貧程度.生物豐度指數(shù)越高，反映區(qū)域內(nèi)風(fēng)險發(fā)生時風(fēng)險受體受影響的程度越大，生態(tài)風(fēng)險越高.參考《生態(tài)環(huán)境狀況評價國家技術(shù)規(guī)范(試行)》[20]，計算公式

(9)

其中，B是生物豐度指數(shù)，S1，S2，S3，S4，S5，S6分別表示區(qū)域內(nèi)林地、耕地、灘涂、水域、建設(shè)用地和未利用地的面積，S為土地總面積.

3)生態(tài)系統(tǒng)服務(wù)價值是指土地利用結(jié)構(gòu)與功能的變化直接影響著生態(tài)系統(tǒng)服務(wù)價值的變化.通過對平潭縣生態(tài)系統(tǒng)服務(wù)價值進(jìn)行估算，能較好地反映當(dāng)?shù)丨h(huán)境狀態(tài)和調(diào)節(jié)能力.因此，生態(tài)系統(tǒng)服務(wù)價值越高，生態(tài)風(fēng)險發(fā)生時風(fēng)險受體可能受到的影響較大，生態(tài)風(fēng)險越高.計算公式[21]

(10)

其中，V為生態(tài)系統(tǒng)服務(wù)價值指數(shù)，Aij為第j個區(qū)域i類土地利用類型面積，Aj為第j個樣方總面積，Vi為第i類土地利用類型的生態(tài)系統(tǒng)服務(wù)價值系數(shù).參照謝高地[22]等以及陳鵬[23]對生態(tài)系統(tǒng)服務(wù)價值的研究確定生態(tài)服務(wù)價值系數(shù).

采用層次分析法構(gòu)建層次結(jié)構(gòu)模型，通過咨詢地理學(xué)、生態(tài)環(huán)境等相關(guān)領(lǐng)域的9位專家意見和評分結(jié)果，采用1～9評判標(biāo)準(zhǔn)對兩兩之間的相對重要程度進(jìn)行比較賦值，構(gòu)建判斷矩陣.通過一致性檢驗來判斷專家打分的可信度，目標(biāo)層-準(zhǔn)則層判斷矩陣A-B的一致性比CR=0.051 6，準(zhǔn)則層-指標(biāo)層的判斷矩陣B1-C1～3，B2-C4～7，B3-C8～10一致性比例CR分別為0.017 6、0.022 7、0.051 6，可以看出，CR均小于0.1時，說明判斷矩陣的結(jié)果令人滿意.最后確定生態(tài)風(fēng)險評價各指標(biāo)的權(quán)重值，如表1所示.

表1 平潭縣生態(tài)風(fēng)險評價指標(biāo)因子權(quán)重值

2.5 平潭縣生態(tài)風(fēng)險評價模型的構(gòu)建 在區(qū)域生態(tài)風(fēng)險評價中，每個指標(biāo)對生態(tài)風(fēng)險值的貢獻(xiàn)率不同，采用綜合指數(shù)法[24]計算S，R，E的值，計算公式

ERI=S+R+E ，

(11)

(12)

(13)

(14)

其中，ERI為生態(tài)風(fēng)險指數(shù)，S為風(fēng)險源指標(biāo)，R為風(fēng)險受體指標(biāo)，E為風(fēng)險效應(yīng)指標(biāo)，Wi為指標(biāo)i的權(quán)重值，Xi為指標(biāo)i的量化值.

3 結(jié)果與分析

3.1 平潭縣土地利用分類及精度檢驗 參考《土地利用分類》中的分類標(biāo)準(zhǔn)，結(jié)合平潭縣實際情況，采用決策樹分類輔以目視解譯對2期遙感影像進(jìn)行分類，得到平潭縣土地利用分類圖，如圖2所示.采用基于誤差矩陣的精度評價方法進(jìn)行精度驗證.經(jīng)檢驗，2005年Landsat 5 TM多光譜遙感影像的總體分類精度為87.65%，Kappa系數(shù)為0.845 0，2014年Landsat 8 OLI多光譜遙感影像總體分類精度為88.57%，Kappa系數(shù)為0.858 7，且各類用地的使用者精度和生產(chǎn)者精度均在85%以上，分類結(jié)果比較理想.

a b圖2 平潭縣土地利用分類圖

3.2 近十年土地利用變化分析

3.2.1 土地利用總體變化特征分析 2005～2014年間平潭縣土地利用變化顯著，最明顯的是建設(shè)用地的增加和耕地的減少，如圖3所示.相比于2005年，2014年平潭縣林地、耕地、水域和灘涂的比例均有所下降，分別下降了3.79%、7.96%、3.56%、0.04%；建設(shè)用地與未利用地比例均增大，分別上升了10.61%、1.41%.從變化面積、幅度和變化速率分析(表2)，相比于2005年，平潭縣建設(shè)用地變化最大，增加3 405.6 hm2，增加幅度達(dá)89.84%，年均變化速度為8.98%；水域減少1 007.91 hm2，變化幅度為36.71%；平潭縣建設(shè)用地主要由耕地、林地、水域等轉(zhuǎn)變而來，由于填海造陸、海島開發(fā)建設(shè)等原因，水域、耕地的面積減少最多.未利用地增加了467.55 hm2，變化幅度較大，為46.88%.在遙感影像分類中，未利用地包括森林中的裸地、火燒跡地及待建設(shè)用地等，由于生態(tài)環(huán)境遭受一定程度的破壞，使得森林變成裸地；大規(guī)模的開發(fā)、征用農(nóng)用地等導(dǎo)致林地、耕地轉(zhuǎn)為待建設(shè)用地等，均導(dǎo)致未利用土地面積上升.

表2 2005～2014年平潭縣土地利用變化面積、幅度及速率

3.2.2 平潭縣土地利用轉(zhuǎn)移矩陣分析 利用2期土地利用分類圖在ERDAS中進(jìn)行空間建模，生成平潭縣2005～2014年各類土地利用類型變化分布圖，如圖4所示. 結(jié)果表明：2005～2014年間林地主要轉(zhuǎn)為耕地、建設(shè)用地和未利用地，面積分別為1 024.11 hm2、720.18 hm2、212.13 hm2，占林地總面積的比例分別為9.74%、6.85%、2.02%.耕地主要轉(zhuǎn)為建設(shè)用地、林地和未利用地，面積分別為2 265.3 hm2、972.81 hm2、236.97 hm2，占耕地總面積的比例分別為17.57%、7.54%、1.84%.水域主要轉(zhuǎn)為建設(shè)用地，面積為811.89 hm2，占水域總面積的比例為29.57%.建設(shè)用地、灘涂、未利用地的面積無明顯變化.近十年來平潭縣處于大規(guī)模的開發(fā)建設(shè)階段，基礎(chǔ)設(shè)施建設(shè)、填海造陸工程等人類活動導(dǎo)致了大量的耕地、林地和水域向建設(shè)用地轉(zhuǎn)移；自然災(zāi)害及濫砍濫伐、圍湖造田等人為因素造成土壤退化、水土流失，導(dǎo)致林地、耕地向未利用地轉(zhuǎn)移，均導(dǎo)致平潭縣土地利用類型發(fā)生巨大變化.

圖4 2005～2014年平潭縣土地利用類型變化分布圖

3.3 近十年平潭縣生態(tài)風(fēng)險評價

3.3.1 生態(tài)風(fēng)險評價指標(biāo)標(biāo)準(zhǔn)化處理 通過ArcGIS中鄰域分析計算3×3矩形的像元范圍內(nèi)中心像元的平均值生成多風(fēng)險源指標(biāo)、多風(fēng)險效應(yīng)指標(biāo)圖層，利用Fragstates 3.3的應(yīng)用移動窗口法生成多風(fēng)險受體指標(biāo)圖層，最后進(jìn)行歸一化處理.

3.3.2 多風(fēng)險源強(qiáng)度空間分布 通過對多風(fēng)險源指標(biāo)進(jìn)行地圖代數(shù)運算，生成2期多風(fēng)險源強(qiáng)度分布圖，如圖5所示.研究發(fā)現(xiàn)：與2005年相比，2014年平潭縣整體的風(fēng)險源強(qiáng)度值有增大的趨勢，主要在北厝鎮(zhèn)的西部、潭城鎮(zhèn)、流水鎮(zhèn)、澳前鎮(zhèn)以及中樓鄉(xiāng).由于近十年平潭縣經(jīng)濟(jì)的快速發(fā)展，基礎(chǔ)設(shè)施、填海造陸工程等導(dǎo)致林地、耕地、水域轉(zhuǎn)為建設(shè)用地，建設(shè)用地比例急劇上升；人口增長造成對土地的需求量增大，促使土地集約化利用，土地利用程度增大，對生態(tài)環(huán)境產(chǎn)生一定的威脅.

a 2005年 b 2014年圖5 平潭縣多風(fēng)險源強(qiáng)度分布圖

a 2005年 b 2014年圖6 平潭縣多風(fēng)險受體暴露度分布圖

3.3.3 多風(fēng)險受體暴露度空間分布 通過對多風(fēng)險受體指標(biāo)進(jìn)行地圖代數(shù)運算，生成2期多風(fēng)險受體暴露度分布圖，如圖6所示.研究發(fā)現(xiàn)：與2005年相比，2014年平潭縣風(fēng)險受體的暴露度整體有所上升.北厝鎮(zhèn)、澳前鎮(zhèn)、敖東鎮(zhèn)、蘇澳鎮(zhèn)、東庠鄉(xiāng)等鄉(xiāng)鎮(zhèn)進(jìn)行基礎(chǔ)設(shè)施建設(shè)、農(nóng)村居民點擴(kuò)建，加劇了區(qū)域景觀結(jié)構(gòu)的破壞，土地利用程度豐富，景觀破碎度高，導(dǎo)致鄉(xiāng)鎮(zhèn)沿海地帶風(fēng)險受體的暴露度增大.

3.3.4 多風(fēng)險效應(yīng)空間分布 通過對多風(fēng)險效應(yīng)指標(biāo)進(jìn)行地圖代數(shù)運算，生成2期多風(fēng)險效應(yīng)分布圖，如圖7所示.研究發(fā)現(xiàn)：與2005年相比，2014年平潭縣風(fēng)險效應(yīng)值整體有所降低.由于建設(shè)用地大量增加，林地、耕地、水域等大面積轉(zhuǎn)為建設(shè)用地，導(dǎo)致其生態(tài)系統(tǒng)服務(wù)價值降低，生物豐度指數(shù)降低，風(fēng)險發(fā)生時風(fēng)險受體不易受到影響，故生態(tài)風(fēng)險效應(yīng)值降低.

a 2005年 b 2014年圖7 平潭縣多風(fēng)險效應(yīng)分布圖

a 2005年 b 2014年圖8 平潭縣生態(tài)風(fēng)險分布圖

3.3.5 綜合生態(tài)風(fēng)險分布 為了綜合分析平潭縣生態(tài)風(fēng)險狀況，對多風(fēng)險源、多風(fēng)險受體和多風(fēng)險效應(yīng)3個指標(biāo)圖層進(jìn)行疊加分析，得出生態(tài)風(fēng)險分布圖，如圖8所示.2014年平均生態(tài)風(fēng)險值為0.403 3，較之于2005年，上升了0.03.從圖9可以看出，2005～2014年間大部分鄉(xiāng)鎮(zhèn)生態(tài)風(fēng)險平均值上升，只有南海鄉(xiāng)的生態(tài)風(fēng)險平均值降低，但其變化值不大.

3.3.6 生態(tài)風(fēng)險等級分布 參考相關(guān)生態(tài)風(fēng)險的分類方法[25]，將平潭縣劃分為5個等級，得到2005年與2014年生態(tài)風(fēng)險等級圖，如圖10所示.從縣域尺度上看，整個區(qū)域的極低生態(tài)風(fēng)險、低生態(tài)風(fēng)險、中等生態(tài)風(fēng)險等級的面積減少，高生態(tài)風(fēng)險面積增加.從鄉(xiāng)鎮(zhèn)尺度上看，與2005年相比，2014年平潭縣潭城鎮(zhèn)(含嵐城鄉(xiāng))、東庠鄉(xiāng)、蘇澳鎮(zhèn)、流水鎮(zhèn)、敖東鎮(zhèn)、嶼頭鄉(xiāng)由極低生態(tài)風(fēng)險、低生態(tài)風(fēng)險、中等生態(tài)風(fēng)險向高生態(tài)風(fēng)險轉(zhuǎn)移.

a 2005年 b 2014年圖10 平潭縣生態(tài)風(fēng)險等級圖

4 結(jié)束語

從多風(fēng)險源、多風(fēng)險受體、多風(fēng)險效應(yīng)3方面構(gòu)建土地利用變化下的平潭縣生態(tài)風(fēng)險評價指標(biāo)體系和模型，采用層次分析法計算指標(biāo)的權(quán)重值，利用綜合指數(shù)法計算生態(tài)風(fēng)險平均值，分析平潭縣2005～2014年的生態(tài)風(fēng)險變化.結(jié)果表明：1)2005年與2014年的土地利用分類精度分別為87.65%、88.57%，Kappa系數(shù)分別為0.845 0、0.858 7；2)2005～2014年平潭縣正處于大力發(fā)展的時期，填海造陸和海島開發(fā)建設(shè)等導(dǎo)致平潭土地利用類型發(fā)生著較大的變化，其中建設(shè)用地增加3 405.6 hm2，變化幅度最大；3)較之于2005年，2014年平潭縣總體平均生態(tài)風(fēng)險值上升0.03，大部分鄉(xiāng)鎮(zhèn)平均生態(tài)風(fēng)險值上升，極低生態(tài)風(fēng)險、低生態(tài)風(fēng)險、中等生態(tài)風(fēng)險等級的面積減少，高生態(tài)風(fēng)險面積增加.綜上所述，由于平潭縣人口增長及社會經(jīng)濟(jì)的快速發(fā)展，基礎(chǔ)設(shè)施建設(shè)和填海造陸工程等導(dǎo)致平潭土地利用發(fā)生較大變化，加之平潭縣作為兩岸交流合作的先行區(qū)、海峽西岸經(jīng)濟(jì)區(qū)先行先試的突破口，政策的導(dǎo)向必然促進(jìn)平潭快速發(fā)展，使得生態(tài)風(fēng)險總體呈上升趨勢，對生態(tài)環(huán)境存在較大威脅.為了保護(hù)平潭的生態(tài)環(huán)境，解決土地利用與生態(tài)環(huán)境之間的關(guān)系，制定合適的生態(tài)風(fēng)險管理對策：加強(qiáng)森林保護(hù)，積極培育幼林，嚴(yán)格控制森林采伐量，提高森林覆蓋率；嚴(yán)格控制建設(shè)用地的擴(kuò)張，防止城市無限制地蔓延；改變不合理的土地利用方式，建立土地利用安全評價指標(biāo)體系；堅持以人為本的區(qū)域生態(tài)風(fēng)險管理理念，建設(shè)可持續(xù)發(fā)展的生態(tài)城市.

區(qū)域生態(tài)風(fēng)險的定量化綜合評價為生態(tài)環(huán)境管理與社會經(jīng)濟(jì)可持續(xù)發(fā)展提供科學(xué)依據(jù)，同時為土地利用生態(tài)風(fēng)險指標(biāo)體系與模型的科學(xué)合理構(gòu)建提供參考，但本研究在構(gòu)建評價指標(biāo)體系時，對經(jīng)濟(jì)社會因素的考慮尚有不足，這將是下一步研究工作方向.

[1] Mooney H A, Duraiappah A, Larigauderie A. Evolution of natural and social science interactions in global change research programs [J]. PNAS, 2013, 110(S1): 3 665-3 672.

[2] Sterling S M, Ducharne A, Poleher J. The impact of global land-cover change on the terrestrial watercycle [J]. Nature Climate Change, 2012, 3(4): 385-390.

[3] Foley J A, DeFries R, Asner G P, et al. Global consequences of land use [J]. Science, 2005, 309(5734): 570-574.

[4] Liu J, Kuang W, Zhang Z, et al. Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s [J]. Journal of Geographical Sciences, 2014, 24(2): 195-210.

[5] 紀(jì)亞洲, 顧和和, 李保杰. 基于多層感知器神經(jīng)網(wǎng)絡(luò)的土地利用數(shù)據(jù)庫更新模型及應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報, 2015(7): 227-237.

[6] Hayes E H, Landis W G. Regional ecological risk assessment of a near shore marine environment: Cherry Point, WA [J]. Human and Ecological Risk Assessment, 2004, 10(2): 299-325.

[7] 傅伯杰, 陳利頂, 馬克明. 黃土丘陵區(qū)小流域土地利用變化對生態(tài)環(huán)境的影響——以延安市羊圈溝流域為例[J]. 地理學(xué)報, 1999, 54(3): 241-10.

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

[9] 周啟剛, 張曉媛, 王兆林. 基于正態(tài)云模型的三峽庫區(qū)土地利用生態(tài)風(fēng)險評價[J]. 農(nóng)業(yè)工程學(xué)報, 2014(23): 289-297.

[10] 趙巖潔, 李陽兵, 邵景安. 基于土地利用變化的三峽庫區(qū)小流域生態(tài)風(fēng)險評價——以草堂溪為例[J]. 自然資源學(xué)報, 2013, 28(6): 944-956.

[11] Obery A M,Landis W G. A regional multiple stressor risk assessment of the Codorus Creek watershed applying the relative risk model [J]. Human & Ecological Risk Assessment An International Journal, 2002, 8(2): 405-428.

[12] 傅麗華, 謝炳庚, 張曄, 等. 長株潭城市群核心區(qū)土地利用生態(tài)風(fēng)險評價[J]. 自然災(zāi)害學(xué)報, 2011, 20(2): 96-101.

[13] 孫洪波, 楊桂山, 蘇偉忠, 等. 沿江地區(qū)土地利用生態(tài)風(fēng)險評價——以長江三角洲南京地區(qū)為例[J]. 生態(tài)學(xué)報, 2010, 30(20): 5 616-5 625.

[14] Fried M A, Brodeley C E. Decision tree classification of land cover from remotely sensed data [J]. Remote Sensing of Environment, 1997, 61(3): 399-409.

[15] 楊子生. 滇東北山區(qū)土地墾殖指數(shù)與退耕還林規(guī)劃[J]. 山地學(xué)報, 1999, 17(S1): 45-48.

[16] 吳金華, 李紀(jì)偉, 梁晶晶. 土地利用程度與效益關(guān)系研究——以延安市為例[J]. 中國土地科學(xué), 2011, 25(8): 54-60.

[17] 劉紀(jì)遠(yuǎn). 中國資源環(huán)境遙感宏觀調(diào)查與動態(tài)研究[M]. 北京: 中國科學(xué)技術(shù)出版社, 1996: 171-179.

[18] McGarigal K, Marks B J. FRAGSTATS, Spatial Pattern Analysis Program for Quantifying Landscape Structure [M]. [S.l.]:[s.n.]: 1994.

[19] O’Neill R V, Krummel J R, Gardner R H, et al. Indices of landscape pattern [J]. Landscape Ecology, 1988, 1(3): 153-162.

[20] 國家環(huán)境保護(hù)部. HJ/T192—2006生態(tài)環(huán)境狀況評價技術(shù)規(guī)范(試行)[S]. 北京: 中國環(huán)境科學(xué)出版社, 2006.

[21] Costanza R, d’Arge R, de Groot R, et al. The value of the world’s ecosystem services and natural capital [J]. Science, 1997, 387: 253-260.

[22] 謝高地, 甄霖, 魯春霞, 等. 一個基于專家知識的生態(tài)系統(tǒng)服務(wù)價值化方法[J]. 自然資源學(xué)報, 2008, 23(5): 911-919.

[23] 陳鵬. 廈門濕地生態(tài)系統(tǒng)服務(wù)功能價值評估[J]. 濕地科學(xué), 2006, 4(2): 101-107.

[24] 朱國宇, 熊偉. 模糊評價法與綜合指數(shù)法在生態(tài)影響后評價中的應(yīng)用比較研究[J]. 東北農(nóng)業(yè)大學(xué)學(xué)報, 2011, 42(2): 54-59.

[25] 龔建周, 夏北成, 劉彥隨. 基于空間統(tǒng)計學(xué)方法的廣州市生態(tài)安全空間異質(zhì)性研究[J]. 生態(tài)學(xué)報, 2010, 30(20): 5 626-5 634.

Land Use Change and Ecological Risk Assessment in Pingtan from 2005 to 2014

Shi Tingting1, Xu Zhanghua1,2

(1. College of Environment and Resources, Fuzhou University, Fuzhou 350116, China；2. Fujian Provincial Key Laboratory of Resources and Environment Monitoring & Sustainable Management and Utilization, Fuzhou 350002, China)

In order to study the effects of land use change on regional ecological environment, the method of decision tree classification with visual interpretation was proposed for the land use change in Pingtan, the remote sensing images of Landsat 5 TM and Landsat 8 OLI were used as the main data sources; the multiple risk sources, multiple risk receptors and multiple risk effects were used to establish the ecological assessment index system and model, and the ecological risk changes from 2005 to 2014 in Pingtan were analyzed. The results indicated that: 1. the area of built-up land has increased by 3 405.6 hm2due to land reclamation and the island development from 2005 to 2014; 2. the overall average of ecological risk valued increase by three percent from 2014 to 2005, with an growing tendency in most of villages and towns; 3. the area with three risk levels, namely, very low ecological risk region, low ecological risk region, and moderate ecological risk region decreased, whereas the level of high ecological risk increased. The quantitative assessment of ecological risk in Pingtan establish a scientific foundation for the sustainable development of regional ecological environment and socio-economy, and provide a meaningful reference for scientific and reasonable construction of the ecological assessment index system and model.

regional ecological risk assessment; land use change; risk management; Pingtan

2016-04-20

國家自然科學(xué)基金青年科學(xué)基金項目(41501361，51408130)；福建省自然科學(xué)基金面上項目(2016J01188)；福建省資源環(huán)境監(jiān)測與可持續(xù)經(jīng)營利用重點實驗室開放基金項目(ZD1403)；福州大學(xué)人才基金項目(XRC-1345)

施婷婷(1993-)，女，福建平潭人，福州大學(xué)2015級碩士研究生，研究方向：資源環(huán)境遙感、城鄉(xiāng)規(guī)劃，E-mail: shitingting93@163.com

許章華(1985-)，男，福建福清人，博士，講師，碩導(dǎo)，研究方向：資源環(huán)境遙感、城鄉(xiāng)規(guī)劃與GIS應(yīng)用，E-mail: fafuxzh@163.com

1004-1729(2016)03-0278-11

P 951；TP 79

A DOl：10.15886/j.cnki.hdxbzkb.2016.0042