基于生態(tài)足跡模型的白洋淀流域生態(tài)承載力評價與時空格局

摘要：以白洋淀流域為研究區(qū)，基于生態(tài)足跡模型，結合統(tǒng)計年鑒數(shù)據(jù)、遙感影像數(shù)據(jù)，測算該流域2000年、2010年、2020年的生態(tài)足跡、生態(tài)承載力以及生態(tài)盈余（赤字）狀況，并從時間和空間上分析其演化特征。結果表明， 2000—2020年白洋淀流域整體以及大部分區(qū)縣的人均生態(tài)足跡呈下降趨勢，且具有南部高、西部低的格局；2000—2020年流域整體的人均生態(tài)承載力波動較小且存在空間分異，西北部人均生態(tài)承載力較高，中部、東南部人均生態(tài)承載力較低；2000—2020年流域整體處于生態(tài)赤字狀態(tài)，但生態(tài)赤字水平在逐漸好轉，西北地區(qū)生態(tài)盈余，中東部生態(tài)赤字，經(jīng)濟發(fā)展水平與生態(tài)盈余（赤字）緊密相連。

關鍵詞：生態(tài)承載力；時空格局；生態(tài)足跡；白洋淀流域

中圖分類號：X826" " " " "文獻標識碼：A

文章編號：0439-8114（2024）12-0040-06

DOI：10.14088/j.cnki.issn0439-8114.2024.12.007 開放科學（資源服務）標識碼（OSID）：

Evaluation and spatial-temporal pattern of ecological carrying capacity in Baiyangdian Basin based on ecological footprint model

LIU Xiao1， GAO Xing1， YUAN Zi-hua1， GAO Xiao-chen2

（1.School of Public Administration， Hebei University of Economics and Business， Shijiazhuang" 050061， China；

2. Hebei Province Real Estate Registration Center， Shijiazhuang" 050031， China）

Abstract： Taking Baiyangdian Basin as the research area， based on the ecological footprint model， combined with statistical yearbook data and remote sensing image data， the ecological footprint， ecological carrying capacity and ecological surplus （deficit） of the basin in 2000， 2010 and 2020 were estimated， and their evolution characteristics were analyzed from time and space. The results showed that， from 2000 to 2020， the ecological footprint of Baiyangdian Basin as a whole and most districts and counties showed a downward trend， with a pattern of high in the south and low in the west. From 2000 to 2020， the ecological carrying capacity of the basin as a whole fluctuated little and had spatial differentiation. The per capita ecological carrying capacity in northwestern areas was relatively high， while the per capita ecological carrying capacity in central and southeastern areas was relatively low. From 2000 to 2020， the whole basin was in a state of ecological deficit， but the level of ecological deficit was gradually improving. The ecological surplus in northwest China and the ecological deficit in central and eastern China were closely related to the economic development level.

Key words： ecological carrying capacity； spatial-temporal pattern； ecological footprint； Baiyangdian Basin

收稿日期：2024-03-12

基金項目：河北省社會科學基金項目（HB22GL030）；河北省研究生創(chuàng)新資助項目

作者簡介：劉 瀟（2002-），女，河北保定人，在讀碩士研究生，研究方向為土地資源利用與評價，（電話）18032064803（電子信箱）18032064803@163.com；通信作者，高 星（1985-），男，河北石家莊人，副教授，博士，主要從事土地資源利用與評價研究，（電話）13552525085（電子信箱）gaoxing_85@126.com。

生態(tài)承載力是指在特定的時空范圍內，在資源環(huán)境和社會經(jīng)濟等協(xié)調可持續(xù)發(fā)展的前提下，自然-經(jīng)濟-社會復合生態(tài)系統(tǒng)對人類活動的最大承受能力［1］。它是評價一個地區(qū)人口、資源、環(huán)境可持續(xù)發(fā)展的重要指標［2］，其強弱受到土地利用變化的深刻影響。土地利用變化是人類活動的直觀表現(xiàn)，產(chǎn)業(yè)構成的變動、新興產(chǎn)業(yè)的涌入、地方政策等都可以使土地利用類型發(fā)生巨大改變［3］。雄安新區(qū)建設旨在疏解北京非首都功能、推進京津冀協(xié)同發(fā)展，該過程中必然伴隨著區(qū)域內外大規(guī)模的土地利用變化。與當?shù)厣鷳B(tài)環(huán)境不相適宜的土地利用變化勢必會影響區(qū)域的生態(tài)承載能力，引發(fā)生態(tài)環(huán)境風險［4］。因此，評價白洋淀流域的生態(tài)足跡、生態(tài)承載力以及生態(tài)赤字演變格局有助于推進白洋淀全流域治理和保護，對守住雄安新區(qū)生態(tài)環(huán)境安全邊界具有重要意義。

