吉林省西部地區(qū)生態(tài)脆弱性的動(dòng)態(tài)評(píng)價(jià)

田 超, 楊夏捷, 鞠園華, 郭福濤

(1.福建農(nóng)林大學(xué) 林學(xué)院, 福州 350002; 2.海峽兩岸紅壤區(qū)水土保持協(xié)同創(chuàng)新中心, 福州 350002)

隨著社會(huì)經(jīng)濟(jì)的快速發(fā)展和人類(lèi)活動(dòng)的不斷加劇，對(duì)生態(tài)環(huán)境的干擾日益頻繁，導(dǎo)致生態(tài)脆弱性問(wèn)題尤為突出。目前，有關(guān)生態(tài)環(huán)境脆弱性的研究已引起了諸多學(xué)者的廣泛關(guān)注，并成為探討全球變化和生態(tài)環(huán)境可持續(xù)發(fā)展的重要切入點(diǎn)[1-3]。目前國(guó)內(nèi)外學(xué)者已開(kāi)展了有關(guān)生態(tài)脆弱性的實(shí)踐性研究[4-6]、生態(tài)脆弱性的動(dòng)態(tài)評(píng)價(jià)和預(yù)測(cè)預(yù)報(bào)[7-9]等，這些研究有助于全面了解區(qū)域生態(tài)環(huán)境質(zhì)量狀況，并為生態(tài)環(huán)境恢復(fù)和重建工作提供科學(xué)依據(jù)[10-11]。

吉林省西部地區(qū)是我國(guó)典型的生態(tài)脆弱區(qū)。區(qū)內(nèi)土地資源的開(kāi)發(fā)歷史雖不長(zhǎng)，但由于長(zhǎng)期的過(guò)度放牧、墾殖，不合理的土地利用，加之惡劣的自然條件的影響，生態(tài)脆弱性十分明顯。2000年以來(lái)吉林省政府生態(tài)環(huán)境保護(hù)政策的陸續(xù)出臺(tái)，加強(qiáng)了對(duì)生態(tài)問(wèn)題嚴(yán)重區(qū)的治理力度，西部地區(qū)的生態(tài)環(huán)境不斷改善[12]。然而，目前還缺少對(duì)吉林省西部地區(qū)，尤其是2000年以后的生態(tài)環(huán)境脆弱性長(zhǎng)時(shí)間尺度動(dòng)態(tài)變化研究。鑒于此，本研究選取降水、溫度、人口密度和土地利用類(lèi)型等10個(gè)因子，基于ArcGIS 10.2，SPSS 19.0等技術(shù)軟件，利用主成分分析法和綜合指數(shù)法，對(duì)吉林省西部地區(qū)2000年、2006年、2012年的生態(tài)脆弱度進(jìn)行定量評(píng)價(jià)，明確生態(tài)脆弱性的主導(dǎo)因素，研究結(jié)論為吉林省生態(tài)建設(shè)和生態(tài)恢復(fù)提供科學(xué)依據(jù)。

1 研究區(qū)概況

吉林省西部地區(qū)位于松嫩平原的西南部，地跨東經(jīng)123°09′—124°22′，北緯44°57′—45°46′，其研究范圍包括白城及松原市所管轄的9個(gè)縣市(區(qū))，國(guó)土總面積為46 000 km2，約占吉林省總面積的1/4。研究區(qū)內(nèi)地形平緩，海拔為120～180 m。氣候類(lèi)型屬半干旱半濕潤(rùn)的大陸性季風(fēng)氣候，四季分明，雨熱同期。夏季炎熱多雨，冬季寒冷漫長(zhǎng)。熱量資源豐富，年平均氣溫5.2℃。平均降水量約為400 mm，年平均蒸發(fā)量1 500～1 900 mm。區(qū)內(nèi)植被類(lèi)型以針葉混交林、落葉松林為主。土壤類(lèi)型主要有黑鈣土、草甸土、風(fēng)沙土等。2012年末，吉林省吉林西部各縣(市)人均國(guó)內(nèi)生產(chǎn)總值為39.911 億元，年末實(shí)有耕地面積為296.157 hm2，人口總數(shù)為624.34萬(wàn)人[13]。

