雷軍成，徐海根，曹銘昌，關(guān)慶偉，丁 暉

（1.南京林業(yè)大學(xué) 森林資源與環(huán)境學(xué)院，南京210037；2.環(huán)境保護部 南京環(huán)境科學(xué)研究所，南京210042）

土地利用／覆被變化（LUCC）能夠引起地表輻射收支、碳循環(huán)及生物多樣性等一系列地表生態(tài)過程的變化，因此土地利用變化研究成為全球變化研究的核心研究領(lǐng)域之一［1－5］。不同尺度、不同地區(qū)、不同時段土地利用變化及其生態(tài)環(huán)境效應(yīng)的典型個案研究既是對全球土地利用變化研究的有益補充，也可以為區(qū)域土地利用規(guī)劃和相關(guān)政策的制定提供科學(xué)的理論依據(jù)［6－8］。國內(nèi)外學(xué)者從不同尺度、針對不同地區(qū)開展了大量的土地利用變化及其生態(tài)環(huán)境效應(yīng)研究工作［9－14］。在土地利用變化及其生態(tài)環(huán)境效應(yīng)研究中，土地利用動態(tài)度、土地利用程度綜合指數(shù)、馬爾科夫轉(zhuǎn)移矩陣、生態(tài)系統(tǒng)服務(wù)價值等分析指標和方法得到了廣泛的應(yīng)用［14－18］。

江蘇省以其獨特的區(qū)位優(yōu)勢和得天獨厚的自然條件，在我國經(jīng)濟發(fā)展過程中充當(dāng)著排頭兵的作用。改革開放以來江蘇省經(jīng)濟發(fā)展迅速，人地之間矛盾突出，但對于江蘇省土地利用時空變化特征的研究尚顯不足，主要表現(xiàn)在：（1）以往對于江蘇省內(nèi)土地利用變化的研究多集中在沿海或環(huán)太湖等熱點地區(qū)；（2）分析多采用較單一指標描述土地利用變化的時空特征；（3）對土地利用變化引起的生態(tài)環(huán)境效應(yīng)，尤其生物多樣性變化特征的研究涉及很少［19－21］。從國內(nèi)其他省域尺度土地利用變化的研究看，研究內(nèi)容多集中在土地利用變化特征分析、驅(qū)動力分析和生態(tài)系統(tǒng)服務(wù)價值評價這三個方面［22－24］。本研究采用多指標分析，全方位揭示1990—2008年江蘇省土地利用的時空變化特征，并分析土地利用變化對江蘇省生物多樣性的影響，一方面希望為今后全省土地利用規(guī)劃及其相關(guān)政策的制定提供更詳實的基礎(chǔ)資料；另一方面希望通過對經(jīng)濟發(fā)達的平原地區(qū)土地利用變化及其生物多樣性變化特征的研究，為其他類似區(qū)域的研究提供一定的借鑒。

1 數(shù)據(jù)來源及處理

1.1 數(shù)據(jù)來源

研究所使用土地利用數(shù)據(jù)來源于中國科學(xué)院地理資源與科學(xué)研究所解譯的1990年、1995年、2000年、2005年、2008年5期全國土地利用柵格數(shù)據(jù)，空間分辨率為1km×1km；1∶400萬全國和江蘇省行政區(qū)劃矢量數(shù)據(jù)來源于國家基礎(chǔ)地理信息系統(tǒng)（http：∥nfgis.nsdi.gov.cn／）。

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

利用江蘇省1∶400萬行政區(qū)劃矢量數(shù)據(jù)，通過ArcGIS的空間分析功能，提取出江蘇省1990—2008年5期土地利用類型數(shù)據(jù)。研究所使用土地利用原始數(shù)據(jù)將土地利用類型劃分耕地、林地、水體等6個一級地類，在一級地類基礎(chǔ)上又進一步將土地利用類型細分為水田、旱地、有林地等25個二級地類。為研究方便和便于橫向?qū)Ρ龋狙芯繉@得的土地利用原始數(shù)據(jù)整理歸并到一級分類中的6類，即耕地、林地、草地、水體、建設(shè)用地和未利用地，并分別編碼分別為1—6。