生態(tài)承載力是從承載力、生態(tài)足跡演變發(fā)展而來的，經(jīng)歷了資源承載力、環(huán)境承載力和生態(tài)承載力等階段［5］。有關生態(tài)承載力的研究已經(jīng)成為多學科的研究熱點。目前生態(tài)承載力的研究集中在時空演變［6，7］和驅動因素分析［8-10］兩個方面，多以水生態(tài)承載力變化及其驅動因素探析等為研究重點［11］，與土地利用變化相結合的文獻占比較小。從研究尺度來看，國家［12］、城市群［13］、流域［14］、省域［15］、市域［16］等尺度均有涉及。從研究方法來看，當前學界對生態(tài)承載力的測算主要采用生態(tài)足跡法［17］、綜合指標評價法［18］和系統(tǒng)動力學方法［19］。其中，生態(tài)足跡模型由加拿大學者Rees［20］首次提出，該模型因自身方法簡明、可操作性強，在流域生態(tài)壓力評價領域得到廣泛應用［21］，被國內外學術界認為是比較生態(tài)承載力供需平衡的一種簡單而有效的方法。因此本研究采用生態(tài)足跡模型測算研究區(qū)的生態(tài)承載力與生態(tài)足跡并對其進行評價，以便更精確地分析生態(tài)承載力演化特征。

近年來，隨著雄安新區(qū)的快速城鎮(zhèn)化建設，白洋淀流域的生態(tài)足跡與生態(tài)承載力勢必受到影響。因此，本研究選取白洋淀流域為研究區(qū)，利用“3S”技術獲得該流域2000年、2010年、2020年3期土地利用狀況，在此基礎上運用生態(tài)足跡模型分析該流域2000—2020年生態(tài)足跡、生態(tài)承載力、生態(tài)盈虧狀況的時空格局演變特征，以期為流域生態(tài)環(huán)境保護與可持續(xù)發(fā)展提供科學依據(jù)。

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

1.1 研究區(qū)概況

本研究中的白洋淀流域是指流入白洋淀的9條河流所流經(jīng)的河北省縣級行政區(qū)域，涵蓋保定市、滄州市、張家口市、石家莊市、衡水市5個市所轄的35個區(qū)縣。它以白洋淀水系為紐帶、以各區(qū)縣為基本單元，兼顧了自然、人文、區(qū)縣的宏觀發(fā)展。白洋淀流域地貌類型多樣，以耕地、林地和建設用地為主，區(qū)域內地形復雜多樣，主要有平原、山脈、丘陵和洼地，地勢西高東低，高差達2 862 m（圖1）。流域總面積約3.4萬km2，占河北省總面積的18.1%，總人口約為1 373.5萬人，占全省人口的18.5%。流域內擁有華北地區(qū)最大的淡水湖泊，土壤肥沃、生態(tài)資源豐富，樹木多為喬木，農作物以玉米和小麥為主，經(jīng)濟作物主要為棉花，主要產(chǎn)業(yè)為輕工業(yè)、水產(chǎn)以及旅游業(yè)。近年來，在加快建設“天藍、地綠、水秀的美麗河北”以及“生態(tài)優(yōu)先、綠色發(fā)展”的背景下，評價流域生態(tài)承載力守住生態(tài)環(huán)境安全邊界尤為重要。

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