2 研究方法

2.1 評(píng)價(jià)指標(biāo)的選擇

指標(biāo)是進(jìn)行生態(tài)環(huán)境評(píng)價(jià)的基本尺度和衡量標(biāo)準(zhǔn)[14]。引發(fā)區(qū)域生態(tài)環(huán)境脆弱性的原因主要包括自然因素和人為因素兩個(gè)方面。在遵循科學(xué)性、全面性及數(shù)據(jù)可獲取性等原則的基礎(chǔ)上，結(jié)合前人研究成果，共選取了10個(gè)指標(biāo)(表1)。具體指標(biāo)選取原則如下：自然因素是生態(tài)系統(tǒng)各種因素長(zhǎng)期作用的結(jié)果，也是生態(tài)環(huán)境構(gòu)成的基礎(chǔ)[15]?；谇叭搜芯?，本文選取了“年平均降水量”、“年平均相對(duì)濕度”、“年平均溫度”反映氣候特征變化[16]，選取“坡度”、“高程”因子反映地形地貌表征[17]，選取“植被覆蓋度”反映植被的生長(zhǎng)狀況[18]。人類(lèi)活動(dòng)因子的選取對(duì)生態(tài)環(huán)境脆弱度有重要影響。前人研究表明，“土地利用類(lèi)型”直接影響人類(lèi)的生存環(huán)境，每種土地利用類(lèi)型在一定程度上承載不同的生態(tài)利用價(jià)值，是決定生態(tài)脆弱程度的關(guān)鍵因素[19]。姚昆等[20]研究表明，“人均GDP”、“人口密度”、“耕地面積”是生態(tài)脆弱性的潛在因子，可對(duì)生態(tài)系統(tǒng)的穩(wěn)定性造成破壞[21,14]。因此，本研究選取以上10個(gè)指標(biāo)對(duì)吉林省西部2000年、2006年、2012年的生態(tài)脆弱性進(jìn)行評(píng)價(jià)。

本研究以柵格為基本評(píng)價(jià)單元，所有數(shù)據(jù)均采用Albers投影，空間數(shù)據(jù)柵格大小為100 m×100 m。研究中高程、坡度數(shù)據(jù)由DEM提取，氣象數(shù)據(jù)包括“年平均降水”、“年平均溫度”和“年平均濕度”，取吉林省西部9個(gè)市縣及周邊地區(qū)氣象站點(diǎn)數(shù)據(jù)。“人口密度”、“人均GDP”、“年末實(shí)有耕地面積”來(lái)源于吉林省2000年、2006年、2012年統(tǒng)計(jì)年鑒，根據(jù)吉林省西部地區(qū)各市縣的社會(huì)經(jīng)濟(jì)數(shù)據(jù)，基于ArcGIS 10.2地統(tǒng)計(jì)分析中的克里金插值法進(jìn)行空間確定性插值，實(shí)現(xiàn)其空間化。

表1 生態(tài)脆弱性評(píng)價(jià)指標(biāo)描述及來(lái)源

2.2 有效性檢驗(yàn)

有效性檢驗(yàn)是所選取的指標(biāo)是否適用于該方法進(jìn)行檢驗(yàn)。因此，在評(píng)價(jià)前，需要驗(yàn)證相關(guān)矩陣中各變量的相關(guān)性，即各個(gè)變量間是否相互獨(dú)立。借助SPSS 19.0的因子分析(Bartlett′s球形檢驗(yàn))對(duì)各指標(biāo)進(jìn)行檢驗(yàn)，若KMO統(tǒng)計(jì)量在0.5以上，說(shuō)明適合應(yīng)用因子分析法。若Bartlett′s球形檢驗(yàn)結(jié)果變量拒絕原假設(shè)，說(shuō)明生態(tài)環(huán)境指標(biāo)間存在一定的相關(guān)性，因子分析有效[22]。由于指標(biāo)的生成和提取需要從研究區(qū)內(nèi)選取一定數(shù)量的隨機(jī)點(diǎn)作為樣點(diǎn)進(jìn)行研究，在參考其他文獻(xiàn)的基礎(chǔ)上[19]，基于ArcGIS 10.2隨機(jī)生成研究區(qū)900個(gè)點(diǎn)，根據(jù)這些點(diǎn)提取各指標(biāo)所需數(shù)據(jù)。

