鄧志位 全斌 謝紅群 張海波

引用格式：鄧志位，全斌，謝紅群，等.基于強(qiáng)度分析與動(dòng)態(tài)度模型的長(zhǎng)株潭都市圈土地利用變化特征分析[J].國(guó)土資源導(dǎo)刊，2023，20（04）：104-113.

Reference format：Deng Zhiwei，Quan Bin，Xie Hongqun，et al.Analysis of land use change characteristics in the Changsha-Zhuzhou-Xiangtan metropolitan area based on intensity analysis and dynamic degree model[J].Land & Resources Herald，2023，20（04）：104-113.

摘 要：土地利用變化信息的測(cè)度是探索土地景觀動(dòng)態(tài)特征及其趨勢(shì)和土地資源管理的重要基礎(chǔ)，但傳統(tǒng)的定量方法與思路尚未對(duì)用地類型動(dòng)態(tài)轉(zhuǎn)換信息進(jìn)行深入挖掘，難以系統(tǒng)提供多層次的過程理解。本研究基于地理信息系統(tǒng)（GIS）工具，采用強(qiáng)度分析與動(dòng)態(tài)度模型，探討長(zhǎng)株潭都市圈土地利用時(shí)空動(dòng)態(tài)特征。結(jié)果表明，1980—2020年間，長(zhǎng)株潭都市圈土地利用整體變化加速，其變化強(qiáng)度擴(kuò)增16倍。建設(shè)用地的轉(zhuǎn)入與耕地的轉(zhuǎn)出表現(xiàn)為穩(wěn)定的活躍。建設(shè)用地持續(xù)增加主要來源于耕地和林地。建設(shè)用地從林地處獲得較大面積，因此建設(shè)用地的面積穩(wěn)定地趨向耕地，而穩(wěn)定地避免林地。近40年來長(zhǎng)株潭都市圈建設(shè)用地增長(zhǎng)突出，致使生態(tài)、生產(chǎn)土地縮減明顯，該演化過程主要受社會(huì)經(jīng)濟(jì)因素的驅(qū)動(dòng)，未來該區(qū)應(yīng)堅(jiān)持生態(tài)優(yōu)先、綠色發(fā)展的理念。此外，本研究闡釋了強(qiáng)度分析與其他量化手段的分析思想，可用于其他區(qū)域土地利用變化探測(cè)，以幫助決策者制定更有針對(duì)性的可持續(xù)土地管理策略。

關(guān)鍵詞：土地利用變化；強(qiáng)度分析；土地利用動(dòng)態(tài)度；長(zhǎng)株潭都市圈

中圖分類號(hào)：F301.24? ? ? ? 文獻(xiàn)標(biāo)志碼：A? ? ? ? ? 文章編號(hào)：1672-5603（2023）01-104-10

Analysis of Land Use Change Characteristics in the Changsha-Zhuzhou-Xiangtan Metropolitan Area Based on Intensity Analysis and Dynamic Degree Model

Deng Zhiwei，Quan Bin，Xie Hongqun，Zhang Haibo

（College of Geography and Tourism， Hengyang Normal University， Hengyang Hunan 421002）

Abstract： The measurement of land use change information is an important foundation for exploring the dynamic characteristics and trends of land landscape and land resource management. However， traditional quantitative methods and approaches have not yet deeply explored the dynamic transformation information of land use types， making it difficult to provide a systematic understanding of multi-level processes. This study is based on geographic information system （GIS） tools and uses intensity analysis and dynamic degree models to explore the spatiotemporal dynamic characteristics of land use in the Changsha Zhuzhou Xiangtan metropolitan area. The results indicate that from 1980 to 2020， the overall land use change in the Changsha Zhuzhou Xiangtan metropolitan area accelerated， with a 16 times increase in intensity. The transfer of construction land and the transfer of cultivated land show stable and active behavior. The continuous increase in construction land mainly comes from cultivated land and forest land. Although construction land obtains a larger area from forest land， the increase in construction land steadily tends towards cultivated land， while steadily avoiding forest land. In the past 40 years， the growth of construction land in the Changsha Zhuzhou Xiangtan metropolitan area has been prominent， resulting in a significant reduction in ecological and production land. This evolutionary process is mainly driven by socio-economic factors. In the future， the area should adhere to ecological priority and green development. In addition， this study elucidates the analytical ideas of intensity analysis and other quantitative methods， which can be used for detecting land use changes in other regions to help decision-makers develop more targeted sustainable land management strategies.

Keywords： land use change; strength analysis; land use dynamic degree; Changsha- Zhuzhou-Xiangtan metropolitan area

0 引言

快速的城鎮(zhèn)化伴隨著社會(huì)經(jīng)濟(jì)的蓬勃發(fā)展，同時(shí)也加快了地球表面土地利用演化進(jìn)程[1]。土地利用變化作為反映人類活動(dòng)對(duì)陸地表層系統(tǒng)作用的重要因素，被認(rèn)為是導(dǎo)致全球環(huán)境變化的主要決定因素之一，對(duì)生態(tài)系統(tǒng)、全球生物地球化學(xué)循環(huán)、氣候變化和生物多樣性等產(chǎn)生重大影響[2]。同時(shí)，城市地區(qū)及其居民由于這些影響而面臨嚴(yán)重的風(fēng)險(xiǎn)?？茖W(xué)評(píng)估土地利用變化時(shí)空動(dòng)態(tài)有助于揭示人類活動(dòng)與自然環(huán)境的交互機(jī)制[3]。

