基于生態(tài)網(wǎng)絡(luò)分析的黃河三角洲碳代謝分析

路 昌,王之語

基于生態(tài)網(wǎng)絡(luò)分析的黃河三角洲碳代謝分析

路 昌*,王之語

(山東建筑大學(xué)管理工程學(xué)院,山東 濟南 250101)

為深入探討黃河三角洲碳代謝生態(tài)關(guān)系與國土空間利用的內(nèi)在聯(lián)系,為國土空間格局調(diào)整和綠色發(fā)展提供參考,構(gòu)建“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間分類體系,利用生態(tài)網(wǎng)絡(luò)分析模型(ENA)對黃河三角洲2000～2020年的碳代謝過程進行分析,并運用標(biāo)準(zhǔn)差橢圓和生態(tài)關(guān)系變化率(ERCR)評價各空間之間的生態(tài)關(guān)系時空分布和變化特征.結(jié)果表明,黃河三角洲2000～2015年生態(tài)空間被農(nóng)業(yè)空間占用等現(xiàn)象頻發(fā),碳代謝呈紊亂局面,2015年后正碳流顯著提高,碳代謝趨于平衡;研究區(qū)主導(dǎo)生態(tài)關(guān)系為競爭關(guān)系,碳代謝生態(tài)沖突較為強烈,2015年后西北地區(qū)互惠共生關(guān)系逐漸增多,國土空間演變重點區(qū)域由南向北轉(zhuǎn)變;研究區(qū)大部分生態(tài)關(guān)系變化率滿足ERCR>0,生態(tài)關(guān)系規(guī)模增減與空間轉(zhuǎn)移面積的變化趨勢一致.說明黃河三角洲國土空間仍有不合理利用現(xiàn)象,碳代謝壓力尚存,需進一步優(yōu)化國土空間格局,促進綠色發(fā)展.

“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間；碳代謝；生態(tài)網(wǎng)絡(luò)分析模型；生態(tài)關(guān)系變化率；黃河三角洲

綠色、低碳的發(fā)展方式既是習(xí)近平生態(tài)文明思想的主要內(nèi)容,也是貫徹新發(fā)展理念、推動高質(zhì)量發(fā)展的必經(jīng)之路[1].21世紀(jì)以來,城鎮(zhèn)化進程快步推進,城市面積不斷擴張,國土空間利用方式的變化帶來了極其活躍的城市碳代謝過程.由于“碳源”的增加,“碳匯”的有限,溫室氣體逐漸增多,全球變暖問題愈演愈烈,生態(tài)環(huán)境面臨著前所未有的挑戰(zhàn),氣候異常甚至伴隨嚴(yán)重的蟲災(zāi),導(dǎo)致全球糧食產(chǎn)量下降[2-3].因此,探尋國土空間利用與碳代謝的聯(lián)系,優(yōu)化國土空間格局,對于生態(tài)環(huán)境的治理維護和發(fā)展方式“綠色轉(zhuǎn)型”目標(biāo)的實現(xiàn)具有重要意義.

“碳代謝”是指碳在生物圈、大氣及城市系統(tǒng)的不同部分之間流動,包括碳的吸收和排放[4].學(xué)術(shù)界通過建立定量模型的方法對碳代謝及相關(guān)問題進行了廣泛的研究.Sovacool等[5]和Karakiewicz等[6]最早將碳作為單一研究對象,并與“代謝”概念結(jié)合研究城市碳代謝過程,但其研究僅局限于基礎(chǔ)的流量分析;直到Patten和Fath等[7-8]提出并推廣的生態(tài)網(wǎng)絡(luò)分析(ENA)模型被應(yīng)用到碳代謝研究領(lǐng)域,為其復(fù)雜的機理過程提供了網(wǎng)絡(luò)視角的分析方法;Zhang和Xia等[9-10]將土地利用引入碳代謝研究,運用ENA模型分析了土地利用/覆被變化與碳代謝的關(guān)聯(lián)性,將碳代謝的理論分析與土地實際結(jié)合,提出了研究碳減排問題的新視角.

過往研究對于碳代謝過程與“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間體系下的國土空間利用相關(guān)研究較少,且對碳代謝生態(tài)關(guān)系與國土空間演變的動態(tài)分析有所欠缺[11].本研究以黃河三角洲為研究區(qū),構(gòu)建“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間分類體系,對其2000～2020年由于國土空間變化而產(chǎn)生的碳代謝過程進行深入分析,對碳代謝生態(tài)關(guān)系的時空分布和變化特征作探討和研究,為黃河三角洲國土空間格局優(yōu)化提出科學(xué)指導(dǎo)和理論依據(jù),為其生產(chǎn)生活方式“綠色轉(zhuǎn)型”目標(biāo)的實現(xiàn)提供參考.

1 數(shù)據(jù)來源與研究方法

1.1 研究區(qū)概況

黃河三角洲(117°31′E～119°18′E,36°55′N～ 38°16′N)(圖1)位于渤海灣南岸和萊州灣西岸,與日本、韓國隔海相望,北接京津唐經(jīng)濟區(qū),南連山東半島藍(lán)色經(jīng)濟區(qū),向西輻射廣大內(nèi)陸地區(qū),是環(huán)渤海經(jīng)濟區(qū)的重要節(jié)點、山東半島城市群的重要組成部分,處于連接中原經(jīng)濟區(qū)與東北經(jīng)濟區(qū)、京津唐經(jīng)濟區(qū)與膠東半島經(jīng)濟區(qū)的樞紐位置[12].本研究選取東營市的東營區(qū)、河口區(qū)、墾利區(qū)和利津縣為研究區(qū),總面積約為6360km2.研究區(qū)屬典型的暖溫帶半濕潤季風(fēng)氣候區(qū),四季分明,雨熱同季,多年平均氣溫均大于10℃,多年平均降水量約為585mm.研究區(qū)地處環(huán)渤海經(jīng)濟圈,社會經(jīng)濟增長動能明顯加強,國土空間供需矛盾日益加劇,社會生態(tài)系統(tǒng)失衡,對區(qū)域人地關(guān)系的調(diào)節(jié)與優(yōu)化刻不容緩[13].

圖1 研究區(qū)位置示意

1.2 數(shù)據(jù)來源

