中國縣域碳匯時空格局及影響因素

張 赫,彭千芮,王 睿,*,強文麗,張建勛

1 天津大學建筑學院, 天津 300072 2 蘭州大學資源環(huán)境學院, 蘭州 730000

土地利用/覆被變化深刻作用于陸地生態(tài)系統(tǒng)的碳循環(huán),對全國和區(qū)域氣候變化產(chǎn)生了舉足輕重的影響[1]。1980年到2005年間,由于土地利用方式的轉(zhuǎn)變,我國人為源排放量由陸地生態(tài)系統(tǒng)存儲的3倍增長到10倍,生態(tài)環(huán)境惡化明顯[2]。同時,作為全球溫室氣體排放的大國,我國2019年全球溫室氣體排放量高居首位。事實上,我國承諾到2020年碳排放強度比2005年下降40%—45%[3]?？h域作為城鄉(xiāng)土地增減掛鉤的直接單元[4],其非建設(shè)用地的土地利用方式對大氣溫室氣體含量的調(diào)節(jié)作用不容忽視,利用其碳貯存量增加而產(chǎn)生減排效果,對推動我國城鄉(xiāng)控碳減排有極大意義[5]。在此背景下,加強對我國縣域非建設(shè)用地碳匯規(guī)模與規(guī)律研究,對構(gòu)建基于低碳導向的國土空間組織體系,優(yōu)化國土空間開發(fā)格局,實施差異化減排,增強生態(tài)系統(tǒng)固碳能力具有重要的理論和現(xiàn)實意義。

隨著碳中和(carbon neutral)概念的提出,“碳匯”指個人或組織等以碳中和方式向從事“碳中和”業(yè)務的中介機構(gòu)支付基金,最后以植樹等形式用于碳減排的過程和機制[6-7]。通過改變非建設(shè)用地的土地利用/覆被方式,對其進行合理的開發(fā)利用,可以降低土地利用碳排放占人為碳排放的比例,將此時的土地利用變化發(fā)揮出碳匯作用[8]。迄今為止,國內(nèi)外學者針對土地利用碳源/匯相關(guān)研究主要集中于以下三個方面：首先,通過構(gòu)建計算模型對碳匯/源進行總量計算和效應評估來探析不同區(qū)域的碳匯時空特征。Naipal等[9]通過建立土壤侵蝕模型,研究了在一定時期內(nèi)土地利用變化下萊茵河流域的碳循環(huán)特征；彭文甫等[10]構(gòu)建了碳排放、碳足跡等模型計算并分析四川省土地利用碳匯的時空變化特征。其次,研究土地利用方式轉(zhuǎn)變對碳匯的影響。裴杰等[11]利用地理信息技術(shù)分析了深圳市土地利用/覆被變化情況及其導致的碳效應變化；張海鳳和崔桂善[12]通過總結(jié)不同中碳收支的特點對土地覆被類型轉(zhuǎn)化下的碳源/匯變化的影響進行了分析,羅谷松和李濤[13]運用Super SBM-undesirable DEA等方法分析了碳排放影響下中國省域土地利用效率的影響因素。第三,探究影響土地利用碳源/匯的相關(guān)因素以及對應的控碳減排策略。李小康等[14]通過分析土地利用結(jié)構(gòu)及碳排放的變化發(fā)現(xiàn)控制建設(shè)用地總量對控碳減排意義重大；Chen等[15]結(jié)合渦度協(xié)方差法分析發(fā)現(xiàn)氣候帶類型和植樹造林是影響碳封存的主要因素。

縱觀現(xiàn)有研究成果,研究尺度主要集中于國家、區(qū)域和省域?qū)用?對縣域的研究相對缺乏,已有研究中僅趙榮欽等[16- 17]對中原經(jīng)濟區(qū)和河南省的縣域碳收支狀況進行了核算和分析?？h域作為中國經(jīng)濟社會發(fā)展的基本單元,開展碳匯在其尺度下的時空變化特征和影響因素的分析,對評估縣級尺度的生態(tài)狀況以及推動區(qū)域低碳協(xié)調(diào)發(fā)展意義重大?；诖?針對全國縣域,基于全國土地利用數(shù)據(jù)(1990—2015),運用標準差橢圓分析方法、空間自相關(guān)分析方法和地理加權(quán)回歸分析方法,著重研究：①中國縣域碳匯隨時間變化的特征；②中國縣域碳匯的空間分異格局；③中國縣域碳匯的時空變化規(guī)律和影響因素。

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