2 研究方法

2.1 綜合土地利用動態(tài)度分析

綜合土地利用動態(tài)度不考慮土地利用變化的內(nèi)在過程，定量地反映研究區(qū)總體上土地利用變化的速度?？紤]土地利用類型轉(zhuǎn)入與轉(zhuǎn)出雙向變化過程的綜合土地利用動態(tài)度計算方法如公式（1）所示：

式中：ΔUi－out——研究時段T內(nèi)土地利用類型i轉(zhuǎn)變?yōu)槠渌恋乩妙愋偷拿娣e之和；ΔUi－in——研究時段T內(nèi)其他土地利用類型轉(zhuǎn)變?yōu)閕類土地利用類型的面積之和；Uia——研究區(qū)土地利用類型i的面積；n——研究區(qū)土地利用類型數(shù)量；T——研究時段時間間隔。

2.2 土地利用轉(zhuǎn)移分析

2.3 土地利用程度變化分析

土地利用程度反映了土地利用的廣度和深度，它不僅能反映土地利用類型本身的自然屬性，還能反映出土地系統(tǒng)中人類因素和自然環(huán)境的綜合效應(yīng)。劉紀遠等從生態(tài)學(xué)的角度出發(fā)，提出了土地利用程度綜合分析方法，即將土地利用程度按照土地自然綜合體在社會因素影響下的自然平衡狀態(tài)分為4個級別：土地未利用級、土地自然再生利用級、土地人為再生利用級和土地非再生利用級，并分別賦值1～4（見表1）。土地利用程度計算公式為［27－29］：

式中：Ai——第i級土地利用程度分級指數(shù)；Ci——第i級土地利用程度分級面積百分比；n——研究區(qū)土地利用類型數(shù)。

其中Yt為產(chǎn)出，Xt=(x1t,x2t,…,xnt)為要素投入向量，xnt為第n種要素投入，A(t)為技術(shù)系數(shù)。假設(shè)函數(shù)為希克斯技術(shù)中性和規(guī)模收益不變，對(1)式兩邊求導(dǎo)，則有：

表1 土地利用程度分級賦值

2.4 土地利用類型重心轉(zhuǎn)移分析

將1990—2008年江蘇省5期土地利用柵格數(shù)據(jù)按不同土地利用類型轉(zhuǎn)換為矢量多邊形數(shù)據(jù)后，利用ArcGIS多邊形重心分析功能，求取1990—2008年江蘇省各土地利用類型的重心。

2.5 土地利用轉(zhuǎn)化頻次分析

根據(jù)地圖代數(shù)原理，計算獲得1990—2008年研究區(qū)土地利用類型轉(zhuǎn)移歷史序列柵格數(shù)據(jù)，公式如下：

式中：Ki×j——1990—2008年土地利用類型轉(zhuǎn)移歷史序列柵格圖；A1990i×j——1990年研究區(qū)土地利用類型柵格圖，其余分別為1995年、2000年、2005年、2008年研究區(qū)土地利用類型柵格圖；

獲得的土地利用類型轉(zhuǎn)移歷史序列柵格數(shù)據(jù)的萬位數(shù)到個位數(shù)依次代表1990年、1995年、2000年、2005年、2008年各柵格點上依次出現(xiàn)的土地利用類型。統(tǒng)計研究時段內(nèi)各柵格點上土地利用類型轉(zhuǎn)換的頻次。

2.6 土地利用類型空間分布