轉(zhuǎn)移矩陣作為表征土地變化信息常用的量化途徑，用于分析某一時(shí)段內(nèi)土地利用類別之間的轉(zhuǎn)換規(guī)模差異[4]。然而，一方面?zhèn)鹘y(tǒng)分析思路通常關(guān)注轉(zhuǎn)換規(guī)模大小，往往闡述規(guī)模較大的轉(zhuǎn)換過程，而忽視了規(guī)模較小的轉(zhuǎn)換信息單元。分析轉(zhuǎn)換規(guī)模思路固然直觀，但難以表達(dá)變化面積在所涉類別自身面積中的比率，致使未能有效描述用地類型動(dòng)態(tài)變化的劇烈程度。另一方面，該方法應(yīng)用多集中在對(duì)單一時(shí)期下不同土地利用類型的轉(zhuǎn)換過程以及轉(zhuǎn)入轉(zhuǎn)出狀況，缺乏對(duì)多個(gè)時(shí)段轉(zhuǎn)移矩陣信息進(jìn)行整合并深入挖掘，導(dǎo)致無法提供更多的數(shù)值信息以供研究者分析連續(xù)多個(gè)時(shí)期形成的土地利用演化規(guī)律。這在一定程度上限制了對(duì)土地利用變化時(shí)序穩(wěn)定性等特征的進(jìn)一步理解。為此，Huang等[5]提出強(qiáng)度分析方法用于探測(cè)土地利用動(dòng)態(tài)中哪些變化、轉(zhuǎn)換過程更為劇烈。目前，該方法已被廣泛應(yīng)用于“三生”用地[6]，流域景觀動(dòng)態(tài)[7]、城鎮(zhèn)擴(kuò)張[8]、草原沙漠化[9]等領(lǐng)域。

在2022年2月，長(zhǎng)株潭都市圈發(fā)展規(guī)劃被國(guó)家批復(fù)同意，這意味著長(zhǎng)株潭都市圈成為我國(guó)第四個(gè)國(guó)家級(jí)都市圈[10]。長(zhǎng)株潭作為全國(guó)“資源節(jié)約型和環(huán)境友好型”社會(huì)典型示范區(qū)，近幾十年來，該區(qū)大量鄉(xiāng)村人口遷入城市，高強(qiáng)度的人類活動(dòng)與土地利用開發(fā)引發(fā)了耕地流失、水資源短缺、空氣污染等系列問題[11]。一些學(xué)者在長(zhǎng)株潭城市群、長(zhǎng)沙中心城區(qū)尺度上開展了富有成效的土地變化相關(guān)研究[12]，主要圍繞在生態(tài)系統(tǒng)服務(wù)評(píng)價(jià)[13]、土地利用模擬[14]、城市增長(zhǎng)及其邊界劃設(shè)[15]、生態(tài)網(wǎng)絡(luò)構(gòu)建[16]等方面。然而，在都市圈尺度上，開展長(zhǎng)時(shí)間序列土地利用變化格局與過程的系統(tǒng)定量分析尚不多見。國(guó)家設(shè)立長(zhǎng)株潭都市圈旨在加快未來長(zhǎng)株潭一體化、同城化的步伐，這不可避免地會(huì)加速城鎮(zhèn)空間的增長(zhǎng)，進(jìn)而對(duì)當(dāng)?shù)厣鷳B(tài)環(huán)境承載帶來巨大壓力。鑒于此，本研究以最新規(guī)劃文本中長(zhǎng)株潭都市圈為研究區(qū)，借助強(qiáng)度分析與動(dòng)態(tài)度模型，探討1980—2020年間長(zhǎng)株潭都市圈土地利用變化規(guī)模、強(qiáng)度以及穩(wěn)定性特征，揭示該區(qū)城鎮(zhèn)擴(kuò)張過程對(duì)土地景觀格局變化的影響機(jī)制，為該都市圈未來實(shí)現(xiàn)高質(zhì)量發(fā)展和城市可持續(xù)發(fā)展提供參考依據(jù)。

1 研究區(qū)概況

長(zhǎng)株潭都市圈地處111°53′32″～114°15′28″ E和27°13′31″～28°39′56″ N之間，位于湖南省中北部，作為長(zhǎng)江中游城市群重要組成部分，是推動(dòng)長(zhǎng)江經(jīng)濟(jì)帶協(xié)同發(fā)展的核心力量之一[10]。本研究以《長(zhǎng)株潭都市圈發(fā)展規(guī)劃》（以下簡(jiǎn)稱《規(guī)劃》）中所劃定的長(zhǎng)株潭都市圈為研究區(qū)，包括長(zhǎng)沙市全域、株洲市中心城區(qū)及醴陵市、湘潭市中心城區(qū)及韶山市和湘潭縣，面積1.89萬km2（圖1）。2021年，該區(qū)常住人口1 484萬，城鎮(zhèn)化率達(dá)到80.9%，經(jīng)濟(jì)總量1.79萬億元，占全省經(jīng)濟(jì)總量的40%[17]。此外，該區(qū)交通區(qū)位優(yōu)勢(shì)凸顯，京廣、滬昆高鐵等多條交通路網(wǎng)交匯，基本構(gòu)成“一縱一橫”交通網(wǎng)格局。

2 數(shù)據(jù)與方法

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