1.1 研究區(qū)域

按照中國現(xiàn)行的行政區(qū)劃,縣級行政區(qū)包括地級市的市轄區(qū)、縣級市、縣、自治縣、旗、自治旗、林區(qū)、特區(qū)等8種。由于地級市的市轄區(qū)、縣級市以及部分特區(qū)的土地利用以城鎮(zhèn)建設(shè)用地為主,其功能主要體現(xiàn)在碳排放上而固碳能力較低,所以不在討論的范圍內(nèi)。綜上,通過篩選最后選定包含縣、自治縣、旗、自治旗、林區(qū)在內(nèi)的1330個縣級行政單位。

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

土地利用數(shù)據(jù)來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心(http://www.resdc.cn),是以1990—2015年6期Landsat TM/ETM遙感影像為主要數(shù)據(jù)源,通過人工目視解譯獲取。通過ArcGIS 10.2,用2018年中國縣域矢量行政邊界裁剪,經(jīng)過重采樣為1 km的全國6期土地利用數(shù)據(jù)；應用ENVI 4.8對土地利用數(shù)據(jù)進行編碼和歸類,將土地利用類型劃分為林地、草地、水域、未利用地4大類、13小類。耕地的農(nóng)作物通過光合作用吸收空氣中的CO2,但絕大多數(shù)農(nóng)作物在短期內(nèi)又通過呼吸作用釋放到空氣中去,農(nóng)作物生物量作為碳匯的效果不明顯[18],因此不在本文進行討論。

全國縣域國民生產(chǎn)總值、產(chǎn)業(yè)結(jié)構(gòu)等數(shù)據(jù)來源于2015年《中國城市統(tǒng)計年鑒》、《區(qū)域經(jīng)濟統(tǒng)計年鑒》、相關(guān)地市級統(tǒng)計年鑒及環(huán)境統(tǒng)計公報等。

2 計算與研究方法

2.1 計算方法

土地利用的碳匯估算主要涉及耕地、林地、草地、水域和未利用地的碳匯。耕地作為碳源[18],在本文中不進行研究；林地、草地、水域和未利用地為碳匯,為主要研究對象。因此,碳匯估算模型可表達如下:

(1)

式中,Cs為碳匯量,si為第i種土地利用類型產(chǎn)生的碳匯量；Ai為第i種土地利用類型(林地、草地、水域、未利用地)面積,αi為第i種土地利用類型的碳匯系數(shù)。

2.1.1碳匯系數(shù)的確定

碳匯系數(shù)根據(jù)不同用地類型,通過跟蹤國內(nèi)外期刊的研究成果進行確定(表1)。將林地按照郁閉度大小分為有林地、灌木林、疏林地和其他林地4種類型。將草地按照覆蓋度大小分為高覆蓋度草地、中覆蓋度草地和低覆蓋度草地3種類型。將我國水域分為河渠、湖泊、水庫坑塘、灘涂和灘地5種類型。未利用地指我國還未利用的土地或難以利用的土地類型,包括沙地、戈壁、鹽堿地、沼澤地、裸土地、裸巖石質(zhì)地和其他未利用地等。我國未利用地有的是碳匯、有的是碳排放,但都較弱,根據(jù)賴力等[23]的研究,取未利用地碳匯系數(shù)為0.0005 t hm-2a-1。

表1 土地利用類型及對應的碳匯系數(shù)

2.1.2碳匯強度

考慮到用地面積大小對碳匯總量的影響,本文構(gòu)建碳匯強度的概念來對單位土地碳匯能力進行衡量,并且就碳匯強度進行局部自相關(guān)分析。碳匯強度計算公式如下：

Cp=Cs/Su

(2)

式中,Cp為碳匯強度；Cs為碳匯總量(×104t),Su為非建設(shè)用地面積。

2.2 研究方法

2.2.1標準差橢圓分析法

基于ArcGIS 10.2的空間分析功能,采用標準差橢圓分析方法[24]對全國整體的碳匯在空間位置的分布情況進行分析。橢圓的長軸表征數(shù)據(jù)集分布的主要方向,其面積表征數(shù)據(jù)集分布的集聚(分散)程度。

2.2.2空間相關(guān)分析法

基于ArcGIS 10.2的局部空間自相關(guān)分析功能,使用空間聯(lián)系局域指標(Local Indicators of Spatial Association, LISA)對區(qū)域整體相關(guān)程度進行衡量[25- 26]：當Moran′sI為正數(shù)時,表示碳匯量在空間上呈集聚分布；當Moran′sI為負數(shù)時,表示碳匯量在空間上呈分散分布。冷點與熱點的識別能進一步討論的空間分布是聚集、分散或是隨機的,目的是識別全國縣域碳匯熱點與冷點的空間分布。在ArcGIS中,采用 Getis-Ord G*統(tǒng)計量測度全國縣域碳匯在局部空間的依賴性及異質(zhì)性,Gi*(d) 值為正表示高值聚類,即“熱點區(qū)”；反之為低值聚類,即“冷點區(qū)”[27- 28]。

2.2.3地理加權(quán)回歸分析法

地理加權(quán)回歸分析(Geographically Weighted Regression,GWR),可以量化反映空間數(shù)據(jù)關(guān)系中的異質(zhì)性特征,與傳統(tǒng)回歸模型相比,可以對不同空間進行局部地加權(quán)回歸模型,并求解[29]?？h域碳匯的空間異質(zhì)性是應用地理加權(quán)回歸模型進行分析的前提條件,采用AIC法與核密度估計求得帶寬,利用高斯函數(shù)來確定權(quán)重,選取Fotheringham等提出的方法作為計算準則,即當AIC值取最小值時,則此時的帶寬為最佳帶寬[30]。以上所有操作均通過ArcGIS 10.2軟件完成。

3 中國縣域碳匯時空分布特征

3.1 1990—2015中國縣域碳匯的時序變化

3.1.1時序變化總體特征

1990—2015年中國碳匯總量呈波動下降趨勢,碳匯量從13307.79×104t下降至13198.27×104t(圖1)。可以將凈碳匯總量的變化劃分為3個階段,即大幅下降-上升-緩慢下降。具體表現(xiàn)為,1990—1995年間,碳匯呈現(xiàn)急劇下降的趨勢,從1990年的13307.79×104t,至1995年達到13003.41×104t的谷值；之后近10年內(nèi)回升,到2005年,碳匯達到13231.46×104t；隨后緩慢下落,到2015年達到13198.27×104t。

圖1 1990—2015年中國碳匯總量及比重Fig.1 Total carbon sink volume and proportion in China, 1990—2015

3.1.2時序變化比例結(jié)構(gòu)

碳匯類型主要以林地碳匯為主,占比約80%,依次是草地(約13.5%)、水域(約4.0%)和未利用地(約2.5%)(圖1)。進一步對不同碳匯類型的比例結(jié)構(gòu)展開分析,雖然比例結(jié)構(gòu)大體沒變,但在演變趨勢上略有不同,高覆蓋度草地、灘地碳匯都在1995年達到最低值,并在之后的年份里回升；灘涂、疏木林、沼澤地的碳匯均呈現(xiàn)先上升再波動下降的趨勢；而其他類型的碳匯則呈現(xiàn)波動維穩(wěn)的狀態(tài)。

3.2 1990—2015中國縣域碳匯的空間格局

3.2.1空間總體格局

圖2 1990—2015年中國縣域碳匯空間分布Fig.2 Spatial distribution of county carbon sinks in China, 1990—2015

從全國縣域碳匯空間分布來看(圖2),1990—2015年間,中國縣域碳呈現(xiàn)差異化分布的空間格局,在區(qū)域上[31]呈現(xiàn)不同的變化趨勢。從碳匯總量來看,全國縣域碳匯大致呈現(xiàn)“西部>東北>南部>中部”的 “西高東低”格局；碳匯低值區(qū)域主要集中在新疆西北部、內(nèi)蒙古西南部、甘肅、陜西、寧夏、山西、河北、河南、山東、安徽以及江蘇北部；而碳匯高值區(qū)域主要集中于內(nèi)蒙古東北部、西藏、浙江、江蘇、福建地區(qū)；四川、云南、貴州廣西等地碳匯高值與低值區(qū)域相間分布。分時段看,1990—1995年間,碳匯空間格局有所變化,內(nèi)蒙古中部、西藏北部碳匯均有所下降；1995—2015年間,內(nèi)蒙古中部、西藏北部碳匯量回升,四川西北部碳匯增加。整體上看,我國內(nèi)蒙古中部、黑龍江、吉林、遼寧、山東、河南、安徽等地碳匯總量下降較多,而新疆南部、陜西、陜西以及我國華中、華南、西南部分地區(qū)的碳匯則呈上升趨勢。

3.2.2空間分異特征

圖3 1990—2015年中國縣域碳匯標準差橢圓 Fig.3 Ellipse standard deviation of county carbon sinks in China, 1990—2015

從標準差橢圓方法揭示的空間分異格局來看(圖3),碳匯整體以東北-西南的方向為主導,空間分布中心存在向西南移動的趨勢,空間分布范圍呈現(xiàn)收縮態(tài)勢。在分異形狀上,1991—2011 年,全國凈碳匯空間分布標準差橢圓短軸與長軸的比值總體呈現(xiàn)出兩階段的規(guī)律,1990—2010 年間,全國凈碳匯空間分布標準差橢圓短軸與長軸的比重呈現(xiàn)出增加的趨勢,隨后2010—2015 年間基本上保持不變,下降趨勢極慢。

在分異方向上,1990—2015年間,全國碳匯空間分布標準差橢圓方位角在波動中呈現(xiàn)總體增長的趨勢(表2)。1990—1995年,方位角不斷縮小,西南部地區(qū)對全國碳匯空間格局的影響作用加強；1995—2000年,方位角不斷增大,東北地區(qū)對全國凈碳匯空間格局的影響作用加強； 2000—2005年,方位角不斷縮小,隨后保持不變,西南部地區(qū)的碳匯量對全國碳匯空間格局的影響作用加強后,趨于平穩(wěn)。

表2 1990—2015年中國縣域碳匯標準差橢圓測算結(jié)果

3.2.3空間聚集特征

從冷熱點的空間分布變化看(圖4),中國縣域碳匯總量的冷熱點空間分布集聚程度呈現(xiàn)波動穩(wěn)定的特征。整體上看,碳匯總量熱點以黑龍江、內(nèi)蒙古北部,西藏南部和廣東、湖南交界處為中心,在其周圍呈連片分布；碳匯總量冷點在山西、陜西、河北、河南、安徽、江蘇等地連片分布,并且在四川、重慶和貴州交界處有少量分布。在時空上,1990—1995年間,我國西部的碳匯總量熱點從西藏藏族自治區(qū)的尼瑪縣、雙湖縣和改則縣消失,演變?yōu)閲@西藏南部的加查縣、那曲縣、班戈縣等地分布的格局；1995—2000年,我國西部的碳匯總量的熱點逐漸變化為1990年的格局特征,并在2000—2015年間保持較為平穩(wěn)的態(tài)勢。