空間洛倫茨曲線和基尼系數(shù)（Gini coefficient，GI）是重要的經(jīng)濟學(xué)分析方法，在非經(jīng)濟領(lǐng)域得到一定的應(yīng)用［30－32］。本研究以江蘇省13個地級市國土面積占全省國土面積累積百分比為橫坐標，以13個地級市各土地利用類型占全省該土地利用類型累積百分比為縱坐標，繪制1990—2008年江蘇省各土地利用類型在13個地級市空間分布的洛倫茨曲線，并計算基尼系數(shù)?；嵯禂?shù)越趨近于0，說明某種地類在13個地級市的空間分布越均勻；越趨近于1，則說明空間分布越不均勻。參照聯(lián)合國有關(guān)組織對于居民收入基尼系數(shù)的規(guī)定，本研究將GI＜0.2表示“分布平均”，0.2≤GI＜0.3表示“相對平均”，0.3≤GI＜0.4表示“相對合理”，0.4≤GI＜0.5表示“差距大”，GI＞0.6表示“差距懸殊”。

2.7 MSA變化分析

為了評估土地利用變化、生境破碎化、人為干擾等因素對不同生態(tài)系統(tǒng)生物多樣性的影響，Alkemade等提出GLOBIO3評估框架。GLOBIO3評估框架通過平均物種豐富度 （Mean Species Abundance，MSA）指數(shù)反應(yīng)評價單元上生物多樣性的時空變化。MSA實際上是衡量生態(tài)系統(tǒng)偏離其原始或自然狀態(tài)的一個指標。以未受干擾自然狀態(tài)下的生態(tài)系統(tǒng)的MSA值為1，受干擾越大其MSA值降低越多。GLOBIO3評估框架首先依據(jù)不同生態(tài)系統(tǒng)受人為干擾的強弱，賦予不同生態(tài)系統(tǒng)一個初始的MSA值，然后進一步考慮生境破碎化、面積效應(yīng)、邊界效應(yīng)等對生態(tài)系統(tǒng)生物多樣性的影響，最終得出評價單元的MSA值。本研究基于GLOBIO3評估框架，各生態(tài)系統(tǒng)類型初始MSA賦值，參考《中國生態(tài)系統(tǒng)服務(wù)與管理戰(zhàn)略》，同時結(jié)合研究區(qū)實際，取值為耕地0.2，林地0.5，草地0.35，河流與湖泊0.7，建設(shè)用地0.05，未利用地0.7［33－34］。

3 結(jié)果與分析

3.1 土地利用結(jié)構(gòu)變化特征

圖1 1990－2008年江蘇省各土地利用類型面積變化

3.2 綜合土地利用動態(tài)度特征

江蘇省1990—2008年綜合土地利用動態(tài)度變化如圖2所示。1990—2008年江蘇省土地利用的平均年變化率為0.45%，但不同時段的變化率有明顯的差別，1995—2000年的變化率最低為0.18%，2000年以后各時段的變化率要明顯高于歷史時期。

3.3 土地利用類型轉(zhuǎn)移特征

1990—2008年江蘇省土地利用類型轉(zhuǎn)移的馬爾科夫轉(zhuǎn)移矩陣如表2所示（僅分析土地利用類型的最終轉(zhuǎn)移結(jié)果，不考慮轉(zhuǎn)移的中間過程）。從不同土地利用類型的轉(zhuǎn)移模式來看，1990—2008年江蘇省建設(shè)用地的增加主要來源于耕地，有5 586km2的耕地最終轉(zhuǎn)化為建設(shè)用地，占耕地流失總面積的81.36%和建設(shè)用地凈增加面積的98.16%；耕地除主要轉(zhuǎn)化為建設(shè)用地外，還轉(zhuǎn)化為水體，達1 081km2，占耕地流失總面積的15.74%；草地的減少主要轉(zhuǎn)化為耕地、水體和建設(shè)用地，三者總面積達449km2，合占草地流失總面積的98.68%；林地的減少部分主要轉(zhuǎn)化為耕地和建設(shè)用地，二者合占林地流失總面積的93.62%。從總體上可以看出建設(shè)用地成為其他土地利用類型流失的主要去向。

