新疆土地利用類型轉變的碳強度變化規律分析*

王志強，李雪薇

(新疆農業大學管理學院，烏魯木齊市，830000)

0 引言

氣候問題已經成為全球關注的焦點，2009年丹麥召開的哥本哈根會議拉開了低碳革命的序幕，作為社會經濟發展的產物，碳排放是導致溫室效應的重要原因[1]，碳循環過程與人類活動，特別是土地利用、土地覆被變化有著密切的關系[2]，國內外學者對于土地利用類型轉變與碳排放的研究已經較為成熟，曲福田等[3]透析了土地利用變化對碳排放的影響機理，Houghton[4]研究發現地類的轉換引起土壤有機碳的變化，從而產生碳源或碳匯的作用，Chi[5]認為土地轉換前植被類型和土地轉換后的時間對碳動態影響程度起決定性作用，林地、草地要優于受人類干擾的[6]，余德貴等[7]提出了土地利用的低碳調控方式，分別有學者研究了重慶[8]、京津冀地區[9]、黑龍江[10]等地區的土地利用和碳排放的時空變化趨勢及土地利用碳排放效應，并進行土地利用與碳排放的相關關系[11]和預測[12]，發現集約利用土地對區域碳排放效率產生正向作用，這種影響隨著碳排放效率的增加而減弱[13]。現有研究成果缺少對于土地資源充裕、經濟欠發達典型區域的論證，本文立足于新疆，進行遙感影像圖的空間分析，深入探析新疆土地利用類型轉變的碳強度變化，以期為新疆今后的土地利用與規劃提供數據支撐，為區域經濟高質量發展提供最優的土地配置方案。

1 研究材料

1.1 研究區域概況

新疆地處東經73°40′～96°18′，北緯34°25′～48°10′之間，是中國陸地面積最大的省級行政區，總面積約為中國陸地面積的1/6，邊界線長度約為中國陸地邊界線總長度的1/4，土地資源豐富，土地利用率僅為39.18%，建設用地占比不足1%，境內有獨具特色的大冰川，共計1.86萬余條，冰川總面積約為2.4×103km2，占全國冰川面積的42%。年平均氣溫10.5 ℃，年均降水量145.7 mm，陽光充裕，日照時間為2 691.9 h。主要土壤類型為干旱土、漠土、鹽堿土和高山土；主要植被類型為草原、草甸、高山植被和荒漠。新疆能源豐富，已探明的煤炭、石油、天然氣儲量均位于全國前列，依靠能源帶動城鎮化和工業化快速發展，土地利用類型轉變過程中帶來的碳排放量巨大，環境問題凸顯，亟需轉向綠色、低碳、可持續的發展模式。

1.2 數據來源及處理

1.2.1 土地利用數據來源及處理

研究采用的100萬植被類型空間分布數據、土壤類型空間分布數據以及1985年、1990年、2000年、2005年、2010年、2015年和2020年7個時間截面的土地利用遙感監測數據均來源于中國科學院資源環境科學數據中心。以GB/T 2010—2017 《土地利用現狀分類》為根本，兼顧其他分類標準，對新疆土地利用類型進行歸并、重編碼，共分為6大類，分別是耕地、林地、草地、水域、建設用地和未利用地，詳見表1。借助ArcGIS的轉換工具將植被類型空間分布數據、土壤類型空間分布數據和土地利用遙感監測數據轉為矢量文件并統一地理坐標系，利用1985年新疆土地利用圖層裁剪得出新疆植被類型空間分布圖，利用1990年新疆土地利用圖層裁剪得出新疆土壤類型空間分布圖。

1.2.2 碳密度數據來源及處理

土地利用類型碳密度包括植被碳密度和土壤碳密度兩部分。以新疆植被類型空間分布圖為基礎，借鑒賴力[14]學者關于不同類型的植被碳密度相關研究成果(表2)，通過ArcGIS的空間統計功能得到新疆植被碳密度分布圖；同理，以新疆土壤類型空間分布圖為基礎，借鑒于東升等[15]關于不同類型的土壤碳密度相關研究成果(表3)，得到新疆土壤碳密度分布圖，不同土地利用類型的綜合碳密度由土地利用類型植被碳密度和土地利用類型土壤碳密度加和而成。