土地利用數(shù)據(jù)來源于中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http：//www.resdc.cn/）。本研究采用5期土地利用數(shù)據(jù)，時(shí)間為1980、1990、2000、2010、2020年。該數(shù)據(jù)基于Landsat TM/ETM+/OLI多時(shí)相影像資料，通過人機(jī)交互解譯而成，空間分辨率30 m，解譯精度均在85%以上[18]。區(qū)域行政邊界、鐵路與公路等交通基礎(chǔ)設(shè)施、政府駐地?cái)?shù)據(jù)來自國(guó)家基礎(chǔ)地理信息中心的全國(guó)地理信息資源目錄服務(wù)系統(tǒng)（http：//www.webmap.cn/）。30 m分辨率DEM數(shù)據(jù)來自地理空間數(shù)據(jù)云中ASTER GDEM數(shù)據(jù)產(chǎn)品（http：//www.gscloud.cn/）。土地利用分類參考最新修訂的《土地利用現(xiàn)狀分類》文件，并根據(jù)研究區(qū)實(shí)地調(diào)查，一級(jí)地類為耕地、林地、草地、水域和建設(shè)用地。

2.2 研究方法

2.2.1 強(qiáng)度分析

強(qiáng)度分析方法以轉(zhuǎn)移矩陣為基礎(chǔ)，由自上而下的間隔層、類別層和轉(zhuǎn)換層構(gòu)成[19]。間隔層也被稱為“時(shí)間間隔層”，用于分析哪一時(shí)間間隔（時(shí)期）土地利用變化得相對(duì)快速或緩慢。用地類別層次是進(jìn)一步解析某一時(shí)間間隔下哪些類別的增加（轉(zhuǎn)入）和減少（轉(zhuǎn)出）是相對(duì)活躍或沉寂的，而用地轉(zhuǎn)換層次則更深入地剖析其他地類轉(zhuǎn)換為某一特定地類是趨向或避免的。本研究采用該方法分析研究區(qū)土地利用動(dòng)態(tài)的變化強(qiáng)度特征。圖2顯示了土地利用轉(zhuǎn)移矩陣的數(shù)學(xué)表達(dá)，其中J為用地類型總數(shù)；Ctij則表示t時(shí)段地類i轉(zhuǎn)換為地類j的面積規(guī)模，符號(hào)釋義下同。

時(shí)間間隔層次研究不同時(shí)段下的總體變化差異，計(jì)算不同時(shí)期的總體變化強(qiáng)度以及整個(gè)研究時(shí)段的平均變化強(qiáng)度，計(jì)算分別見公式（1）和公式（2）[12]。通過比較St與U的大小判斷不同時(shí)期總體變化快慢。若St > U，判斷該時(shí)期的變化是快速的；若St < U，判斷在t時(shí)段內(nèi)的變化是緩慢的。

式中，St為t時(shí)段土地利用總體變化強(qiáng)度；U為t時(shí)段內(nèi)St的平均強(qiáng)度；t為時(shí)段步長(zhǎng)；T為整個(gè)研究時(shí)段長(zhǎng)度，符號(hào)釋義下同。

類別層次研究各地類變化的差異，計(jì)算地類的轉(zhuǎn)出強(qiáng)度與轉(zhuǎn)入強(qiáng)度。通過比較Lti、Gtj和St的大小判斷地類的轉(zhuǎn)入、轉(zhuǎn)出是活躍或沉寂的，計(jì)算見公式（3）和公式（4）[12]。若Gtj > St，判斷類別j的轉(zhuǎn)入是活躍的；若Gtj < St，判斷類別j的轉(zhuǎn)入是沉寂的；對(duì)地類轉(zhuǎn)出的判斷亦同。

式中，Gtj為t時(shí)段地類j的轉(zhuǎn)入強(qiáng)度；Lti為t時(shí)段內(nèi)地類i的轉(zhuǎn)出強(qiáng)度。

轉(zhuǎn)換層次研究其他地類轉(zhuǎn)換為某一地類時(shí)的轉(zhuǎn)換差異，計(jì)算其他地類轉(zhuǎn)換為某一地類的強(qiáng)度以及其他地類轉(zhuǎn)換為該地類的平均轉(zhuǎn)換強(qiáng)度，計(jì)算見公式（5）和公式（6）[12]。若Rtin > Wtn，判斷地類i趨向轉(zhuǎn)換為地類n，換言之，地類n是趨向從一種地類i處獲得增加；同理，若Rtin < Wtn，判斷地類i避免轉(zhuǎn)換為地類n。在該層次中，本研究的地類n選取了建設(shè)用地，旨在更進(jìn)一步理解建設(shè)用地動(dòng)態(tài)轉(zhuǎn)換過程。

式中，Rtin為t時(shí)段地類i轉(zhuǎn)換為特定地類n的轉(zhuǎn)換強(qiáng)度；Wtn為t時(shí)段其他地類轉(zhuǎn)換為特定地類n的平均轉(zhuǎn)換強(qiáng)度；地類n代表本研究中建設(shè)用地；Ctin為t時(shí)段地類i轉(zhuǎn)換為特定地類n的轉(zhuǎn)換規(guī)模面積。

2.2.2 動(dòng)態(tài)度模型

動(dòng)態(tài)度模型是測(cè)算土地利用/覆被變化劇烈程度的常用指標(biāo)，可分為單一動(dòng)態(tài)度與綜合動(dòng)態(tài)度。單一動(dòng)態(tài)度通常被用于描述不同地類在某一時(shí)段內(nèi)的變化劇烈程度，見公式（7）[20]。綜合動(dòng)態(tài)度旨在計(jì)算所有地類的總體變化劇烈程度，見公式（8）[20]。