圖2 1990－2008年江蘇省綜合土地利用動態(tài)度變化

3.4 土地利用程度變化特征

1990—2008年江蘇省土地利用程度綜合指數(shù)呈逐漸升高的趨勢（見圖3），從1990年的289.2上升到2008年的295.1。從不同時期土地利用程度變化量來看，各時期土地利用程度綜合指數(shù)變化量均大于0，其中2005—2008年3a間江蘇省土地利用程度綜合指數(shù)升高2.77，較1990—2005年15a間研究區(qū)土地利用程度綜合指數(shù)升高的量高出近12.11%。

圖3 江蘇省1990－2008年土地利用程度綜合指數(shù)變化

表2 1990－2008年江蘇省土地利用類型轉(zhuǎn)移矩陣km2

3.5 土地利用重心轉(zhuǎn)移特征

1990—2008年，江蘇省各土地利用類型重心發(fā)生了不同程度的轉(zhuǎn)移（見表3），以建設(shè)用地重心轉(zhuǎn)移規(guī)律最為明顯。如圖4所示，1990—2008年江蘇省建設(shè)用地重心明顯向南偏東方向移動。1990—2008年18a間江蘇省建設(shè)用地重心的位移距離達80km，年均向南偏東21.42°方向移動4.44km。

表3 1990－2008年江蘇省各土地利用類型重心轉(zhuǎn)移（1°）

圖4 1990－2008年江蘇省建設(shè)用地重心變化

3.6 土地利用轉(zhuǎn)化頻次特征

江蘇省1990—2008年各柵格點上土地利用類型轉(zhuǎn)換頻次如圖5所示。可以看出江蘇省1990—2008年各柵格點上發(fā)生土地利用類型轉(zhuǎn)換的地區(qū)以1～2次轉(zhuǎn)換為主，極少數(shù)地區(qū)發(fā)生3次以上土地利用類型的轉(zhuǎn)換。從發(fā)生土地利用類型轉(zhuǎn)換地區(qū)的空間分布上來講，蘇南地區(qū)發(fā)生土地利用類型轉(zhuǎn)換柵格的數(shù)量和轉(zhuǎn)換的頻次都要明顯高于蘇北地區(qū)。

3.7 土地利用類型空間分布變化特征

江蘇省1990—2008年各土地利用類型在13個地級市空間分布的洛倫茨曲線和基尼系數(shù)變化如圖6所示。如果某種地類基尼系數(shù)變大，說明該地類在13個地級市的空間分布趨于越來越不均勻，反之則表示趨于越來越均勻。1990—2008年江蘇省耕地和水體在13個地級市空間分布的基尼系數(shù)有增大的趨勢，但增大趨勢并不明顯，各年份基尼系數(shù)均小于0.2，表明研究時段內(nèi)江蘇省耕地和水體在13個地級市的空間分布“分布平均”；建設(shè)用地基尼系數(shù)的變化呈現(xiàn)以2000年為時間節(jié)點，1990—2000年建設(shè)用地基尼系數(shù)減小，2000—2008年建設(shè)用地基尼系數(shù)變大，但總體上研究時段內(nèi)建設(shè)用地在13個地級市的空間分布仍可以認為是“分布平均”；2005年以前，林地基尼系數(shù)變化不大，2008年較2005年有明顯的上升，總體上研究時段內(nèi)林地在13個地級市的空間分布“相對合理”；草地的基尼系數(shù)有明顯下降局勢，表明草地在研究時段內(nèi)在13個地級市的空間分布有趨向越來越均勻的趨勢；未利用地面積小，基尼系數(shù)變化較大。

3.8 MSA變化特征

1990—2008年江蘇省各柵格點上MSA指數(shù)平均值的變化狀況如圖7所示?？傮w上，1990—2008年，江蘇省MSA指數(shù)呈持續(xù)下降趨勢，但是各時間段MSA指數(shù)的下降幅度不盡相同。1990—2005年MSA指數(shù)下降比較緩慢，2005—2008年MSA指數(shù)的下降幅度較大，甚至超過了1990—2005年15a間總的下降幅度。

