關于湖南省“雙碳”監測工作的技術探析*

謝優平，肖祥紅

（湖南省第二測繪院，湖南長沙 410119）

根據近幾年IPCC 氣象清單報告，因人類活動和自然變化釋放的溫室氣體（主要為CO2、N2O、CH4）脅迫地球表層溫度正以每十年0.2 ℃的幅度上升，導致地表冰土面積減少，海平面上升，極端災害（海嘯、地震、森林火災）增多，逐漸危害人類的生存安全。早在上世紀80年代，為應對改善人類生存環境的需要，不少國家或機構已經開始對碳源、碳匯進行了研究工作。2021 年9月22 日，在“中共中央國務院關于完整準確全面貫徹新發展理念作為碳達峰碳中和工作的意見(中發36 號文件)”中指出，“要依托和拓展自然資源調查監測體系，建立生態系統碳匯監測核算體系”；
2021 年我國在世界氣候峰會上向全世界承諾“中國力爭2030 年前CO2排放達到峰值，2060年底實現碳中和?！敝链?，“雙碳”監測工作正式提上國家層面，“雙碳”監測思路、任務及工作目標逐步明晰，包括前期的“雙碳”監測，中期的碳中和、碳達峰、碳交易，后期的碳規劃、清潔能源節點循環等內容。近年來，基于生態學的碳氮水循環研究及“雙碳”監測工作已取得一定進展，但由于自然生態的多樣性及繁雜性、人類活動干擾的不確定性以及技術手段的局限性，使得研究成果內容不充分、不全面、不徹底，用于碳收支的各種生物計量模型結果各有差異，研究技術路線不完善、不成熟。目前，大氣圈中基于衛星遙感的碳衛星監測技術和基于化學反應的碳捕集、碳存儲技術尚在初始研究階段，技術不成熟。因此，基于陸地生態循環的“雙碳”監測研究在今后相當長的時期內仍為工作重點。筆者就陸地生態系統中“雙碳”監測工作的若干個人觀點闡述如下，并對這些工作的開展進行技術探析。

在生態系統中，碳氮水循環主要通過光合作用、呼吸作用和水分蒸騰作用等生物物理方式實現土壤、植被、大氣三者之間的物質與能量交換。因此，本文認為碳水循環監測主要采取如下兩種方式。

一是通過監測碳氮水在大氣圈和生物圈之間的循環流向，遵守物質守恒定律。碳既是生命體最基本的元素，也是自然界中的較穩定的元素，無論是以有機碳形式存在于植物、微生物、藻類、土壤中，還是以CO2、CO 等形式存在于大氣中，生態系統的所含碳的總量不變，陸地生態系統物質循環和能量循環的基本流程如圖1[1-3]。因此，開展碳源碳匯監測時，既可以根據碳元素的流向來監測，也可以根據光合作用所形成的生物量來監測[4]。針對某一碳匯生態系統時，生物量的主要指標計算方式如下（式1~式3）：

圖1 陸地生態系統物質循環和能量循環流程Fig.1 Material cycle and energy cycle process of terrestrial ecosystem

若針對某一周期的區域生態系統，還需要將該時間段內人類活動及自然變化的擾動因素綜合納入計算中，如火燒、采伐、風蝕、水泄等對碳的去向干擾（式4）。

其中，NPP為植被凈初級生產力，APAR 為吸收的光合有效輻射，C 為實際光合利用率，RSG為太陽總輻射，FPAR 為入射太陽有效輻射的吸收率，GPP 為總初級生產力，NEP 為凈生態系統生產力，NBP 為凈生物群系生產力，R 自為植物自身呼吸所排放的CO2，R 異為同一生態鏈上微生物群呼吸所排放的CO2，NRP 為凈區域生產力，R 生物群為生態系統的其他動物和昆蟲呼吸所排放的CO2，R 泄是人為采伐、焚燒、風蝕、水泄等所產生的碳泄漏。

二是監測在某一特定時間內通過大氣中垂直界面的CO2濃度和脈動量，采用微氣象學中的空氣動力學通過通量觀測來監測，遵守動量和能量守恒定律[5-7]。

所有生物量監測都要考慮季節溫度、降水、土壤水分蒸騰蒸發、土壤肥力、近地表CO2濃度、物候條件等環境因子對生物生長的干擾和協迫作用。

