旅游型海岸區(qū)域土地利用變化?生態(tài)效應及驅(qū)動力分析

余中元 溫文柱

摘要 以三亞市海棠區(qū)為例，運用ArcGIS和ENVI軟件對2002—2022年遙感影像進行處理，探索20年來土地利用變化趨勢；計算人均生態(tài)足跡和生態(tài)承載力，分析土地利用變化產(chǎn)生的生態(tài)環(huán)境效應。結(jié)果表明：海棠區(qū)的土地利用狀況發(fā)生了巨大的變化，其中變化較為顯著的是耕地、林地、建設(shè)用地；耕地和林地的流向主要是建設(shè)用地。這一時期海棠區(qū)高速發(fā)展，人類活動對土地利用的影響加劇。海棠區(qū)土地利用生態(tài)效應不宜樂觀，海棠區(qū)生態(tài)足跡與生態(tài)承載力供給呈反方向發(fā)展趨勢，海棠區(qū)人口對自然資源的利用已逐漸接近自然生態(tài)系統(tǒng)的生態(tài)承載力飽和范圍。

關(guān)鍵詞 旅游型海岸區(qū)；土地利用變化；生態(tài)效應；驅(qū)動力

中圖分類號 F301.2? 文獻標識碼 A

文章編號 0517-6611（2024）10-0089-06

doi：10.3969/j.issn.0517-6611.2024.10.020

Analysis of Land Use Change， Ecological Effect and Driving Force in Tourism Coastal Region—Taking Haitang District of Sanya City as an Example

YU Zhong-yuan1，2，WEN Wen-zhu1

（1.School of Geography and Tourism， Huizhou University， Huizhou，Guangdong 516007;2.College of Geography and Environmental Science， Hainan Normal University， Haikou，Hainan 571158）

Abstract Taking Haitang District of Sanya City as an example，ArcGIS and ENVI software were used to process remote sensing images from 2002 to 2022， explore the change trend of land use over the past 20 years，calculate per capita ecological footprint and ecological carrying capacity， and analyze the ecological environmental effects generated by land use changes.The results showed that the land use status in Haitang District had undergone great changes， among which the most significant changes were cultivated land， forest land and construction land.The decrease of cultivated land and forest land was mainly for construction land.The rapid development of Haitang District in this period reflected the intensification of human activities on land use.The ecological effect of land use in Haitang District should not be optimistic.The ecological footprint and ecological carrying capacity supply in Haitang District developed in the opposite direction.The utilization of natural resources by population in Haitang District gradually approached the saturation range of ecological carrying capacity of natural ecosystem.

Key words Tourism coastal region;Land use change;Ecological effects;Driving force

基金項目 國家自然科學基金項目（41961041）;廣東省駐幫鎮(zhèn)扶村農(nóng)村科技特派員重點派駐項目（KTP20210412）；海南省自然科學基金重點項目（ZDYF2016148）。

作者簡介 余中元（1969—）， 男， 四川廣安人，教授，博士，從事土地利用與生態(tài)經(jīng)濟、旅游規(guī)劃與資源利用研究。

收稿日期 2023-06-25

經(jīng)濟全球化、新型城市化以及對土地的過度利用，使得區(qū)域土地景觀的結(jié)構(gòu)、功能都受到一定程度的破壞，導致嚴重的生態(tài)環(huán)境問題，比如水土流失、空氣污染、生物群落失衡以及自然災害的頻繁發(fā)生［1-3］。自十八大以來，我國加強國土資源的有效管理，實施節(jié)約、合理、有序的開發(fā)政策，努力實現(xiàn)綠色、有序、可持續(xù)的國土空間開發(fā)，促進經(jīng)濟社會的健康增長［4］。

海岸帶是陸地和海洋的過渡地帶，是人類活動和自然環(huán)境之間的重要界面，是海洋生態(tài)系統(tǒng)重要的棲息地和繁殖地。隨著市場經(jīng)濟的繁榮和人們收入的提高，居民旅游消費需求和消費能力增強。旅游開發(fā)等人類活動引起的土地利用改變對海岸生態(tài)環(huán)境造成重要影響。旅游產(chǎn)業(yè)的快速發(fā)展和過快的城鎮(zhèn)化，給海岸生態(tài)環(huán)境帶來了潛在的威脅。海灘沙灘清理、海洋開發(fā)、濱海城市擴展等一定程度上破壞了原有的生態(tài)環(huán)境。大量游客的涌入和游艇、潛水、浮潛等活動，可能破壞海洋的生態(tài)平衡，給海洋動植物帶來很大的損害和威脅，引發(fā)潛在的生態(tài)災難。此外，旅游業(yè)發(fā)展還會帶來垃圾量的增加、水體的污染和噪音等問題。隨著城市化的加速和人們海濱旅游的增加，海岸帶面臨的環(huán)境挑戰(zhàn)不斷增多。