圖5 1990－2008年江蘇省土地利用類型轉(zhuǎn)換頻次

圖6 1990－2008年江蘇省13個地級市土地利用空間分布洛倫茨曲線及基尼系數(shù)變化

圖7 1990－2008年江蘇省MSA指數(shù)變化

4 結(jié)論與討論

本文利用1990—2008年江蘇省土地利用類型遙感解譯數(shù)據(jù)，對江蘇省土地利用和生物多樣變化性特征進行分析。研究結(jié)果表明，1990—2008年江蘇省土地利用及其生物多樣性狀況都發(fā)生了比較明顯的變化。

（1）從總體上看，1990—2008年，江蘇省土地利用變化的主要特點是：建設(shè)用地、水體面積的增加和耕地、草地面積的減少。耕地成為土地利用流失的主要“源”而建設(shè)用地成為土地利用流失的主要“匯”。究其原因可能主要包含如下幾個方面：①1990—2008年江蘇省總?cè)丝趶? 766.9萬增加至7 676.5萬，增長13.44%，國 內(nèi) 生 產(chǎn) 總 值 從 1 416.5 億 元 增 加 至30 981.98億元，增長20.87倍，人口增長和經(jīng)濟快速發(fā)展導(dǎo)致對居住用地、交通建設(shè)用地、企業(yè)建設(shè)用地等需求量增加；② 農(nóng)村由于缺乏統(tǒng)一的規(guī)劃和管理不嚴等眾多原因，非法占地、超標占地、一戶多宅的現(xiàn)象比較普遍，造成宅基地不斷向村外擴展，直接造成農(nóng)村建設(shè)用地的擴張和耕地面積的減少［35］；③耕地基本位于平原地區(qū)和居民區(qū)附近，轉(zhuǎn)化為居住用地、企業(yè)用地等建設(shè)用地的成本相對較低。

（2）1990—2008年，江蘇省土地利用處于較快速發(fā)展時期，但土地利用類型變化表現(xiàn)出比較明顯的時間和空間不均衡性。2000年以后，土地利用變化的數(shù)量和強度明顯高于2000年以前；蘇南地區(qū)土地利用變化的轉(zhuǎn)換頻次明顯高于蘇北地區(qū)。導(dǎo)致這一現(xiàn)象產(chǎn)生的主要原因可能是經(jīng)濟發(fā)展的時空差異。1990—2000年，江蘇省GDP總量增長503.86%，但由于基數(shù)低，10a間GDP增長的總量僅為7 137.19億元；2000—2008年，江蘇省GDP總量增幅（262.21%），雖然低于1990—2000年的增幅，但GDP增長的總量卻達22 428.29億元，增長量為1990—2000年GDP增長量的3倍多［36］。江蘇省經(jīng)濟發(fā)展的空間差異明顯，蘇南地區(qū)地處“長三角”，受上海市經(jīng)濟輻射的影響和地區(qū)工業(yè)發(fā)展的規(guī)模效應(yīng)，經(jīng)濟發(fā)展迅速，大量土地流向居住用地、企業(yè)用地和交通用地等建設(shè)用地，而蘇北地區(qū)由于經(jīng)濟發(fā)展相對遲緩，第一產(chǎn)業(yè)生產(chǎn)總值占地區(qū)生產(chǎn)總值比重較高，土地利用類型變動需求小，土地利用類型轉(zhuǎn)換的頻次和發(fā)生的面積也相對較低［37－38］。

（3）研究時段內(nèi)，江蘇省各土地利用類型在13個地級市的空間分布除未利用地外，林地和草地空間分布的基尼系數(shù)也相對較高，究其原因可能是由于江蘇省自然條件差異所致。林地和草地主要分布于低山丘陵等人工轉(zhuǎn)化難度較大的地區(qū)，而江蘇省低山丘陵主要位于全省的北部和西南部。

