基于土地利用時空格局演變的生態(tài)足跡分析

摘要：自20世紀70年代至今，可持續(xù)發(fā)展已經(jīng)成為大部分國家共同達成的認識。中國將繼續(xù)走可持續(xù)發(fā)展道路，努力實現(xiàn)碳達峰和碳中和。本文以生態(tài)足跡分析法為主，結(jié)合土地利用時空格局演變，衡量濟南市2010—2020年可持續(xù)發(fā)展水平。其間利用ArcGIS 10.8軟件和TerrSet 2020軟件對濟南市2010—2020年的土地利用/覆被數(shù)據(jù)進行處理，得到土地利用轉(zhuǎn)移矩陣及土地利用類型增勢和減勢圖譜，進而分析研究期內(nèi)濟南市土地利用時空格局演變。根據(jù)濟南市2011年及2021年統(tǒng)計年鑒數(shù)據(jù)，計算2010年、2020年濟南市人均生態(tài)足跡和人均生態(tài)承載力。結(jié)果發(fā)現(xiàn)，2010—2020年，濟南市耕地出現(xiàn)大面積轉(zhuǎn)出，轉(zhuǎn)出面積為153.61 km2；濟南市2010年和2020年都出現(xiàn)生態(tài)赤字，2020年生態(tài)赤字相比2010年更大，濟南市可持續(xù)發(fā)展狀況亟待改善。

關鍵詞：土地利用變化；生態(tài)足跡；新舊動能轉(zhuǎn)換；濟南市

中圖分類號：X22 文獻標識碼：A 文章編號：1008-9500（2023）05-00-09

DOI：10.3969/j.issn.1008-9500.2023.05.017

Abstract： Since the 1970s， sustainable development has become a common understanding among most countries. China will continue to follow the path of sustainable development and strive to achieve carbon peak and carbon neutrality. Based on the ecological footprint analysis method， this paper combines the spatiotemporal pattern evolution of land use to measure the sustainable development level of Jinan City from 2010 to 2020. During this period， ArcGIS 10.8 software and TerraSet 2020 software are used to process the land use/cover data of Jinan City from 2010 to 2020， obtaining the land use transfer matrix and the map of land use type increase and decrease， and then analyzing the spatiotemporal pattern evolution of land use in Jinan City during the research period. Based on the 2011 and 2021 statistical yearbook data of Jinan City， the per capita ecological footprint and per capita ecological carrying capacity of Jinan City in 2010 and 2020 are calculated. The results show that from 2010 to 2020， a large area of arable land was transferred out of Jinan City， with an area of 153.61 km2; Jinan City experienced ecological deficits in both 2010 and 2020， compared to 2010， the ecological deficit in 2020 was even larger， and the sustainable development situation of Jinan City urgently needs improvement.

Keywords： land use change; ecological footprint; conversion of old and new kinetic energy; Jinan City

可持續(xù)發(fā)展是指既能滿足當代人的需求，又不對后代人滿足其需求的能力構(gòu)成威脅的發(fā)展[1]。從20世紀90年代可持續(xù)發(fā)展理論的誕生，到現(xiàn)如今的21世紀，可持續(xù)發(fā)展的內(nèi)涵被全球大多數(shù)國家認可，其被納入21世紀國家發(fā)展戰(zhàn)略規(guī)劃[2]。可持續(xù)發(fā)展正逐步滲透到全球各個國家的不同產(chǎn)業(yè)，學者將原本用來評價國家發(fā)展是否和諧的理論應用到產(chǎn)業(yè)發(fā)展是否和諧的問題解決上，都表明可持續(xù)發(fā)展的大勢所趨。我國對于可持續(xù)發(fā)展極為重視，已經(jīng)將其寫入國家戰(zhàn)略[1]。