根據(jù)檢驗(yàn)結(jié)果顯示，2000年、2006年、2012年Bartlett′s的卡方值為分別為5 742.838，42 313.103，6 494.726，且顯著性均Sig.<0.05(即p值<0.05)，故應(yīng)拒絕原假設(shè)，認(rèn)為相關(guān)系數(shù)矩陣與單位陣有顯著差異。KMO>0.5，根據(jù)Kaiser給出的KMO度量標(biāo)準(zhǔn)可知原有變量適合進(jìn)行因子分析。由此可見(jiàn)，所選取的生態(tài)環(huán)境指標(biāo)適宜作因子分析。

2.3 數(shù)據(jù)處理

2.3.1 評(píng)價(jià)指標(biāo)標(biāo)準(zhǔn)化 生態(tài)評(píng)價(jià)指標(biāo)具有不同的含義，因?yàn)樵跀?shù)量級(jí)和量綱上都不同[23]，為保證數(shù)據(jù)的客觀性和科學(xué)性，需要按照一定的標(biāo)準(zhǔn)對(duì)參評(píng)因子進(jìn)行標(biāo)準(zhǔn)化處理。由于評(píng)價(jià)指標(biāo)與生態(tài)環(huán)境脆弱性關(guān)系有正逆兩種，其算法也不同。標(biāo)準(zhǔn)化采用以下兩種公式：

(1)

(2)

式中：Xi′為第i個(gè)指標(biāo)的標(biāo)準(zhǔn)化值；Xi表示第i個(gè)指標(biāo)的初始值；Xmax，Xmin分別表示第i個(gè)指標(biāo)的最大值和最小值。Xi′越大，表明生態(tài)環(huán)境越脆弱，越容易受到外界干擾而破壞。其中，定量正向指標(biāo)包括：“人口密度”、“高程”、“坡度”、“年均溫度”、“年末實(shí)有耕地面積”。定量負(fù)向指標(biāo)包括：“人均GDP”、“年均降水”、“年均濕度”、“土壤類(lèi)型”、“植被覆蓋度”。當(dāng)指標(biāo)因子值與生態(tài)環(huán)境呈正相關(guān)時(shí)使用公式(1)；當(dāng)指標(biāo)因子值與生態(tài)環(huán)境呈負(fù)相關(guān)時(shí)使用公式(2)。

定性指標(biāo)中的土地利用類(lèi)型通過(guò)分等級(jí)賦值法實(shí)現(xiàn)指標(biāo)的定量處理，最后進(jìn)行極差標(biāo)準(zhǔn)化計(jì)算(表2)[4]。

表2 定性指標(biāo)賦值標(biāo)準(zhǔn)

2.3.2 貢獻(xiàn)率及累計(jì)貢獻(xiàn)率 第h個(gè)主成分累計(jì)貢獻(xiàn)率達(dá)到70%～90%時(shí)[24]，取前q個(gè)主成分(表3)。貢獻(xiàn)率公式為：

(3)

2.3.3 特征值和特征向量 主成分因子負(fù)荷是主成分因子與原始變量因子間的相關(guān)系數(shù)。主成分的特征向量，又稱(chēng)系數(shù)向量，等于對(duì)應(yīng)的載荷系數(shù)除以特征值的平方根[14]。

2.3.4 綜合得分模型 各主成分的系數(shù)分別乘上所對(duì)應(yīng)的貢獻(xiàn)率再除以4個(gè)主成分貢獻(xiàn)率之和，最后所有主成分每個(gè)系數(shù)前的值相加，即可得到綜合得分模型：