式中，Ki為地類i的單一動(dòng)態(tài)度值；Uib為地類i在t時(shí)段末的面積；Uia為地類i在t時(shí)段初的面積；LCt是t時(shí)段內(nèi)所有地類的綜合動(dòng)態(tài)度值；[LUij]為該時(shí)段內(nèi)地類i的轉(zhuǎn)出面積；LUi代表地類i在該時(shí)段內(nèi)初時(shí)的面積。

3 結(jié)果與分析

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

圖3顯示長(zhǎng)株潭都市圈土地利用結(jié)構(gòu)及其空間格局。從圖3可見，在1980—2020年間，長(zhǎng)株潭都市圈土地利用類型面積最大為林地，耕地次之，兩者約占研究區(qū)域90%。這表明該區(qū)生態(tài)資源分布廣泛，擁有較好的資源稟賦。建設(shè)用地面積增長(zhǎng)持續(xù)上升，占比從1980年的2.20%增加到7.52%；與之相反，耕地和林地從1980年的33.3%、61.5%縮減至2020年30.0%、59.3%。通過鷹眼放大觀察（圖4），可見湘潭市與醴陵市兩處建設(shè)用地迅速向外擴(kuò)張與蔓延，不同程度地侵占周邊耕地與林地等生態(tài)資源，說明盡管建設(shè)用地初始時(shí)比重小，但其數(shù)量增長(zhǎng)是該區(qū)最明顯的特征。

3.2 土地利用變化強(qiáng)度分析

3.2.1 土地利用變化時(shí)序分析

表1顯示了研究區(qū)四個(gè)時(shí)期土地利用總體變化的規(guī)模、強(qiáng)度以及動(dòng)態(tài)度值。從表1可見，無論是變化規(guī)模、變化強(qiáng)度還是綜合動(dòng)態(tài)度，都表現(xiàn)出上升趨勢(shì)。從變化強(qiáng)度可見，由于四個(gè)時(shí)期的變化強(qiáng)度都不等于整個(gè)研究時(shí)期的平均強(qiáng)度（0.26%），長(zhǎng)株潭都市圈在1980—2020年間土地利用總體變化具有不穩(wěn)定特征。值得注意的是，自2000年后，該區(qū)變化規(guī)模與強(qiáng)度出現(xiàn)突增，從1980—1990年的年均變化規(guī)模6.32 km2增加到2010—2020年的103.90 km2，說明前一個(gè)時(shí)期的總體變化相對(duì)于后一個(gè)時(shí)期是相對(duì)緩慢的。綜合動(dòng)態(tài)度突增形勢(shì)基本與變化強(qiáng)度保持一致。

3.2.2 土地利用類別變化分析

表2和表3分別顯示1980—2000年與2000—2020年長(zhǎng)株潭都市圈地類的年均變化規(guī)模、強(qiáng)度、單一動(dòng)態(tài)度。從表2與表3可見，在1980—2020年間的4個(gè)時(shí)期，建設(shè)用地年均轉(zhuǎn)入面積分別為3.53 km2、7.55 km2、51.55 km2和46.78 km2，其轉(zhuǎn)入強(qiáng)度遠(yuǎn)超平均強(qiáng)度，且在前3個(gè)時(shí)期呈現(xiàn)出遞增趨勢(shì)。此外，建設(shè)用地的轉(zhuǎn)出規(guī)模與強(qiáng)度均遠(yuǎn)低于其轉(zhuǎn)入，進(jìn)一步反映了建設(shè)用地處于相對(duì)活躍的數(shù)量增長(zhǎng)狀態(tài)。與之相反，林地的轉(zhuǎn)入與轉(zhuǎn)出量較大，但由于林地的轉(zhuǎn)入與轉(zhuǎn)出強(qiáng)度在4個(gè)時(shí)期里均小于平均強(qiáng)度，因此表現(xiàn)為穩(wěn)定沉寂。林地在整個(gè)研究期間轉(zhuǎn)出大于轉(zhuǎn)入，表現(xiàn)為凈減少。耕地的轉(zhuǎn)出強(qiáng)度始終大于轉(zhuǎn)入強(qiáng)度，并且由于耕地轉(zhuǎn)出強(qiáng)度始終大于平均強(qiáng)度，表現(xiàn)為穩(wěn)定活躍，相反的是轉(zhuǎn)入強(qiáng)度一直小于平均強(qiáng)度，因此耕地轉(zhuǎn)入表現(xiàn)為穩(wěn)定沉寂。耕地轉(zhuǎn)入規(guī)模出現(xiàn)逐漸增加的趨勢(shì)，這與自21世紀(jì)以來長(zhǎng)株潭都市圈耕地占補(bǔ)平衡政策的有力執(zhí)行和土地可持續(xù)發(fā)展的逐漸普及有關(guān)。水域表現(xiàn)為凈增加，其轉(zhuǎn)入與轉(zhuǎn)出強(qiáng)度均表現(xiàn)為穩(wěn)定活躍。對(duì)于草地而言，除在2010—2020年間轉(zhuǎn)入強(qiáng)度小于平均強(qiáng)度以外，其余為大于平均強(qiáng)度，表現(xiàn)為活躍。