所使用數(shù)據(jù)包括DEM高程數(shù)據(jù)以及2000年、2010年、2020年遙感影像數(shù)據(jù)、土地利用數(shù)據(jù)和各區(qū)縣社會經(jīng)濟數(shù)據(jù)。其中，DEM高程數(shù)據(jù)和2000年、2010年、2020年遙感影像數(shù)據(jù)來源于地理空間數(shù)據(jù)云。利用ENVI5.3軟件，結合研究區(qū)實際情況，對遙感影像進行監(jiān)督分類和人工目視解譯，將其分為6個一級地類，用于研究使用。土地利用數(shù)據(jù)及社會經(jīng)濟數(shù)據(jù)來源于《第三次國土調查主要數(shù)據(jù)公報》《中國統(tǒng)計年鑒》（2000—2020年）、《河北省統(tǒng)計年鑒》《保定市經(jīng)濟統(tǒng)計年鑒》《河北農村統(tǒng)計年鑒》《中國縣域統(tǒng)計年鑒》和《全國漁業(yè)發(fā)展統(tǒng)計公報》。由于生產(chǎn)與消費的不對等性、貿易的跨區(qū)域流動和不同類型土地生物資源的消耗總量數(shù)據(jù)收集存在困難，在計算生態(tài)足跡時，采用各類生物資源的生產(chǎn)總量近似于其消耗總量。

2 研究方法

1992年生態(tài)足跡模型首次提出，起初該模型的核心是一定區(qū)域內人類對生態(tài)生產(chǎn)性土地面積的需求；1997年二維生態(tài)足跡模型引入生態(tài)承載力的概念，對比生態(tài)承載力與生態(tài)足跡獲得區(qū)域生態(tài)盈余或生態(tài)赤字情況［22］。生態(tài)足跡模型的常用單位是“全球公頃”，但其無法精確反映國家尺度下的實際生產(chǎn)力狀況［23］。針對傳統(tǒng)生態(tài)足跡模型的不足，本研究用“國家公頃”代替“全球公頃”來計算區(qū)域生態(tài)足跡與生態(tài)承載力。另外依據(jù)流域內的實際生產(chǎn)，耕地消耗產(chǎn)品為糧食、棉花、油料、蔬菜，林地為園林水果，水域為水產(chǎn)品，建設用地為電力，草地為肉類和禽蛋（表1）。

2.1 生態(tài)足跡測算

個人或區(qū)域生態(tài)足跡是指生產(chǎn)這些人口所消費的資源和吸納這些人口所產(chǎn)生的廢棄物所需要的生物生產(chǎn)性土地總面積和水資源總量，其計算過程如下［24］。

首先，計算各類消費項目占用的土地面積。

[Si=Ci/Yii=1，2，…，6] （1）

式中，i為消費項目類型；Si為各消費項目消耗的生物生產(chǎn)面積；Ci為i項消費總量；Yi為i項消費項目的國家年平均產(chǎn)量。

其次，計算占用的某一類型土地的總面積。

[Sj=Sij] （2）

式中，j為土地類型；Sj為折算的各類土地利用類型的總面積；Sij為第i種消費項目占用的第j類土地利用類型的面積。

最后，計算人均生態(tài)足跡。

[ef=（Sj×fj）/p] （3）

式中，ef為人均生態(tài)足跡；fj為均衡因子，是指生態(tài)系統(tǒng)中各要素之間的平衡狀態(tài)，本研究使用河北省均衡因子和產(chǎn)量因子［25］，詳見表2；p為總人口。

2.2 生態(tài)承載力測算

生態(tài)承載力是指區(qū)域所提供給人類生物生產(chǎn)型土地（包括水域）的面積［21］，其計算式如下。

[ec=ai×fj×yj] （4）

式中，ec為人均生態(tài)承載力；aj為人均第i種資源的生產(chǎn)性土地面積；yj為產(chǎn)量因子，即某個地區(qū)某類生產(chǎn)性土地的平均生產(chǎn)力與世界同類土地平均生產(chǎn)力的比率。

2.3 生態(tài)盈余（赤字）

在三維模型中存在以下關系：

[ed=ec-ef] （5）

式中，ed為人均生態(tài)盈余（赤字），正值表示盈余，負值表示赤字。

3 結果與分析

3.1 生態(tài)足跡

