游巍斌，何東進，林 立，王 韌，蔡金標，王 鵬，張中瑞，肖石紅，鄭曉燕

(1.福建農(nóng)林大學(xué)林學(xué)院，福建福州350002;2.福鼎市林業(yè)局，福建福鼎355200;3.霞浦縣林業(yè)局，福建 霞浦355100;4.寧德市林業(yè)局，福建 寧德352100)

濕地、森林、海洋并稱為全球三大生態(tài)系統(tǒng)，它們在凈化環(huán)境、調(diào)蓄洪水、抵御自然災(zāi)害、保護生物多樣性、穩(wěn)定海岸等方面發(fā)揮著重要作用［1］.然而，全球人口的爆炸式增長、快速城市化和工農(nóng)業(yè)經(jīng)濟發(fā)展導(dǎo)致的盲目開墾、污染排放等人類行為給濕地生態(tài)系統(tǒng)健康造成巨大壓力，濕地退化問題愈發(fā)嚴重，探索濕地保護理論和方法已成為人類可持續(xù)發(fā)展的重大課題［2］.其中濕地生態(tài)脆弱性問題更是當(dāng)前國際濕地科學(xué)研究前沿領(lǐng)域關(guān)鍵問題之一.濱海濕地處于海陸生態(tài)系統(tǒng)頻繁相互作用的邊緣地帶，是典型的脆弱生態(tài)系統(tǒng)或生態(tài)敏感區(qū);城市邊緣帶上的濕地生態(tài)系統(tǒng)受城市化、工業(yè)化等人為干擾或破壞，生態(tài)系統(tǒng)脆弱性問題尤為突出，濱海濕地生態(tài)脆弱性的評價成為備受關(guān)注的研究熱點［3－6］.濕地脆弱性是指在自然環(huán)境和人為壓力下濕地退化的程度和可能性［7］.濕地脆弱性的定量評價是濕地脆弱性研究的核心.通過分析濕地生態(tài)脆弱性的特征及其影響因子評價濕地生態(tài)系統(tǒng)的生態(tài)脆弱性［6，8－13］，多屬于靜態(tài)評價.對濱海濕地生態(tài)系統(tǒng)脆弱性動態(tài)評價研究鮮見報道.

閩東地區(qū)兼具山地與濱海資源特色，且為我國秋茄紅樹林天然分布的北緣區(qū)，濱海濕地生態(tài)環(huán)境狀況特殊而復(fù)雜.在人類城市化過程以及臺風(fēng)、風(fēng)暴潮等諸多自然災(zāi)害的共同作用下，閩東天然濕地面積不斷萎縮，污染加劇，濕地退化問題日趨嚴重，濱海濕地生態(tài)系統(tǒng)日漸脆弱［14－16］.因此，通過構(gòu)建閩東濱海濕地生態(tài)系統(tǒng)脆弱性評價體系，對研究區(qū)2001、2010年生態(tài)系統(tǒng)脆弱性開展動態(tài)評價，以期為閩東濱海濕地的保護、恢復(fù)與重建提供科學(xué)依據(jù).

1 研究區(qū)概況

福建省閩東(寧德市)位于福建省東北部，南連省會福州，北接浙江溫州，西鄰武夷名山，東與臺灣島隔海相望.土地面積1.34萬km2，地形特點為丘陵山地和沿海小平原，年平均氣溫13.4－20.2℃，年平均降雨量1250－2350 mm，無霜期235－300 d.東臨臺灣海峽，海岸線長878 km，占全省海岸線總長的28.35%;海域面積4.45萬km2，占全省海洋漁場總面積的35.63%，區(qū)域內(nèi)有島、礁、沙、灘、岬角、水道、河口共1215個，濱海濕地資源豐富，水產(chǎn)豐饒，擁有海洋生物600多種.近年來，閩東濱海城市快速發(fā)展，大型工程修建和其他人為活動正急劇地改變著原始的濕地景觀，使得濱海濕地生態(tài)系統(tǒng)脆弱性不斷增大.

2 研究方法

2.1 濱海濕地分類體系的建立