從單一動(dòng)態(tài)度來看，在1980—2020年間，耕地與林地出現(xiàn)凈流失，呈現(xiàn)為負(fù)值，而建設(shè)用地與水域則表現(xiàn)出正值，反映建設(shè)用地與水域面積凈增加。建設(shè)用地第四個(gè)時(shí)期時(shí)動(dòng)態(tài)度值相較于第一個(gè)時(shí)期擴(kuò)大約5倍。值得注意的是，建設(shè)用地與耕地的動(dòng)態(tài)度絕對(duì)值呈現(xiàn)“先增后降”趨勢(shì)，在2000—2010年達(dá)最大值9.17%、0.47%。在整個(gè)研究時(shí)段，草地的動(dòng)態(tài)度值正負(fù)反復(fù)。

3.2.3基于建設(shè)用地增長(zhǎng)的轉(zhuǎn)換過程分析

圖4顯示長(zhǎng)株潭都市圈4個(gè)時(shí)期下其他地類轉(zhuǎn)換為建設(shè)用地的轉(zhuǎn)入規(guī)模與強(qiáng)度。從圖4左側(cè)的年均轉(zhuǎn)換規(guī)模可見，在這4個(gè)時(shí)期耕地年均轉(zhuǎn)入建設(shè)用地面積分別為1.94 km2、4.00 km2、26.73 km2、24.93 km2，表現(xiàn)出“先升后降”的趨勢(shì)。由于耕地轉(zhuǎn)換強(qiáng)度大于均衡強(qiáng)度，耕地穩(wěn)定地轉(zhuǎn)換為建設(shè)用地。相反，由于林地的轉(zhuǎn)換強(qiáng)度在4個(gè)時(shí)期都小于均衡強(qiáng)度，林地穩(wěn)定地避免轉(zhuǎn)換為建設(shè)用地。但是，林地的轉(zhuǎn)換面積規(guī)模是巨大的，僅次于耕地，在這4個(gè)時(shí)期年均轉(zhuǎn)入建設(shè)用地面積分別為1.59 km2、3.11 km2、23.10 km2、20.84 km2。此外，草地穩(wěn)定地避免轉(zhuǎn)換為建設(shè)用地。水域的轉(zhuǎn)換強(qiáng)度表明水域經(jīng)歷了“輕微變化-劇烈變化-輕微變化”的轉(zhuǎn)換過程。某一地類轉(zhuǎn)入的均衡轉(zhuǎn)換強(qiáng)度可描述該地類擴(kuò)張的劇烈程度。因此，該區(qū)建設(shè)用地的擴(kuò)張均衡強(qiáng)度變化表現(xiàn)為先增加、后減少的狀態(tài)。

4 討論

4.1 量化土地利用/覆被動(dòng)態(tài)信息

土地變化劇烈程度作為土地變化檢測(cè)和評(píng)價(jià)的重要信號(hào)之一，可以刻畫土地類別空間范圍動(dòng)態(tài)過程和類別轉(zhuǎn)換的激烈程度[6]。多數(shù)關(guān)于遙感影像分類的土地利用變化探測(cè)的方法與指數(shù)，都能在不同程度上反映土地變化劇烈程度，譬如轉(zhuǎn)移矩陣、統(tǒng)計(jì)地類面積變化、單一與綜合土地利用動(dòng)態(tài)度以及強(qiáng)度分析等。

轉(zhuǎn)移矩陣實(shí)質(zhì)上是通過不同時(shí)期的土地利用地圖交互疊置，進(jìn)而統(tǒng)計(jì)類別相互轉(zhuǎn)換面積數(shù)據(jù)而形成的特殊表格，可反映研究區(qū)內(nèi)類別間轉(zhuǎn)換規(guī)模以及土地流向特征，為國(guó)內(nèi)外廣泛采用[19]。此外，它也是馬爾科夫（Markov）模型計(jì)算地類轉(zhuǎn)移概率的重要基礎(chǔ)，有助于歷史趨勢(shì)情景下的土地變化預(yù)測(cè)分析。然而，一方面，轉(zhuǎn)移矩陣是以某一時(shí)段的表達(dá)形式出現(xiàn)，并不利于分析連續(xù)的不同時(shí)期下地類變化與轉(zhuǎn)移所形成的時(shí)序規(guī)律。另一方面，常見的一種思路是將類別間轉(zhuǎn)換規(guī)模作為變化劇烈程度的信號(hào)。例如，Zhang等[21]分析了衡陽市2010—2015年和2015—2018年土地利用轉(zhuǎn)移特征，發(fā)現(xiàn)林地流向建設(shè)用地規(guī)模最大。傳統(tǒng)的分析思路是聚焦轉(zhuǎn)換規(guī)模較大的動(dòng)態(tài)信息，而對(duì)規(guī)模較小的轉(zhuǎn)換過程重視不足。直接利用轉(zhuǎn)換規(guī)模進(jìn)行分析盡管直觀，但未能體現(xiàn)變化面積在所涉類別自身面積中的比率。相較于變化規(guī)模，變化強(qiáng)度（比率）更能體現(xiàn)變化劇烈程度[22]。若采用強(qiáng)度分析方法挖掘該研究中兩個(gè)時(shí)段的轉(zhuǎn)移矩陣結(jié)果，可發(fā)現(xiàn)耕地轉(zhuǎn)換為建設(shè)用地的劇烈程度遠(yuǎn)比林地轉(zhuǎn)換為建設(shè)用地的大，且建設(shè)用地的增加穩(wěn)定地避免林地。因此，大規(guī)模林地流向建設(shè)用地可解釋為林地的初始規(guī)模大，相比其他地類具備更大的可獲得性。