Fi=β1X1+β2X2+…+βmXm

(4)

綜合得分模型中的各指標(biāo)前的系數(shù)(β1—βm)即為各指標(biāo)因子得分值。

表3 各成分的特征值、貢獻(xiàn)率和累計(jì)貢獻(xiàn)率

2.3.5 確定權(quán)重 各指標(biāo)的權(quán)重等于公式(4)中每個(gè)指標(biāo)的得分值除以指標(biāo)的總得分。

2.3.6 生態(tài)環(huán)境脆弱度 生態(tài)環(huán)境脆弱度是用于研究研究區(qū)內(nèi)生態(tài)環(huán)境狀況和表征生態(tài)環(huán)境脆弱度一個(gè)重要標(biāo)準(zhǔn)[11]。通過(guò)多因子綜合分析，結(jié)合研究區(qū)內(nèi)自然因素、社會(huì)經(jīng)濟(jì)因子等方面構(gòu)建生態(tài)脆弱綜合指數(shù)(Ecological Vulnerable Degree)。用公式(5)表達(dá)吉林省西部地區(qū)生態(tài)脆弱度：

(5)

式中：EVDi為第i個(gè)評(píng)價(jià)區(qū)域生態(tài)脆弱度，其值為(0，1)；Pij為第i個(gè)區(qū)域j個(gè)指標(biāo)；wj為各指標(biāo)權(quán)重值。EVDi值越高，脆弱程度越嚴(yán)重。

研究區(qū)生態(tài)脆弱等級(jí)的劃分，對(duì)全面認(rèn)識(shí)區(qū)內(nèi)生態(tài)環(huán)境狀況具有重要意義[25]。結(jié)合國(guó)內(nèi)外研究現(xiàn)狀[7,14]，并綜合吉林省西部地區(qū)自然環(huán)境和生態(tài)脆弱度的表現(xiàn)特征，按照生態(tài)脆弱度的大小將研究區(qū)生態(tài)環(huán)境脆弱度劃分為5個(gè)等級(jí)(表4)。

2.3.7 生態(tài)脆弱度變化趨勢(shì) 使用定量的生態(tài)環(huán)境脆弱性綜合性指數(shù)(EVSI)能夠更加直觀全面地反映研究區(qū)生態(tài)環(huán)境質(zhì)量，而對(duì)研究區(qū)生態(tài)脆弱度變化趨勢(shì)可根據(jù)變化率進(jìn)行描述[18]，計(jì)算公式如下：

(6)

(7)

式中：EVSI為生態(tài)脆弱綜合指數(shù)；Pi為第i類(lèi)脆弱性等級(jí)值；Ai為第i類(lèi)脆弱性面積；S為研究區(qū)總面積(km2)；R為綜合指數(shù)變化率；Aim和Ain分別為m，n時(shí)間，第i類(lèi)脆弱性面積。若R<0，則生態(tài)環(huán)境呈恢復(fù)趨勢(shì)，若R>0，則生態(tài)環(huán)境呈惡化趨勢(shì)，不容樂(lè)觀。

3 結(jié)果與分析

3.1 研究區(qū)生態(tài)脆弱性評(píng)價(jià)

表5顯示，2000年各指標(biāo)權(quán)重排序?yàn)椋耗昶骄鶜鉁?年平均降水>年平均濕度>人口密度>年末實(shí)有耕地面積>植被覆蓋度>坡度>土地利用類(lèi)型>人均GDP>高程，2006年指標(biāo)權(quán)重排序：年平均降水>土地利用類(lèi)型>年平均濕度>人均GDP>年末實(shí)有耕地面積>年平均氣溫>植被覆蓋度>坡度>人口密度>高程，2012年各指標(biāo)權(quán)重排序?yàn)椋和恋乩妙?lèi)型>年平均降水>年平均濕度>年平均氣溫>人口密度>年末實(shí)有耕地面積>植被覆蓋度>人均GDP>坡度>高程。