土地利用變化是研究可持續(xù)發(fā)展和環(huán)境變化的方法之一，土地變化研究源于20世紀80年代，當時的學者認為土地利用是全球發(fā)生碳循環(huán)的關鍵因素，土地利用被學者看作碳源與碳匯[3]。20世紀90年代初，不少研究發(fā)現(xiàn)土地利用變化造成土地退化，人類環(huán)境耦合系統(tǒng)變得脆弱。經(jīng)過研究者的不斷努力，一門具有實際應用意義的土地變化科學誕生[4]。土地利用側(cè)重于兩方面，一是地球表面生長的植被，二是在自然生長過程中人類活動參與后所導致的土地利用/土地覆被變化[5]。20世紀90年代起，中國學者緊跟國際研究動態(tài)，對區(qū)域土地利用引起的變化也進行廣而深的研究[6]。王秀蘭等[7]通過收集國外研究土地利用動態(tài)變化的模型方法并加以匯總，將其分為3種類型，即系統(tǒng)診斷模型、土地利用動態(tài)變化模型、土地利用綜合評估模型。我國土地利用研究稍顯滯后，國外研究方法相對成熟，我國主要借鑒國外研究方法和模型進行土地利用研究，未來需要充分結(jié)合國情進行創(chuàng)新，服務土地利用規(guī)劃管理。

生態(tài)足跡是指能夠持續(xù)提供資源或消納廢物、具有生物生產(chǎn)力的地域空間[8]，生態(tài)足跡計算相對簡單，結(jié)果直觀全面，數(shù)據(jù)獲取途徑方便。相比國外，國內(nèi)生態(tài)足跡研究起步較晚，其概念于20世紀末被引入我國[9]。進入21世紀后，我國開始積極開展生態(tài)足跡研究。生態(tài)足跡模型的主要作用是對生態(tài)足跡與生態(tài)承載力進行比較，衡量一個地區(qū)的可持續(xù)發(fā)展狀況[10]。該模型運算所需要的數(shù)據(jù)一般有兩種獲取途徑。一是利用研究區(qū)發(fā)布的統(tǒng)計資料中各土地類型對應的消費產(chǎn)品數(shù)量進行折算[11]；二是借助調(diào)查問卷進行直接獲取[12]。陳薇等[13]改進了傳統(tǒng)的生態(tài)足跡模型，利用各種消耗品數(shù)據(jù)計算了研究區(qū)的生態(tài)足跡和可持續(xù)性。張海鷗[14]利用改進的生態(tài)足跡模型對2007—2009年張家界市旅游業(yè)的可持續(xù)發(fā)展狀況做出評價。崔維軍等[15]分析重化工業(yè)生態(tài)足跡的形成機制，將其生態(tài)足跡的形成分成兩個階段：一是資源消耗，二是環(huán)境排放。我國生態(tài)足跡研究較國外稍晚，但仍然取得一定成果，生態(tài)足跡的應用范圍從大尺度的國家、省域?qū)哟无D(zhuǎn)變?yōu)樾〕叨鹊氖杏蚝彤a(chǎn)業(yè)發(fā)展層次，研究尺度有所改變，研究的創(chuàng)新程度還落后于國外，我國需要不斷對生態(tài)足跡模型進行研究創(chuàng)新與應用實證。

改革開放后，我國社會經(jīng)濟快速發(fā)展，土地利用逐年發(fā)生變化，除受人類活動影響，它還受到政治、經(jīng)濟、科技和自然等客觀因素的影響。土地利用方式發(fā)生變化，致使當?shù)厣鷳B(tài)環(huán)境發(fā)生一定的變化，進而影響研究區(qū)的生態(tài)足跡和生態(tài)承載力[16]。隨著經(jīng)濟發(fā)展進入新常態(tài)，我國經(jīng)濟已由高速增長階段轉(zhuǎn)向高質(zhì)量發(fā)展階段。2015年，我國提出“新舊動能轉(zhuǎn)換”的概念。2017年，國務院辦公廳印發(fā)了《關于創(chuàng)新管理優(yōu)化服務 培育壯大經(jīng)濟發(fā)展新動能 加快新舊動能接續(xù)轉(zhuǎn)換的意見》，這是我國培育新動能、加速舊動能接續(xù)轉(zhuǎn)換的第一份文件，拉開新舊動能轉(zhuǎn)換的序幕。進入21世紀以來，濟南市經(jīng)濟發(fā)展迅速，對土地利用和環(huán)境產(chǎn)生重要影響，并帶來相應的變化，其涉及當?shù)卣摺a(chǎn)業(yè)類型結(jié)構(gòu)以及人類活動等方面[17]。濟南市作為新舊動能轉(zhuǎn)換綜合試驗區(qū)之一，擁有全國第一個新舊動能轉(zhuǎn)換先行區(qū)，新舊動能轉(zhuǎn)換背景下，濟南市土地利用研究和可持續(xù)發(fā)展研究具有極其重要的價值，對于經(jīng)濟發(fā)展階段轉(zhuǎn)變形勢下把握自然環(huán)境變化的規(guī)律具有重要的探索意義。因此，濟南市土地利用變化及其生態(tài)可持續(xù)發(fā)展狀況亟待研究。本文從土地利用時空格局演變和生態(tài)足跡兩方面進行分析，對濟南市2010—2020年的可持續(xù)發(fā)展水平做出評價。