湖南省域尺度陸地生態系統碳匯監測要綜合考慮自然因素、國家尺度的人為活動等多方面影響。監測對象除了常見的耕地、林地、草地、濕地四個一級碳匯種類外，還要包括水域碳匯類以及工業、建設用地碳源類。

不同于土地利用現狀分類，“雙碳”監測根據不同植被、不同時令應有更細的二級分類和影響因素考慮（表1）。如林地包含的闊葉、針葉、灌木，其生長季節、茂盛程度及表露特征不一樣，并且還要兼顧生物群系的林下植被；
又如耕地類作物可分為水稻（一季或兩季）、玉米、油菜、大豆、紅薯、棉花等，它們的生長季節、成熟周期不一樣，果實莖塊所處的位置（地上、地下）也不一樣；
再如水域中淡水湖泊有水下藻類的光合作用，海水有水下珊瑚礁的鈣化，此外，還要充分考慮水對CO2的溶解吸收，以及水的流速影響，流速越快，吸收能力越強。碳源監測不僅要考慮鋼鐵、水泥、化石燃料等污染工業的CO2排放統計，也要考慮建筑區、商業體育等人類活動區域CO2的排放監測[8]。

表1 湖南省“雙碳”監測陸地生態系統分類Table 1 Classification of terrestrial ecosystem for double carbon monitoring in Hunan Province

單一生物類的“源”與“匯”隨時間、季節動態變化，甚至可以互相轉換。因此，源與匯的表現應依據該生物類逐月、逐年或單一生長周期所固定或消耗的干物質量來判定。

在陸地生態系統的碳循環監測中，因包含的生物種類繁多，且各自生長季節不同，導致同一區域不同季節的碳固速率和碳固潛力迥然不同。要較為準確求取某一區域尺度生態系統的碳儲增量，按照季節或某一周期準確劃分每一植被地類的邊界尤為重要[9]。

植物地類的邊界識別以及智能提取主要分兩種方式。一是采用多光譜數據，通過對植被覆蓋長時間序列的生長周期觀測，計算生物生長過程的植被指數，如NDVI、RVI、EVI 等，根據各指數的閾值變化趨勢及收斂情況，判讀提取不同作物的邊界[14]。二是采用分辨率較高的衛星影像或航片影像進行提取。根據作物在影像上不同的紋理特征，結合作物生長周期特性，參考農作物樣本數據庫進行提取。兩種方法的提取結果需實地核查或比對三調數據庫地類進行調整和修正，以確保邊界準確度[10]。

陸地生態系統的生物計量及碳收支方法主要有如下兩類：（1）傳統的生物計量及碳收支方法。包括樣點連續清查法（薄記模型）、氣象數據統計模型、生物生態過程模擬、渦度相關系統通量觀測等。其中，樣點連續清查法具有極大的生態破壞性，而氣象數據統計模型及生物生態過程模擬普遍存在參數太多且求取困難、模型復雜且模擬周期序列長、工作量大等缺點。（2）遙感反演方法。包括基于衛星遙感的多光譜數據模型反演、碳衛星探測的溫室氣體（CO2）濃度反演、航測法計算活立木蓄積量反演碳儲量等。衛星遙感數據模型反演方法切合了我們地理信息行業的工作數據優勢和技術優勢，采用的衛星遙感數據時間序列長，覆被信息豐富，其中常用的光能利用模型驅動變量少，參數求取簡單，模型通俗易懂并可反復模擬[11]。

4.1 傳統或常用的生物計量方法

4.1.1 樣點連續清查法

耕地農作物、林地、草地最傳統的生物計量方式主要采用連續清查法加統計法。以月、季度、年為計量周期，以耕地農作物收支、森林、草地及土壤資源清查數據為統計基礎，通過生物量轉換因子分別將森林資源材積量轉化為森林生物量和碳儲量，利用草地資源清查資料及地下生物量數據，估算草地生態系統碳儲量，利用全國土壤普查資料和季節農作物收成估算農作物的地上、地下生物量，地下土壤碳儲量主要采用重鉻酸鉀稱重法[12]。

其中，森林材積轉化為生物量表示為（式5）：

V 為木材蓄積量，D 為木材密度，BEF是蓄積量與生物量之間的轉換因子，R為地上、地下生物量之比，CF為物質含碳量。

其它耕地農作物及草地收支亦可利用上式進行類推。樣點連續清查法主要特征是每個作物類、每個格網點都要大量設置樣點，工作量巨大繁重，破壞性大，優點是精度相對較高。