旅游行業(yè)急需建立旅游用地生態(tài)效應評估和監(jiān)管制度［5-6］。開展土地利用景觀格局演變及預測研究，研究海岸生態(tài)效應和生態(tài)環(huán)境變化的機理、特點和規(guī)律，對于保護海洋生態(tài)系統(tǒng)、減輕環(huán)境破壞、提升生態(tài)旅游價值、對區(qū)域土地合理開發(fā)保護、提高社會發(fā)展的可持續(xù)均具有重要意義。該研究以旅游型海岸區(qū)三亞市海棠區(qū)為例，研究區(qū)景區(qū)土地利用變化趨勢、通過計算生態(tài)足跡和生態(tài)承載分析土地利用變化產(chǎn)生的生態(tài)環(huán)境效應，為政府制定合理的管理政策和合理利用方案提供科學依據(jù)。

1 資料與方法

1.1 研究區(qū)概況

三亞市海棠區(qū)位于我國海南省南部，地理位置為109°52′30″～110°18′00″E、18°12′30″～18°37′30″N，東南臨瓊海，西接三亞市崖州區(qū)和天涯區(qū)，北鄰樂東黎族自治縣，南瀕南海（圖1）。海棠區(qū)總面積為166 km2，其中陸地面積為38.13 km2，水域面積為127.87 km2，海岸線長達54.6 km。海棠區(qū)地質(zhì)基底主要由古生代基巖和新生代沉積巖構(gòu)成的，屬于海島型火山巖山地和沿海平原地貌。

1.2 數(shù)據(jù)來源 該研究所需的遙感數(shù)據(jù)及土地利用數(shù)據(jù)遙感數(shù)據(jù)來源于中國地理監(jiān)測云平臺，由中國科學院提供，包含海棠區(qū)2002年9月21日的Landsat7衛(wèi)星影像、2012年12月27日的Landsat7衛(wèi)星影像、2022年4月26日的Landsat8衛(wèi)星影像。生態(tài)足跡和環(huán)境承載力數(shù)據(jù)源選自2002、2012、2022年《三亞市統(tǒng)計年鑒》，其中包括耕地（谷類、豆類、薯類、油料等）、林地（原始林、次生林和人工林等）、草地（豬肉、牛肉、羊肉、禽肉、奶類、禽蛋等用地）、水域（河流、湖泊、水庫）、建設(shè)用地（居民區(qū)、商業(yè)區(qū)、行政機構(gòu)、休閑設(shè)施等）。

1.3 研究方法 該研究定性和定量分析海棠區(qū)生態(tài)足跡與生態(tài)承載力供給的現(xiàn)狀與發(fā)展趨勢，為管理者制定合理的管理政策和土地利用方案提供科學依據(jù)。具體策略：構(gòu)建土地利用變化模型方法和合適的指標體系；進行土地利用類型變化分析和空間轉(zhuǎn)移分析，探索土地利用變化趨勢；計算人均生態(tài)足跡和生態(tài)承載力，分析土地利用變化產(chǎn)生的生態(tài)環(huán)境效應；提出土地利用及生態(tài)效應優(yōu)化的措施（圖2）。

1.3.1 土地利用變化分析。

運用ArcGIS和ENVI軟件對遙感影像進行處理，根據(jù)遙感解譯結(jié)果，對土地利用數(shù)據(jù)進行采集、處理、空間綜合分析。將不同時期的三亞市海棠區(qū)海岸帶土地利用分類柵格圖像按照國家標準的土地利用分類進行矢量提取、輸出，并運用ArcMap對所選地物信息進行重分類、修整等，獲得2002、2012、2022年的土地利用分類專題圖。統(tǒng)計2002、2012和 2022年三亞市海棠區(qū)土地利用分類專題圖面積信息，運用ArcGIS軟件進行土地利用情況分析。對土地利用現(xiàn)狀進行統(tǒng)計分析，了解該地的土地使用情況。