本研究所使用的空間數(shù)據(jù)為中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(https://www.resdc.cn/)獲取的中國多時期土地利用土地覆被遙感監(jiān)測數(shù)據(jù)集(CNLUCC)[14],空間分辨率為30m×30m.研究區(qū)人口、農(nóng)業(yè)、能源等數(shù)據(jù)來源于2001年、2006年、2011年、2016年及2021年的《東營統(tǒng)計年鑒》、《中國縣域統(tǒng)計年鑒(縣市篇)》,工業(yè)、交通運輸?shù)葦?shù)據(jù)來源于東營市各年度各區(qū)縣的國民經(jīng)濟和社會發(fā)展統(tǒng)計公報以及中華人民共和國交通運輸部網(wǎng)站(https://www.mot.gov.cn/).對于部分缺失和異常的能源數(shù)據(jù)和港口吞吐量數(shù)據(jù)采用插值的方法進行了補齊和修正.

1.3 研究方法

選取黃河三角洲作為研究區(qū),基于已有研究構(gòu)建“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間分類體系,對研究區(qū)2000～2020年的國土空間變化產(chǎn)生的碳代謝過程利用ENA模型進行分析,研究不同國土空間之間的生態(tài)關(guān)系,評價不同國土空間相互轉(zhuǎn)移對研究區(qū)碳代謝的綜合作用情況,并利用標(biāo)準(zhǔn)差橢圓和生態(tài)關(guān)系變化率分析各類生態(tài)關(guān)系的時空分布和變化特征(圖2).

圖2 研究框架

1.3.1 構(gòu)建“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間分類體系 對國土空間依據(jù)其利用方式進行分類,可以反映人類活動對國土利用的影響作用,有助于國土空間利用管理和科研工作的開展.對于某一地域空間來說,其在功能上具有多宜性,但總歸有一種用途發(fā)揮著主導(dǎo)作用,這種主導(dǎo)作用一般是依據(jù)空間利用主體的主觀意愿所確定的.因此以人類活動對于國土空間功能的依賴性和差異性為基礎(chǔ),結(jié)合十八大提出的國土空間“三生”(生態(tài)、生產(chǎn)、生活)發(fā)展目標(biāo)和中共中央、國務(wù)院提出的“三區(qū)”(生態(tài)、農(nóng)業(yè)、城鎮(zhèn)區(qū))劃定要求[15],并考慮數(shù)據(jù)的可獲得性,劃分出黃河三角洲的“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間分類體系(表1).考慮到研究區(qū)分布大量鹽田,根據(jù)2020年自然資源部辦公廳印發(fā)的《國土空間調(diào)查、規(guī)劃、用途管制用地用海分類指南(試行)》將其歸為工礦用地,屬于其他建設(shè)用地,故在本研究的“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間分類體系中劃分為城鎮(zhèn)建設(shè)空間.

表1 黃河三角洲“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間分類體系

1.3.2 碳代謝計算 碳代謝計算主要包括對國土空間利用變化過程中產(chǎn)生的碳吸收量、碳排放量及碳流量進行計算.目前學(xué)術(shù)界尚無統(tǒng)一核算方法,本研究采用碳轉(zhuǎn)化系數(shù)法[16]進行計算.

(1)碳吸收量.對國土空間碳吸收量的計算就是對“碳匯”的碳吸收量進行計算.基于黃河三角洲的國土空間分類體系,需要考慮其碳吸收量的國土空間類型包括:林業(yè)生態(tài)空間、草地生態(tài)空間、水域生態(tài)空間、其他生態(tài)空間和農(nóng)業(yè)生產(chǎn)空間.計算公式為:

式中:為碳吸收量;為國土空間面積;為碳吸收系數(shù)[17-18](表2).

表2 各類空間碳吸收系數(shù)(t/hm2)

(2)碳排放量.碳排放量的計算是研究碳代謝過程的重要環(huán)節(jié),本研究的碳排放量主要考慮對具有顯著碳排放過程的國土空間類型進行計算,包括:水域生態(tài)空間、農(nóng)業(yè)生產(chǎn)空間、農(nóng)村生活空間和城鎮(zhèn)建設(shè)空間.其中水域生態(tài)空間參考舒心等[19]的研究成果,選取東營港的港口吞吐量為計算指標(biāo);農(nóng)業(yè)生產(chǎn)空間考慮各種農(nóng)業(yè)活動中產(chǎn)生的機械能耗和能源消耗,選取農(nóng)機能耗量、灌溉面積、農(nóng)用柴油量、農(nóng)業(yè)用電量、化肥施用量、農(nóng)膜使用量、農(nóng)藥使用量和水稻種植面積作為計算指標(biāo),考慮研究區(qū)的主要作物為小麥和玉米,相關(guān)研究[20]表明,小麥和玉米光合作用所消耗的CO2遠(yuǎn)大于呼吸作用排出的CO2,一般不會引起CO2濃度增加,因此本研究不考慮小麥和玉米種植的碳排放量;農(nóng)村生活空間考慮農(nóng)村居民的生活能耗和呼吸帶來的碳排放,以及農(nóng)村居民在飼養(yǎng)家畜過程中產(chǎn)生的代謝形成的碳排放,選取農(nóng)村居民的生活用電量、家用天然氣量、家用液化石油氣量、農(nóng)村居民人口數(shù)以及豬、牛、羊、家禽的飼養(yǎng)量為計算指標(biāo);城鎮(zhèn)建設(shè)空間考慮城鎮(zhèn)居民的日常生活能耗和呼吸帶來的碳排放,以及工業(yè)生產(chǎn)消耗能源所產(chǎn)生的碳排放量,選取計算指標(biāo)為城鎮(zhèn)居民的生活用電量、家用天然氣量、家用液化石油氣量、城鎮(zhèn)居民人口數(shù)、汽車年行駛里程[21]以及各種工業(yè)能源消耗量.計算公式為:

式中:WE為水域生態(tài)空間碳排放量;PT為東營港港口吞吐量;PT為港口吞吐的碳排放系數(shù);AP為農(nóng)業(yè)生產(chǎn)空間碳排放量;K(=AM, IA, AD, AE, CF, AF, AC)為對應(yīng)的碳排放系數(shù);AM為農(nóng)機能耗量;IA為灌溉面積;AD為農(nóng)用柴油量;AE為農(nóng)業(yè)用電量;CF為化肥施用量;AF為農(nóng)膜使用量;AC為農(nóng)藥使用量;RC為水稻種植面積;RL為農(nóng)村生活空間碳排放量;K(=CE, CN, CG, CP)為對應(yīng)的碳排放系數(shù);CE為農(nóng)村居民生活用電量;CN為農(nóng)村居民家庭用天然氣量;CG為農(nóng)村居民家庭用液化石油氣量;CP為農(nóng)村居民人口;1,2,3,4分別為豬、牛、羊、家禽的飼養(yǎng)量;B和E為對應(yīng)家畜的呼吸和排泄的碳排放系數(shù)(=1, 2, 3, 4);UB為城鎮(zhèn)建設(shè)空間碳排放量;K(=UE, UN, UG, UP, SV, OV)為對應(yīng)的碳排放系數(shù);UE為城鎮(zhèn)居民生活用電量;UN為城鎮(zhèn)居民家庭用天然氣量;UG為城鎮(zhèn)居民家庭用液化石油氣量;UP為城鎮(zhèn)居民人口;SV和OV分別為為小型車輛和大型車輛年行駛里程;P為各種工業(yè)能源消耗量;K為對應(yīng)的工業(yè)能源的碳排放系數(shù)(=1, 2, ..., 13).

以上所用碳排放系數(shù)(表3)等數(shù)據(jù)主要參考了文獻[4,19,22-25]的結(jié)果,以及IPCC (https://www. ipcc-nggip.iges.or.jp/)和中華人民共和國生態(tài)環(huán)境部(https://www.mee.gov.cn/)所提供的數(shù)據(jù).

表3 各類空間相關(guān)指標(biāo)碳排系數(shù)

續(xù)表3

(3)碳流量.參考李晶等[26]的研究方法,定義單位面積的凈碳排放量為凈碳代謝密度,以此量化不同空間的碳代謝能力,借助碳代謝密度差和國土空間轉(zhuǎn)移面積衡量國土空間變化過程中產(chǎn)生的碳流量及其方向.計算公式為:

式中:Δ是碳代謝密度差;W和W分別表示空間和的凈碳代謝密度;V和V分別表示空間和的凈碳排放量;a和a分別表示空間和的面積;f表示從空間向空間流動的碳流量;Δa表示空間向空間轉(zhuǎn)移的面積.一般地,f>0表示正碳流,說明碳代謝趨于平衡;f<0表示負(fù)碳流,說明碳代謝趨于紊亂.

1.3.3 生態(tài)網(wǎng)絡(luò)分析模型 生態(tài)網(wǎng)絡(luò)分析(ENA)模型廣泛應(yīng)用于環(huán)境科學(xué)與資源利用領(lǐng)域研究[27-28],它定量地描述了生態(tài)系統(tǒng)各個分室之間的物質(zhì)能量流動以及相互關(guān)系[10],包括流量分析和效用分析[29-30].本研究利用ENA模型描述國土空間碳代謝生態(tài)關(guān)系,反映兩類空間之間的碳代謝空間分布特征和相互作用關(guān)系.

(1)流量和效用分析.通過構(gòu)造直接效用矩陣和整體效用矩陣來完整反映空間之間直接和間接的碳代謝關(guān)系[9].計算公式為:

式中:為直接效用矩陣;f表示從空間向空間流動的碳流量;T表示空間的碳通量;為整體效用矩陣;為單位矩陣.

(2)生態(tài)關(guān)系判定.將整體效用矩陣的元素正負(fù)號作為元素,定義符號矩陣sgn()來定性判斷空間之間的相互作用關(guān)系,通常有互惠共生(+, +),競爭(-, -),掠奪(+, -)和限制(-, +)這4種關(guān)系.其中,互惠共生關(guān)系表示兩空間相互作用過程中均增加了效用;競爭關(guān)系則表示兩空間相互作用過程中均損失了效用;掠奪關(guān)系與限制關(guān)系總是“捆綁”出現(xiàn),當(dāng)一類空間利用另一類空間增加了效用,使得另一空間效用損失,則這一空間對于另一空間表現(xiàn)為掠奪關(guān)系,相應(yīng)的另一類空間對其則表現(xiàn)為限制關(guān)系[31].

1.3.4 標(biāo)準(zhǔn)差橢圓 標(biāo)準(zhǔn)差橢圓(SDE)在描述空間數(shù)據(jù)分布特征方面具有廣泛應(yīng)用[32].本研究利用標(biāo)準(zhǔn)差橢圓扁率進行全局分析,以反映研究區(qū)各生態(tài)關(guān)系空間分布的離散程度和方向性[33].標(biāo)準(zhǔn)差橢圓扁率的計算公式為:

式中:為標(biāo)準(zhǔn)差橢圓扁率;和分別為標(biāo)準(zhǔn)差橢圓長軸和短軸.通常越大,空間分布軸線性越明顯,方向性就越顯著;反之,則分布特征越離散.

1.3.5 生態(tài)關(guān)系變化分析 通過對比不同研究期及時段間兩空間演變的各生態(tài)關(guān)系數(shù)量和空間分布,直觀反映不同時段兩空間演變的主導(dǎo)生態(tài)關(guān)系類型及生態(tài)關(guān)系變化情況;同時參考城市擴展速率[34]的概念,構(gòu)建生態(tài)關(guān)系變化率(ERCR),用以描述不同時段兩空間演變的生態(tài)關(guān)系變化情況,反映生態(tài)關(guān)系與國土空間演變的動態(tài)關(guān)系.生態(tài)關(guān)系變化率ERCR的計算公式為:

式中:ΔERA為第(=1, 2, 3, 4)種生態(tài)關(guān)系分布面積的變化量;Δ2S為空間與空間(和分別為生態(tài)空間、農(nóng)業(yè)空間、城鎮(zhèn)空間中不同的兩類)相互轉(zhuǎn)移面積在某一時段的變化量;ERA(t)和ERA(t+1)表示第種生態(tài)關(guān)系分別在研究期和+1的分布面積;ΔS(t)和ΔS(t+1)表示空間與空間在研究期和+1的相互轉(zhuǎn)移面積(即空間向空間轉(zhuǎn)移的面積與空間向空間轉(zhuǎn)移的面積之和);Δa(t)表示研究期空間向空間轉(zhuǎn)移的面積(Δa(t)同理).

由于空間之間的生態(tài)關(guān)系是基于國土空間變化引起的碳流而存在的,即每一研究期四種生態(tài)關(guān)系的分布面積總和與這一時期國土空間相互轉(zhuǎn)移的總面積是相等的(公式12),這使得每兩類國土空間的四種生態(tài)關(guān)系的ERCR在每一時段和為1.一般地,ERCR<0表示該生態(tài)關(guān)系分布面積變化方向與國土空間相互轉(zhuǎn)移面積的變化方向不一致, ERCR>0表示該生態(tài)關(guān)系分布面積變化方向與國土空間相互轉(zhuǎn)移面積的變化方向一致.一種生態(tài)關(guān)系ERCR的正負(fù)往往受到其他生態(tài)關(guān)系ERCR的影響,這與國土空間變化方式及其引起的碳代謝特征密切相關(guān).

2 結(jié)果與分析

2.1 碳代謝時空演變

利用ArcMap 10.6獲取研究區(qū)4個研究時段的碳代謝強度及分布情況.2000～2005年(圖3(a)),大部分空間轉(zhuǎn)移所帶來的均是較強的負(fù)碳流,主要集中在研究區(qū)西北部和東部區(qū)域,來源于城鎮(zhèn)建設(shè)空間的緩慢擴張和草地生態(tài)空間向農(nóng)業(yè)生產(chǎn)空間的轉(zhuǎn)移,研究區(qū)整體正負(fù)碳流比約為1:35.97.2005～2010年(圖3(b)),負(fù)碳流的區(qū)域和強度均有所增加,增加的區(qū)域主要分布在研究區(qū)北部、東部和中南部,其中北部和中南部的新增負(fù)碳流主要來源于河口區(qū)和東營區(qū)城鎮(zhèn)建設(shè)空間的快速擴張,東部地區(qū)的新增負(fù)碳流主要集中在黃河三角洲自然保護區(qū),該區(qū)域草地生態(tài)空間向水域生態(tài)空間轉(zhuǎn)移,國土空間碳吸收能力變?nèi)?正碳流明顯增加,主要源于“退耕還濕,退農(nóng)還灘”政策,整體正負(fù)碳流比約為1:3.72.2010～ 2015年(圖3(c)),碳流分布零星且正負(fù)碳流都較弱,負(fù)碳流強度下降顯著,僅為上一研究期的5.38%,整體正負(fù)碳流比約為1:12.08.2015～2020年(圖3(d)),正碳流顯著增多,主要分布在研究區(qū)北部的河口區(qū),該區(qū)域?qū)嵤竦乇Ｗo修復(fù)工作,大量城鎮(zhèn)建設(shè)空間向水域生態(tài)空間轉(zhuǎn)移,形成的正碳流為研究期總正碳流的98.73%,整體正負(fù)碳流比約為1:0.18.