4.1.2 基于生物生長過程模擬的葉面積指數法（林木異速生長方程）

在生態系統中，葉面積指數是生態系統的一個重要結構參數，它反映植物葉面數量、冠層結構變化、植物群落生命活力及其環境效應，為植物冠層物質和能量交換提供結構化的定量信息。通過葉面積指數可以鑒定其它組織如葉莖、枝、軀干、根莖所對應的比重，即計算林分冠層的導度，通過建立林木生物的生長過程模型，確定某一時期段生物量的增長[13]（式6~式7）。

式中，LAI為葉面積指數，LE為實際有效葉面積指數，由LAI-2000 冠層分析儀測量，VE為葉總面積與族面積之比，對于闊葉VE=1，ΩE為集聚效應[18]。

式中，k為透光率，L為葉面積指數，Y 為NPP與GPP 間的轉化率，δ 為樹冠覆蓋率，ARPR 為太陽輻射，△t為時間。

采用樣點調查法的碳固量增長值對上式生物量模型求得NPP 進行校正，進而求得某一區域尺度的生態系統的固碳速率和固碳潛力。

利用葉面積指數法計算生物量時，其葉面積實際指數測定和異速生長方程求取困難，且干擾因素多，計算結果準確性較樣點清查法差。

4.1.3 通量觀測法

碳通量觀測是監測某一時刻大氣與植被之間、大氣與土壤之間、土壤與根系之間碳氮水的交換過程，是大氣圈、生態圈、土壤圈3 個界面之間相互作用的結果，是碳流向監測最為直接的一種方式[14]。

碳通量觀測包括生態系統的能量輸入、能量輸出（土壤-大氣界面、植被-大氣界面的輻射通量、顯熱通量和潛熱通量）、動量傳輸通量及氣體交換通量幾個方面的監測（式8）。

因切應力位于大氣層之下與地表層之上，風速儀與CO2脈動裝置必須架設在監測對象植被的上方。而林地植被普遍較高，導致通量觀測實施比較困難，因此，通量觀測方式主要適用于開闊建設用地上的碳源觀測和耕地碳匯觀測。

微氣象學中碳氮水通量觀測主要采取空氣動力學法、熱平衡法和渦度相關法。其中渦度相關法是通過測量大氣中湍流運動所產生的風速脈動和CO2濃度脈動，求算能量和物質通量的直接測定法，該方法是目前世界上CO2和水通量測量的標準方法，也是被認定為誤差來源最小的方法。

4.2 衛星航測遙感技術在生物計量測定中的運用

隨著遙感技術的日趨成熟，各種多光譜衛星數據源在時間空間獲取方面有了很大的便利性和適時性。在模擬植被蒸騰、土壤水分子蒸發及植被光合作用過程時，植被覆蓋是生物控制大氣、溫度、雨水、土壤侵蝕等方面的關鍵因素，可以貼切地反演林木生態過程和生物量的增長。而遙感數據的影像光譜、空間分辨率以及時相特征能夠反映不同空間尺度上的植被覆蓋信息及變化趨勢。因此利用遙感的光譜數據與植被覆蓋之間密切的對應關系，建立了40多種反映生物生長態勢的植被指數，并通過植被碳利用率CUE、植被光能利用率LUE、植被水分利用率WUE 幾種方式來求取植被凈生產力NPP。其中光能利用率模型因理論基礎性強、參數少且易獲取、數據覆蓋時空分辨率高、計算方法簡單等優點，成為當下求取植被凈生產力的主流方式[15]。

（1）植被光能利用率模型主要有四類驅動變量，①太陽光照輻射因子：包括有逐日、逐月、逐年光照輻射和散射數據；
②植被指數因子：主要有歸一植被指數NDVI、增強指數EVI、葉面積指數LAI 等；
③溫度脅迫因子：日平均溫度、日最低溫度、日最高溫度、白天地表溫度、夜間地表溫度；
④水分脅迫因子：包括水氣壓差VBD、蒸發比EF和地表水分指數LSWI[16]。

通過光能利用率模型，利用四類驅動因子來研究碳氮水在大氣和生物圈中的相互作用與循環。因不同植被類型在不同區域、不同時間、不同物候條件下光能利用率均有較大差異，需利用同期通量觀測數據來校正模型準確性，進行模型參數本地化（式9-10）。

光能利用率通用模型[17]（式9-10）：