1.3.2 生態(tài)效應評價。

生態(tài)足跡是衡量生態(tài)效應的一種方法，最早由加拿大人威廉·里斯特·理查德森（William Rees）和馬思·麥克斯·溫厄姆（Mathis Wackernagel）于1992年提出并完善，又稱為生態(tài)占用、生態(tài)痕跡、生態(tài)腳印。已被廣泛應用于生態(tài)學、環(huán)境經(jīng)濟學、社會學等領(lǐng)域，成為評估和促進生態(tài)可持續(xù)發(fā)展的重要方法［7］。

生態(tài)足跡的基本思想是將人類生態(tài)足跡與生態(tài)承載力轉(zhuǎn)化為可以共同比較的土地面積，通過比較二者的差值來分析人類對生態(tài)環(huán)境的利用情況［8］。計算給定人口和經(jīng)濟條件下，通過將個人消耗的資源以及吸納人類所排放的廢物折合成全球統(tǒng)一的具有生產(chǎn)力的地域面積，來計算區(qū)域生態(tài)的供給與需求差值，即生態(tài)赤字或生態(tài)盈余，以測度不同區(qū)域內(nèi)人類活動對生態(tài)環(huán)境現(xiàn)狀的影響［9］。

生態(tài)足跡計算公式如下：

EF=N×ef=N×（rj×aj）（1）

式中：EF為總的生態(tài)足跡;N為總?cè)丝跀?shù);ef為人均生態(tài)足跡;rj為均衡因子;aj為人均擁有生物生產(chǎn)性土地面積。

生態(tài)承載力計算公式如下：

EC=N×ec=N×（aj×rj×yj）（2）

式中：EC為總的生態(tài)承載力;N為總?cè)丝跀?shù)；ec為人均生態(tài)承載力；aj為人均擁有生物生產(chǎn)性土地面積；rj為均衡因子；yj為產(chǎn)量因子。

生態(tài)赤字（ED）計算公式如下：

ED=EF-EC（3）

若ED>0，則表示生態(tài)赤字，若ED<0，則表示生態(tài)盈余。為方便計算，用人均數(shù)值來計算生態(tài)赤字或生態(tài)盈余的情況。

在計算生態(tài)足跡時，根據(jù)海棠區(qū)海岸帶生物生產(chǎn)性特點，土地面積計算選擇考慮5種土地利用的類型，分別是耕地、林地、草地、水域和建設(shè)用地。為了將這5種具有不同生物生產(chǎn)力的土地面積轉(zhuǎn)化為可加總的具有比較性的土地面積，需要對各生物生產(chǎn)性土地面積乘以均衡因子［10-11］。某類生物生產(chǎn)性土地面積的均衡因子等于全球該類生物生產(chǎn)性土地面積的平均生態(tài)生產(chǎn)力除以全球所有生物生產(chǎn)性土地面積的平均生態(tài)生產(chǎn)力?？紤]到各市（縣）商品和資源的流通與交易量難以統(tǒng)計以及生態(tài)系統(tǒng)的完整性，該研究選用海南省產(chǎn)量因子來計算三亞市海棠區(qū)的生態(tài)承載力。海南省產(chǎn)量因子可以根據(jù)全國產(chǎn)品平均產(chǎn)量與海南省平均產(chǎn)量之間的比值確定［12-13］。

2 結(jié)果與分析

2.1 三亞市海棠區(qū)的土地利用變化

2.1.1 土地利用類型變化。

從圖3和表1可以看出，2002—2012年用地類型主要是林地和耕地，草地、水域和建設(shè)用地面積占比較小。草地面積近20年有大幅減少，2022年較2002年減少1.93 km2，減少了31.54%。水域面積波動幅度較大，2002—2012年增加2.49 km2；2012—2022年減少了0.77 km2。建設(shè)用地面積不斷增長，2002—2012年增加1.43 km2；2012—2022年又增加14.27 km2，2022年較2002年增加了266.10%。耕地面積基本呈減少趨勢，2002—2012年減少1.45 km2，2012—2022年又減少了4.65 km2。

2.1.2 土地利用類型轉(zhuǎn)換。

