1986-2013年南京市土地利用/覆被景觀格局時空變化及驅動力因素分析

佟光臣, 林 杰, 陳 杭, 顧哲衍, 唐 鵬, 張金池

(1.江蘇省南方現代林業協同創新中心, 江蘇省水土保持與生態修復重點實驗室, 南京 210037;2.江蘇省水利勘測設計研究院有限公司, 江蘇 揚州 225127)

1986-2013年南京市土地利用/覆被景觀格局時空變化及驅動力因素分析

佟光臣1, 林 杰1, 陳 杭2, 顧哲衍2, 唐 鵬1, 張金池1

(1.江蘇省南方現代林業協同創新中心, 江蘇省水土保持與生態修復重點實驗室, 南京 210037;2.江蘇省水利勘測設計研究院有限公司, 江蘇 揚州 225127)

以南京市為研究對象，采用南京市1986年、1996年、2002年、2013年的四期土地利用/覆被數據，運用景觀生態學原理，借助ArcGIS軟件和Fragstats軟件，從土地利用類型結構、變化速度、轉化關系和景觀格局等方面分析南京市近30年土地利用/覆被和景觀格局的時空變化狀況及其驅動因素。結果表明：(1) 南京市近30 a中土地利用/覆被發生了劇烈變化，耕地在整個研究區起到主導作用，其面積比例均在48%以上，但是其優勢度逐漸降低，大量自然景觀向人為景觀轉變是研究區土地利用變化最顯著的特征；(2) 整體景觀類型變化速度都在增大，而草地和其他用地變化最為活躍，林地相對最穩定，土地利用結構處于快速調整階段；(3) 耕地的轉出最大，主要轉向為居民點及建設用地；(4) 景觀格局變化明顯，破碎化程度越來越高，景觀格局趨向復雜化，景觀結構不穩定性增加；(5) 人口增長、經濟發展、區域政策和城市規劃等共同推動了研究區景觀格局的變化。該研究結果對于南京市的科學規劃及可持續發展提供重要參考。

土地利用; 景觀格局; 驅動力因素; 南京市

當下，許多發展中國家正經歷著城市化浪潮，城市化發展勢頭迅猛，成為根本上改變全球土地利用/覆被景觀格局的最強大因素，引起土地利用/覆被景觀格局和動態演變的復雜化，進而對生態系統的結構、功能和演變造成了深刻的影響[1]，產生諸如城市熱島效應[2-3]、地表徑流與土壤侵蝕[4-5]、土壤性質變異[6-8]、水資源污染[9]、生物多樣性減少[10]等問題。目前，土地利用/覆被變化及景觀格局變化是地理學和景觀生態學領域長期關注的熱點問題[11-12]，也已成為全球變化研究的前沿和熱點課題[13-15]。賈保全[16]等選擇南京市1995年、2000年、2008年的土地利用數據，從景觀轉移矩陣與景觀斑塊的空間穩定性進行了探討，王慧[17]選擇南京市1989年、2000年、2009年土地利用數據，從南京市土地利用現狀及變化情況、景觀格局變化特征、馬爾科夫模型預測三方面對南京市景觀格局演變及發展趨勢作了分析。南京作為長三角地區重要的中心城市，自改革開放后經濟發展突飛猛進，在十二五時期(2011—2015)，通過撤縣設區調整行政區劃進入了“無縣時代”，同時為舉辦2014年青奧會而進行一系列大規模的城市建設活動，城市土地利用/覆被變化劇烈，因此有必要增加時間序列來分析從改革開放到現在的南京市景觀格局演變過程及驅動因素。

本文以南京市為研究對象，從土地利用變化速度、轉移方向、景觀指數和驅動力因素等方面探索城市土地利用及景觀格局時空動態變化對快速城市化進程的響應，以期為南京市土地資源可持續開發和利用、生態環境保護等提供重要參考。

1 研究區與數據源

1.1 研究區概況

南京市，簡稱寧，江蘇省省會，地處中國東部地區，是長江三角洲重要的中心城市。地理坐標為31°14″—32°37″N,118°22″—119°14″E。全市總面積為6 582 km2?，F轄玄武區、秦淮區、建鄴區、鼓樓區、雨花臺區、棲霞區、浦口區、江寧區、六合區、溧水區和高淳區11個區。氣候屬于北亞熱帶季風氣候，年降水1 200 mm，年平均溫度15.4°C。改革開放后南京市社會經濟快速發展，截至2013年末全市戶籍總人口為643.09萬人，常住人口為818.78萬人，GDP為8 011.7億元。南京市第一產業、第二產業和第三產業的比重由1986年的14.1∶60.9∶25.0，轉變為2013年的2.55∶43.07∶54.38。