光能利用率模型還可加入渦度相關數據、植被生長呼吸、物候條件、葉片陰陽兩面對光照吸收的特性差異等因素影響項，對模型進行優化和改進[18]（式11-12）。

兩葉光能利用率模型（式13）：

上述各種驅動變量數據可從MODIS 遙感數據集、CERN數據共享系統、中國陸地通量觀測網絡和中科院資源環境科學數據中心下載。

其中，RSG 為太陽總輻射，PAR 為光合有效輻射，FPAR 為光合輻射吸收比例，C 為光能利用率，T為溫度脅迫因子，W 為水分脅迫因子，Cmax為最大光能利用率，PAR0為光合有效輻射半飽和參數，C陽及PAR陽為陽葉最大光能利用率和陽葉吸收的太陽輻射，C陰及PAR陰為陰葉最大光能利用率和陰葉吸收的太陽輻射，P 為物候條件，a、b為呼吸線性參數，NEE夜為夜間通觀測值[19]。

（2）利用激光雷達和激光掃描儀進行生物量的測定（主要針對森林碳匯）

激光雷達數據穿透力強，受水霧干擾因素小，可以準確地測定森林樹木、植物的高度。又因雷達數據還可以測定地物的平面位置，因此也可求取林木的胸徑。DSM 與DEM 之差即樹木或植物的高度（圖2）[20]，高度值乘以面積得到森林體積。與連續清查法類似，上述數據模型再匹配森林生物量與蓄積量的轉換因子、木材密度、碳密度，可求取某一特定生態系統的碳儲量[21]（式14~16）。

圖2 森林或植被高度計算原理圖Fig.2 Schematic diagram of forest or vegetation height calculation

其中，CP 是森林碳儲量，V 是木材蓄積量，S 表示森林面積，D 是木材密度，BEF 是生物量與蓄積量的轉換因子，CF是碳密度。

采用生物量模型與遙感反演技術相結合的方式，可以部分替代傳統連續清查法實地調查工作內容，具有速度快，周期短，時效性強，經濟效益高，破壞性低等特點。但相對于傳統清查法來講，精度低，準確性不夠，需要大量的調查數據進行校正。

在區域尺度的碳匯、碳源監測中，森林、耕地、草地生物量的計算主要依靠生物生態過程模型來進行。在此推算過程中，要針對不同地域、不同類型，按一定格網設置樣點，進行生物量實測。再采用最小二乘法，通過迭代的方式對模型進行分類校正和平滑處理，處理后的模型值還要用實測值進行驗證，以確保均方差∑?S2最小或擬合優度（決定系數）R2最大（圖3）[22]。

圖3 陸地生態系統中計量生物碳儲量校正流程Fig.3 Calibration process of measuring biological carbon storage in terrestrial ecosystem

區域各類用地在某一周期的碳儲增量即為該區域尺度生態系統同周期的碳匯或者碳源[23]（式17）。

CER：某一時間周期的凈碳排放速率，CSR：某一時間周期的固碳速率。

其中耕地要針對不同的農作物或季節性的輪種農作物進行統計，且林地、耕地生物量要用歷史清查數據進行驗證，并計算年度碳匯能力和最大碳匯潛力。水域、雪地、巖石地、裸露地的碳匯碳源可以忽略不計。海洋、湖泊應分別監測統計。建設用地及工業區可以采取通量觀測，用該區域年度的化石燃料的排放統計量來檢驗碳源總量。

區域尺度的陸地生態系統碳源、碳匯監測，涉及數學、微氣象學、林學、農學、攝影遙感、地信GIS等基礎學科、邊緣學科和交叉學科，研究對象包括自然進展和人類活動的干擾，不確定因素太多。陸地生態系統的碳氮水循環其物理化學過程尚不明晰，沒有規范的理論體系和技術路線。迄今研究的對象僅限于林地、耕地、草地和海洋幾個領域。而人類活動對碳源、碳匯的擾動因素研究成果積累很少，建設用地、水域碳源、碳匯研究幾乎沒有。且各研究機構所取得的監測成果差異較大，無法作為準確的監測數據進行參考，各種生物計量模型亟待改進完善。同時，分布于國內耕地、林地的通量觀測網點少，觀測周期不長。種種因素極大制約了“雙碳”監測的工作發展。在總結、繼承前人研究成果的基礎上，亟需探索制定一套貼合實踐的監測方法和體系。