根據(jù)生態(tài)足跡測算的公式，得到白洋淀流域3期人均生態(tài)足跡為0.421～0.785 hm2，在時間序列上呈逐年下降的趨勢。相較于2000年，2010年白洋淀流域整體人均生態(tài)足跡由0.785 hm2下降至0.586 hm2，下降了25.35%；2020年流域整體人均生態(tài)足跡進一步降至0.421 hm2，相較于2010年下降了28.16%。流域整體生態(tài)足跡的降低得益于各類土地人均生態(tài)足跡的降低（圖2）。耕地人均生態(tài)足跡呈先減少后上升的趨勢，由2000年的0.19 hm2下降至2010年的0.17 hm2，后又上升至2020年的0.19" hm2，整體波動不大；林地、草地、水域和建設用地的人均生態(tài)足跡呈持續(xù)下降趨勢，其中草地變化較大，由2000年的0.57 hm2下降至2010年的0.41 hm2，再到2020年的0.22 hm2，但在6個地類的人均生態(tài)足跡中仍居首位。

各區(qū)縣的生態(tài)足跡整體呈持續(xù)下降趨勢，僅有任丘市、保定市清苑區(qū)、安平縣、深澤縣的人均生態(tài)足跡整體呈上升趨勢（圖3）。同時，還有部分區(qū)縣，如蔚縣、涿鹿縣、定興縣的人均生態(tài)足跡出現(xiàn)波動，呈先上升后下降的趨勢。流域整體生態(tài)足跡的下降是好的發(fā)展趨勢，這代表著流域內居民對自然資源的消耗降低，人類活動對土地生態(tài)系統(tǒng)的消耗也在逐漸減少。

白洋淀流域各區(qū)縣人均生態(tài)足跡呈南部高、西部低的格局，流域整體人均生態(tài)足跡下降趨勢明顯。2000—2020年中東部地區(qū)生態(tài)足跡變化明顯，如易縣人均生態(tài)足跡由2000年的0.59 hm2下降至2020年的0.29 hm2。西部地區(qū)始終保持較低的人均生態(tài)足跡。

3.2 生態(tài)承載力

通過計算得到流域整體的人均生態(tài)承載力呈下降趨勢，但是整體波動較小，3期人均生態(tài)承載力為0.22～0.26 hm2。同時，大部分區(qū)縣的人均生態(tài)承載力在時間序列上都是下降的，僅有淶源縣、定興縣、滿城區(qū)、保定市清苑區(qū)、深澤縣、無極縣的人均生態(tài)承載力是上升的，下降較明顯的為雄縣，這與雄安新區(qū)建設有一定關系。每類土地的生態(tài)承載力變動都會對流域整體的生態(tài)承載力產(chǎn)生影響，建設用地和草地的人均生態(tài)承載力呈逐年上升趨勢，建設用地的人均生態(tài)承載力由0.037 hm2上升到0.063 hm2，草地的人均生態(tài)承載力由0.006 hm2上升到0.014 hm2（圖4）。耕地、林地的人均生態(tài)承載力均下降，耕地下降幅度較大，林地下降幅度非常??；水域的人均生態(tài)承載力先下降后上升，這與白洋淀流域河流生態(tài)改善有很大的關系。隨著經(jīng)濟社會的進步與發(fā)展，人們對于自然的開發(fā)與利用逐漸增加，對建設用地的需求也在增加，因此流域建設用地的生態(tài)承載力一直保持上升。由于流域內土地數(shù)量有限，建設用地面積的大幅增加勢必會導致其他類型土地面積的減少，從而也會對生態(tài)承載力產(chǎn)生影響。在上述5種土地類型中，流域內耕地生態(tài)承載力占比仍為最大，因此，要想提高流域內的整體生態(tài)承載力，勢必要加大對流域耕地的保護。

從空間格局上看，白洋淀流域各區(qū)縣人均生態(tài)承載力呈西北高、東南低的分布格局，且生態(tài)承載力在雄縣、安新縣、任丘市、望都縣變化較大。各區(qū)縣的人均生態(tài)承載力大致可以分為兩個等級（圖4）。淶源縣、蔚縣、阜平縣、涿鹿縣的人均生態(tài)承載力較高，介于0.58～0.76 hm2。這些縣域均位于白洋淀流域上游、太行山一帶，海拔相對較高，多為山地，生態(tài)盈余較高。剩下的縣域人均生態(tài)承載力較低，介于0.08～0.39 hm2。

3.3 生態(tài)盈余（赤字）

白洋淀流域整體生態(tài)赤字，但是生態(tài)赤字在逐漸好轉，人均生態(tài)赤字由2000年的0.52 hm2下降到0.20 hm2，說明近年來流域的生態(tài)治理工作取得了一定的成效，但仍有很大的發(fā)展空間。