另一種量化思路是按類別統(tǒng)計(jì)凈變化面積以及凈變化大小占初始時(shí)間點(diǎn)所涉類別的大小，后者即為單一動(dòng)態(tài)度的計(jì)算邏輯。單一動(dòng)態(tài)度能表征地類凈變化的正負(fù)之分，即凈減少或凈增加，該指數(shù)未能同時(shí)體現(xiàn)在空間格局上同時(shí)發(fā)生的地類轉(zhuǎn)入（增加）與轉(zhuǎn)出（減少）。從圖5可見，在某一時(shí)段內(nèi)，地類增加的部分與其初期的面積大小和空間分布并無直接聯(lián)系，卻構(gòu)成其末期面積；換而言之，類別的減少面積必然是其初期面積的一部分。單一動(dòng)態(tài)度在刻畫地類動(dòng)態(tài)程度上存在對(duì)地類的分配變化的忽視。交換變化是指地類與其他地類在空間上發(fā)生區(qū)域的錯(cuò)位，即相互轉(zhuǎn)化。因此，地類的總變化為區(qū)域A與B面積之和，凈變化為B與A面積之差，而總變化減去凈變化即為交換變化部分。

研究人員運(yùn)用綜合動(dòng)態(tài)度旨在比較分析區(qū)域內(nèi)不同時(shí)期土地利用總體變化劇烈程度[20]。對(duì)比強(qiáng)度分析時(shí)間間隔層次公式（1）和綜合動(dòng)態(tài)度公式（8）以及分析表1，不難發(fā)現(xiàn)公式（8）為公式（1）的一半。若觀察轉(zhuǎn)移矩陣，對(duì)角線上的元素即為某一時(shí)期內(nèi)的不變區(qū)域，而其余部分應(yīng)為動(dòng)態(tài)區(qū)域。然而，綜合動(dòng)態(tài)度在分母上為研究區(qū)面積的2倍，嘗試將地類的轉(zhuǎn)出轉(zhuǎn)入進(jìn)行平均，這在一定程度上導(dǎo)致對(duì)區(qū)域內(nèi)土地總體變化的低估。

強(qiáng)度分析作為一種定量分析框架，具有分層遞進(jìn)式特征，可彌補(bǔ)傳統(tǒng)量化途徑的局限性，在間隔、類別和轉(zhuǎn)換層次上耦合土地變化的規(guī)模、強(qiáng)度以及穩(wěn)定性[7]。強(qiáng)度分析是從早期的測(cè)度單一間隔的方法演化而來的，并且那些早期版本已經(jīng)變得流行，早期方法成為系統(tǒng)性轉(zhuǎn)換分析方法[23]。然而，強(qiáng)度分析方法更易理解，鑒于它的層次清晰、圖形表達(dá)，特別是關(guān)于土地時(shí)間變化潛在原因有了更深刻的解釋[24]。例如，為什么地類X轉(zhuǎn)換為地類Y的規(guī)模大于其他地類轉(zhuǎn)換為地類Y的規(guī)模?？赡茉蛴卸浩湟?，若在初始時(shí)間點(diǎn)地類X的面積比其他地類的大，那么自然擁有更多面積用于轉(zhuǎn)換成Y。倘若地類X只以均衡轉(zhuǎn)換強(qiáng)度從其他地類獲得轉(zhuǎn)換土地時(shí)，地類Y也會(huì)從地類X處獲得更多的面積轉(zhuǎn)換；其二，若在初始時(shí)間點(diǎn)，地類Y從地類X處獲得轉(zhuǎn)換面積，相比從其他地類處擁有更大的強(qiáng)度來獲得轉(zhuǎn)換面積。有時(shí)這兩種情況同時(shí)存在。強(qiáng)度分析提供了這樣的一種途徑，方便學(xué)者從關(guān)注土地變化格局再到分析相關(guān)的轉(zhuǎn)換過程。為了更全面且直觀地反映轉(zhuǎn)換特征，有研究在強(qiáng)度分析的轉(zhuǎn)換層次結(jié)果基礎(chǔ)上相繼提出了“轉(zhuǎn)換圖譜”，以豐富強(qiáng)度分析方法論[7-8，22]。

4.2長(zhǎng)株潭都市圈土地利用變化及其過程

土地利用變化強(qiáng)度分析結(jié)果表明，2000年后的兩個(gè)時(shí)段內(nèi)的總體變化速率明顯加快（表1），這可能與后期人口增長(zhǎng)和經(jīng)濟(jì)發(fā)展具有密切聯(lián)系。對(duì)于類別層次（表2和表3），建設(shè)用地是最大的轉(zhuǎn)入地類，相反耕地是最大的轉(zhuǎn)出地類。林地減少面積較大，同時(shí)林地又經(jīng)歷著沉寂的減少。這一表征被描述為“地類沉寂現(xiàn)象”[25]。在轉(zhuǎn)換層次，建設(shè)用地增長(zhǎng)主要來源于耕地與林地。同時(shí)，耕地、林地的景觀破碎化加劇，間接反映耕地與林地減少且其斑塊趨于縮小趨勢(shì)，進(jìn)而加劇了耕地被侵占或撂荒的可能性。因?yàn)槊娣e規(guī)模較小的斑塊更易被領(lǐng)域內(nèi)的優(yōu)勢(shì)斑塊侵吞占用，也容易被土地監(jiān)管部門所忽視。耕地轉(zhuǎn)換為建設(shè)用地表現(xiàn)為穩(wěn)定的趨向（圖5），這一轉(zhuǎn)換過程的變化與該地區(qū)非農(nóng)經(jīng)濟(jì)活動(dòng)增長(zhǎng)可能有關(guān)。