基于對閩東濱海濕地資料的收集與濕地勘察，將1993、2001和2010年3個時期閩東沿海濕地的TM和ETM圖像進行幾何校正、圖片配準，再結(jié)合研究區(qū)1∶10000地形圖和實地考察資料，對遙感圖像進行監(jiān)督分類和人機交互解譯.使用ENVI軟件的主要和次要分析功能對解譯完成的圖像進行去除椒鹽噪聲等分類后處理，像元個數(shù)不足16個且明顯少于周圍濕地類型的類別歸為周圍濕地.處理完得到最終分類圖像，1993年、2001年、2010年的分類精度分別為90.00%、90.12%和89.58%.閩東濱海濕地劃分為自然濕地(水域、潮間裸灘)、人工濕地(農(nóng)田、人工庫塘、養(yǎng)殖場)、非濕地景觀(建筑用地、其他非濕地景觀)三大類和水域、潮間裸灘、農(nóng)田、人工庫塘、養(yǎng)殖場、建筑用地、其他非濕地景觀7個子類［14］.數(shù)據(jù)來源于3個時期閩東社會經(jīng)濟統(tǒng)計年鑒.

2.2 指標體系的建立

基于生態(tài)系統(tǒng)健康理論，采用由聯(lián)合國合作開發(fā)署建立的“壓力一狀態(tài)一響應(yīng)(PSR)”模型框架，對閩東濱海濕地生態(tài)系統(tǒng)的脆弱性進行評價.一般而言，生態(tài)脆弱性影響因子包括自然因素和人為因素兩大類，自然因素主要包括氣候、地形地貌、植被、土壤等;人為因素包括過度墾殖、采伐、放牧等生產(chǎn)活動和工農(nóng)業(yè)污染等方面.基于指標選取的目的性、整體性、數(shù)據(jù)易獲取性、代表性等基本原則［17，18］，利用PSR模型建立閩東濱海濕地生態(tài)脆弱性評價指標體系(表1)，采用層次分析法確定評價指標權(quán)重［19］.

2.3 單項指標的計算

2.3.1 壓力層指標 壓力層主要描述生態(tài)系統(tǒng)所承受的壓力大小，目前閩東濱海濕地生態(tài)系統(tǒng)的壓力主要來源于人口增長和對自然濕地的開發(fā)，故從土地壓力與人口壓力兩方面來反映濕地生態(tài)環(huán)境面臨的壓力［22］.選取土地墾殖率作為土地壓力的評價指標，該指標可以反映土地供人類生存所需的潛力.土地墾殖率與生態(tài)系統(tǒng)健康成負相關(guān)［22］，土地墾殖率=農(nóng)用地面積/研究區(qū)總面積.選擇人口密度和人口干擾指數(shù)作為人口壓力指標.人類干擾度指數(shù)采用不同期間與人為干擾密切相關(guān)的景觀類型(農(nóng)田、人工庫塘、養(yǎng)殖場、建筑用地)占研究區(qū)濕地總面積的比例來表示.

2.3.2 狀態(tài)層指標 參考前人研究［18，20］，采用景觀多樣性、平均斑塊面積、初級生產(chǎn)力、濕地平均彈性度和污染負荷指數(shù)表征濕地生態(tài)系統(tǒng)的現(xiàn)狀.系統(tǒng)活力可由初級生產(chǎn)力反映，初級生產(chǎn)力用植被歸一化指數(shù)(NDVI)度量.系統(tǒng)組織主要反映濕地生態(tài)系統(tǒng)的復(fù)雜性，選擇景觀多樣性指數(shù)(H)和斑塊平均面積指數(shù)(MPS)2個指標來度量.不同濕地景觀類型的彈性指標計算公式表示如下:

式中，ECOmax為各種景觀類型的平均彈性度，Pi為景觀類型i所占的比例，m為景觀類型的數(shù)目，Bi為第i類景觀類型彈性度分值.ECOmax越小，表明系統(tǒng)的生態(tài)彈性度越低，生態(tài)系統(tǒng)越脆弱.閩東濱海濕地不同景觀類型生態(tài)彈性度賦值［21］如下:水域 0.8，潮間裸灘 0.6，養(yǎng)殖場 0.4，建筑用地 0.3，農(nóng)田 0.5，人工庫塘 0.3，其他非濕地類型景觀 0.7.

采用自然濕地蓄水量和污染負荷指數(shù)2個指標表征濕地生態(tài)服務(wù)功能.自然濕地面積包括水域和潮間裸灘兩部分.污染負荷指數(shù)是以海灣的重金屬Pb含量來反映濕地凈化水質(zhì)的能力［16，21］.當(dāng)某地區(qū)污染物越多，污染物負荷指數(shù)越大，則濕地凈化污染物數(shù)量愈少，保持水質(zhì)能力愈差.