1 研究區(qū)概況與數(shù)據(jù)處理

1.1 研究區(qū)概況

濟南市是山東省省會，位于華北平原東南部、山東省中西部，北接京津冀，南連長江三角洲，區(qū)域面積近1.03萬km2[18]。根據(jù)第七次全國人口普查數(shù)據(jù)，截至2020年11月1日零點，濟南市常住人口有近920.3萬人[19]，它是我國14個特大城市之一。2020年，濟南市地區(qū)生產(chǎn)總值達到10 140.9億元，成功晉升為“萬億俱樂部”的成員之一，其中第二產(chǎn)業(yè)的經(jīng)濟增長率最為顯著。2018年，濟南市被列為山東新舊動能轉(zhuǎn)換綜合試驗區(qū)之一，其在2021年經(jīng)國務院批準設立，濟南新舊動能轉(zhuǎn)換起步區(qū)成為繼雄安新區(qū)之后的全國第二個起步區(qū)。

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

1.2.1 數(shù)據(jù)來源與分類體系

土地利用/土地覆被數(shù)據(jù)來源[20]于中國科學院地理科學與資源研究所的資源環(huán)境科學與數(shù)據(jù)中心（http：//www.resdc.cn）。2010年、2020年土地利用/土地覆被數(shù)據(jù)作為基礎數(shù)據(jù)源，將濟南市行政矢量數(shù)據(jù)轉(zhuǎn)化為柵格數(shù)據(jù)作為掩模文件，利用ArcGIS 10.8軟件提取濟南市2010年和2020年土地利用數(shù)據(jù)，將數(shù)據(jù)轉(zhuǎn)換為ASCII格式后利用TerrSet 2020軟件將2010年和2020年濟南市土地利用數(shù)據(jù)分別進行重分類（見表1），最終通過CrossTab模塊得出土地利用轉(zhuǎn)移矩陣。

1.2.2 增勢圖譜和減勢圖譜

利用ArcGIS 10.8軟件，轉(zhuǎn)換濟南市2010年、2020年兩期土地利用的矢量數(shù)據(jù)圖，使矢量數(shù)據(jù)轉(zhuǎn)換為柵格數(shù)據(jù)，然后將2010年和2020年濟南市土地利用的柵格數(shù)據(jù)圖做疊加分析，并將其命名為zstp（增勢圖譜），在其屬性表中添加code字段，依據(jù)表2進行歸類編碼，得到濟南市2010—2020年土地利用類型的增勢圖譜。減勢圖譜的制作方法同增勢圖譜操作類似，依據(jù)表3進行歸類編碼后，得到2010—2020年濟南市土地利用變化的減勢圖譜。

1.2.3 生態(tài)足跡計算

生態(tài)足跡本質(zhì)就是某一區(qū)域人口消耗的資源和產(chǎn)出的廢棄物所需要的生物生產(chǎn)土地的總面積和水資源量[10]，主要將其分為耕地、林地、草地、水域、建設用地、化石能源用地6大類[21-22]。鑒于一個地區(qū)各地類的差異性，將其進一步劃分成生物資源賬戶、建設用地賬戶、化石能源賬戶三類[23]。在一定條件下，可以用生態(tài)足跡來衡量研究區(qū)所生活的人對于資源和環(huán)境的需要，若研究區(qū)的生態(tài)足跡數(shù)值增大，則說明該研究區(qū)對資源和環(huán)境的需求增大，反之，則減小。本文均利用各生態(tài)生產(chǎn)性土地的人均面積進行計算[24]。