圖6表明，自1984年以來，長(zhǎng)株潭地區(qū)的總?cè)丝诜€(wěn)步上升。2001—2019年期間城市人口的增長(zhǎng)速度是1984—2000年間的3倍。自21世紀(jì)以來，該地區(qū)呈現(xiàn)出更快速的經(jīng)濟(jì)發(fā)展。城市人口的急劇增加導(dǎo)致需要更多的住房、交通和其他基礎(chǔ)設(shè)施來滿足人類社會(huì)經(jīng)濟(jì)發(fā)展。此外，本研究發(fā)現(xiàn)其他用地轉(zhuǎn)入建設(shè)用地的強(qiáng)度在前三個(gè)時(shí)段均呈上升趨勢(shì)，但在第四個(gè)時(shí)段內(nèi)轉(zhuǎn)換強(qiáng)度稍有下降（表2、表3）。歐陽曉等[15]研究長(zhǎng)株潭城市群城市用地發(fā)現(xiàn)其擴(kuò)張強(qiáng)度呈現(xiàn)出先加速后減緩的趨勢(shì)，本研究結(jié)果與之類似。圖4顯示，在2000—2010年、2010—2020年間的建設(shè)用地是趨向轉(zhuǎn)換為耕地。在城市化背景下，建設(shè)用地幾乎難以轉(zhuǎn)換為耕地，故而發(fā)生這一過程通常被視為數(shù)據(jù)質(zhì)量問題。然而在我國(guó)，這與現(xiàn)行的耕地占補(bǔ)平衡制度息息相關(guān)。近20年來，耕地占補(bǔ)平衡制度歷經(jīng)了由國(guó)土部門注重?cái)?shù)量保護(hù)的起步階段發(fā)展到數(shù)量與質(zhì)量相耦合的發(fā)展階段，再到新型城鎮(zhèn)化背景下，數(shù)量、質(zhì)量、生態(tài)三位一體階段[26]。

長(zhǎng)株潭都市圈的設(shè)立旨在充分釋放長(zhǎng)沙、株洲、湘潭三個(gè)地區(qū)的潛能，優(yōu)勢(shì)互補(bǔ)[10]。相較于深圳、廣州、南京都市圈以及武漢都市圈的社會(huì)經(jīng)濟(jì)，長(zhǎng)株潭都市圈建設(shè)用地的適當(dāng)蔓延仍然是該區(qū)社會(huì)經(jīng)濟(jì)發(fā)展的內(nèi)在需求，但其增長(zhǎng)需符合生態(tài)文明建設(shè)與糧食安全保障的要求。城鎮(zhèn)擴(kuò)張不僅造成城市周邊耕地與林地急劇減少以及潛在負(fù)面環(huán)境影響，也威脅了地區(qū)生態(tài)安全。城市的可持續(xù)能力依賴于城市區(qū)域中各種生態(tài)要素和生態(tài)過程所組成的景觀整體的穩(wěn)定性與可持續(xù)性[16]。因此，未來研究綜合構(gòu)建都市圈生態(tài)網(wǎng)絡(luò)，耦合斑塊生成土地利用模擬模型，優(yōu)化該區(qū)城鎮(zhèn)增長(zhǎng)邊界，提出城鎮(zhèn)增長(zhǎng)管控策略。

5 結(jié)論

（1）長(zhǎng)株潭都市圈土地利用主要類型為林地、耕地，兩者共占研究區(qū)比重約90%。通過研究1980—2020年土地利用時(shí)序總體變化，長(zhǎng)株潭都市圈的變化區(qū)域持續(xù)上升，其變化強(qiáng)度增加16倍。這與幾十年的中國(guó)經(jīng)濟(jì)高速發(fā)展和快速城市化具有密切聯(lián)系。

（2）在類別層次，耕地的減少、建設(shè)用地與水域的轉(zhuǎn)入表現(xiàn)出穩(wěn)定性趨勢(shì)。相反，林地減少呈現(xiàn)為沉寂的，盡管其流失規(guī)模較大。在1980—2020年，建設(shè)用地的增加主要來源于耕地、林地，且建設(shè)用地的增加穩(wěn)定地趨向耕地，但穩(wěn)定地避免林地、草地。自2000年來，受“耕地占補(bǔ)平衡”政策影響，耕地的增加趨向建設(shè)用地。

（3）本研究以中部地區(qū)典型快速城市化區(qū)域?yàn)槔\(yùn)用強(qiáng)度分析與動(dòng)態(tài)度模型深入分析近40年土地利用演化信息與轉(zhuǎn)換規(guī)律，試圖揭示各模型在土地變化信息定量化表達(dá)上的能力與局限性，為進(jìn)一步挖掘土地變化信息提供有益思考。

參考文獻(xiàn)/References

[1]Aspinall R J， Staiano M， Pearson D M. Data， time， change and land-system dynamics[J]. Journal of Land

[2]劉紀(jì)遠(yuǎn)， 寧佳， 匡文慧， 等. 2010—2015年中國(guó)土地利用變化的時(shí)空格局與新特征[J]. 地理學(xué)報(bào)， 2018， 73（5）： 789-802.

[3]匡文慧， 張樹文， 杜國(guó)明， 等. 2015—2020年中國(guó)土地利用變化遙感制圖及時(shí)空特征分析[J]. 地理學(xué)報(bào)， 2022， 77（5）： 1056-1071.