1.2 數據源及其預處理

部分數據來源于中國科學院資源環境科學數據中心提供的1986年、1996年、2002年三期土地利用覆被數據。另外2013年南京市土地利用覆被數據是基于南京市2013年8月11日30 m分辨率Landsat8 OLI影像數據矢量化得到，影像可在美國地質調查局(USGS)官網免費下載。首先應用ENVI 5.1軟件對該影像進行輻射定標、大氣糾正、圖像融合(與第八波段15 m分辨率的全色波段融合)、幾何裁剪等處理，然后導入ArcGIS 9.3軟件進行人機交互解譯。經實地采樣驗證，解譯精度在88%以上，滿足本研究需要，至此，共建立4期LUCC數據庫。

2 研究方法

2.1 土地利用類型劃分

2013年數據處理的土地利用分類標準為了與前三期數據相一致，采用全國農業區劃委員會1984年9月發布的《土地利用現狀調查技術規程》制定的《土地利用現狀分類及含義》?？紤]到研究區情況，最終將土地利用景觀類型劃分為6類：耕地、林地、草地、居民點及建設用地、水域、其他用地。

2.2 土地利用變化分析

2.2.1 南京市土地利用變化分析 為充分反映南京市土地利用在時間和空間上的動態演變過程，分別采用單一土地利用類型動態度、綜合土地利用動態度[18-19]和土地利用轉移矩陣，從不同景觀類型各研究時段內的變化差異、劇烈程度以及土地利用類型變化的主導過程3方面對南京市土地利用變化進行解析，主要計算公式如下：

(1) 單一土地利用類型動態度

(2) 綜合土地利用動態度

式中：LC為綜合土地利用類型動態度；LCi為研究初期第i類土地利用類型的面積(hm2)；ΔLUi-j為研究時段內第i類土地利用類型轉為其他土地利用類型面積的絕對值(hm2)；T2-T1為研究時段(a)。

(3) 土地利用轉移矩陣。轉移概率矩陣在景觀動態分析中被廣泛使用，該方法可以全面顯示區域土地利用變化的細節結構特征。研究中采用土地利用面積作為土地利用狀態轉移矩陣中的向量，制作了1986—1996，1996—2002，2002—2013年三個時期的轉移概率矩陣。

2.2.2 南京市土地利用景觀格局變化分析 景觀指數能夠定量描述景觀空間結構和格局的變化，本文分別從景觀水平和斑塊類型2個尺度方面選用被國內外廣泛使用且公式意義明確的指標進行景觀格局分析。首先在斑塊類型上主要選擇斑塊數(NP)、斑塊類型比例(PLAND)、平均斑塊面積(AREA_MN)、周長—面積分維數(PAFRAC)；景觀水平上主要選擇多樣性指數(SHDI)、均勻度指數(SHEI)、景觀分離度(DⅣISION)、斑塊密度(PD)、景觀形狀指數(LSI)、蔓延度指數(CONTAG)、景觀連通度指數(COHESION)。各類指標的計算公式及生態學意義已經常見，在此不再闡述，所有指標的計算采用Fragstas3.3軟件。

2.3 景觀格局演變驅動力因素分析

城市景觀的形成過程中，人類活動對土地利用/覆被景觀格局的變化具有很大的影響。本研究主要從經濟發展、區域政策、人口因素三方面選取固定資產投資額、城鎮居民人均可支配收入、財政收入、總人口、農民人均純收入、GDP等指標來分析這些因素。

3 結果與分析

3.1 土地利用/覆被變化分析

在1986年，南京市的優勢景觀主要為耕地、居民點及建設用地和林地，分別占土地利用/覆被總面積的67.29%，11.85%，10.9%。到2013年，優勢景觀排序沒變，但是相應比重變成了48.58%，27.27%，12.60%，耕地作為本底景觀的優勢在不斷下降(圖1)。由圖2可知，近30a來南京市耕地面積一直處于下降趨勢，與之相反，居民點及建設用地的面積處于增長態勢，表現出自然景觀快速向人為景觀轉變，這與南京市城市化快速發展有關，大量農田轉化成城市建設用地。

圖1 1986-2013年南京市土地覆被結構變化

圖2 1986-2013年各種土地覆被類型面積變化