通過《濟南統(tǒng)計年鑒2011》[25]與《濟南統(tǒng)計年鑒2021》[19]，分別搜集濟南市2010年、2020年的生態(tài)足跡相關數(shù)據(jù)。鑒于無法直接從統(tǒng)計年鑒找出相關的土地利用面積數(shù)據(jù)，可將當?shù)氐囊恍┵Y源生產(chǎn)量和消費量折算成具有相應生產(chǎn)能力的生態(tài)生產(chǎn)性土地[23]，如表4所示。

生物資源消費量主要受研究區(qū)的產(chǎn)量影響，研究區(qū)該生物資源的進口量影響較小，因此生物資源賬戶采用耕地、林地、草地、水域四類土地中能夠為研究區(qū)人民提供生物資源的產(chǎn)量數(shù)據(jù)來進行計算[26]，利用該資源的全球平均產(chǎn)量[27]將不同種類的生物資源生產(chǎn)量轉(zhuǎn)換為統(tǒng)一的單位公頃（hm2）后進行運算，各種生物資源的全球平均產(chǎn)量對照如表5所示。化石能源消費量則與研究區(qū)環(huán)境能夠處理的由人類活動所排放的CO2以及其他溫室氣體有關，因此研究區(qū)居民所消耗的化石能源數(shù)據(jù)參與該計算。建設用地賬戶的數(shù)據(jù)來源主要是依據(jù)電力消費數(shù)據(jù)進行計算[28]。根據(jù)折算系數(shù)與全球平均能源足跡[29]，將化石能源賬戶與建設用地賬戶的相關數(shù)據(jù)（見表6）轉(zhuǎn)換為研究區(qū)電力消耗所需要的土地利用面積。為便于計算，本文最終計算結(jié)果為人均層次，人口數(shù)據(jù)則采用第六次全國人口普查、第七次全國人口普查中濟南市全域的常住人口數(shù)量，研究中計算所需要的均衡因子和產(chǎn)量因子[30]如表7所示。

2 結(jié)果與分析

首先開展土地利用時空格局演變分析，涉及土地利用轉(zhuǎn)移矩陣和土地利用變化圖譜，然后開展生態(tài)足跡分析。

2.1 土地利用時空格局演變分析

2.1.1 土地利用轉(zhuǎn)移矩陣分析

2010—2020年，濟南市土地利用轉(zhuǎn)移矩陣如表8所示。在研究期內(nèi)，轉(zhuǎn)出最大的為耕地，貢獻率達76.63%，主要向建設用地轉(zhuǎn)化；其次為建設用地，貢獻率達10.71%，主要向耕地轉(zhuǎn)化。在濟南市土地利用的轉(zhuǎn)入面積中，最大為建設用地，有近146.28 km2，貢獻率高達72.97%；其次為耕地，貢獻率為13.76%。同時，林地的轉(zhuǎn)出面積中，過半數(shù)向草地轉(zhuǎn)化，草地轉(zhuǎn)出面積中近50%都向耕地轉(zhuǎn)化，水域主要向耕地轉(zhuǎn)化。除此之外，在林地轉(zhuǎn)入面積中，大部分來源于草地，水域和草地的轉(zhuǎn)入面積則大部分源自耕地。

2.1.2 土地利用變化圖譜分析

2010—2020年，濟南市土地利用類型變化圖譜如圖1所示。由土地利用類型增勢圖譜可見，大部分區(qū)域土地利用類型基本保持不變，不過仍有處于新增狀態(tài)的土地利用類型，其中新增建設用地最為明顯，主要集中在濟南市的中部地區(qū)，北部的商河縣、濟陽區(qū)，東南部的萊蕪區(qū)也有少量新增建設用地存在；其次為耕地。由土地利用類型減勢圖譜可見，大部分區(qū)域的土地利用類型保持不變，減勢最為明顯的就是耕地，濟南市各個區(qū)縣基本上都有耕地利用類型的減少；其次為建設用地的減少，主要集中在濟南市西部的槐蔭區(qū)，同時還有少量草地利用類型的減少。經(jīng)圖譜對比分析，新增的建設用地所在的區(qū)域位置和減少的耕地有高度重合，因此可以推斷濟南市2010—2020年新增建設用地主要由耕地轉(zhuǎn)化；新增水域也與耕地減少有著密不可分的關系；在減勢圖譜中，仍有部分未知的土地利用類型轉(zhuǎn)化方向尚不明朗。