（4）研究時段內(nèi)，江蘇省生物多樣性MSA指數(shù)呈現(xiàn)持續(xù)下降趨勢，其直接原因主要是由于MSA指數(shù)較高的土地利用類型轉(zhuǎn)化為MSA指數(shù)較低的土地利用類型，如草地和林地轉(zhuǎn)化為耕地，耕地轉(zhuǎn)化為MSA指數(shù)更低的建設(shè)用地；土地利用類型發(fā)生轉(zhuǎn)化的同時，伴隨著生態(tài)系統(tǒng)的破碎化。

本研究利用地圖代數(shù)原理，充分發(fā)揮柵格數(shù)據(jù)的優(yōu)勢，構(gòu)建土地利用類型流轉(zhuǎn)歷史序列數(shù)據(jù)，統(tǒng)計并空間化表達每個柵格點上土地利用類型轉(zhuǎn)換的頻次，較傳統(tǒng)的馬爾科夫轉(zhuǎn)移矩陣分析方法能更直觀、更詳盡地展示土地利用類型轉(zhuǎn)換的強度，可以為以后相關(guān)分析提供一定的參考。

需要指出的是，由于不同的研究需要及研究人員對于各土地利用類型的理解和認知不同，同一土地利用方式在不同的研究中可能被劃分為不同的土地利用類型，尤其是兼具多種土地利用類型特點的土地利用方式，如本研究及有些相關(guān)研究將鹽田歸為建設(shè)用地，有些研究則將其歸為水體，有些研究甚至直接將其歸為單獨一類［39－42］。因此，對于土地利用類型的劃分應(yīng)主要立足于研究需要，在此基礎(chǔ)上進一步兼顧不同研究之間的可比性。

本研究也存在一定的不足之處，主要表現(xiàn)在：（1）所使用遙感解譯數(shù)據(jù)的空間分辨率仍不夠高，土地利用類型在更細微空間上的變化無法得到充分體現(xiàn)；（2）研究所使用土地利用數(shù)據(jù)的時間分辨率不夠高，無法反應(yīng)土地利用年季的甚至更高時間分辨率的變化特征；（3）本研究將研究區(qū)土地利用類型分為6類，掩蓋了更低級別土地利用類型的時空轉(zhuǎn)換特征；（4）由于沒有研究時段內(nèi)研究區(qū)詳細的道路數(shù)據(jù)，對生態(tài)系統(tǒng)破碎化導(dǎo)致的MSA指數(shù)的降低可能存在過低估計。因此，采用高時空分辨率土地利用數(shù)據(jù)，基于更精細的土地類型分類，是未來土地利用變化時空特征研究需要優(yōu)先考慮的方面。

［1］ Tyson P D，F(xiàn)uchs R，F(xiàn)u C，et al.The earth system：regional－global linkages［J］.Regional Environmental Change，2001，2（3）：128－140.

［2］ Verburg P H，Veldkamp W S A，Espaldon R L V.Modeling the spatial dynamics of regional land use：the CLUE－S model［J］.Environmental Management，2002，30（3）：391－405.

［3］ Yu W H，Zang S Y，Wu C S，et al.Analyzing and modeling land use land cover change（LUCC）in the Daqing City，China ［J］.Applied Geography，2011，31（2）：600－608.

［4］ 盛晟，劉茂松，徐馳，等.CLUE－S模型在南京市土地利用變化研究中的應(yīng)用［J］.生態(tài)學(xué)雜志，2008，27（2）：235－239.

［5］ 張躍紅，安裕倫，馬良瑞，等.1960—2010年貴州省喀斯特山區(qū)陡坡土地利用變化［J］.地理科學(xué)進展，2012，31（7）：878－884.