整體來看,2000～2015年,黃河三角洲主要注重經(jīng)濟發(fā)展,著力于城鎮(zhèn)建設(shè)空間的擴張,國土空間演變以消極方式為主,對黃河流域的生態(tài)治理不足以抵消其帶來的負(fù)碳流,研究區(qū)整體碳代謝紊亂; 2015～2020年,研究區(qū)更加注重生態(tài)安全問題,改變了之前消極的國土空間轉(zhuǎn)移方式,提高了水域生態(tài)空間的碳匯作用,研究區(qū)逐漸趨于碳代謝平衡(圖4).

圖3 2000～2020年黃河三角洲碳代謝空間分布

圖4 2000～2020年黃河三角洲碳代謝變化情況

2.2 生態(tài)關(guān)系判定及空間分布

2.2.1 生態(tài)關(guān)系判定 通過ENA模型對黃河三角洲4個研究期的國土空間碳代謝進行流量及效用核算,借助笛卡爾坐標(biāo)系4個象限的符號組合進行描述,得到其國土空間碳代謝生態(tài)關(guān)系判定結(jié)果(圖5).從數(shù)量上看,在4個研究期中占據(jù)主導(dǎo)地位的生態(tài)關(guān)系是競爭關(guān)系,占生態(tài)關(guān)系總數(shù)量的29.77%,其中2000～2005年占比最大,達到38.10%(表4),表明黃河三角洲在研究期內(nèi)不同空間類型之間的碳代謝沖突比較強烈;而互惠共生關(guān)系、掠奪關(guān)系和限制關(guān)系在整個研究期內(nèi)總體占比平衡,分別占23.81%、23.21%和23.21%.