2.2 生態(tài)足跡分析

按照數(shù)據(jù)處理公式，結(jié)合濟南市2011年和2021年統(tǒng)計年鑒中所提取的2010年和2020年相關數(shù)據(jù)進行計算，結(jié)果如表9、表10所示，然后計算出這兩個時間刻度的生態(tài)盈余/赤字。在計算生物資源的生態(tài)承載力時，要扣除12%的生物多樣性保護面積[33]。除此之外，還要注意化石能源用地在現(xiàn)實生活中沒有相對應的生態(tài)承載力[31]。由表9、表10可見，濟南市2010年處于生態(tài)赤字，為0.437 hm2；2020年同樣處于生態(tài)赤字，為0.623 hm2，可持續(xù)發(fā)展水平下降。

除去沒有對應生態(tài)承載力的化石能源用地，2010年濟南市耕地和建設用地生態(tài)承載力高于生態(tài)足跡，即耕地和建設用地處于生態(tài)盈余；相比之下，水域、草地和林地都處于生態(tài)赤字，如圖2所示。

由圖3可見，2020年草地、林地和水域表現(xiàn)為生態(tài)赤字，草地赤字貢獻率顯著，林地、水域分別占據(jù)的人均生態(tài)足跡也有一定程度的減少，除此之外，雖然耕地的人均生態(tài)承載力相對于人均生態(tài)足跡有所下降，但耕地和建設用地仍保持生態(tài)盈余。

由圖4可見，2020年生態(tài)足跡有所改善，同時生態(tài)承載力水平發(fā)生顯著下降，導致2020年生態(tài)赤字比2010年生態(tài)赤字高，表明2020年濟南市可持續(xù)發(fā)展水平較2010年有所下降。

由圖5可見，2010年濟南市耕地、草地、水域的人均生態(tài)足跡明顯高于2020年，化石能源用地人均生態(tài)足跡則低于2020年。

3 討論

研究發(fā)現(xiàn)，濟南市土地利用變化十分迅速，顯著表現(xiàn)為大量耕地轉(zhuǎn)為建設用地，這與董梁蘊等[32]和劉鳳蓮等[34]得出的結(jié)論相似，即減少的耕地大部分轉(zhuǎn)化為建設用地，耕地面積優(yōu)勢減弱。2010—2020年，濟南市減少的耕地大部分位于黃河沖積平原、山前傾斜平原等平坦地區(qū)，地勢條件良好，水資源較為豐富，適合建設城鎮(zhèn)，這表明濟南市城鎮(zhèn)化可能以犧牲耕地為代價[35]。新增耕地與減少的建設用地也存在一定的空間重合，原因可能是當?shù)卣吲c其他因素共同作用，需要做進一步討論。

我國生態(tài)足跡研究稍晚，創(chuàng)新程度不高，本文使用的生態(tài)足跡模型屬于傳統(tǒng)模式。從數(shù)據(jù)來源來看，對濟南市公布的統(tǒng)計數(shù)據(jù)中各土地類型相對應的產(chǎn)品數(shù)量進行折算，轉(zhuǎn)變?yōu)樯鷳B(tài)生產(chǎn)性土地進行計算。分析發(fā)現(xiàn)，2010—2020年，濟南市化石能源用地生態(tài)足跡增長0.124 hm2，原因可能是濟南市城鎮(zhèn)化的不斷推進和常住人口數(shù)量的增長使得研究區(qū)對能源、建筑材料等高耗能產(chǎn)業(yè)的需求增加。此外，耕地生態(tài)承載力占很大比重，但2020年人均耕地生態(tài)承載力也低于2010年，減少0.228 hm2，這可能與耕地面積減少、農(nóng)業(yè)科技化水平有關。新舊動能轉(zhuǎn)換背景下，相比2010年生態(tài)足跡，濟南市2020年生態(tài)足跡有所改善，但生態(tài)承載力水平顯著降低，生態(tài)赤字擴大。當前，濟南市要發(fā)展新動能、轉(zhuǎn)換舊動能，保護耕地、林地及水域，從土地利用角度改善生態(tài)平衡，提高生態(tài)承載力，減少生態(tài)赤字，助力可持續(xù)發(fā)展。