表1 閩東濱海濕地生態(tài)脆弱性評價指標體系Table 1 Evaluation index system for coastal wetland ecological vulnerability in eastern Fujian Province

2.3.3 響應(yīng)層指標 將濕地生態(tài)系統(tǒng)響應(yīng)層分為自然系統(tǒng)與人類系統(tǒng)兩部分.采用自然和人工濕地面積之和的變化反映自然對生態(tài)系統(tǒng)健康的響應(yīng).考慮到在正常情況下相對惡劣的生態(tài)環(huán)境無法供給、維持一個良好的人類社會，并且經(jīng)濟相對較弱，因而采用人均GDP值來反映人類系統(tǒng)對生態(tài)系統(tǒng)變化的響應(yīng).

2.4 評價模型的建立

采用各單項指標因子和濕地生態(tài)脆弱性綜合評估閩東濱海濕地生態(tài)脆弱性.

2.4.1 單因子評價模型 由于濕地生態(tài)系統(tǒng)并非一個平衡態(tài)的系統(tǒng)，濕地內(nèi)部各種反饋機制也不是單純的線性關(guān)系，而是在不斷地與外界進行能量、物質(zhì)交換的有序耗散結(jié)構(gòu)，因此，簡單的線性關(guān)系不能反映濕地生態(tài)系統(tǒng)中各方面和層次的狀態(tài).這里選用邏輯斯蒂增長曲線模型來對各個單項指標進行評價［20］.該模型可表示為:

式中，P表示單項指標的生態(tài)環(huán)境質(zhì)量指數(shù)評價值(無量綱值);R表示單項指標測度值(%).

對于常數(shù) a、b，采用以下方法確定:當(dāng) R=0.01時，P值近似為0.001;當(dāng) R=0.99時，P值近似為0.999，則此時方程中的a、b分別為4.595和9.19.鑒于指標變化對生態(tài)脆弱性影響一致性問題，評價指標量值增加與生態(tài)系統(tǒng)健康的增加方向相同時的單項指標計算公式為:

當(dāng)評價指標量值增加與生態(tài)系統(tǒng)健康的增加方向相反時，單項指標計算公式為:

式中，各指標含義同式(2).

2.4.2 多因子綜合評價模型 在各個單項指標值的基礎(chǔ)上，采用加權(quán)平均計算生態(tài)脆弱性綜合評價指數(shù)，其公式為:

式中:V為濕地生態(tài)脆弱性綜合評價指數(shù);Wi表示i個指標的權(quán)重值;Pi表示第i個單項指標的評價值.

2.5 脆弱性等級劃分

根據(jù)已有研究結(jié)果［17，18，20，21］按照綜合評價得分高低，從高到低排序，反映研究區(qū)的生態(tài)脆弱狀況.生態(tài)脆弱性程度分級和生態(tài)脆弱狀態(tài)含義表示如下.(1)一級(0≤V＜0.2)，潛在脆弱.濕地景觀結(jié)構(gòu)完整，狀態(tài)自然，活力強，功能水平正常，恢復(fù)力很強.(2)二級(0.2≤V＜0.4)，輕度脆弱.濕地景觀結(jié)構(gòu)較完整，狀態(tài)較自然，活力較強，功能水平正常，恢復(fù)力較強.(3)三級(0.4≤V＜0.6)，中度脆弱.濕地景觀結(jié)構(gòu)還算完整，自然狀態(tài)已受一定影響，活力衰退，功能水平有一定退化，對外界干擾的恢復(fù)力減弱.(4)四級(0.6≤V＜0.8)，重度脆弱.濕地景觀的結(jié)構(gòu)破碎，自然狀態(tài)受到極大破壞，活力很差，功能水平很大程度上退化，對外界干擾響應(yīng)迅速.(5)五級(0.8≤V≤1.0)，極度脆弱.濕地景觀的結(jié)構(gòu)完全破碎，自然狀態(tài)已經(jīng)徹底破壞，功能水平極低，對外界干擾響應(yīng)迅速.

3 結(jié)果與分析

3.1 閩東濱海濕地生態(tài)脆弱性單因子指標評價