圖5 2000～2020年黃河三角洲碳代謝生態(tài)關(guān)系判定

表4 碳代謝生態(tài)關(guān)系數(shù)量結(jié)構(gòu)(%)

2.2.2 生態(tài)關(guān)系空間分布 利用ArcMap 10.6將黃河三角洲4期的國土空間碳代謝生態(tài)關(guān)系空間可視化.2000～2005年(圖6(a)),掠奪關(guān)系主要分布于研究區(qū)北部和東部地區(qū),1.96萬hm2的其他生態(tài)空間向水域生態(tài)空間轉(zhuǎn)移,而同期水域生態(tài)空間向其他生態(tài)空間轉(zhuǎn)移僅為170.46hm2,由于兩空間碳匯能力相差不大,因此水域生態(tài)空間獲得了更多的碳儲量,而其他生態(tài)空間失去碳儲量,從而形成了掠奪、限制關(guān)系.掠奪關(guān)系標(biāo)準(zhǔn)差橢圓扁率為0.21(表5),軸線特征較弱,分布較為均衡.2005～2010年(圖6(b)),互惠共生關(guān)系居主導(dǎo)地位,主要分布在墾利區(qū)北部和南部、東營區(qū)北部以及河口區(qū)中部.互惠共生關(guān)系主要存在于農(nóng)業(yè)生產(chǎn)空間和其他生態(tài)空間之間,農(nóng)業(yè)生產(chǎn)空間產(chǎn)生的碳排放被其他生態(tài)空間所吸收,碳排放壓力的緩解有利于農(nóng)業(yè)生產(chǎn),二者形成互惠共生關(guān)系.互惠共生關(guān)系標(biāo)準(zhǔn)差橢圓扁率為0.60,軸線特征較為明顯,主要為“西南-東北”分布.2010～2015年(圖6(c)),出現(xiàn)了國土空間“單向轉(zhuǎn)移”現(xiàn)象,如林業(yè)生態(tài)空間向城鎮(zhèn)建設(shè)空間單向轉(zhuǎn)移,碳排放量大增,碳匯能力銳減,使得掠奪關(guān)系占主導(dǎo)地位.掠奪關(guān)系標(biāo)準(zhǔn)差橢圓扁率較大,呈“東南-西北”分布.2015～2020年(圖6(d)),互惠共生關(guān)系分布最廣,主要分布于河口區(qū)北部,這源于河口區(qū)進行的濕地修復(fù)以及鹽田養(yǎng)殖產(chǎn)業(yè)的調(diào)整,大規(guī)模城鎮(zhèn)建設(shè)空間轉(zhuǎn)移為水域生態(tài)空間,碳排放量相對減少,新增的水域生態(tài)空間也成為了能力顯著的“碳匯”.互惠共生關(guān)系標(biāo)準(zhǔn)差橢圓扁率為0.44,大致呈“西南-東北”分布.