圖4 1990—2015年中國縣域碳匯的熱點分析Fig.4 The of hot spot analysis of county carbon sinks in China, 1990—2015

1990—2015年,我國東南部碳匯總量熱點呈現(xiàn)先擴大再逐漸縮小的過程；1990—1995年,我國東南部碳匯熱點區(qū)域呈現(xiàn)向湖南桂陽縣、嘉禾縣、臨武縣等地聚集特征,熱點分布范圍有擴大趨勢；1995—2015,我國西南部碳匯總量熱點分布范圍逐年縮小。1990—2000年間,我國東北部碳匯總量熱點區(qū)域不斷縮小,體現(xiàn)在內(nèi)蒙古阿魯科爾、巴林右旗、克什克騰、扎魯特旗等地逐漸退出碳匯總量熱點區(qū)；2000—2015年,我國東北部碳匯總量熱點區(qū)域逐漸擴大。從我國縣域碳匯總量冷熱點的時空分布可以看出,我國碳匯總量整體格局是西南、東南、東北多,中部少,空間集聚性和整體分布波動維穩(wěn)。

圖5 1990—2015年中國縣域碳匯強度的LISA集聚圖Fig.5 LISA cluster map of per capita county carbon sinks in China, 1990—2015

從局部空間自相關(guān)來看(圖5),中國縣域碳匯強度的空間集聚程度變化不明顯,具有高值與低值聚集,高-低區(qū)域零星分布的特點。其中,碳匯強度高值區(qū)域集聚在湖南、江西、福建、浙江及云南、四川、黑龍江、吉林、遼寧的部分地區(qū)；碳匯強度低值區(qū)域集聚在新疆西北部、內(nèi)蒙古西部,青海東部,以及寧夏、甘肅、山西、河北、山東、河南、安徽等地；高-低區(qū)域較稀少,零星分布于陜西、山西、河南及京津冀地區(qū)。