當(dāng)外界壓力超過濕地生態(tài)系統(tǒng)所能承受范圍時，便會對濕地生態(tài)系統(tǒng)的結(jié)構(gòu)和功能造成損害，甚至導(dǎo)致濕地生態(tài)系統(tǒng)功能退化.閩東濱海濕地不同時期各指標值見表2，生態(tài)脆弱性單因子評價值與測度值見表3.1993－2010年，閩東濱海濕地生態(tài)脆弱性各單因子指標變化趨勢各異，土地墾殖率由3.66%上升至4.5%，上升了近0.9%;人口密度上升了0.89%;人類干擾面積由1993年占研究區(qū)面積的9.9%增至13.7%.景觀多樣性指數(shù)增加0.0240，表明近二十多年來閩東濱海濕地多樣性程度有所提高.平均斑塊面積由81.11 hm2增加至87.34 hm2.然而，濕地平均彈性度減少 0.0133，地表蓄水量減少 4.64%;污染負荷指數(shù)有所增加.1993－2010年人均GDP增長9.5倍的同時，濕地總面積減少了364.23 hm2，這與人口和社會經(jīng)濟的發(fā)展帶來的土地利用需求關(guān)系密切.

表2 不同時期閩東濱海濕地生態(tài)環(huán)境脆弱性評價指標值Table 2 Assessment index value of ecological environment vulnerability in coastal wetland of eastern Fujian Province in different periods

從項目層各指標變化來看，1993－2010年閩東濱海濕地壓力普遍增加，土地墾殖率不斷提高、人口密度增加，以及伴隨而來的人類干擾強度增強是閩東濱海濕地生態(tài)系統(tǒng)壓力的主要來源.閩東地區(qū)二十多年來共增加近34萬人，從人口的地區(qū)分布上看，人口匯聚在沿?？h.人口大量增加造成人類生存所需資源以及對環(huán)境條件的需求提高，對環(huán)境資源索取能力提高與生產(chǎn)生活相關(guān)的人為干擾活動(灘涂養(yǎng)殖，圍墾養(yǎng)殖，圍海造陸等)使得濕地生態(tài)系統(tǒng)面臨較大壓力，脆弱性程度提高.人類活動雖然小幅度提高了閩東濱海的景觀多樣性，但也導(dǎo)致自然濕地景觀破碎度增大;平均斑塊面積增大表明濕地生態(tài)系統(tǒng)有趨于集中分布的趨勢.

一般而言，生態(tài)系統(tǒng)的彈性表征系統(tǒng)對外界干擾的抵抗能力，彈性越高，生態(tài)系統(tǒng)越穩(wěn)定.閩東濱海濕地彈性度降低，表明濕地對外界干擾的抵抗能力逐漸下降.濕地蓄水量的降低和污染負荷指數(shù)的提高主要是由于研究區(qū)自然濕地逐漸向人工景觀轉(zhuǎn)移所致［14］.

3.2 閩東濱海濕地生態(tài)脆弱性綜合評價

基于濕地生態(tài)系統(tǒng)總體生態(tài)脆弱性，將閩東濱海濕地壓力、活力、組織、彈性、服務(wù)、自然系統(tǒng)響應(yīng)、人類系統(tǒng)響應(yīng)等生態(tài)脆弱性單因子指標加權(quán)后得到研究區(qū)生態(tài)脆弱性綜合評價指數(shù)(V).如表3所示，2010年和2001年閩東濱海濕地的生態(tài)脆弱性綜合評價指數(shù)(V)分別為0.5374和0.5384，均屬中度脆弱等級.此等級濕地生態(tài)特征表現(xiàn)為:濕地景觀結(jié)構(gòu)還算完整，但自然狀態(tài)已受一定影響，活力衰退，功能水平有一定退化，對外界干擾的恢復(fù)力減弱.進一步比較2個時期生態(tài)脆弱性綜合評價指數(shù)可知，2001－2010年研究區(qū)生態(tài)脆弱性指數(shù)值僅下降了0.001，脆弱性程度變化微小，說明此期間，雖然閩東濱海濕地壓力層各單因子指標均有所增加，而研究區(qū)脆弱性程度卻略有下降，這也從一個側(cè)面反映壓力因素雖然是濕地脆弱性增大的重要因素，然而，濕地總體脆弱性程度還需綜合考慮濕地系統(tǒng)固有結(jié)構(gòu)所維持的狀態(tài)及其相應(yīng)的自然和人類響應(yīng)因素.1993－2001年和1993－2010年各單因子指標變化方向和程度不同，而研究區(qū)綜合脆弱性程度變化程度小，這可能是由于濕地生態(tài)系統(tǒng)非線性作用機制的存在，各單因子指標的線性變化經(jīng)過濕地非線性系統(tǒng)內(nèi)多重反饋機制綜合作用所致.