由圖6可知，2000—2020年各區(qū)縣的生態(tài)盈余（赤字）情況呈現(xiàn)空間分化，流域西北部的區(qū)縣處于生態(tài)盈余，其中2020年處于生態(tài)盈余的縣有蔚縣、涿鹿縣、淶源縣、唐縣、易縣、淶水縣、阜平縣、曲陽縣、高陽縣，約占全部區(qū)縣的25%，另外75%的區(qū)縣全部處于生態(tài)赤字。其主要原因是白洋淀流域上游海拔較高，多為山地，經(jīng)濟發(fā)展較緩慢，人口稀疏，生態(tài)環(huán)境質量較高；而流域下游地處華北平原，人口聚集，城鎮(zhèn)建設用地較多，產(chǎn)業(yè)以工業(yè)為主，因此生態(tài)赤字水平較高。從時間動態(tài)變化來看，2020年大部分區(qū)縣的生態(tài)環(huán)境都有一定程度的好轉。生態(tài)盈余的區(qū)縣數(shù)量呈上升趨勢，2000年生態(tài)盈余的縣只有5個，2020年上升為9個，這與當?shù)卣陙淼纳鷳B(tài)保護政策密不可分，但是上升最大的石家莊市藁城區(qū)和正定縣仍處于生態(tài)赤字。但也有部分區(qū)縣生態(tài)赤字狀況逐漸加深，占比為11.1%，下降較明顯的有3個，分別是深澤縣、任丘市和安平縣。因此，各區(qū)縣在發(fā)展過程中應注重經(jīng)濟發(fā)展與生態(tài)保護相結合，走可持續(xù)發(fā)展的道路。

4 小結

生態(tài)環(huán)境是社會經(jīng)濟發(fā)展的根基，生態(tài)承載力評價是了解區(qū)域生態(tài)環(huán)境、制定后續(xù)發(fā)展規(guī)劃的重要依據(jù)。本研究基于生態(tài)足跡模型對白洋淀流域的生態(tài)承載力進行時空評價，得到以下結論。

1）2000—2020年白洋淀流域整體以及大部分區(qū)縣的人均生態(tài)足跡都呈下降趨勢，各地類（除耕地外）的生態(tài)足跡也呈下降趨勢，流域內居民對自然資源的消耗逐漸減少；從空間上看，人均生態(tài)足跡呈南部高、西部低的格局。

2）2000—2020年白洋淀流域整體的人均生態(tài)承載力波動較小，大部分區(qū)縣的人均生態(tài)承載力都呈下降趨勢。耕地、林地的人均生態(tài)承載力都在下降，應加強對這2種地類的生態(tài)保護；從空間上看，流域人均生態(tài)承載力可以分為2個等級，且存在空間分異，西北太行山一帶人均生態(tài)承載力較高，中部、東南部地區(qū)人均生態(tài)承載力較低。

3）2000—2020年流域整體一直處于生態(tài)赤字，但人均生態(tài)赤字在逐漸好轉。從空間上看，西北地區(qū)生態(tài)盈余，中東部生態(tài)赤字，經(jīng)濟發(fā)展水平與生態(tài)盈余（赤字）緊密相連。

綜合來看，白洋淀流域的生態(tài)環(huán)境雖受到一定破壞，但是經(jīng)過政府的治理和保護以及居民生態(tài)保護意識的加強，流域生態(tài)環(huán)境有了很大的改善，處于可持續(xù)發(fā)展狀態(tài)。

參考文獻：

［1］ 李小勇，賈海發(fā)，趙萬梅. 基于 DPSIR-TOPSIS 模型的青海省生態(tài)承載力評價［J］. 水土保持通報，2022，42（2）： 166-173，183.

［2］ LIU D， FENG Z M， YANG Y Z， et al. Spatial patterns of ecological carrying capacity supply-demand balance in China at county level［J］. Journal of geographical sciences， 2011， 21（5）： 833-844.

［3］ 胡唐松，吳鑫淼，王曉夢，等.近十年雄安新區(qū)景觀格局動態(tài)變化及生態(tài)系統(tǒng)服務功能評價［J］.林業(yè)與生態(tài)科學，2024，39（1）：66-75.