表1 新疆土地利用類型分類Tab. 1 Land use classification in Xinjiang

表2 不同植被類型的碳密度Tab. 2 Carbon density of different vegetation types

表3 不同土壤類型的碳密度Tab. 3 Carbon density of different soil types

2 研究方法

2.1 土地利用動態度

土地利用類型的轉換導致地類碳源碳匯狀況改變，同時影響人類能源活動的方式和強度，進一步影響整個區域的碳過程[16]。基于ArcGIS的相交分析功能和Excel的數據透視功能，以5年為間隔，分別建立了新疆2000—2005年、2005—2010年、2010—2015年及2015—2020年的土地利用轉移矩陣，由此不僅可以清晰看到不同土地利用類型之間的轉化情況，還可以揭示不同土地利用類型間的轉移速率，詳見表4。

表4 新疆2000—2020年各階段土地利用轉移矩陣Tab. 4 Land use transfer matrix of Xinjiang in 2000—2020 km2

土地利用動態度模型可以模擬分析新疆各類型土地在一定時期內的變化速度，計算如式(1)所示。

(1)

式中：K——某一類型土地的變化速率；

Ua——某一類型土地在研究初期時的面積；

Ub——某一類型土地在研究末期時的面積；

T——研究時長，本文選取2000—2020年數據，因此T=21。

綜合土地利用動態度模型可以模擬分析新疆土地利用類型的整體變化趨勢，計算如式(2)所示。

(2)

式中：S——某一時段綜合土地利用動態度；

t——土地類型總數；

Ui——第i類土地在研究初期的面積；

ΔUi-j——研究時段內第i類土地轉為其他非i類土地面積的絕對值。

土地利用程度綜合指數模型可以測算土地利用效果，反映了人類活動與生態環境的交互作用，了解土地利用是否合理、土地利用與生態環境是否協調，具體公式如式(3)所示。

(3)

式中：L——土地利用程度綜合指數；

n——土地利用等級數；

Mi——第i類土地的利用程度分級指數；

Ci——第i類土地利用類型面積比例。

參照土地利用分級原則，得到建設用地的利用程度分級指數為4、耕地的利用程度分級指數為3、林地、草地的利用程度分級指數為2、未利用地的利用程度分級指數為1。

2.2 土地利用類型轉變的碳強度核算

為消除社會經濟的影響，更加直觀、充分地分析新疆土地利用類型轉變導致的碳強度變化，本文參考張梅等[17]的研究方法，運用ArcGIS的疊加分析和分區統計得到不同土地類型的綜合碳密度，并通過數學模型運算得出土地利用類型轉變的碳強度變化，流程如圖1所示。

圖1 土地利用類型轉變的碳強度計算流程Fig. 1 Calculation process of carbon intensity of land use type transformation

土地利用類型決定區域碳源碳匯能力的強弱，參考張婷[18]研究成果，發現林地及耕地的碳匯能力較強，建設用地表現為強碳源性，草地、水域及未利用地的碳源碳匯性不顯著，具體數值如表5所示。

不同土地利用類型的綜合碳密度差異顯著，因此可以用地類轉變前后的碳密度差估算由于土地利用類型轉變而導致的碳排放量變化。轉變后地類的碳匯能力決定了區域碳收支狀況。同時，考慮到土地利用的長期性以及土地利用過程中的生態響應的周期性，一般將生態系統對碳排放響應的時間周期定為20年。綜上，不同土地利用類型轉變的碳強度計算見式(4)。

(4)

式中：Pi——第i種土地利用類型轉變的碳強度；