［6］ Long H L，Wu X Q，Wang W J，et al.Analysis of urban－rural land－use change during 1995—2006and its policy dimensional driving forces in Chongqing，China［J］.Sensors，2008，8（2）：681－699.

［7］ 蒲雯，張超，王淑華，等.基于元胞自動機復(fù)合模型的村鎮(zhèn)土地利用變化模擬［J］.中國農(nóng)業(yè)大學(xué)學(xué)報，2012，17（5）：182－189.

［8］ TurnerⅡB L，Lambin E F，Reenberg A.The emergence of land change science for global environmental change and sustainability［J］.Proceedings of the National Academy of Sciences，2007，104（52）：20666－20671.

［9］ Liu J Y，Zhan J Y，Deng X Z.Spatio－temporal patterns and driving forces of urban land expansion in China during the economic reform era［J］.Royal Swedish Acade－my of Sciences，2005，34（6）：450－455.

［10］ Reidsma P，Tekelenburg T，Van Den Berg M，et al.Impact of land－use change on biodiversity：An assessment of agricultural biodiversity in the European Union［J］.Agriculture，Ecosystems and Environment，2006，114（1）：86－102.

［11］ Temme A J A M，Verburg P H.Mapping and modeling of changes in agricultural intensity in Europe［J］.Agriculture，Ecosystems and Environment，2011，140（1）：45－56.

［12］ 張寧，曾光建，關(guān)國鋒.近25年來黑龍江省生態(tài)系統(tǒng)服務(wù)價值對土地利用變化的響應(yīng)［J］.國土與自然資源研究，2011（5）：16－18.

［13］ 侯鵬，王橋，王昌佐，等.流域土地利用／土地覆被變化的生態(tài)效應(yīng)［J］.地理研究，2011，30（11）：2092－2098.

［14］ 錢江波，王恩，章銀柯.杭州西湖湖西景區(qū)土地利用／覆被變化與生態(tài)服務(wù)價值評估［J］.西北林學(xué)院學(xué)報，2011，26（6）：208－212.

［15］ 柯新利，韓冰華，劉蓉霞，等.1990年以來武漢城市圈土地利用變化時空特征研究［J］.水土保持研究，2012，19（1）：76－81.

［16］ 王秀蘭，包玉梅.土地利用動態(tài)變化研究方法探討［J］.地理科學(xué)進展，1999，18（1）：81－87.

［17］ 張楊，劉艷芳，丁慶，等.1996—2006年武漢市土地利用變化研究［J］.生態(tài)環(huán)境學(xué)報，2010，19（11）：2534－2539.

［18］ 胡英敏，高瓊，蘭玉芳，等.太仆寺旗2000—2008年EVI對氣候及土地利用變化的響應(yīng)［J］.自然資源學(xué)報，2012，27（7）：1200－1213.

［19］ 李閩，孫在宏，楊素靜，等.江蘇省環(huán)太湖地區(qū)土地利用結(jié)構(gòu)演變與驅(qū)動力分析［J］.地球信息科學(xué)，2008，10（2）：136－141.

［20］ 王千，金曉斌，周寅康.江蘇沿海地區(qū)耕地景觀生態(tài)安全格局變化與驅(qū)動機制［J］.生態(tài)學(xué)報，2011，31（20）：5903－5909.

［21］ 嚴長清，孫偉，陸效平，等.江蘇省沿海灘涂土地利用與生態(tài)保護研究［J］.生態(tài)科學(xué)，2007，26（3）：263－268.

［22］ 趙哲遠，馬奇，華元春，等.浙江省1996—2005年土地利用變化分析［J］.中國土地科學(xué)，2009，23（11）：54－60.

［23］ 王剛，謝世友.基于土地利用變化的重慶市生態(tài)系統(tǒng)服務(wù)功能價值評估［J］.西南師范大學(xué)學(xué)報：自然科學(xué)版，2011，36（5）：177－181.