利用ArcGIS軟件對海棠區(qū)土地利用轉(zhuǎn)換進行空間疊加分析，得出2002—2012和2012—2022年的轉(zhuǎn)移矩陣圖（圖4）和轉(zhuǎn)移矩陣表（表2～3）。結(jié)果顯示：耕地2002—2012年主要轉(zhuǎn)化為林地、水域和建設(shè)用地，2012—2022年主要轉(zhuǎn)化為建設(shè)用地、水域和林地；林地2002—2012年主要轉(zhuǎn)化為建設(shè)用地、耕地和水域，2012—2022年主要轉(zhuǎn)化為建設(shè)用地、耕地和水域。草地2002—2012年主要轉(zhuǎn)化為林地和耕地；2012—2022 年主要轉(zhuǎn)化為水域，其次是耕地和建設(shè)用地。水域2002—2012年主要轉(zhuǎn)化為林地和耕地；2012—2022年主要轉(zhuǎn)化為林地，其次轉(zhuǎn)化成耕地和建設(shè)用地。建設(shè)用地2002—2012 年主要轉(zhuǎn)化為耕地和水域，2012—2022 年主要轉(zhuǎn)化為林地。

2.1.3 土地利用變化的驅(qū)動力。

由于自然因素對于該研究的時間尺度而言基本可以忽略不計，因此海棠區(qū)的土地利用類型變化驅(qū)動力分析中，自然驅(qū)動力不占主要原因，在此不做過多討論。在Excel工作表中運用CORREL函數(shù)計算研究區(qū)2009—2022年不同土地利用類型與各因子之間的相關(guān)系數(shù)（表4），得出耕地、林地、草地、建設(shè)用地面積與人口、經(jīng)濟的變化關(guān)系（圖5）。由表4可知，耕地、林地、草地、建設(shè)用地面積與人口、人均可支配收入的相關(guān)系數(shù)的絕對值都大于0.85，其中耕地、林地、草地面積與人口、人均可支配收入呈強負相關(guān)性，而建設(shè)用地面積與人口、人均可支配收入呈強正相關(guān)性，這說明人口、經(jīng)濟對耕地、林地、草地、建設(shè)用地面積的變化有重要影響。而水域面積與各影響因子的相關(guān)系數(shù)的絕對值都比較小，說明人口、經(jīng)濟在水域面積變化的影響因素中不占主要地位。

（1）人口變化是導致土地利用變化的主要原因。

海棠區(qū)總?cè)丝趶?002年的44 944人增至2012年的47 637人，到了2022年人口更是達到了85 086人。隨著海棠區(qū)總?cè)丝诘募眲≡龆?，海棠區(qū)建設(shè)用地面積也必然會迅速增加，從2002年的5.90 km2增至2012年的7.33 km2，2022年又大幅增至21.60 km2。對比人口變化趨勢和建設(shè)用地面積變化趨勢，可以發(fā)現(xiàn)2002—2022年人口數(shù)增加且建設(shè)用地面積也相應增加，建設(shè)用地面積增加必然導致耕地面積減少。2002—2022年耕地面積的變化趨勢與建設(shè)用地面積的變化趨勢有很好的對應關(guān)系。

（2）經(jīng)濟增長是土地利用變化的最重要驅(qū)動力。

近年來，海棠區(qū)的旅游業(yè)極大地推動了全區(qū)經(jīng)濟社會和城市建設(shè)快速發(fā)展。經(jīng)濟增長和社會發(fā)展使農(nóng)業(yè)生產(chǎn)部門從傳統(tǒng)農(nóng)業(yè)向現(xiàn)代農(nóng)業(yè)轉(zhuǎn)變，改變了農(nóng)民的耕種方式，也影響了土地利用類型的變化。城市化與經(jīng)濟增長和社會發(fā)展之間存在正相關(guān)關(guān)系。經(jīng)濟增長越快，人口增加越多，城市土地利用變化越大；反之，人口增長也會促進城市土地利用變化。

2.2 生態(tài)足跡與承載力的動態(tài)變化

2.2.1 生態(tài)足跡的變化。

土地利用類型的模式及其轉(zhuǎn)化趨勢能夠在很大程度上反映城市發(fā)展的結(jié)構(gòu)和規(guī)模，生態(tài)足跡可以反映出海棠區(qū)土地利用變化的生態(tài)效應［14-15］。世界環(huán)境發(fā)展委員會的研究表明，在評估生態(tài)承載能力的過程中，應考慮扣除12%的生物多樣性的保護面積［16］。