瑞士研究人員最新發現，長期使用手機產生的射頻電磁場可能對特定腦區的記憶能力產生不良影響。研究結果顯示，在經常使用手機接聽電話者的右側頭部，射頻電磁場對圖形記憶功能產生了負面影響。不過，研究發現，與手機通話相比，用手機發送短信、玩游戲或瀏覽互聯網則對青少年記憶力發展并沒有明顯影響，這可能是因為這類操作產生的射頻電磁場相對微小。

C1i——轉變前地類的綜合碳密度；

C2i——轉變后地類的綜合碳密度；

N——生態系統對碳排放響應的時間周期，即N=20；

Ii——轉變后地類的碳匯值。

土地利用類型的轉變影響著區域生態系統碳排放的變化，隨著時間變化其碳匯能力也在變化，直至趨于穩定。陸地生態系統碳排放計算見式(5)。

(5)

式中：LC——陸地生態系統碳排放；

Ai——第i種地類轉變的面積；

m——土地利用類型轉變的種類數量，本研究中共有30種，即m=30。

表5 不同土地利用類型的碳匯能力Tab. 5 Carbon sink capacity of different land use types

3 結果與分析

3.1 土地利用動態度分析

不同土地利用類型的碳源碳匯能力不同，土地利用類型轉變影響區域的碳源碳匯格局，這種影響在耕地、林地、草地和建設用地4大地類之間的轉換尤為明顯。經研究林地、耕地表現為較強的碳匯能力，建設用地則為主要的碳源。林地轉變為耕地或草地產生高強度碳排，退耕還林、還草可以增加20%～40%的碳匯量[14]。根據2000年、2005年、2010年、2015年和2020年5期的新疆土地利用數據，得到不同時期的土地利用類型面積占比變化圖如圖2所示。其中，耕地面積占比由3.70%穩步上升至5.56%，增加了1.86%；林地面積由2.19%降至1.49%，下降了0.7%，其中，2015—2020年，林地面積顯著減少，由34 725.67 km2減少至24 392.73 km2，極大降低了固碳能力，增加了資源和環境壓力；建設用地面積占比由0.23%上升至0.45%，城鎮化的推進和工業崛起，增加了二氧化碳的排放，引起區域碳匯壓力增強；草地面積占比維持在28.6%～29.6%的穩定狀態；2000—2020年，水域面積減少了35.36%，水資源狀況不容樂觀。

圖2 新疆2000—2020年各土地利用類型面積占比變化圖Fig. 2 Change chart of area proportion of different land use types in Xinjiang from 2000 to 2020

土地利用的速率變化影響地表植被類型的改變，進而影響區域土壤的碳循環，改變土壤碳密度[19]，根據式(1)和式(2)，結合表4數據，得到新疆土地利用動態度，如表6所示。新疆綜合土地利用動態度在2005—2010年下降，但總體呈現上升的趨勢，特別是2015—2020年間，上升明顯。其中，水域、林地和建設用地變化較快，然而水域和林地的動態度是負向增長，阻礙了新疆的生態發展和綠色發展。2005—2010年間，建設用地的利用動態度僅為0.664 7%，在2010—2015年間躍升至6.430 0%，建設用地作為重要的碳源地類，其面積擴增與利用程度加強不斷增強區域的碳強度，對生態環境和區域綠色發展產生負向作用。

表6 新疆土地利用動態度Tab. 6 Land use dynamics in Xinjiang

人類活動施加于土地的強度反饋至土地的利用效果和土壤碳通量變化，土壤有機碳的吸收和排放影響碳循環速率。根據式(3)并結合表4數據，得到新疆土地利用程度，見表7。較于2000年，2020年新疆綜合土地利用程度提升了0.238 3%，其中未利用地的利用程度最高，其次是草地；未利用地主要開發為草地或種植林木，整體上增強了新疆的碳匯能力。耕地、建設用地的利用程度增加明顯，分別提升了5.584 4%和0.910 8%。林地、水域和未利用地的利用程度出現不同幅度的降低，其中，林地的利用程度由4.369 2%降為2.975 5%，水域的利用程度由6.239 7%降為4.061 3%，各自降低了1.393 7%和2.178 4%。