此外，根據(jù)表5可知，2000年氣象因子所占比重較大，說(shuō)明2000年吉林省西部地區(qū)的生態(tài)脆弱度主要受氣象因素的影響。2006年和2012年，年平均降水和土地利用類(lèi)型在2個(gè)年份中權(quán)重最大，分別為0.216，0.204。說(shuō)明2006—2012年吉林省西部生態(tài)環(huán)境脆弱度主要受氣象因素和人為活動(dòng)共同影響。

表5 評(píng)價(jià)指標(biāo)權(quán)重

3.2 生態(tài)脆弱性空間分析

根據(jù)生態(tài)脆弱度分類(lèi)標(biāo)準(zhǔn)，基于ArcGIS 10.2軟件中的重分類(lèi)函數(shù)對(duì)生態(tài)脆弱圖進(jìn)行分類(lèi)，從而得到吉林省西部地區(qū)生態(tài)脆弱性分級(jí)表(表6)。并根據(jù)2000年、2006年、2012年吉林省西部地區(qū)生態(tài)脆弱度可對(duì)其進(jìn)行半方差函數(shù)結(jié)構(gòu)和克里格差值，基于ArcGIS 10.2區(qū)域分析中的分區(qū)統(tǒng)計(jì)功能可得到不同地區(qū)生態(tài)脆弱度的空間平均分布(圖1)。

表6 吉林省西部生態(tài)脆弱性分級(jí)

根據(jù)表6可知，2000年、2006年、2012年生態(tài)脆弱性指數(shù)平均值分別為：0.600±0.129，0.531 9±0.111，0.506±0.118，整體處于中度脆弱水平。2006年與2000年相比，吉林省西部地區(qū)輕度脆弱的土地面積增長(zhǎng)較快，增加了15.362%，中度脆弱的面積增長(zhǎng)不是很大，增加了0.621%，而“微度脆弱”、“強(qiáng)度脆弱”和“極強(qiáng)度脆弱”都有所減少，分別降低了4.764%，0.142%，11.075%。2012年較2006年相比，吉林省西部地區(qū)“輕度脆弱”土地面積增長(zhǎng)最快，增長(zhǎng)了18.281%，“強(qiáng)度脆弱”面積有較小幅度的增加，增長(zhǎng)了0.508%，而“輕度脆弱”、“中度脆弱”、“極強(qiáng)度脆弱”都有不同程度的減少，分別降低了10.892%，2.52%，5.376%。吉林省西部地區(qū)生態(tài)環(huán)境脆弱度由2000年以“輕度脆弱”和“極強(qiáng)度脆弱”為主變?yōu)?006年以“輕度脆弱”為主，再變?yōu)?012年以“微度脆弱”和“輕度脆弱”為主的趨勢(shì)。由圖1可知，2000—2012年研究區(qū)生態(tài)脆弱性空間分布特征明顯，呈現(xiàn)自西向東降低的變化特征。整個(gè)時(shí)間段里，“強(qiáng)度脆弱”和“極強(qiáng)度脆弱”主要分布在研究區(qū)西部包括洮南市、白城市和通榆縣等。該地區(qū)為與研究區(qū)西部，人口密度大、氣候干燥、鹽堿化較為嚴(yán)重，自然環(huán)境較為惡劣。這些自然條件和人為因素特征導(dǎo)致該地區(qū)的生態(tài)自我修復(fù)能力較差，生態(tài)問(wèn)題顯著?！爸卸却嗳酢焙汀拜p度脆弱”主要分布在研究區(qū)中部，隨著年份的增長(zhǎng)逐漸向西部擴(kuò)展?！拔⒍却嗳酢敝饕植荚谘芯繀^(qū)東部包括前郭縣、長(zhǎng)嶺縣等地區(qū)。該區(qū)域降水量相對(duì)較多，植被覆蓋度較高，土地鹽堿化程度較輕，經(jīng)濟(jì)水平較高，加之近些年來(lái)的生態(tài)環(huán)境保護(hù)措施的實(shí)施和治理力度的加大，對(duì)生態(tài)環(huán)境起到了良好的作用，沙漠化得到基本控制。這些地區(qū)抗干擾和生態(tài)恢復(fù)能力較強(qiáng)，生態(tài)環(huán)境狀況比較樂(lè)觀。