自改革開放至今，我國城鎮(zhèn)化從未中斷，一座城市的崛起不僅代表經(jīng)濟騰飛，也是耕地面積縮減的表征。進入21世紀，城鎮(zhèn)化的弊端開始浮現(xiàn)，首當其沖的就是糧食安全。目前，我國耕地保護力度依然不夠，農(nóng)業(yè)技術(shù)水平有待提高。經(jīng)土地利用變化生態(tài)足跡分析，濟南市耕地面積減少，大部分耕地向建設用地轉(zhuǎn)化，耕地生態(tài)承載力降低。作為山東省省會，在新舊動能轉(zhuǎn)換中，濟南市要堅持新發(fā)展理念，加強生態(tài)保護，減少生態(tài)赤字，建設綠色宜居共享的國際化智慧新城。未來，可對不同城市開展類似研究，提高耕地生態(tài)承載力，以切實保護耕地，改善農(nóng)業(yè)生態(tài)環(huán)境。

4 結(jié)論

本文以濟南市為研究區(qū)，結(jié)合土地利用時空格局演變，分析其生態(tài)足跡。2010—2020年，濟南市土地利用變化從未終止，顯著表現(xiàn)為轉(zhuǎn)出的耕地主要向建設用地轉(zhuǎn)化，新增建設用地與減少耕地在地域空間上高度重合，這表明研究期內(nèi)濟南市城鎮(zhèn)化可能以耕地轉(zhuǎn)出為代價。濟南市2020年耕地生態(tài)足跡低于2010年，生態(tài)承載力也較2010年降低0.228 hm2。

原因可能是濟南市城鎮(zhèn)化速度過快，常住人口數(shù)量不斷增加，資源需求量增大。濟南市可持續(xù)發(fā)展狀況有待改善。化石能源用地的生態(tài)足跡趨于增加，它是2010年和2020年濟南市生態(tài)足跡的重要組成部分。

參考文獻

1 戴云菲.可持續(xù)發(fā)展理論文獻綜述[J].商，2016（13）：111.

2 陳月英.可持續(xù)發(fā)展理論綜述[J].長春師范學院學報，2000（5）：50-52.

3 Ii B L T，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 of the United States of America，2007（52）：20666-20671.

4 Rindfuss R R，Walsh S J，Turner B L，et al.Developing a science of land change：Challenges and methodological issues[J].Proceedings of the National Academy of Sciences，2004（39）：13976-13981.

5 袁子坤.土地利用/土地覆被變化研究綜述[J].甘肅農(nóng)業(yè)科技，2016（9）：73-77.

6 路云閣，蔡運龍，許月卿.走向土地變化科學：土地利用/覆被變化研究的新進展[J].中國土地科學，2006（1）：55-61.

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

8 錢 易.環(huán)境保護與可持續(xù)發(fā)展[M].北京：高等教育出版社，2010：182.

9 龍愛華，張志強，蘇志勇.生態(tài)足跡評介及國際研究前沿[J].地球科學進展，2004（6）：971-981.

10 范 瑋，張 麗，李惠芳，等.生態(tài)足跡研究現(xiàn)狀及展望[J].四川林業(yè)科技，2019（1）：87-91.

11 蔣依依，王仰麟，卜心國，等.國內(nèi)外生態(tài)足跡模型應用的回顧與展望[J].地理科學進展，2005（2）：13-23.

12 周國忠.國內(nèi)外生態(tài)足跡研究進展[J].浙江學刊，2010（6）：167-171.