表7 新疆土地利用程度Tab. 7 Land use degree in Xinjiang

3.2 土地利用類型轉變碳強度變化分析

3.2.1 不同土地利用類型的綜合碳密度

以新疆植被類型空間分布圖為基礎，結合表2數據，通過ArcGIS的統計分析功能，得到新疆不同土地利用類型的植被碳密度，如表8所示。不同地類的植被碳密度差距較大，其中林地的植被碳密度最高，達148.45 kg/m2，其次是耕地、建設用地和草地。建設用地植被碳密度為39.01 kg/m2，與耕地相差7.96 kg/m2，說明城鎮綠化水平較高，具有一定的碳吸收水平。

表8 新疆不同土地利用類型的植被碳密度Tab. 8 Vegetation carbon density of different land use types in Xinjiang

以新疆土壤類型空間分布圖為基礎，結合表3數據，通過ArcGIS的統計分析功能，得到新疆不同土地利用類型的土壤碳密度，如表9所示。土壤碳密度最高的是林地，13.13 kg/m2，耕地和草地的土壤碳密度相差0.38 kg/m2，同時，耕地和草地間的轉換面積和速率處于較高水平；由于氣候、降水等原因，新疆水域大面積縮減，干涸的水域土壤有機質含量較低，土壤活性較差，難以開發利用；未利用地主要由戈壁、鹽堿地、沼澤地、裸土地、裸巖石質地等組成，植被覆蓋度低，土壤養分低，開發為農、牧、園、林難度較大，碳源碳匯性差。

表9 新疆不同土地利用類型的土壤碳密度Tab. 9 Soil carbon density of different land use types in Xinjiang

碳密度不僅可以用來衡量固碳能力，還可以準確估算區域碳儲量，對分析、研究、預測區域碳強度變化有重要作用。不同土地利用類型的綜合碳密度由土地利用類型植被碳密度和土地利用類型土壤碳密度加和而成，見表10。森林具有巨大的碳儲潛力，對于平衡區域碳排放，緩解碳壓力至關重要，林地的綜合碳密度最高，為161.58 kg/m2。新疆耕地的碳密度高于草地碳密度，差值為17.2 kg/m2，水域和未利用地的綜合碳密度較低，分別為22.04 kg/m2和13.11 kg/m2，不同土地利用類型的綜合碳密度差值明顯。

表10 新疆不同土地利用類型的綜合碳密度Tab. 10 Comprehensive carbon density of different land use types in Xinjiang

3.2.2 土地利用類型轉變的碳強度變化分析

根據式(4)、表4、表5和表10計算得出新疆土地利用類型轉變的碳強度變化，由表11可以看出：耕地轉變為林地可降低10.258 6 kg/m2的碳強度，增加區域碳匯效應，轉變為草地、水域和建設用地則會增加碳強度，而轉變為建設用地的碳源效應增加最明顯(3.638 1 kg/m2)，耕地棄耕將增加2.011 4 kg/m2的碳強度；林地轉換為其他地類均會增加碳強度，其中轉為建設用地、未利用地和水域的碳強度明顯高于轉化為耕地和草地的碳強度，分別為8.696 7 kg/m2，7.070 0 kg/m2和6.844 8 kg/m2、轉變為耕地、草地分別為3.758 6 kg/m2和5.377 6 kg/m2；草地轉為耕地、林地可以提升區域碳匯值，但草場退化、轉為水域、或進行城鎮建設等活動不同程度的增加了碳強度；未利用地除開發為建設用地增加碳強度，用于作物耕種、造林植樹、發展畜牧、建設水域等均會增加區域碳匯能力，其中開發為林地碳匯值增加最高(12.27 kg/m2)。

表11 新疆土地利用類型轉變的碳強度Tab. 11 Carbon emission intensity of land use change in Xinjiang kg/m2