圖1 吉林省西部地區(qū)生態(tài)環(huán)境脆弱性評(píng)價(jià)

3.3 生態(tài)脆弱度變化趨勢(shì)分析

根據(jù)公式(6)—(7)，可得到及吉林省西部地區(qū)的生態(tài)環(huán)境脆弱性綜合指數(shù)和變化率(表7)。由表7可知，吉林省2000年、2006年、2012年生態(tài)環(huán)境脆弱指數(shù)分別為2.891，2.610，2.251。在由此可見(jiàn)，吉林省西部地區(qū)生態(tài)脆弱指數(shù)在整個(gè)時(shí)段內(nèi)呈下降趨勢(shì)，可能由于“強(qiáng)度脆弱”和“極強(qiáng)度脆弱”的綜合指數(shù)均呈不斷下降趨勢(shì)，在一定程度上影響了吉林省西部地區(qū)生態(tài)環(huán)境狀況的變化。

表7 吉林省西部生態(tài)環(huán)境脆弱性綜合指數(shù)變化

由吉林省西部生態(tài)脆弱度的變化率可知，2000—2006年、2006—2012年、2000—2012年的R值均小于0，說(shuō)明生態(tài)環(huán)境呈不斷改善趨勢(shì)。在一定程度上，自然條件的改善和人為活動(dòng)的限制對(duì)研究區(qū)生態(tài)環(huán)境產(chǎn)生了良好影響，加之2000年開(kāi)始吉林省政府加大對(duì)西部地區(qū)生態(tài)環(huán)境建設(shè)和保護(hù)，在對(duì)江河洪水調(diào)蓄區(qū)、防風(fēng)固沙區(qū)、水土流失保護(hù)區(qū)等實(shí)施了生態(tài)治理和保護(hù)項(xiàng)目[12]，退耕還林、還草等措施已初見(jiàn)成效。在一定程度上減少了生態(tài)環(huán)境壓力緩解了生態(tài)環(huán)境的脆弱程度。綜上所述，吉林省西部十多年來(lái)生態(tài)環(huán)境狀態(tài)發(fā)展較好，但生態(tài)環(huán)境保護(hù)措施和力度仍需提高和加強(qiáng)。

4 討論與結(jié)論

(1) 吉林省西部地區(qū)生態(tài)環(huán)境脆弱性的空間分布差異明顯，自西向東逐漸降低。整體上，研究區(qū)的東部和中部地區(qū)的生態(tài)環(huán)境質(zhì)量較好，生態(tài)環(huán)境脆弱性較輕。而西部地區(qū)由于人類(lèi)干擾頻繁導(dǎo)致生態(tài)脆弱性較高。生態(tài)脆弱性綜合指數(shù)(EVSI)表明，2000—2012年吉林省西部地區(qū)生態(tài)環(huán)境脆弱性不斷減弱，生態(tài)環(huán)境呈不斷改善趨勢(shì)。這與麻素挺等[26]研究相一致。

(2) 生態(tài)指標(biāo)的權(quán)重排名表明，2000—2012年吉林省西部地區(qū)生態(tài)脆弱變化的影響因子由以氣象因素為主變?yōu)橐酝恋乩脿顩r為主。近10 a間，由于吉林省西部地區(qū)的耕地面積大幅度增加，草地面積大幅度減少，人工建設(shè)用地規(guī)模不斷增長(zhǎng)，加劇了土地利用變化[27]。