表3 閩東濱海濕地生態(tài)環(huán)境脆弱性評價結(jié)果1)Table 3 Assessment results of coastal wetland ecological vulnerability in eastern Fujian Province

4 小結(jié)與討論

基于生態(tài)系統(tǒng)健康理論與景觀生態(tài)學(xué)理論，運用邏輯斯蒂模型對閩東濱海濕地2001年和2010年2個時期濕地生態(tài)系統(tǒng)脆弱性進行了定量評價.結(jié)果表明:2001年和2010年閩東濱海濕地生態(tài)系統(tǒng)生態(tài)脆弱性綜合評價指數(shù)(V)分別為0.5384和0.5374，均處中度脆弱水平，閩東濱海濕地景觀結(jié)構(gòu)還算完整，但自然狀態(tài)已受一定影響，活力衰退，功能水平有一定退化，對外界干擾的恢復(fù)力減弱.1993－2010年各指標測度值有所差異，但由于濕地生態(tài)系統(tǒng)非線性反饋機制的綜合作用，此期間濕地脆弱性變化程度較小.

當(dāng)前濕地生態(tài)系統(tǒng)脆弱性評價主要集中在關(guān)于“時點”的靜態(tài)研究方面，這種基于脆弱性現(xiàn)狀特征的靜態(tài)評價，難以反映濕地脆弱性的動態(tài)變化趨勢.為此，本文以1993年的濕地生態(tài)系統(tǒng)狀態(tài)為基期(此時期閩東城市開發(fā)水平較低)，通過評價2001年與2010年不同時期的脆弱性狀態(tài)，以濕地脆弱性“時段”特征實現(xiàn)動態(tài)評估.“構(gòu)建指標體系—確定評價單元—確定評價方法—劃分濕地脆弱性等級—綜合分析”濕地脆弱性評價程序基本得到認可和統(tǒng)一［7］，但由于研究資料、研究技術(shù)、評價對象的時空尺度性、濕地系統(tǒng)自身復(fù)雜性及自然和人文影響因素相互作用的復(fù)雜性，濕地評價中的不確定性普遍存在.這種不確定性問題在本文中主要表現(xiàn)在兩方面.一方面，生態(tài)脆弱性等級采用等距劃分，主觀性較強;同一級別內(nèi)的不同類別在評分的時候沒有顯示出區(qū)別而產(chǎn)生的不確定性，或不同級別的分界處的評分過分拉大也會產(chǎn)生不確定性.目前降低不確定性的方法有隨機數(shù)學(xué)方法、模糊數(shù)學(xué)及灰色系統(tǒng)方法［23］，但現(xiàn)有的有限的零散濕地監(jiān)測數(shù)據(jù)，使得上述方法受到一定限制.另一方面，在計算各單因子評價指標時，需要根據(jù)研究區(qū)的實際情況，確定測度值增量對生態(tài)系統(tǒng)脆弱性或健康程度的影響方向.如1993－2010年單位面積NDVI值上升了1%，實地調(diào)查發(fā)現(xiàn)，這種上升主要是因為研究區(qū)原本裸露的潮間灘涂和淺海水域圍墾成為農(nóng)田.而對于濱海濕地生態(tài)系統(tǒng)而言，潮間裸灘的影響大于農(nóng)田，故這種變化對閩東濱海濕穩(wěn)定性的維持反而不利.同時，寧德的飛鸞、漳灣和霞浦等地，均有不同面積的互花米草入侵，導(dǎo)致NDVI值提高，然而，互花米草的生物入侵對于濱海濕地生態(tài)系統(tǒng)來說是極為不利的因素.因而，濕地脆弱性評價過程中，不論是指標體系構(gòu)建還是評價模型計算都應(yīng)保持對不確定性問題的關(guān)注，這也將是今后濕地脆弱性評價研究的難點和重要突破口.

［1］崔保山，楊志峰.濕地生態(tài)系統(tǒng)健康評價指標體系Ⅱ.方法與案例［J］.生態(tài)學(xué)報，2002，22(8):1231－1239.