圖6 2000～2020年黃河三角洲生態(tài)關(guān)系空間分布

表5 生態(tài)關(guān)系空間分布SDE扁率(α)

2.3 生態(tài)關(guān)系變化分析

從“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”國土空間彼此的轉(zhuǎn)移過程看,生態(tài)空間與農(nóng)業(yè)空間在2000～2005年為互惠共生關(guān)系主導(dǎo)(圖7(a)),隨后互惠共生關(guān)系占比持續(xù)下降,競爭關(guān)系逐漸居于主導(dǎo)地位.通過生態(tài)關(guān)系變化率的計算,生態(tài)空間與農(nóng)業(yè)空間的4類生態(tài)關(guān)系除競爭關(guān)系外在三個時段(Ⅰ、Ⅱ、Ⅲ)大多呈現(xiàn)ERCR>0的特點(表6),這說明互惠共生關(guān)系、掠奪關(guān)系和限制關(guān)系的分布規(guī)模變化與兩類空間的轉(zhuǎn)移面積變化基本呈現(xiàn)正向關(guān)系;而競爭關(guān)系在時段Ⅰ、Ⅱ的ERCR值為負(fù),這主要因為農(nóng)業(yè)空間向生態(tài)空間轉(zhuǎn)移的規(guī)模較大,兩空間的掠奪和限制關(guān)系顯著,但對碳儲量的爭奪未達到抗衡狀態(tài),反映為競爭關(guān)系分布減少.

生態(tài)空間與城鎮(zhèn)空間前期以競爭關(guān)系為主,在2010～2015年為掠奪關(guān)系主導(dǎo)(圖7(b)),這是由于二者的碳代謝“角色”特征鮮明,生態(tài)空間向城鎮(zhèn)空間的轉(zhuǎn)移實質(zhì)就是城鎮(zhèn)空間對生態(tài)空間碳儲量的掠奪.2010年后,城鎮(zhèn)空間發(fā)展占用生態(tài)空間的現(xiàn)象有所改善,互惠共生關(guān)系開始出現(xiàn),到2015～2020年,隨著河口區(qū)的濕地修復(fù)使大面積城鎮(zhèn)空間轉(zhuǎn)移為生態(tài)空間,互惠共生關(guān)系也發(fā)展為主導(dǎo)生態(tài)關(guān)系,這也是時段Ⅲ掠奪關(guān)系和限制關(guān)系ERCR<0的主要原因.

農(nóng)業(yè)空間與城鎮(zhèn)空間都具有一定的負(fù)碳流性質(zhì),因此在2000～2005年為競爭關(guān)系為主(圖7(c)); 2005～2010年,兩類空間均向生態(tài)空間有較大面積的轉(zhuǎn)移,使得二者的碳排放量均有所減少,互惠共生關(guān)系居于主要地位,這也使得時段Ⅰ除互惠共生關(guān)系外,其余生態(tài)關(guān)系的ERCR值均為負(fù);2010～2015年農(nóng)業(yè)空間轉(zhuǎn)移為城鎮(zhèn)空間的規(guī)模遠(yuǎn)大于后者向前者的轉(zhuǎn)移,因此出現(xiàn)大面積掠奪關(guān)系,時段Ⅱ掠奪關(guān)系的ERCR值為負(fù),這是由于二者面積轉(zhuǎn)移量下降,但掠奪關(guān)系分布面積大量增加;2015～2020年隨著生態(tài)治理的推進,二者生態(tài)關(guān)系趨于互惠共生,由于在時段Ⅲ這兩類空間的轉(zhuǎn)移面積微弱減少,但互惠共生關(guān)系分布規(guī)模激增,而同時段競爭關(guān)系和限制關(guān)系的變化不大,這使得掠奪關(guān)系相對呈現(xiàn)銳減的態(tài)勢,因此出現(xiàn)了互惠共生關(guān)系ERCR=-10.60094和掠奪關(guān)系ERCR=10.90047的現(xiàn)象.

圖7 2000～2020年“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間生態(tài)關(guān)系變化

表6 “生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間生態(tài)關(guān)系變化率(ERCR)

3 討論

2021年中共中央、國務(wù)院印發(fā)了《黃河流域生態(tài)保護和高質(zhì)量發(fā)展規(guī)劃綱要》,明確指出要“建設(shè)黃河下游綠色生態(tài)走廊,加大黃河三角洲濕地生態(tài)系統(tǒng)保護修復(fù)力度”.黃河三角洲由于其獨特的地理位置和生態(tài)特征,在如今的經(jīng)濟和生態(tài)發(fā)展中愈發(fā)受到重視,本文對研究區(qū)近20a的國土空間變化及其碳代謝特征進行了研究,發(fā)現(xiàn)存在以下幾點問題:(1)生態(tài)空間轉(zhuǎn)移為農(nóng)用現(xiàn)象依然存在,等質(zhì)等量新增生態(tài)空間不到位,導(dǎo)致生態(tài)空間萎縮,增加區(qū)域碳代謝壓力;(2)生態(tài)修復(fù)多集中在北部河口區(qū),東部和南部相對較少,生態(tài)關(guān)系改善區(qū)主要位于西北地區(qū),其他地區(qū)碳代謝壓力猶存,碳代謝平衡區(qū)域較為集中;(3)“單向轉(zhuǎn)移”現(xiàn)象頻發(fā),掠奪和限制關(guān)系仍殘存于部分空間之間,對綠色經(jīng)濟發(fā)展存在一定的阻力影響.基于這些問題提出以下建議:(1)科學(xué)合理地修改和完善國土空間規(guī)劃,堅持對黃河三角洲生態(tài)空間有效利用和持續(xù)修復(fù),盡量避免農(nóng)業(yè)空間與生態(tài)空間進出不平衡現(xiàn)象,調(diào)節(jié)生態(tài)空間與其他空間的協(xié)調(diào)發(fā)展,改善其生態(tài)脆弱問題,提升生態(tài)系統(tǒng)服務(wù)功能;(2)堅持城鎮(zhèn)建設(shè)空間綠色低碳發(fā)展,實行生產(chǎn)生活方式“綠色轉(zhuǎn)型”,調(diào)和鹽田養(yǎng)殖等產(chǎn)業(yè)與濕地生態(tài)系統(tǒng)的矛盾,加強區(qū)域開發(fā)強度與資源環(huán)境承載力之間的耦合[35],以實現(xiàn)“生態(tài)-經(jīng)濟-社會”效益最大化,達到人地關(guān)系協(xié)調(diào)發(fā)展;(3)建立并完善生態(tài)發(fā)展監(jiān)管體系和制度,劃定生態(tài)保護相關(guān)區(qū)域,制定城鎮(zhèn)碳減排政策,對區(qū)域綠色發(fā)展提供保障.