任何地類轉換為林地均會減少碳強度，增強碳匯能力，任何地類轉為建設用地或荒廢則會增加碳排放，提高碳強度。耕地轉變為建設用地、未利用地、水域和草地均會增加碳排放；林地碳匯能力遠高于其他地類，當其轉變為其他地類時均會增加區域碳排放，碳強度由大到小依次是建設用地、未利用地、水域、草地、耕地；草地和水域之間的轉換對于碳排碳匯產生較為微弱的影響，碳強度變化較小；建設用地碳強度最高，建設用地轉為林地、耕地時碳匯作用比較明顯；通過人類活動，改變未利用地地表覆被，產生不同程度的碳匯能力，當未利用地用作建設時則產生碳排放。

3.2.3 土地利用類型轉變的碳強度變化時序分析

根據式(5)及表4計算得出新疆2000—2020年各階段土地利用類型轉變的碳強度變化及2000—2020年生態系統碳排放。2000—2020年新疆整體表現為碳匯作用，雖然2000—2015年這一階段是新疆大力推進城鎮化和工業化的時期，但是建設用地面積占比仍舊很小，林、耕、草占比超過35%，能夠化解生態系統碳排放量，2010—2015年新疆碳匯作用明顯加強，但效果仍不顯著，2015—2020年相較于前一階段碳匯作用降低了0.517 Tg。

表12 新疆2000—2020年土地利用類型轉變的碳強度變化Tab. 12 Carbon intensity change of land use change in Xinjiang from 2000 to 2020 Tg

新疆2000—2020年土地利用類型轉變導致的生態系統碳排放為-0.22 Tg，長期來看，總體表現為碳匯作用。

4 結論與對策

1) 從土地利用變化速率、土地利用程度可以明顯看出新疆2000—2020年的土地利用類型變化趨勢：其中耕地、建設用地占比逐年增加，2000年耕地面積占比3.7%，2020年占比為5.56%，建設用地面積2000年占比0.23%，2020年占比為0.45%；林地和水域面積卻逐年縮減，20年間林地面積縮減了11 488.87 km2，水域面積縮減了18 108.07 km2，生態問題較為嚴重。同時，2000年新疆土地利用率為38.53%，截至2020年土地利用率也僅為39.18%，新疆未利用地開發難的問題仍未有效解決。

2) 土地利用類型轉變通過影響植被覆蓋變化和土壤碳通量，改變區域的碳強度。不同地類的綜合碳密度決定了地類轉變前后區域的碳平衡狀況，新疆林地的綜合碳密度最高，達161.58 kg/m2，其次是耕地，綜合碳密度為55.35 kg/m2。通過計算得出：林地是主要的碳匯場所，其地類的轉變引起碳強度增加；建設用地是主要的碳源場所，當其他地類轉換為建設用地，碳強度顯著增強。通過不同方式開發未利用地，區域碳匯效果不同，大面積造林對于區域碳匯起到正向作用。

綜上，結合社會經濟發展需要，新疆的土地利用規劃應著重考慮以下幾方面。

1) 林地固碳聚碳能力顯著，同時可以保持水土，調節氣候，具有良好的生態效益和社會效益。2000—2020年，新疆的林地面積縮減，林地利用程度降低，利用動態度變化大，今后的發展中需落實保護林地政策。建成林時選擇適應新疆氣候的樹種，提高樹種成活率，充分挖掘林地潛力。

2) 未利用地轉化為其他地類時應以低碳為約束條件，結合新疆實際，應優先開發未利用地為林地，不僅可以穩碳固碳，還可以防風固沙、涵養水源、改良土壤；其次，新疆風力強勁、日照充裕，依靠資源稟賦優勢，優先在山地、戈壁等難利用地類設置風車、太陽能等可再生資源設備，通過技術創新優化能源消費結構，進一步節能減排，帶動區域經濟高質量發展；第三，未利用地開發為建設用地可以有效刺激經濟，帶動產業發展，但是新疆能源儲量豐富，資源驅動型經濟發展模式極易導致“資源詛咒”現象，需要更加嚴格的環境規制；第四，推進用地城鎮化和工業化過程中應合理規劃土地，形成居民集中、產業集聚，提高能源利用效率，實現綠色發展。