對比中國縣域碳匯總量的冷熱點和中國縣域碳匯強度的局部空間自相關(guān)時空分布,整體上看,我國東北部鄂倫春自治旗、漠河縣、塔城縣等地,東南部桂陽縣、嘉禾縣等地的碳匯總量與碳匯強度都出現(xiàn)了高值聚集的特征；而西藏加查縣、那曲縣等地的碳匯總量雖然呈現(xiàn)熱點分布,但是碳匯強度卻在低值聚類區(qū)域；而我國中部和東部沿海地區(qū)的陜西、江蘇、河北等地的縣域碳匯總量和碳匯強度等處于冷點和低值聚類區(qū)域。

4 影響因素分析

在對碳匯影響因素的研究中,大多數(shù)學者從氣候條件、植被類型角度對碳匯能力進行分析。但本文在變量選取時主要考慮了以下兩個方面: 首先,本文在非建設(shè)用地的土地利用視角下對縣域碳匯的空間分異進行研究,同時想要考察城鎮(zhèn)化過程以及產(chǎn)業(yè)結(jié)構(gòu)調(diào)整對我國的碳匯是否真正起到了促進作用; 其次,地理加權(quán)回歸方法對各個變量是否存在多重共線性要求較高,選取多個變量容易產(chǎn)生多重共線性,基于以上原因,本文重點選取了經(jīng)濟發(fā)展、產(chǎn)業(yè)結(jié)構(gòu)和土地利用程度作為研究的解釋變量。