南京市近30a來各類用地的變化見圖3。在1986—1996年圖形變化并不顯著，而1996—2002，2002—2013年則變化非常明顯，最直觀的變化就是居民點及建設用地的大幅擴張與耕地的快速減少，這與近20a間城市快速城市化進程密切相關；林地、水域、其他用地面積的比例相對穩定，其中林地略有增長，各時段林地面積均占到總面積的10%以上，這與南京建設現代化國際性人文綠都、為舉辦2014年青奧會而實施的一系列生態綠化工程有關。

圖3 南京市1986-2013年土地利用類型空間分布

3.2 土地利用/覆被變化動態度分析

土地利用動態度是指區域一定時間范圍內，某種土地利用類型變化速度的情況，可定量描述土地利用變化速度[20]。單一動態度是正值時，表示該類型面積增加，反之則減??；具體衡量時使用絕對值，絕對值越低，表明該土地利用類型轉換數量少，狀態則相對穩定，反之則變化相對劇烈。從同一時段內來看，第一時段內所有景觀類型變化速度均較低，居民點及建設用地變化相對最快；后兩個時段內所有景觀類型變化速度明顯增大，且都是草地的變化最大。從不同時間段變化來看，除水域變化速度先變大后變小之外，其他景觀類型變化速度均呈上升趨勢，而草地和其他用地變化最為活躍，林地相對最穩定。

綜合動態度值可對整個南京市土地利用類型動態度進行綜合描述。參考劉紀遠[21]按土地利用動態度將中國土地利用時間動態特征劃分為4種類型：0%～3%(極緩慢變化型)；4%～12%(慢速變化型)；13%～20%(快速變化型)；21%～24%(急劇變化型)。由表1可知，南京市在第一個研究時段內屬于極緩慢變化性，之后兩個時段均大于24%，已超過急劇變化型，由此可見南京市近20 a來土地利用/覆被變化是非常劇烈的。

表1 南京市不同時段土地利用類型動態度指數

3.3 土地利用類型轉化

轉移矩陣能夠揭示地類變化過程中的細節特點，可全面、具體地刻畫區域土地利用變化的結構特征與各用地類型變化方向[22]。由表2—4可知：

1986—2013年土地利用類型相互轉換情況：耕地轉出面積最大，主要轉向為居民點及建設用地，其次分別為水域和林地，這與南京城市快速化發展、高淳區擴大水產養殖、居民點搬遷建立苗圃生態旅游等政策密切相關；居面點及建設用地的增加主要靠占用耕地，其次是林地；林地的主要轉出類型是居民點及建設用地，其次是耕地；變化最小的是草地和其他用地，原因是在整個土地利用面積中所占份額最少。

表2 南京市1986-1996年土地利用景觀轉移矩陣 hm2

表3 南京市1996-2002年土地利用景觀轉移矩陣 hm2

3.4 景觀格局演變分析

3.4.1 景觀水平上的景觀格局變化分析 由表5可知，在景觀尺度上從1986—2013年南京市的多樣性指數、均勻度指數、景觀分離度、斑塊密度、景觀形狀指數一直呈上升趨勢，而蔓延度指數、景觀連通度指數則呈下降趨勢。這說明了南京市近30 a的時間里斑塊數目不斷增加，主要景觀優勢度不斷下降，尤其是以耕地為主導的景觀格局正在發生深刻變化，破碎化程度越來越高，斑塊形狀趨于簡單化，以上種種情形表明人類對南京市景觀格局的干擾不斷增強，景觀格局趨向復雜化，景觀結構不穩定性增加，自我調節能力不斷減弱。目前來看城市化過程依然較快，這對于城市可持續發展和生物多樣性保護應引起足夠的重視。

表4 南京市2002-2013年土地利用景觀轉移矩陣 hm2

表5 南京市不同時期景觀指數

3.4.2 斑塊類型上的景觀格局變化分析 在斑塊類型尺度上，從斑塊數來看，1986—2013年耕地、林地和草地一直呈增長趨勢，居民點及建設用地則不斷下降，斑塊數常用來描述整個景觀的異質性，與景觀的破碎度具有較好的正相關性，由此可見耕地、林地和草地的破碎度在增大，而居民點及建設用地則不斷聚集；從平均斑塊面積來看，耕地、林地和草地一直呈減小趨勢，居民點及建設用地則不斷上升，與斑塊數結合起來看，耕地、林地和草地在破碎的同時平均斑塊面積減小，而居民點及建設用地通過不斷擴張聚集膨脹。周長—面積分維數可以反映景觀形狀的復雜程度，體現人類活動的大小，其值越接近于1表示形狀越簡單，越接近于2形狀越復雜，從表6中可以看到，耕地和草地一直呈上升趨勢，林地呈“V”形變化，居民點及建設用地呈倒“V”形變化，水域和其他用地呈波動增加，由此說明人類對耕地和草地的干預越來越大，而林地的變化反映了南京市打造國際性人文綠都等政策規劃的影響。