本研究應(yīng)用了生態(tài)網(wǎng)絡(luò)分析(ENA)模型、標(biāo)準(zhǔn)差橢圓及生態(tài)關(guān)系變化率,對黃河三角洲近20a的國土空間碳代謝及其生態(tài)關(guān)系進行了分析,ENA模型很好地反映了不同類型空間之間直接和間接的碳代謝關(guān)系,使得分析結(jié)果更為科學(xué)合理;同時ENA模型將不同國土空間類型的相互作用關(guān)系具體化為4類,直觀反映了國土空間變化帶來的碳代謝效應(yīng),結(jié)合ArcGIS從時空上清晰地獲取研究區(qū)碳代謝分布信息,最后再輔以標(biāo)準(zhǔn)差橢圓和生態(tài)關(guān)系變化率,對生態(tài)關(guān)系分布從空間和數(shù)量上進一步細(xì)化,對其動態(tài)變化進行分析,明確了黃河三角洲碳代謝生態(tài)關(guān)系與國土空間利用的內(nèi)在聯(lián)系及發(fā)展趨勢.未來可對黃河三角洲的國土空間利用進行模擬,預(yù)測多情景下的碳代謝分布,利用ENA模型對模擬結(jié)果進行分析,以生態(tài)沖突程度為依據(jù)對黃河三角洲的國土空間規(guī)劃提供科學(xué)的參考資料.

4 結(jié)論

4.1 2000～2015年,黃河三角洲整體為負(fù)碳流主導(dǎo),研究區(qū)負(fù)碳流量遠(yuǎn)高于同期的正碳流量,其中2005～2010年負(fù)碳流量達到峰值,整體上該15a研究區(qū)呈現(xiàn)碳代謝紊亂的局面;2015～2020年,正碳流量顯著提高,正負(fù)碳流比達到了1:0.18,且分布大量集中在研究區(qū)北部位置,這表明研究區(qū)在后期有逐漸趨向于碳代謝平衡的態(tài)勢.

4.2 數(shù)量結(jié)構(gòu)上,黃河三角洲20a間占據(jù)主導(dǎo)地位的生態(tài)關(guān)系是競爭關(guān)系,占生態(tài)關(guān)系總數(shù)量的29.77%,其中在2000～2005年占比最大,達38.10%,這表明研究期內(nèi)該地區(qū)不同空間類型之間的生態(tài)沖突較為強烈;空間上,2000～2010年掠奪和限制關(guān)系分布于研究區(qū)東部和北部,該區(qū)域存在較為劇烈的生態(tài)沖突,2015年以后因生態(tài)修復(fù)、產(chǎn)業(yè)結(jié)構(gòu)調(diào)整,研究區(qū)生態(tài)關(guān)系顯著改善,互惠共生關(guān)系占據(jù)主導(dǎo)地位.總體來看,國土空間演變重點區(qū)域由南向北轉(zhuǎn)變.

4.3 對于“生態(tài)-農(nóng)業(yè)-城鎮(zhèn)”空間兩兩相互之間的轉(zhuǎn)移,生態(tài)空間與農(nóng)業(yè)空間的主導(dǎo)碳代謝生態(tài)關(guān)系由互惠共生關(guān)系向競爭關(guān)系轉(zhuǎn)變,而生態(tài)空間與農(nóng)業(yè)空間的主導(dǎo)生態(tài)關(guān)系由競爭關(guān)系逐漸變?yōu)榛セ莨采P(guān)系,農(nóng)業(yè)空間與城鎮(zhèn)空間的主導(dǎo)生態(tài)關(guān)系無漸變趨勢,但以互惠共生關(guān)系為主;整體上大部分生態(tài)關(guān)系變化率滿足ERCR>0,即生態(tài)關(guān)系分布規(guī)模的增減與空間轉(zhuǎn)移面積的變化趨勢一致.

[1] 武漢大學(xué)國家發(fā)展戰(zhàn)略研究院課題組.中國實施綠色低碳轉(zhuǎn)型和實現(xiàn)碳中和目標(biāo)的路徑選擇[J]. 中國軟科學(xué), 2022,(10):1-12. Research Group of National Institute of Development Strategy at Wuhan University. Path choice for China to implement green low- carbon transformation and achieve carbon neutrality [J]. China Soft Science, 2022,(10):1-12.

[2] Deutsch C A, Tewksbury J T, Tigchelaar M, et al. Increase in crop losses to insect pests in a warming climate [J]. Science, 2018,361(6405):916-919.

[3] 張正峰,張 棟.基于社會網(wǎng)絡(luò)分析的京津冀地區(qū)碳排放空間關(guān)聯(lián)與碳平衡分區(qū)[J]. 中國環(huán)境科學(xué), 2023,43(4):2057-2068. Zhang Z F, Zhang D. Spatial relatedness of CO2emission and carbon balance zoning in Beijing Tianjin Hebei counties [J]. China Environmental Science, 2023,43(04):2057-2068.

[4] 杜金霜,付晶瑩,郝蒙蒙.基于生態(tài)網(wǎng)絡(luò)效用的昭通市“三生空間”碳代謝分析[J]. 自然資源學(xué)報, 2021,36(5):1208-1223. Du J S, Fu J Y, Hao M M. Analyzing the carbon metabolism of “Production-LivingEcological” space based on ecological network utility in Zhaotong [J]. Journal of Natural Resources, 2021,36(5): 1208-1223.

[5] Sovacool B K, Brown M A. Twelve metropolitan carbon footprints: A preliminary comparative global assessment [J]. Energy Policy, 2009, 38(9):4856-4869.

[6] Karakiewicz J. Urban metabolism of low carbon cities [C]//The 47th ISOCARP Congress. International Society of City and Regional Planners (ISOCARP) and Urban Planning Society of China (UPSC). Wuhan, China: UPSC, 2011.

[7] Patten B C. Energy, emergy and environs [J]. Ecological Modelling, 1992,62(1-3):29-69.

[8] Fath B D. Network mutualism: Positive community-level relations in ecosystems [J]. Ecological Modelling, 2007,208(1):56-67.