本研究主要基于主成分分析方法進(jìn)行對(duì)生態(tài)脆弱性進(jìn)行評(píng)價(jià)?；谥鞒煞址治龇椒軌?qū)崿F(xiàn)對(duì)原始數(shù)據(jù)的降維，并考慮了時(shí)間序列的變化[28]，從而實(shí)現(xiàn)了生態(tài)脆弱度動(dòng)態(tài)評(píng)價(jià)，保證了數(shù)據(jù)的統(tǒng)一性和可比性，從而為區(qū)域生態(tài)環(huán)境評(píng)價(jià)提供了科學(xué)依據(jù)。但在構(gòu)建指標(biāo)體系時(shí)，較多考慮數(shù)據(jù)資料的可獲得性，對(duì)其他影響生態(tài)環(huán)境的重要因子如土壤農(nóng)藥用量、水資源開(kāi)發(fā)度等較難獲取的指標(biāo)并未納入生態(tài)環(huán)境評(píng)價(jià)體系中，評(píng)價(jià)指標(biāo)如人類(lèi)社會(huì)經(jīng)濟(jì)活動(dòng)的量化比較困難，有一定的局限性，無(wú)疑會(huì)對(duì)評(píng)價(jià)結(jié)果造成一定的影響，對(duì)今后的評(píng)價(jià)指標(biāo)的選擇和等級(jí)的劃分具有一定的主觀性。因此，對(duì)于指標(biāo)的合理選取和評(píng)價(jià)體系的制定，需后續(xù)研究進(jìn)一步完善。

參考文獻(xiàn)：

[1] Ippolito A, Sala S, Faber J H, et al. Ecological vulnerability analysis: A river basin case study[J]. Science of the Total Environment, 2010,408(18):3880-3890.

[2] 毛曉曦,郭云繼,崔江慧,等.濱海生態(tài)脆弱區(qū)土地生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)變化分析:以黃驊市為例[J].水土保持研究,2016,23(2):249-254.

[3] 靳毅,蒙吉軍.生態(tài)脆弱性評(píng)價(jià)與預(yù)測(cè)研究進(jìn)展[J].生態(tài)學(xué)雜志,2011,30(11):2646-2652.

[4] 周永娟,仇江嘯,王姣,等.三峽庫(kù)區(qū)消落帶生態(tài)環(huán)境脆弱性評(píng)價(jià)[J].生態(tài)學(xué)報(bào),2010,30(24):6726-6733.

[5] 張德君,高航,楊俊,等.基于GIS的南四湖濕地生態(tài)脆弱性評(píng)價(jià)[J].資源科學(xué),2014,36(4):874-882.

[6] Metzger M J, Rounsevell M D A, Acosta-Michlik L, et al. The vulnerability of ecosystem services to land use change[J]. Agriculture, Ecosystems & Environment, 2006,114(1):69-85.

[7] 鐘曉娟,孫保平,趙巖,等.基于主成分分析的云南省生態(tài)脆弱性評(píng)價(jià)[J]. 生態(tài)環(huán)境學(xué)報(bào),2011,20(1):109-113.

[8] 陳佳,楊新軍,尹莎,等.基于VSD框架的半干旱地區(qū)社會(huì)—生態(tài)系統(tǒng)脆弱性演化與模擬[J].地理學(xué)報(bào),2016,71(7):1172-1188.

[9] 廖煒,李璐,吳宜進(jìn),等.丹江口庫(kù)區(qū)土地利用變化與生態(tài)環(huán)境脆弱性評(píng)價(jià)[J].自然資源學(xué)報(bào),2011,26(11):1879-1889.

[10] 廖雪琴,李巍,侯錦湘.生態(tài)脆弱性評(píng)價(jià)在礦區(qū)規(guī)劃環(huán)評(píng)中的應(yīng)用研究:以阜新礦區(qū)為例[J].中國(guó)環(huán)境科學(xué),2013,33(10):1891-1896.

[11] 喬青,高吉喜,王維,等.生態(tài)脆弱性綜合評(píng)價(jià)方法與應(yīng)用[J].環(huán)境科學(xué)研究,2008,21(5):117-123.