由表5可知，海棠區(qū)耕地的人均生態(tài)足跡2002—2012年減少了0.005 km2，到2022年略微增加至0.023 km2；林地的人均生態(tài)足跡在2002—2012年維持在0.001 km2，2022年增至0.002 km2；草地的人均生態(tài)足跡在2002—2012年維持在0.004? km2，2022年增加至0.008? km2；水域的人均生態(tài)足跡2002—2012年維持在0.008? km2，2022年增加至0.023? km2；建設(shè)用地的人均生態(tài)足跡從2002年的0.004 km2增至2012年的0.007? km2，2022年增加至0.086? km2。

2.2.2 承載力的變化。

由表6可知，海棠區(qū)耕地的人均生態(tài)承載力2002—2012年增加了0.002 km2，到2022年驟減至0.090 km2，較2002年減少了48.6%；林地的人均生態(tài)承載力從2002年的0.714 km2減少至2012年的0.474 km2，2022年又減少至0.254 km2；草地的人均生態(tài)承載力從2002年的0.043 km2減至2012年的0.030? km2，2022年又減少至0.011? km2；水域的人均生態(tài)承載力從2002年0.132? km2減至2012年的0.098? km2，2022年又減少至0.053? km2；建設(shè)用地的人均生態(tài)承載力從2002年的0.131 km2增至2012年的0.243? km2，2022年又減少至0.040? km2。

2.2.3 綜合生態(tài)效益的變化。

從圖6可以看出，海棠區(qū)總的人均生態(tài)足跡從2002年的0.044 km2略微減少至2012年的0.042 km2，而到了2022年驟增至0.119 km2，但是人均生態(tài)承載力卻從2002年的1.195 km2降至2012年的1.022 km2，至2022年更是降至0.448 km2。海棠區(qū)的可利用生態(tài)承載力在穩(wěn)定持續(xù)降低，土地利用的人均生態(tài)盈余從2002年的1.008 km2 降至2012年的0.857 km2，2022年進一步降至0.265 km2。海棠區(qū)的生態(tài)足跡與生態(tài)承載力供給呈相反的方向發(fā)展，對當?shù)厣鷳B(tài)環(huán)境的利用雖然還未超出自然生態(tài)系統(tǒng)的生態(tài)承載力范圍，但生態(tài)系統(tǒng)已經(jīng)有所退化，人地關(guān)系趨于緊張。

3 結(jié)論與建議

3.1 結(jié)論

（1）人口和經(jīng)濟發(fā)展是土地利用變化的主要驅(qū)動因素。

土地利用分析的結(jié)果表明海棠區(qū)2002—2022年的土地利用方式有較大的變化，其中變化幅度最大的是耕地、林地和建設(shè)用地。由于人口的增長和城鎮(zhèn)化的發(fā)展，建設(shè)用地面積大幅增加，而建設(shè)用地的增加導致耕地和林地的減少。這說明在這20年間海棠區(qū)的經(jīng)濟正在高速發(fā)展，人口和經(jīng)濟發(fā)展是主要的土地利用變化驅(qū)動因素。

（2）海棠區(qū)的生態(tài)環(huán)境即將進入生態(tài)赤字，土地利用的負面生態(tài)效應顯現(xiàn)。研究區(qū)2002—2022年人均生態(tài)足跡、人均生態(tài)承載力及人均生態(tài)盈余變化情況表明，由于人口的增加和土地利用類型的轉(zhuǎn)變，海棠區(qū)人均生態(tài)足跡在不斷增加，而自然環(huán)境承載力卻在不斷減少，二者呈負相關(guān)的發(fā)展方向，且對自然資源的利用已逐漸接近自然生態(tài)系統(tǒng)的生態(tài)承載力飽和范圍。這極不利于維持海棠區(qū)自然環(huán)境的可持續(xù)發(fā)展，若按照這個趨勢繼續(xù)發(fā)展，海棠區(qū)的生態(tài)環(huán)境即將進入生態(tài)赤字。

3.2 建議

（1）加強土地利用規(guī)劃，明確土地功能，調(diào)整產(chǎn)業(yè)結(jié)構(gòu)。加強土地利用規(guī)劃時，劃分不同的用途和區(qū)域，明確不同區(qū)域的功能，避免重復占用土地，合理布局。調(diào)整海棠區(qū)的產(chǎn)業(yè)結(jié)構(gòu)，適量增加農(nóng)業(yè)用地，維持現(xiàn)存的林地面積，高效利用現(xiàn)今的建設(shè)用地。