［2］竇勇，唐學(xué)璽，王悠.濱海濕地生態(tài)修復(fù)研究進展［J］.海洋環(huán)境科學(xué)，2012，31(4):616－620.

［3］劉世棟，高峻.旅游開發(fā)對上海濱海濕地植被的影響［J］.生態(tài)學(xué)報，2012，32(10):2992－3000.

［4］NHUAN M T，NGOC N T M，HUONG N Q，et al.Assessment of vietnam coastal wetland vulnerability for sustainable use(case study in Xuanthuy Ramsar Site，Vietnam)［J］.Journal of Wetlands Ecology，2009，2(1－2):1－16.

［5］MOSER S C，JEFFRESS WILLIAMS S，BOESCH D F.Wicked challenges at land's end:managing coastal vulnerability under climate change［J］.Annual Review of Environment and Resources，2012，37:51－78.

［6］付博，姜琦剛，任春穎.扎龍濕地生態(tài)脆弱性評價與分析［J］.干旱區(qū)資源與環(huán)境，2011，25(1):49－52.

［7］擺萬奇，尚二萍.濕地脆弱性評價研究進展［J］.濕地科學(xué)，2012，10(3):378－384.

［8］萬忠娟，于少鵬，王海霞，等.松嫩平原典型濕地脆弱性分析與評價［J］.東北師大學(xué)報:自然科學(xué)版，2003，35(2):93－99.

［9］陳久和.城市邊緣濕地生態(tài)環(huán)境脆弱性研究——以杭州西溪濕地為例［J］.科技通報，2003，19(5):395－398.

［10］郭躍東，何巖，鄧偉，等.扎龍河濱濕地水系統(tǒng)脆弱性特征及影響因素分析［J］.濕地科學(xué)，2004，2(1):47－53.

［11］廖富強，劉影，葉慕亞，等.鄱陽湖典型濕地生態(tài)環(huán)境脆弱性評價及壓力分析［J］.長江流域資源與環(huán)境，2008，17(1):133－137.

［12］楊志榮，吳次芳，劉勇.快速城市化地區(qū)濕地脆弱性評價——以杭州灣濕地為例［J］.資源科學(xué)，2009，31(2):243－249.

［13］劉振乾，劉紅玉，呂憲國.三江平原濕地生態(tài)脆弱性研究［J］.應(yīng)用生態(tài)學(xué)報，2001，12(2):241－244.

［14］林立，何東進，王韌，等.閩東濱海濕地景觀分類體系與格局特征［J］.西南林學(xué)院學(xué)報，2012，32(2):62－67.

［15］何東進，鄭開基，王韌，等.閩東濱海濕地不同起源秋茄林重金屬元素Zn、Cd、Cu的累積與分布特征比較［J］.福建農(nóng)林大學(xué)學(xué)報:自然科學(xué)版，2012，41(2):187－192.

［16］何東進，鄭開基，王韌，等.閩東濕地重金屬污染特征與生態(tài)風(fēng)險評價［J］.武夷科學(xué)，2009，25:44－49.

［17］左偉，周慧珍，王橋.區(qū)域生態(tài)安全評價指標體系選取的概念框架研究［J］.土壤，2003，35(1):2－7.

［18］吳姍姍.萊州灣南岸濱海濕地的景觀格局變化及其生態(tài)脆弱性評價［D］.濟南:山東師范大學(xué)，2009.

［19］朱衛(wèi)紅，郭艷麗，孫鵬，等.圖們江下游濕地生態(tài)系統(tǒng)健康評價［J］.生態(tài)學(xué)報，2012，32(21):6609－6618.

［20］李永建.拉魯濕地生態(tài)環(huán)境質(zhì)量評價的景觀生態(tài)學(xué)方法應(yīng)用研究［D］.成都:四川大學(xué)，2002.

［21］高吉喜.可持續(xù)發(fā)展理論探索——生態(tài)承載力理論、方法與應(yīng)用［M］.北京:中國環(huán)境科學(xué)出版社，2001.

［22］黃方，劉湘南，張養(yǎng)貞.GIS支持下的吉林省西部生態(tài)環(huán)境脆弱態(tài)勢評價研究［J］.地理科學(xué)，2003，23(1):95－100.

［23］周建飛，曾光明，黃國和，等.基于不確定性的城市擴展用地生態(tài)適宜性評價［J］.生態(tài)學(xué)報，2007，2(2):774－782.