[4]喬偉峰， 盛業(yè)華， 方斌， 等. 基于轉(zhuǎn)移矩陣的高度城市化區(qū)域土地利用演變信息挖掘——以江蘇省蘇州市為例[J]. 地理研究， 2013， 32（8）： 1497-1507.

[5]Huang J， Pontius Jr R G， Li Q， et al. Use of intensity analysis to link patterns with processes of land change from 1986 to 2007 in a coastal watershed of southeast China[J]. Applied Geography， 2012， 34： 371-384.

[6]鄧志位， 全斌， 祁劍青， 等. “三生”用地格局演變的強(qiáng)度分析及其驅(qū)動(dòng)力研究——以衡陽市為例[J]. 土壤通報(bào)， 2023， 54（4）： 768-780.

[7]王穎慧， 丁建麗， 李曉航， 等. 伊犁河流域LUCC對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響——基于強(qiáng)度分析模型[J]. 生態(tài)學(xué)報(bào)， 2022（8）： 1-13.

[8]楊建新， 龔健， 高靜， 等. 國(guó)家中心城市土地利用變化穩(wěn)定性和系統(tǒng)性特征——以武漢市為例[J]. 資源科學(xué)， 2019， 41（4）： 701-716.

[9]Meng X， Gao X， Li S， et al. Monitoring desertification in Mongolia based on Landsat images and Google Earth Engine from 1990 to 2020[J]. Ecological Indicators， 2021， 129： 107908.

[10]周國(guó)華， 李秋泓. 從長(zhǎng)株潭城市群到長(zhǎng)株潭都市圈——解讀《長(zhǎng)株潭都市圈發(fā)展規(guī)劃》[J]. 中國(guó)投資（中英文）， 2022（Z6）： 66-70.

[11]魏紫玲， 唐常春. 長(zhǎng)株潭都市圈國(guó)土空間格局演變特征研究[J]. 國(guó)土資源導(dǎo)刊， 2023， 20（1）： 14-21.

[12]Quan B， Ren H， Pontius Jr R G， et al. Quantifying spatiotemporal patterns concerning land change in Changsha， China[J]. Landscape and Ecological Engineering， 2018， 14（2）： 257-267.

[13]歐陽曉， 賀清云， 朱翔. 多情景下模擬城市群土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響——以長(zhǎng)株潭城市群為例[J]. 經(jīng)濟(jì)地理， 2020， 40（1）： 93-102.

[14]李欣， 王志遠(yuǎn)， 劉丹丹， 等. 生態(tài)安全格局約束下長(zhǎng)株潭都市圈建設(shè)用地演變模擬與管控[J]. 環(huán)境工程技術(shù)學(xué)報(bào)， 2023， 13（1）： 348-358.

[15]歐陽曉， 朱翔， 賀清云. 基于空間區(qū)位條件的長(zhǎng)株潭城市群城市用地?cái)U(kuò)張研究[J]. 經(jīng)濟(jì)地理， 2020， 40（5）： 183-190.

[16]Wang W， Li B， Su F， et al. Identifying ecological security patterns meeting future urban expansion in Changsha–Zhuzhou–Xiangtan urban agglomeration， China[J]. Remote Sensing， 2023， 15（12）： 3141-3161.

[17]湖南省統(tǒng)計(jì)局.湖南省統(tǒng)計(jì)年鑒2021[M]. 北京： 中國(guó)統(tǒng)計(jì)出版社， 2021.

[18]徐新良， 劉紀(jì)遠(yuǎn)， 張樹文， 等. 中國(guó)多時(shí)期土地利用土地覆被遙感監(jiān)測(cè)數(shù)據(jù)集（CNLUCC）[Z]. 北京： 中國(guó)科學(xué)院資源環(huán)境數(shù)據(jù)中心注冊(cè)與出版系統(tǒng)（http：//www.resdc.cn/DOI）.

[19]Quan B， Pontius Jr R G， Song H. Intensity analysis to communicate land change during three time intervals in two regions of Quanzhou City， China[J]. GIScience & Remote Sensing， 2020， 57（1）： 21-36.

[20]王秀蘭， 包玉海. 土地利用動(dòng)態(tài)變化研究方法探討[J]. 地理科學(xué)進(jìn)展， 1999（1）： 83-89.

[21]Zhang C， Wang P， Xiong P， et al. Spatial pattern simulation of land use based on FLUS model under ecological protection： A case study of Hengyang City[J]. Sustainability， 2021， 13（18）： 10458-10476.

[22]Deng Z， Quan B. Intensity Characteristics and multi-scenario projection of land use and land cover change in Hengyang， China[J]. International Journal of Environmental Research and Public Health， 2022， 19（14）： 8491-8508.

[23]Alo C A， Pontius Jr R G. Identifying systematic land-cover transitions using remote sensing and GIS： The fate of forests inside and outside protected areas of southwestern Ghana[J]. Environment and Planning B： Planning and Design， 2008， 35（2）： 280-295.

[24]Pontius Jr R G. Metrics that make a difference： How to analyze change and error[M]. Cham： Springer Nature Switzerland AG， 2022.

[25]Pontius Jr R G， Gao Y， Giner N， et al. Design and interpretation of intensity analysis illustrated by land change in central Kalimantan， Indonesia[J]. Land， 2013， 2（3）： 351-369.

[26]蔡漢， 朱權(quán)， 羅云建， 等. 快速城鎮(zhèn)化地區(qū)耕地景觀生態(tài)安全格局演變特征及其驅(qū)動(dòng)機(jī)制[J]. 南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2020， 44（5）： 181-188.