3.5 景觀格局驅動力因素分析

與自然因素不同，人類活動是影響城市土地利用/覆被景觀格局變化的重要因素，隨著城市空間不斷擴張，越來越多的自然景觀轉為人為景觀。從1986到2013年，南京市總人口由471.61萬增長到643.09萬，凈增171.48萬人，平均年增加6.35萬人；GDP從95.46億元增長到8 011.78億元，增幅達8 292%，此外，固定資產投資、財政收入、農民人均純收入、城鎮居民人均可支配收入增幅分別為18 743%，10 344%，2 678%和3 800%。除了人口及經濟發展外，區域政策和城市規劃同樣起到了重要作用。南京市2007—2020年規劃當中對城市性質最高定位為國家重要的區域中心城市，城市中心范圍不斷擴大，同時，不斷撤縣并區，提高城市化水平。同時，“綠色南京”工程、生態旅游建設、高淳擴大水產養殖等政策對于景觀格局變化起到了很大的推動作用。

表6 南京市不同時期不同土地利用類型景觀指數

4 結 論

從土地利用/覆被結構變化來看，南京市近30 a土地覆被類型以耕地為主，雖然面積不斷減少，但仍然是該區域的優勢景觀；林地在一系列綠色生態工程推動下小幅增長；居民點及建設用地呈快速增長趨勢，大量自然景觀向人為景觀轉變是南京市景觀格局變化最顯著的特征。

從土地利用/覆被動態變化度來看，整體景觀類型變化速度都在增大，而草地和其他用地變化最為活躍，林地相對最穩定；除第一時段內變化極緩慢之外，之后兩個時段綜合動態度均大于24%，已超過急劇變化型，說明南京市近20 a來土地利用/覆被變化是非常劇烈的。從土地利用/覆被景觀格局變化來看，在斑塊類型上，耕地和草地的平均斑塊面積減少，破碎度在增大，受人類干預強烈，而居民點及建設用地通過不斷擴張聚集膨脹；在景觀水平上，受人類活動干擾作用，斑塊數目不斷增加，主要景觀優勢度不斷下降，破碎化程度越來越高，斑塊形狀趨于簡單化，景觀格局趨向復雜化，景觀結構不穩定性增加。

從景觀格局變化的驅動力因素來看，人口增長、經濟發展、區域政策和城市規劃對于景觀格局變化起到了很大的推動作用。

[1] Deng J S, Wang K, Hong Y, et al. Spatio-temporal dynamics and evolution of land use change and landscape pattern in response to rapid urbanization[J]. Landscape and Urban Planning, 2009,92(3):187-198.

[2] 盛莉.快速城市化背景下城市熱島對土地覆蓋及其變化的響應關系研究[D].杭州:浙江大學,2013.

[3] 彭保發,石憶邵,王賀封,等.城市熱島效應的影響機理及其作用規律:以上海市為例[J].地理學報,2013,68(11):1461-1471.

[4] Li Z, Liu W, Zhang X, et al. Impacts of land use change and climate variability on hydrology in an agricultural catchment on the Loess Plateau of China[J]. Journal of Hydrology, 2009,377(1):35-42.

[5] 李志,劉文兆,張勛昌,等.未來氣候變化對黃土高原黑河流域水資源的影響[J].生態學報,2009,29(7):3456-3464.

[6] 李志,劉文兆,王秋賢.黃土塬區不同地形部位和土地利用方式對土壤物理性質的影響[J].應用生態學報,2008,19(6):1303-1308.

[7] 吳文斌,楊鵬,唐華俊,等.土地利用對土壤性質影響的區域差異研究[J].中國農業科學,2007,40(8):1697-1702.

[8] 高中貴,彭補拙,喻建華,等.經濟發達區土地利用變化對土壤性質的影響:以江蘇省昆山市為例[J].自然資源學報,2005,20(1):44-51.

[9] 彭澄瑤.城市水資源可持續規劃與水生態環境修復[D].北京:北京工業大學,2011.