各解釋變量的定義和選取理由：①經(jīng)濟發(fā)展：不同地區(qū)處于不同經(jīng)濟發(fā)展階段,其經(jīng)濟發(fā)展水平也不同,經(jīng)濟發(fā)展帶來大量的城市開發(fā),碳匯也隨之減少。用國民生產(chǎn)總值(GDP)來表征地區(qū)經(jīng)濟發(fā)展。②產(chǎn)業(yè)結(jié)構(gòu)：工業(yè)高排放與高污染的特征對環(huán)境的脅迫效應明顯,要促進我國生態(tài)文明建設(shè),產(chǎn)業(yè)結(jié)構(gòu)升級和增匯同樣重要。用工業(yè)產(chǎn)值占GDP總量的比例表征產(chǎn)業(yè)結(jié)構(gòu)。③土地利用程度：土地利用程度越低,碳匯通常更高,隨著城市建設(shè)用地的擴張,碳匯量也會減少。用縣域建設(shè)用地面積和縣域總面積的比值表征土地利用程度。

由于各個變量的量綱不同,因此對各個數(shù)據(jù)進行最大-最小值標準化處理,通過ArcGIS10.2對我國1330個縣域數(shù)據(jù)進行地理加權(quán)回歸,回歸結(jié)果為回歸的整體結(jié)果輸出,包括回歸的擬合系數(shù)、最優(yōu)帶寬和AIC值(表3)。

表3 GWR模型的整體估計結(jié)果

通過地理加權(quán)回歸的方法,可以得到每個縣的回歸方程,而這些結(jié)果可以考慮到各省市在空間地理位置的相互影響,因此結(jié)果可以充分考察各地區(qū)的空間異質(zhì)性。從各因素的影響程度來看(圖6),2015年經(jīng)濟發(fā)展和土地利用程度對中國縣域碳匯的影響遠高于產(chǎn)業(yè)結(jié)構(gòu)產(chǎn)生的影響。經(jīng)濟發(fā)展和產(chǎn)業(yè)結(jié)構(gòu)因素除了在個別地區(qū)對碳匯產(chǎn)生負相關(guān)外,大部分區(qū)域均呈正相關(guān)；土地利用方式則在全部區(qū)域與碳匯均呈負相關(guān)(圖7、圖8)。