［24］ Tong S Y，Lu Y Y.A Study on Vertical Spatial Transfer of Cultivated Land in Yunnan Province from 1986to 2008［J］.Agricultural Science ＆ Technology，2011，12（1）：127－132.

［25］ 李吉鵬，張瑞娟，李川，等.基于GIS和 Markov模型的孟達自然保護區(qū)土地利用變化研究［J］.陜西林業(yè)科技，2012（5）：31－33.

［26］ 徐嘉興，李鋼，張海榮.太湖流域土地利用變化及其驅(qū)動力定量分析［J］.國土與自然資源研究，2012（3）：19－21.

［27］ 劉紀遠.中國資源環(huán)境遙感宏觀調(diào)查與動態(tài)研究［M］.北京：中國科學(xué)技術(shù)出版社，1996：158－188.

［28］ 吳明發(fā)，歐名豪，廖榮浩.經(jīng)濟發(fā)達地區(qū)土地利用變化及其驅(qū)動力分析：以廣東省為例［J］.水土保持研究，2012，19（1）：179－183.

［29］ 李陽兵，謝靜，邵景安，等.川東平行嶺谷區(qū)土地利用變化時空分異特征［J］.長江流域資源與環(huán)境，2011，20（4）：422－427.

［30］ 劉景輝，李立軍，王志敏.中國糧食安全指標的探討［J］.中國農(nóng)業(yè)科技導(dǎo)報，2004，6（4）：10－16.

［31］ 劉洋，金鳳君，甘紅.區(qū)域水資源空間匹配分析［J］.遼寧工程技術(shù)大學(xué)學(xué)報，2005，24（5）：657－660.

［32］ 陳軍偉，孔祥斌，張鳳榮，等.基于空間洛倫茨曲線的北京山區(qū)土地利用結(jié)構(gòu)變化［J］.中國農(nóng)業(yè)大學(xué)學(xué)報，2006，11（4）：71－74.

［33］ Alkemade R，van Oorschot M，Miles L，et al.GLOBIO3：A framework to investigate options for reducing global terrestrial biodiversity loss［J］.Ecosystem，2009，12（3）：374－390.

［34］ 陳宜瑜.中國生態(tài)系統(tǒng)服務(wù)與管理戰(zhàn)略［M］.中國環(huán)境科學(xué)出版社，2011.

［35］ 呂亞生，張小林.江蘇省農(nóng)村城鎮(zhèn)化進程中的村莊用地集約化研究［J］.中國土地科學(xué)，2007，21（2）：65－69.

［36］ 江蘇省統(tǒng)計局.江蘇省統(tǒng)計鑒［EB／OL］.［2013－05－27］.http： ∥ www.jssb.gov.cn／jstj／jsnj／2009／nj02／nj0203.htm.

［37］ 曹芳東，吳江，徐敏.基于空間計量經(jīng)濟模型的縣域經(jīng)濟發(fā)展差異研究：以江蘇省為例［J］.地域研究與開發(fā)，2010，29（6）：23－28.

［38］ 熊薇，徐逸倫，王迎英.江蘇省縣域經(jīng)濟差異時空演變［J］.地理科學(xué)進展，2011，30（2）：224－230.

［39］ 劉艷芬，張杰，馬毅，等.融合地學(xué)知識的海岸帶遙感圖像土地利用／覆被分類研究［J］.海洋科學(xué)進展，2010，28（2）：193－202.

［40］ 謝菲，舒曉波，廖富強，等.浮梁縣土地利用變化及驅(qū)動力分析［J］.水土保持研究，2011，18（2）：213－221.

［41］ 孫賢斌，劉紅玉.江蘇鹽城市海濱土地利用對景觀生態(tài)風(fēng)險的影響［J］.國土資源遙感，2011，23（3）：140－145.

［42］ 李振，黃海軍.膠州灣海岸帶土地利用／覆被變化研究［J］.國土資源遙感，2010，22（4）：71－76.