（2）合理農(nóng)業(yè)種植，推廣多種經(jīng)營，提高土地利用效益。在農(nóng)業(yè)生產(chǎn)中，采用科學的耕作技術(shù)、高效肥料和農(nóng)藥管理等手段，保證土地的產(chǎn)量最大化，提高土地的利用效率。農(nóng)業(yè)、畜牧業(yè)、林業(yè)、漁業(yè)等多種經(jīng)營方式有機結(jié)合，實現(xiàn)資源共享，提高土地利用效益。

（3）多層次開發(fā)，提高建筑密度。城市用地可以進行多層次的開發(fā)，例如在原有建筑的基礎(chǔ)上進行拓建或者細化利用，或者通過增加建筑密度可以實現(xiàn)用地的合理利用，提高建筑的容積率，以此提高土地的使用效率。

（4）加強信息化建設(shè)，實現(xiàn)共享監(jiān)督機制。可以通過遙感技術(shù)、地理信息系統(tǒng)等現(xiàn)代技術(shù)手段，加強土地利用信息化建設(shè)和推廣，對大范圍的區(qū)域進行監(jiān)測和分析，因地制宜利用土地，保護水源、土壤等環(huán)境，保障土地資源的可持續(xù)利用。

（5）加強宣傳，發(fā)動公眾參與。加大土地利用監(jiān)測宣傳力度，讓大眾了解其背景、原理和應用的重要性，以推動社會各方面人員更多參與到海岸生態(tài)保護工作中來。

參考文獻

［1］ 彭建，王仰麟，張源，等.土地利用分類對景觀格局指數(shù)的影響［J］.地理學報，2006，61（2）：157-168.

［2］ 王秀茹，韓興，朱國平，等.關(guān)于土地開發(fā)整理與生態(tài)環(huán)境問題的分析［J］.水土保持研究，2004，11（3）：151-153.

［3］馮麗媛，米文寶，馬國慶.基于CA-Markov模型的寧夏沿黃生態(tài)經(jīng)濟帶土地利用變化及模擬研究［J］.水土保持通報，2019，39（5）：218-222，230.

［4］ 劉家明.國內(nèi)外海島旅游開發(fā)研究［J］.華中師范大學學報（自然科學版），2000，34（3）：349-352.

［5］ 余中元，趙志忠.旅游用地生態(tài)效應及調(diào)控研究框架設(shè)計［J］.生態(tài)科學，2013，32（4）：523-528.

［6］ 張生舉，楊海江.生態(tài)系統(tǒng)服務視角下的旅游服務研究進展［J］.旅游縱覽，2023（7）：82-84.

［7］ KOK A L，BARENDREGT W.Understanding the adoption， use， and effects of ecological footprint calculators among Dutch citizens［J］.Journal of cleaner production，2021，326：1-14.

［8］ 蘇美娜.基于生態(tài)系統(tǒng)的廈門灣海洋生物多樣性優(yōu)先保護格局研究［D］.廈門：廈門大學，2019.

［9］ 吳孟孟，賈培宏，潘少明，等.連云港海岸帶土地利用變化生態(tài)效應量化研究［J］.海洋通報，2015，34（5）：530-539.

［10］ 郭湘，聶篤聰，易紅超.湖南省2005年土地消費行為研究［J］.國土資源科技管理，2008（1）：117-120.

［11］ 劉某承，李文華.基于凈初級生產(chǎn)力的中國各地生態(tài)足跡均衡因子測算［J］.生態(tài)與農(nóng)村環(huán)境學報，2010，26（5）：401-406.

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

［13］ 郭秀銳.城市生態(tài)系統(tǒng)健康評價：以廣州市為例［D］.北京：北京師范大學， 2003.

［14］ 熊鷹，艾贛雄，周晨，等.基于改進三維生態(tài)足跡模型的洞庭湖區(qū)生態(tài)可持續(xù)時空演化研究［J］.生態(tài)學報，2022，42（3）：1165-1179.

［15］? 徐中民，程國棟，張志強.生態(tài)足跡方法的理論解析［J］.中國人口·資源與環(huán)境，2006，16（6）：69-78.

［16］ 梁雪石，鄭福云，郭文棟，等.黑龍江省生態(tài)足跡分析及ARIMA模型預測［J］.國土與自然資源研究，2020（2）：14-19.