4.1 經(jīng)濟發(fā)展對碳匯的影響

GDP的回歸系數(shù)波動較大(-11.216305—5.433819),說明經(jīng)濟發(fā)展對全國縣域碳匯影響的地區(qū)差異較大(圖6),但是仍然存在空間異質(zhì)性的現(xiàn)象。整體上看,2015年經(jīng)濟發(fā)展對碳匯影響能力較小,呈現(xiàn)出在西藏東南部和四川西部的回歸系數(shù)高絕對值區(qū)域。2015年林芝市米林縣、墨脫縣、朗縣,昌都市八宿縣、洛隆縣、邊壩縣以及山南市加查縣、錯那縣和嘉黎縣,說明經(jīng)濟發(fā)展對當?shù)氐奶紖R產(chǎn)生了較強的負影響,在這類區(qū)域應該轉(zhuǎn)變經(jīng)濟發(fā)展方式,促進增匯減排。而在甘孜阿壩自治州的德格縣、白玉縣等地回歸系數(shù)為正值且高于全國其它縣域,說明經(jīng)濟發(fā)展對當?shù)氐奶紖R產(chǎn)生了較強的正影響,說明經(jīng)濟發(fā)展促進了當?shù)氐沫h(huán)保業(yè)發(fā)展。

4.2 產(chǎn)業(yè)結(jié)構(gòu)對碳匯的影響

從產(chǎn)業(yè)結(jié)構(gòu)來看(圖7),產(chǎn)業(yè)結(jié)構(gòu)代表第二產(chǎn)業(yè)產(chǎn)值占整個國民生產(chǎn)總值的比重,從這一指標考察產(chǎn)業(yè)結(jié)構(gòu)對碳匯的影響,若碳匯隨著產(chǎn)業(yè)結(jié)構(gòu)的增加而增加,則說明技術(shù)進步在其中起到了積極作用。2015年,產(chǎn)業(yè)結(jié)構(gòu)與碳匯在內(nèi)蒙古赤峰市阿魯科爾自治旗、林西縣、翁牛特旗等地,內(nèi)蒙古昌都市貢覺縣、邊壩縣、查雅縣等地,西藏林芝市米林縣、墨脫縣等地呈正相關(guān),說明技術(shù)水平進步在這些地區(qū)對增匯產(chǎn)生了正面影響。然而在內(nèi)蒙古呼倫貝爾市阿榮縣、興安盟,黑龍江龍江縣、漠河縣、嫩江縣等地均呈現(xiàn)負相關(guān),表明產(chǎn)業(yè)結(jié)構(gòu)對碳匯產(chǎn)生負影響,即二產(chǎn)比值越高,碳匯越低,這代表在這些區(qū)域應該促進產(chǎn)業(yè)結(jié)構(gòu)升級或者技術(shù)進步以達到增匯減排的目的。

圖6 2015年GWR模型經(jīng)濟發(fā)展回歸系數(shù)的空間分布 Fig.6 Spatial distribution of GDP regression coefficients in the GWR model, 2015

圖7 2015年GWR模型產(chǎn)業(yè)結(jié)構(gòu)回歸系數(shù)的空間分布 Fig.7 Spatial distribution of industrial structure regression coefficients in the GWR model, 2015

圖8 2015年GWR模型土地利用程度回歸系數(shù)的空間分布 Fig.8 Spatial distribution of land-use degree regression coefficients in the GWR model, 2015

4.3 土地利用程度對碳匯的影響

土地利用程度越低,碳匯通常更高,隨著城市建設(shè)用地的擴張,碳匯量也會減少。土地利用程度與全國縣域碳匯呈現(xiàn)明顯的負相關(guān)性,影響程度在地域上存在差異性。整體上看(圖8),2015年土地利用程度對碳匯的影響程度沿胡煥庸線呈現(xiàn)鮮明的差異,即在西北半壁的影響程度高于東南半壁的影響程度,說明在我國西北半壁的縣域應該更加謹慎增加建設(shè)用地面積,合理利用和開發(fā)土地；在東南半壁的縣域應注重調(diào)整非建設(shè)用地的土地使用格局,增加生態(tài)用地面積,植樹造林,保護和增加林地、濕地等碳匯系數(shù)高的用地類型。