［4］ 高 星，劉澤偉，李晨曦，等.基于“三生空間”的雄安新區(qū)土地利用功能轉型與生態(tài)環(huán)境效應研究［J］.生態(tài)學報，2020，40（20）：7113-7122.

［5］ 趙東升，郭彩贇，鄭 度，等. 生態(tài)承載力研究進展［J］. 生態(tài)學報，2019，39（2）： 399-410.

［6］ 熊建新，王鑫濱，趙 迪，等.城鎮(zhèn)化進程中生態(tài)承載力系統(tǒng)耦合時空格局及影響因素——以洞庭湖區(qū)為例［J］.經(jīng)濟地理，2022，42（10）：83-91.

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

［8］ 徐智超，劉華民，韓 鵬，等.內蒙古生態(tài)安全時空演變特征及驅動力［J］.生態(tài)學報，2021，41（11）： 4354-4366.

［9］ 楊雪荻，白永平，車 磊，等.甘肅省生態(tài)安全時空演變特征及影響因素解析［J］.生態(tài)學報，2020，40（14）：4785-4793.

［10］ 董建紅，張志斌，張文斌. 基于三維生態(tài)足跡的甘肅省自然資本利用動態(tài)變化及驅動力［J］.生態(tài)學雜志，2019，38（10）： 3075-3085.

［11］ 段海煦，曾維華，陳家軍. 基于可視化分析的“生態(tài)、環(huán)境、資源”承載力研究［J］.中國環(huán)境科學， 2023， 43（9）： 5031-5040.

［12］ 李 旭，古麗·加帕爾，于 濤，等.中亞5國生態(tài)承載力時空動態(tài)監(jiān)測與驅動因素分析［J］.中國科學院大學學報，2023，" " " " 40（4）：474-485.

［13］ 王世豪，黃 麟，徐新良，等. 特大城市群生態(tài)空間及其生態(tài)承載狀態(tài)的時空分異［J］.地理學報，2022，77（1）： 164-181.

［14］ 張永芳，周 斌，李 蛟，等. 黃河流域生態(tài)承載力及人地關系可持續(xù)發(fā)展——基于異速生長模型的研究［J］. 西北人口，2024，45（4）：112-126.

［15］ 柯志成，連海峰，陳 奕，等. 福建省2009—2018年人均生態(tài)足跡和人均生態(tài)承載力的時空動態(tài)［J］.福建農林大學學報（自然科學版），2021，50（5）： 677-685.

［16］ 樂 暉，金文龍，陳亢利. 長三角地區(qū)生態(tài)承載力評價——以蘇州市為例［J］.環(huán)境工程技術學報， 2023， 13（2）： 725-732.

［17］ 王 錦，韓陽媚，叢 林，等. 基于模型參數(shù)本地化修正的河北省生態(tài)足跡演化特征及可持續(xù)發(fā)展分析［J］.干旱區(qū)資源與環(huán)境，2023，37（2）： 42-49.

［18］ 鄒 游，王夢琪，李桂媛.生態(tài)視角下旅游活動的環(huán)境承載力影響及評價研究——以湖北省隨州市為例［J］.生態(tài)經(jīng)濟，2023，39（12）： 139-146.

［19］ 秦歡歡. 華北平原水資源生態(tài)足跡與生態(tài)承載力評估［J］.科技導報，2023，41（11）： 41-51.

［20］ REES W E. Ecological footprints and appropriated carrying capacity： What urban economics leaves out［J］. Environment and urbanization，1992， 4（2）：121-130.

［21］ 葉進紅，李鮮苗，杜柏鋒. 黃河源區(qū)生態(tài)壓力分析：以甘南為例［J］.河北環(huán)境工程學院學報，2023，33（2）：67-72.

［22］ 徐先斌，劉 險，彭 穎.基于生態(tài)足跡模型評價長江經(jīng)濟帶耕地承載力［J］.環(huán)境科學與技術，2023，46（S2）： 210-216.

［23］ 吳開亞，王玲杰.基于全球公頃和國家公頃的生態(tài)足跡核算差異分析［J］.中國人口·資源與環(huán)境，2007（5）：80-83.

［24］ 姚林杰，張佳寧，張恒瑞，等.基于生態(tài)足跡的內蒙古自然資源可持續(xù)利用分析［J］.生態(tài)環(huán)境學報，2023，32（8）： 1525-1536.

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