[9] Zhang Y, Xia L L, Fath B D, et al. Development of a spatially explicit network model of urban metabolism and analysis of the distribution of ecological relationships: case study of Beijing, China [J]. Journal of Cleaner Production, 2016,112:4304-4317.

[10] Xia C Y, Chen B. Urban land-carbon nexus based on ecological network analysis [J]. Applied Energy, 2020,276.

[11] 夏琳琳,張 妍,李名鏡.城市碳代謝過程研究進展 [J]. 生態(tài)學(xué)報, 2017,37(12):4268-4277. Xia L L, Zhang Y, Li M J. A review of urban carbon metabolism [J]. Acta Ecologica Sinica, 2017,37(12):4268-4277.

[12] 歐陽少虎,黃盼盼,陶宗鑫,等.石油污染土壤對蚯蚓氧化應(yīng)激響應(yīng)研究[J]. 中國環(huán)境科學(xué), 2023,43(8):4353-4359. Ouyang S H, Huang P P, Tao Z X, et al. Oxidative stress responseresponse of petroleum-contaminated soil to earthworm [J]. China Environmental Science, 2023,43(8):4353-4359.

[13] 曲衍波,王世磊,朱偉亞,等.黃河三角洲國土空間演變的時空分異特征與驅(qū)動力分析[J]. 農(nóng)業(yè)工程學(xué)報, 2021,37(6):252-263,309. Qu Y B, Wang S L, Zhu W Y, et al. Spatial-temporal differentiation characteristics and driving force of territorial space evolution in the Yellow River Delta [J]. Transactions of the Chinese Society of Agricultural Engineering, 2021,37(6):252-263,309.

[14] 徐新良,龐治國,于信芳.土地利用/覆被變化時空信息分析方法及應(yīng)用 [M]. 北京:科學(xué)技術(shù)文獻出版社, 2014:90-108. Xu X L, Pang Z G, Yu X F. Spatial-temporal pattern analysis of land use/cover change: methods & applications [M]. Beijing: Science and Technical Documentation Press, 2014:90-108.

[15] 歐定華,張 琪,秦 景,等.基于土地利用與其功能動態(tài)耦合性的縣域國土空間分類體系構(gòu)建[J]. 農(nóng)業(yè)工程學(xué)報, 2021,37(24):284-296. Ou D H, Zhang Q, Qin J, et al. Classification system for county-level territorial space using spatiotemporal heterogeneity and dynamic coupling of land use and functionality [J]. Transactions of the Chinese Society of Agricultural Engineering, 2021,37(24):284-296.

[16] 范建雙,周 琳.中國城鄉(xiāng)居民生活消費碳排放變化的比較研究[J]. 中國環(huán)境科學(xué), 2018,38(11):4369-4383. Fan J S, Zhou L. A comparative study on the changes of residential living consumption carbon emissions in urban and rural China [J]. China Environmental Science, 2018,38(11):4369-4383.

[17] 艾曉艷,趙 源,王大川.基于碳效應(yīng)的滎經(jīng)縣土地系統(tǒng)脆弱性分析[J]. 中國農(nóng)學(xué)通報, 2015,31(29):113-122. Ai X Y, Zhao Y, Wang D C. Land system vulnerability analysis based on carbon effect in Yingjing [J]. Chinese Agricultural Science Bulletin, 2015,31(29):113-122.

[18] 夏楚瑜,李 艷,葉艷妹,等.基于生態(tài)網(wǎng)絡(luò)效用的城市碳代謝空間分析——以杭州為例[J]. 生態(tài)學(xué)報, 2018,38(1):73-85. Xia C Y, Li Y, Ye Y M, et al. Analyzing urban carbon metabolism based on ecological network utility: A case study of Hangzhou City [J]. Acta Ecologica Sinica, 2018,38(1):73-85.

[19] 舒 心,夏楚瑜,李 艷,等.長三角城市群碳排放與城市用地增長及形態(tài)的關(guān)系[J]. 生態(tài)學(xué)報, 2018,38(17):6302-6313. Shu X, Xia C Y, Li Y, et al. Relationships between carbon emission, urban growth, and urban forms of urban agglomeration in the Yangtze River Delta [J]. Acta Ecologica Sinica, 2018,38(17):6302-6313.

[20] 田 云,尹忞昊.中國農(nóng)業(yè)碳排放再測算:基本現(xiàn)狀、動態(tài)演進及空間溢出效應(yīng)[J]. 中國農(nóng)村經(jīng)濟, 2022,(3):104-127. Tian Y, Yin M H. Re-evaluation of China’ s agricultural carbon emissions: Basic status, dynamic evolution and spatial spillover effects [J]. Chinese Rural Economy, 2022,(3):104-127.

[21] 李文翔,李 曄,蔡近近.政府-企業(yè)-居民協(xié)同共治的道路交通碳交易機制[J]. 中國環(huán)境科學(xué), 2021,41(9):4426-4438. Li W X, Li Y, Cai J J. The CO2emission trading system for road transport collaboratively governed by the government, enterprises, and residents [J]. China Environmental Science, 2021,41(9):4426-4438.

[22] 夏楚瑜.基于土地利用視角的多尺度城市碳代謝及“減排”情景模擬研究[D]. 杭州:浙江大學(xué), 2019. Xia C Y. Multi-scale studies on urban carbon metabolism from the perspective of land use and Scenario analysis of emission reduction [D]. Hangzhou: Zhejiang University, 2019.