13 陳 薇，孫 靜.基于改進的生態(tài)足跡模型對循環(huán)經(jīng)濟發(fā)展水平的測度研究[J].安徽農(nóng)業(yè)科學，2017（24）：236-238.

14 張海鷗.基于生態(tài)足跡的旅游城鎮(zhèn)可持續(xù)發(fā)展評價研究[J].安徽農(nóng)業(yè)科學，2010（5）：74-75.

15 崔維軍，周飛雪，徐常萍.中國重化工業(yè)生態(tài)足跡估算方法研究[J].中國人口·資源與環(huán)境，2010（8）：137-141.

16 Gascoigne W R，Hoag D，Koontz L，et al.Valuing ecosystem and economic services across land-use scenarios in the Prairie Pothole Region of the Dakotas，USA[J].Ecological Economics，2011（10）：1715-1725.

17 Ghosh S.A city growth and land-use/land-cover change：a case study of Bhopal，India[J].Modeling Earth Systems and Environment，2019（4）：1569-1578.

18 高 登.大都市區(qū)土地利用變化及其生態(tài)系統(tǒng)服務價值響應：以濟南市為例[J].中國資源綜合利用，2021（12）：72-74.

19 濟南市統(tǒng)計局，國家統(tǒng)計局濟南調(diào)查隊.濟南統(tǒng)計年鑒2021[M].北京：中國統(tǒng)計出版社，2021：49-50.

20 徐新良，劉紀遠，張樹文，等.中國多時期土地利用遙感監(jiān)測數(shù)據(jù)集（CNLUCC）[EB/OL].（2018-07-02）[2023-01-12].https：//www.resdc.cn/DOI/DOI.aspx？DOIid=54.

21 張志強，徐中民，程國棟.生態(tài)足跡的概念及計算模型[J].生態(tài)經(jīng)濟，2000（10）：8-10.

22 辛 宇.基于生態(tài)足跡分析法對陜西省可持續(xù)發(fā)展定量分析[D].西安：西北大學，2007：18-19.

23 關青青.長三角地區(qū)三維生態(tài)足跡動態(tài)分析及影響因素研究[D].上海：上海師范大學，2021：15-16.

24 Wackernagel M，Rees W.Our ecological footprint：reducing human impact on the earth[M].Vancouver：New Society Publishers，1996：32-33.

25 濟南市統(tǒng)計局，國家統(tǒng)計局濟南調(diào)查隊.濟南統(tǒng)計年鑒2011[M].北京：中國統(tǒng)計出版社，2011：40-41.

26 張小雪.基于3維生態(tài)足跡的湖南省可持續(xù)發(fā)展研究[D].長沙：中南大學，2013：11-12.

27 謝高地，曹淑艷，曹春霞.中國生態(tài)資源承載力研究[M].北京：科學出版社，2011：1-374.

28 尋芳芳.長江經(jīng)濟帶資源環(huán)境承載力研究[D].北京：中國地質(zhì)大學，2019：14-15.

29 沈 杉.云南省文山州石漠化問題研究[D].昆明：云南財經(jīng)大學，2014：17-18.

30 劉某承，李文華，謝高地.基于凈初級生產(chǎn)力的中國生態(tài)足跡產(chǎn)量因子測算[J].生態(tài)學雜志，2010（3）：592-597.

31 王 鵬，田亞平.衡陽市可持續(xù)發(fā)展狀況的生態(tài)足跡模型定量評估[J].衡陽師范學院學報，2004（6）：89-93.

32 董梁蘊，董好剛，牛瑞卿，等.基于遙感的鄭州市近十年土地利用時空變化及驅(qū)動力分析[J].華南地質(zhì)，2021（4）：398-405.

33 張曉昱，王儷靜.資源型城市生態(tài)環(huán)境修復研究：以河南省鶴壁市為例[J].商丘師范學院學報，2018（8）：94-97.

34 劉鳳蓮，楊人懿，楊 磊.基于遙感影像的武漢市土地利用時空變化研究[J].忻州師范學院學報，2021（5）：49-57.

35 王 懇，李新舉.城鎮(zhèn)化驅(qū)動下濟南市土地利用變化特征[J].山東農(nóng)業(yè)大學學報（自然科學版），2018（2）：352-358.