[12] 吉林省地方志編簨委員會(huì).吉林省志：自然地理卷[M].長(zhǎng)春:吉林人民出版社,1992.

[13] 吉林省統(tǒng)計(jì)局.吉林統(tǒng)計(jì)年鑒[Z].長(zhǎng)春:吉林大學(xué)音像出版社,2001-2013.

[14] 于國(guó)強(qiáng),杜保佳,劉吉平.1986—2012年吉林省西部濕地生態(tài)系統(tǒng)脆弱性評(píng)價(jià)[J].濕地科學(xué),2016,14(3):439-445.

[15] 趙冰,張杰,孫希華.基于GIS的淮河流域桐柏—大別山區(qū)生態(tài)脆弱性評(píng)價(jià)[J].水土保持研究,2009,16(3):135-138.

[16] Xia X S, Zhu X F, Yue-Chen L I, et al. Evaluation for vulnerability of agroecological environment in Three Gorges Reservoir area(Chongqing section)based on AHP-PCA entropy combination weight mode[J]. Journal of Southern Agriculture, 2016,47(4):548-556.

[17] 張麗謙,韓海榮,劉利,等.山地森林生態(tài)脆弱性評(píng)價(jià)指標(biāo)體系構(gòu)建與應(yīng)用:以北京百花山自然保護(hù)區(qū)為例[J].林業(yè)資源管理,2011(1):67-71.

[18] 王延平,劉霞,姚孝友,等.淮河流域沂蒙山區(qū)水土保持生態(tài)脆弱性的AHP分析[J].中國(guó)水土保持科學(xué),2010,8(3):20-27.

[19] 姚雄,余坤勇,劉健,等.南方水土流失嚴(yán)重區(qū)的生態(tài)脆弱性時(shí)空演變[J].應(yīng)用生態(tài)學(xué)報(bào),2016,27(3):735-745.

[20] 姚坤,吳亮,相恒星,等.涼山地區(qū)生態(tài)脆弱性的時(shí)空變化[J].水土保持通報(bào),2017,37(1):329-334.

[21] 向楠,李俊清,于琳倩.內(nèi)蒙古鄂溫克族自治旗生態(tài)脆弱性動(dòng)態(tài)研究[J].西北林學(xué)院學(xué)報(bào),2014,29(4):93-97,150.

[22] 茶增芬,張翊,李銀富,等.基于全局主成分分析的羅平縣資源環(huán)境承載力動(dòng)態(tài)評(píng)價(jià)[J].安徽農(nóng)業(yè)科學(xué),2016,44(32):43-46,52.

[23] 吳健生,宗敏麗,彭建.基于景觀格局的礦區(qū)生態(tài)脆弱性評(píng)價(jià):以吉林省遼源市為例[J].生態(tài)學(xué)雜志,2012,31(12):3213-3220.

[24] 王敏,熊麗君,黃沈發(fā),等.崇明生態(tài)島建設(shè)的生態(tài)環(huán)境指標(biāo)體系研究[J].中國(guó)人口·資源與環(huán)境,2010,20(3):341-344.

[25] 黃潤(rùn),朱誠(chéng),葛向東,等.皖西大別山西北坡水土流失與生態(tài)修復(fù)[J].水土保持通報(bào),2004,24(6):90-93.

[26] 麻素挺,湯潔,林年豐.基于GIS和RS多元空間信息的吉林西部生態(tài)環(huán)境綜合評(píng)價(jià)[J].資源科學(xué),2004,26(4):140-145.

[27] 黃方,劉湘南,張養(yǎng)貞.GIS支持下的吉林省西部生態(tài)環(huán)境脆弱態(tài)勢(shì)評(píng)價(jià)研究[J].地理科學(xué),2003,23(1):95-100.

[28] 陳濤.中部六省經(jīng)濟(jì)發(fā)展?fàn)顩r分析:基于SPSS主成分分析方法[J].湖南工業(yè)職業(yè)技術(shù)學(xué)院學(xué)報(bào),2008,8(4):30-42.