5 結(jié)論與建議

5.1 結(jié)論

通過對1990—2015年間中國縣域碳匯總量時空格局進行分析,基于標準差橢圓、局部空間自相關(guān)、冷熱點分析和地理加權(quán)回歸分析方法,探析中國縣域碳匯的時空分異特征及影響因素,其主要結(jié)論如下：

(1)從時空變化規(guī)律上看,1990—2015年中國碳匯總量呈波動下降趨勢,從13307.79×104t下降至13198.27×104t；但不同碳匯類型的比例結(jié)構(gòu)保持不變,碳匯類型主要以林低碳匯為主,占比約80%,依次是草地(約13.5%)、水域(約4.0%)和未利用地(約2.5%)。在空間分布上,中國縣域碳匯總體呈現(xiàn)差異化分布的空間格局,在區(qū)域上呈現(xiàn)不同的變化趨勢,大致呈現(xiàn)“西部>東北>南部>中部”的 “西高東低”格局。

(2)從空間分異和空間聚集格局來看,碳匯空間分布中心整體存在向西南移動的趨勢,空間分布范圍呈現(xiàn)收縮態(tài)勢,西南地區(qū)對全國縣域碳匯空間格局的影響作用加強；1990—2015年中國縣域碳匯總量的冷熱點空間分布集聚程度呈現(xiàn)波動穩(wěn)定的特征,而中國縣域碳匯強度的空間集聚程度變化不明顯,具有高值與低值聚集,高-低區(qū)域零星分布的特點。

(3)從影響因素來看,2015年經(jīng)濟發(fā)展對全國縣域碳匯產(chǎn)生了較弱的影響,產(chǎn)業(yè)結(jié)構(gòu)對全國縣域碳匯產(chǎn)生了較強的影響,土地利用程度對全國縣域碳匯產(chǎn)生了負影響,且都存在空間異質(zhì)性的現(xiàn)象。

5.2 建議

以全國縣域作為研究對象,以1990—2015每5年的土地利用數(shù)據(jù)為基礎(chǔ),從碳匯類型和碳匯總量等方面討論了中國縣域碳匯的時空變化規(guī)律、空間分布差異性以及其影響因素?；谘芯?在未來發(fā)展過程中可考慮以下低碳策略：①合理規(guī)劃縣域非建設(shè)用地的土地利用方式。林地碳匯作用主要碳匯類型的作用不容忽視,而濕地和高覆蓋度草地碳匯由于強大的碳吸收能力也應得到充分保護。因此,在縣域非建設(shè)用地區(qū)域應采取擴大林地面積, 保護低承載力草地和濕地,恢復退化土地,優(yōu)先開發(fā)未利用地等措施來提高其碳貯量,實現(xiàn)生態(tài)系統(tǒng)增匯減排的目的。②根據(jù)不同縣域的碳匯時空發(fā)展現(xiàn)狀,針對性制定增匯減排目標和方案。例如在我國中部以及部分東部存在碳匯低值集聚的地區(qū),其首要任務是促進產(chǎn)業(yè)的轉(zhuǎn)型升級,加強生態(tài)文明建設(shè),加大生態(tài)保護力度等；而在而在黑龍江、內(nèi)蒙古北部,西藏南部、廣東、湖南等高值集聚的區(qū)縣,則可采用碳交易的方式“以東帶西”,拉動西部經(jīng)濟發(fā)展,彌補東部的“生態(tài)赤字”。③嚴格控制建設(shè)用地規(guī)模。研究期內(nèi)建設(shè)用地的擴張對全國范圍內(nèi)的碳匯均產(chǎn)生負影響,因此應控制建設(shè)用地的低效擴張,提高土地利用效率,集約節(jié)約用地。

由于數(shù)據(jù)限制,計算采用精度1km的土地利用數(shù)據(jù)進行碳匯計算,數(shù)據(jù)的柵格較大,而采用的系數(shù)法對系數(shù)的精度要求較高,會一定程度上影響計算結(jié)果的精度。本文僅評價分析了經(jīng)濟發(fā)展、產(chǎn)業(yè)結(jié)構(gòu)和土地利用程度對碳匯的影響,對于格局形成的其他驅(qū)動因素和內(nèi)在機制有待深入研究。