[23] 鄒 康,舒予晴,李桂娥,等.基于生態(tài)網(wǎng)絡(luò)分析(ENA)的城市地-碳框架構(gòu)建及其時空演變研究[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2022,38(8): 972-982. Zou K, Shu Y Q, Li G E, et al. Urban land-carbon framework construction based on ecological network analysis and its space-time evolution research [J]. Journal of Ecology and Rural Environment, 2022,38(8):972-982.

[24] 胡婉玲,張金鑫,王紅玲.中國農(nóng)業(yè)碳排放特征及影響因素研究[J]. 統(tǒng)計與決策, 2020,36(5):56-62. Hu W L, Zhang J X, Wang H L. Characteristics and influencing factors of agricultural carbon emission in China [J]. Statistics & Decision, 2020,36(5):56-62.

[25] 曾賢剛,余 暢.中國農(nóng)業(yè)農(nóng)村碳排放結(jié)構(gòu)與碳達峰分析[J]. 中國環(huán)境科學(xué), 2023,43(4):1906-1918. Zeng X G, Yu C. Studies on carbon emission structure and carbon peak of agriculture and rural areas in China [J]. China Environmental Science, 2023,43(4):1906-1918.

[26] 李 晶,陳松林,李晨欣,等.基于生態(tài)網(wǎng)絡(luò)分析的廈漳泉地區(qū)土地碳代謝[J]. 應(yīng)用生態(tài)學(xué)報, 2023,34(5):1375-1383. Li J, Chen S L, Li C X, et al. Land carbon metabolism in Xiamen- Zhangzhou-Quanzhou region based on ecological network analysis [J]. Chinese Journal of Applied Ecology, 2023,34(5):1375-1383.

[27] 程 蕾,郝增超,曾維華.基于生態(tài)網(wǎng)絡(luò)分析法的煤化工水網(wǎng)絡(luò)結(jié)構(gòu)解析[J]. 中國環(huán)境科學(xué), 2017,37(6):2363-2368. Cheng L, Hao Z C, Zeng W H. Analysis of the water network structure of coal chemical industry based on ecological network analysis [J]. China Environmental Science, 2017,37(6):2363-2368.

[28] 王悅斌,林常青,李 媛,等.基于多流域生態(tài)網(wǎng)絡(luò)模型的用水系統(tǒng)結(jié)構(gòu)特征演化分析[J]. 中國環(huán)境科學(xué), 2018,38(2):755-765.Wang Y B, Lin C Q, Li Y, et al. Analysis of the network structure evolution of water use systems based on multi-basin ecological network model [J]. China Environmental Science, 2018,38(2):755-765.

[29] Coskun H. Static ecological system analysis [J]. Theoretical Ecology, 2020,13(1):17-52.

[30] 王 菲,曹永強,范帥邦.“雙碳”目標(biāo)下東北三省水-能源紐帶關(guān)系及網(wǎng)絡(luò)特征分析[J]. 生態(tài)學(xué)報, 2022,42(14):5692-5707.Wang F, Cao Y Q, Fan S B. Analysis of water-energy ties and network characteristics in the three provinces of northeast china under the “double carbon” target [J]. Acta Ecologica Sinica, 2022,42(14):5692-5707.

[31] Meng B, Liu J L, Bao K, et al. Water fluxes of Nenjiang River Basin with ecological network analysis: Conflict and coordination between agricultural development and wetland restoration [J]. Journal of Cleaner Production, 2019,213:933-943.

[32] 李云燕,張 碩.中國城市碳排放強度時空演變與影響因素的時空異質(zhì)性[J]. 中國環(huán)境科學(xué), 2023,43(6):3244-3254. Li Y Y, Zhang S. Spatio-temporal evolution of urban carbon emission intensity and spatiotemporal heterogeneity of influencing factors in China [J]. China Environmental Science, 2023,43(6):3244-3254.

[33] 徐 旳,徐 彥.1928～2017年南京城市戶籍人口時空格局演變[J]. 地理學(xué)報, 2022,77(10):2439-2456. Xu D, Xu Y. Spatio-temporal pattern of registered population in Nanjing from 1928 to 2017 [J]. Acta Geographica Sinica, 2022,77(10): 2439-2456.

[34] 楊仕琪,王 冀,竇銀銀,等.1916～2020年北京城市變遷及其與區(qū)域氣候演化的關(guān)系[J]. 地理學(xué)報, 2023,78(3):620-639. Yang S Q, Wang J, Dou Y Y, et al. Spatio-temporal evolution of urbanization and its relationship with regional climate change in Beijing over the past century [J]. Acta Geographica Sinica, 2023,78(3): 620-639.

[35] 盛科榮,樊 杰.地域功能的生成機理:基于人地關(guān)系地域系統(tǒng)理論的解析[J]. 經(jīng)濟地理, 2018,38(5):11-19. Sheng K R, Fan J. The formation mechanism of regional function: An analysis based on the theory of man-earth areal system [J]. Economic Geography, 2018,38(5):11-19.

ENA-model based carbon metabolism analysis of the Yellow River Delta.

LU Chang*, WANG Zhi-yu

(School of Management Engineering, Shandong Jianzhu University, Jinan 250101, China)., 2023,43(11)：6090～9101

The purpose of the study is to explore in depth the intrinsic linkage between carbon metabolism ecological relationships and national spatial utilization in the Yellow River Delta, and to provide a reference for national spatial pattern adjustment and green development. The research was based on the spatial classification system of “ecology-agriculture-urban”, using the ENA-model to analyze the carbon metabolism process in the Yellow River Delta from 2000 to 2020, and using the standard deviation ellipse and the ecological relationship change rate(ERCR) to evaluate the spatial and temporal distribution and change characteristics of ecological relationships between spaces. The results of the analysis: In the Yellow River Delta, the carbon metabolism was in disorder from 2000 to 2015 due to frequent phenomena such as the occupation of ecological space by agricultural space, and the positive carbon flow increased significantly after 2015, resulting in a balanced carbon metabolism; The dominant ecological relationships in the study area were competitive, with stronger ecological conflicts over carbon metabolism, and mutualistic symbiotic relationships gradually increased in the northwest after 2015, with the priority area of national land space evolution shifted from south to north; Most of the ecological relationship change rate in the study area satisfy ERCR>0, and the scale of ecological relationship increase or decrease was consistent with the change trend of spatial transfer area. The conclusion of the study: There are still unreasonable utilization of the national space in the Yellow River Delta, and the pressure of carbon metabolism still exists, so it is necessary to further optimize the national spatial pattern and promote green development.

“ecology-agriculture-urban” space；carbon metabolism；ENA-model；the ecological relationships change rate；the Yellow River Delta

F301.2;X321

A

1000-6923(2023)11-6090-12

路 昌(1990-),男,山東臨沂人,講師,博士,主要從事土地利用轉(zhuǎn)型及生態(tài)環(huán)境效應(yīng)研究.發(fā)表論文20余篇.chang20081990@126. com.

路 昌,王之語.基于生態(tài)網(wǎng)絡(luò)分析的黃河三角洲碳代謝分析 [J]. 中國環(huán)境科學(xué), 2023,43(11):6090-6101.

Lu C, Wang Z Y. ENA-model based carbon metabolism analysis of the Yellow River Delta [J]. China Environmental Science, 2023,43(11):6090-6101.

2023-04-03

國家自然科學(xué)基金青年基金資助項目(42301300);廣東省哲學(xué)社會科學(xué)規(guī)劃項目(GD21YGL05);山東建筑大學(xué)博士科研基金資助項目(XNBS1826)

* 責(zé)任作者, 講師, chang20081990@126.com