[10] Meffert P J, Dziock F. The influence of urbanisation on diversity and trait composition of birds[J]. Landscape Ecology, 2013,28(5):943-957.

[11] 劉吉平,趙丹丹,田學智,等.1954—2010年三江平原土地利用景觀格局動態變化及驅動力[J].生態學報,2014,34(12):3234-3244.

[12] Wood R, Handley J. Landscape dynamics and the management of change[J]. Landscape Research, 2001,26(1):45-54.

[13] 趙陽,余新曉,賈劍波,等.紅門川流域土地利用景觀動態演變及驅動力分析[J].農業工程學報,2013,29(9):239-248.

[14] 李保杰,顧和和.礦區土地復墾景觀格局變化和生態效應[J].農業工程學報,2012,28(3):251-256.

[15] 謝花林.典型農牧交錯區土地利用變化驅動力分析[J].農業工程學報,2008,24(10):56-62.

[16] 賈寶全,王成,邱爾發.南京市景觀時空動態變化及其驅動力[J].生態學報,2013,33(18):5848-5857.

[17] 王慧.南京市土地利用景觀格局變化及預測[D].南京:南京農業大學,2010.

[18] 王秀蘭,包玉海.土地利用動態變化研究方法探討[J].地理科學進展,1999,18(1):81-87.

[19] 劉紀遠,布和敖斯爾.中國土地利用變化現代過程時空特征的研究:基于衛星遙感數據[J].第四紀研究,2000,20(3):229-239.

[20] 張忠輝,楊雨春,謝朋,等.松原市近20年土地利用景觀格局動態變化[J].中國農學通報,2014,30(2):222-226.

[21] 劉紀遠,張增祥,莊大方.20世紀90年代中國土地利用變化的遙感時空信息研究[M].北京：科學出版社,2005.

[22] 古琳,劉波,龔固堂,等.成都市近20年林地景觀變化特征[J].應用生態學報,2010,21(5):1081-1089.

LandUseandLandscapePatternChangesandtheDrivingForceFactorsinNanjingFrom1986to2013

TONG Guangchen1, LIN Jie1, CHEN Hang2, GU Zheyan2, TANG Peng1, ZHANG Jinchi1

(1.CollegeofForestry,CollaborativeInnovationCenterofSustainableForestryinSouthernChinaofJiangsuProvince,NanjingForestryUniversity,KeyLaboratoryofSoilandWaterConservationandEcologicalRestorationofJiangsuProvince,Nanjing210037,China; 2.JiangsuSurveyingandDesignofWaterResourcesCo.,Ltd.,Yangzhou，Jiangsu225127,China)

Taking Nanjing City as the research object, we used the land use/land cover data of 1986, 1996 and 2002 and 2013 as the data source. Based on the theory of landscape ecology, the support of ArcGIS and Fragstats software, land use/land cover and landscape pattern spatial and temporal variation and driving factors in the recent 30 years of Nanjing were analyzed in terms of the land use type structure, speed change, transformation relation and landscape pattern. The results show that: (1) the land use/cover has changed dramatically in Nanjing City in recent 30 years, cultivated land plays the leading role in the whole study area, the area ratio was above 48%, but its advantage gradually reduces, a large number of natural landscape converted to man-made landscape features is the most significant change in the research area; (2) the change speed of overall landscape types increased, while grassland and other land use changes were the most active and forest land remained the most stable, the land use structure was at the rapid adjustment stage; (3) the transformation of cultivated land mainly to residential and construction land was significant; (4) landscape pattern changed significantly, the degree of fragmentation was more and more high, landscape pattern tended to be complicated, and landscape structure instability increased; (5) population growth, economic development, regional policy and city planning promoted the change of landscape pattern together. These results of the study can provide an important reference for the scientific planning and sustainable development of Nanjing City.

land use; landscape; driving force factor; Nanjing City

2016-02-25

：2016-03-25

國家自然科學基金(31200534);江蘇高校優勢學科建設工程資助項目(PAPD)和2014年江蘇省普通高校學術學位研究生科研創新計劃項目(KYLX_0903)

佟光臣(1987—),男,江蘇邳州人,碩士研究生,研究方向:土壤侵蝕遙感監測。E-mail:969714589@qq.com

林杰(1976—),女,遼寧丹東人,副教授,碩士生導師,研究方向:土壤侵蝕遙感監測。E-mail:jielin@njfu.edu.cn

F301； P901

：A

：1005-3409(2017)02-0240-06