遼河口濕地自然植被碳儲量研究

張婷婷 石昊 蘆曉峰 楊國范

摘?要：濕地可以維持自然生態環境穩定和碳循環穩定，隨著城市化步伐的加快、海產養殖業的發展，遼河口濕地生態系統遭到了一定程度的破壞。利用遼河口濕地自然植被實測生物量和遙感影像提取的歸一化植被指數，建立了基于遙感影像的生物量回歸模型，經過對比分析得出二項式函數為最優模型。通過該模型反演出的生物量計算得到遼河口濕地自然植被的碳儲量為969 323.100 t，與實測值1 077 146.018 t相比，相對誤差為10.0%，表明利用遙感影像反演植被生物量，可以得到較精確的遼河口濕地自然植被的碳儲量。

關鍵詞：生物量;遙感影像;NDVI;回歸模型;碳儲量;遼河口濕地

中圖分類號：TM938.84;TV213.2?文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2020.10.019

Study on Carbon Reserves of Natural Vegetation in Liaohe River Estuary Wetland

ZHANG Tingting， SHI Hao， LU Xiaofeng， YANG Guofan

（Shenyang Agricultural University， Shenyang 110866， China）

Abstract：Wetland can maintain the stability of natural ecological environment and carbon cycle. With the acceleration of urbanization and the development of marine aquaculture， the wetland ecosystem of Liaohe River estuary has been damaged to a certain extent. Based on the analysis of measured biomass of natural vegetation and normalized vegetation index extracted from remote sensing images， the biomass regression model based on remote sensing image was established. Through comparative analysis， the binomial function was the optimal model. The carbon reserves of natural vegetation in Liaohe River estuary wetland is 969 323.100 t， which is 10% compared with the measured value of 1 077 146.018 t. Using remote sensing image to retrieve vegetation biomass， we can get more accurate carbon reserves of natural vegetation in Liaohe River estuary wetland.

Key words： biomass; remote sensing image; NDVI; regression model; carbon reserves; Liaohe River estuary wetland

濕地占地球陸地面積的5%～7%，而濕地包含的有機碳含量占陸地碳儲量的15%～30%[1]，濕地在維持自然生態環境穩定和碳循環穩定中具有舉足輕重的作用，研究濕地生態系統的植被碳儲量具有重大意義。

遼河口濕地位于渤海北部，在遼河入海口處，總面積約223 km2[2]，以蘆葦沼澤、草甸、河流水域、淺海灘和海域為主，是世界上生態系統保存最完整的濕地之一，也是東亞—澳大利亞水禽遷徙的中轉站，在濕地和生物多樣性研究中有著重要地位。近些年，隨著城市化步伐的加快、海產養殖業的發展[3]，遼河口濕地生態系統遭到了一定程度的破壞，僅1998—2013年該地區濕地面積減少近10.5%[4]。因此，開展遼河口濕地植被調查，研究遼河口濕地植被的生物量和碳儲量及其變化趨勢，對濕地保護具有重要意義。

1?研究方法

從宏觀角度出發，利用遙感影像計算植被指數，結合野外實測數據，通過相關分析建立遼河口濕地植被生物量模型并計算濕地植被生物量，進而推算出該濕地自然植被的碳儲量。

1.1?野外數據采集與處理

野外數據采集時間為2017年7—8月，該時段是遼河口濕地植被生長期，并且與用于建模的遙感影像的成像時間基本一致。在濕地內共實測樣方24個（見圖1），每個樣方設置1 m×1 m的采樣點兩個，采樣點內植被地上部分齊地面收割，標記后獨立裝袋帶回實驗室。

1.2?影像選取與處理

遼河口濕地植被開始生長的季節為春季，夏季和秋初生長旺盛。借助歸一化植被指數（NDVI）提取植被生物量，NDVI對植被蓋度的敏感度在植被生長初期最高，在植被生長后期最低[5]，而植被生物量在植被生長旺盛期最大，因此，基于研究遼河口濕地植被碳儲量的目的，選取2017年8月份的Landsat-8 OLI遙感影像提取歸一化植被指數。

原始影像中存在各種誤差，因此對遙感影像進行幾何校正—輻射定標—大氣校正—影像裁剪，從而得到遼河口濕地遙感影像圖。

2?建立生物量回歸模型

2.1?實測生物量

在實驗室將采集樣本放入烘箱中設置75 ℃烘干至恒重，取出稱重。取一個樣方中2個采樣點實測數據的平均值作為該樣方的生物量，實測生物量見表1。

2.2?歸一化植被指數提取

歸一化植被指數是植被生長狀態及植被覆蓋度的最佳指示因子[6]，是進行生物量建模的主要參數[7]。NDVI的計算公式：

NDVI=ρNIR-ρRρNIR+ρR（1）

式中：ρNIR為近紅外波段（波長0.85～0.88 μm）的反射值;ρR為紅光波段（波長0.64～0.67 μm）的反射值。

通過軟件對遙感影像進行波段運算，提取各樣方的NDVI平均值，見表1。

2.3?生物量回歸模型的建立

2.3.1?基本原理

通過遙感影像提取各樣方的NDVI平均值和實測生物量進行匹配分析，建立遼河口濕地生物量的半經驗模型，半經驗模型種類很多，其中最具有代表性的有線性函數模型、指數函數模型、二項式函數模型、對數函數模型4種。線性函數模型形式為

y=ax-b（2）

指數函數模型形式為

y=aebx（3）

二項式函數模型形式為

y=ax2+bx+c（4）

對數函數模型形式為

y=aln x+b（5）

式中：y為采樣樣方生物量;x為采樣樣方的NDVI平均值;a、b、c為常數，可以通過實測生物量和相應遙感影像的NDVI分析得出。

2.3.2?生物量回歸模型

由于13、16、24號樣方植被為翅堿蓬，其生物量和NDVI與其他樣方存在顯著差異，遼河口濕地翅堿蓬面積僅占整個保護區自然植被面積的4%[8]，因此在進行生物量回歸模型建立時把這3個樣方的數據略去。

建立的4種回歸模型見表2。對建立的4種回歸模型進行擬合分析，得到生物量回歸模型擬合結果，見圖2。擬合后得到線性函數模型、指數函數模型、二項式函數模型、對數函數模型的確定性系數R2分別為0.735 7、0.726 6、0.834 7、0.716 3，即二項式函數模型的確定性系數最大，擬合的精確度最高，因此采用二項式函數模型計算濕地生物量。

3?遼河口濕地自然植被碳儲量計算

遼河口濕地內由北向南隨機選取10株蘆葦、10株翅堿蓬，以及每種雜草3株，對樣本進行碳元素含量測定，得到蘆葦平均含碳率為0.551，翅堿蓬平均含碳率為0.450，雜草平均含碳率為0.510。保護區內蘆葦面積為377.019 km2，翅堿蓬面積為16.213 km2，混合草地面積為11.503 km2[8]。采用分類計算法、碳密度均值法、平均歸一化植被指數法3種方法對遼河口濕地自然植被碳儲量進行計算。

（1）分類計算法。該方法以實測生物量和植被面積為基礎進行植被碳儲量計算，自然植被碳儲量計算公式：

GC=∑ni=1CiSi（6）

Ci=BiomassiVCi（7）

式中：GC為濕地自然植被碳儲量;Ci為第i種植被的碳密度;Si為第i種植被的面積;Biomassi為第i種植被的生物量（實測值）;VCi為第i種植被的含碳率。

計算結果見表3。

（2）碳密度均值法。該方法根據研究區內所有植被實測生物量的平均值，計算出平均碳密度，進而計算出研究區植被的含碳量。

GC=C∑S （8）

C=Biomass（∑ni=1VCiSi∑S）（9）

式中：C為平均碳密度;Biomass為平均生物量;∑S為自然植被總面積。

計算結果見表4。

（3）平均歸一化植被指數法。根據遙感影像提取研究區內各樣方的NDVI平均值，根據生物量回歸方程y=247 881x2-209 583x+48 495，計算平均生物量Biomass后，根據式（8）、式（9）計算遼河口濕地自然植被碳儲量，計算結果見表4。

4?結果分析

（1）生物量模型分析。二項式函數建立的生物量模型的精確度比另外3種模型的高，R2值為0.834 7。生物量模型精確度見表5（*表示翅堿蓬樣方）。由表5可以看出，該生物量模型不適用于翅堿蓬。

（2）碳儲量計算分析。采用分類計算法計算出遼河口濕地植被碳儲量為1 077 146.017 t，該方法采用實測數據計算植被碳儲量，得出的數據最接近真值，因此以該數據為標準進行分析。

碳密度均值法計算出的碳儲量為976 660.400 t，比分類計算法少100 485.617 t，誤差為-9.3%。通過遙感影像建立的平均歸一化植被指數法計算的碳儲量為926 604.356 t，比分類計算法少150 541.662 t，誤差為-13.9%。進一步分析可見，碳密度均值法以平均生物量為前提求得平均碳密度，該方法受采樣點位置和植被種類影響較大。研究區翅堿蓬面積占自然植被總面積的4%，24個采樣點中翅堿蓬樣方不超過1個較為合適，但研究中翅堿蓬樣方選取了3個，導致該方法計算得到的平均生物量偏低，因此計算出的碳儲量低于實測值。翅堿蓬樣方數量也影響了NDVI的平均值，由此反演得到的Biomass值低于實測值，影響最終計算結果的精確度。將多余的翅堿蓬樣方去掉（NDVI為0.117、0.106的兩個樣方去掉，保留中間值0.112），再取研究區植被的NDVI平均值推算平均生物量為4 387.880 g/m2，進而計算得到碳儲量為969 323.100 t（見表4），誤差為-10.0%，由此可見選擇合適的采樣點提取NDVI反演生物量，能得到較精確的碳儲量。

5?結?語

基于野外實測生物量數據和遙感影像提取的NDVI數據，建立了遼河口濕地自然植被地上部分的生物量回歸模型，通過該模型可以較精確地計算出該濕地的生物量。通過該模型反演出的生物量計算得到遼河口濕地自然植被的碳儲量為969 323.100 t，與實測值1 077 146.017 t相比較，相對誤差為10.0%，表明利用遙感影像反演植被生物量，可以得到較精確的遼河口濕地自然植被的碳儲量。

參考文獻：

[1]?BERNAL B， MITSCH W J. A Comparison of Soil Carbon Pools and Profiles in Wetlands in Costa Rica and Ohio[J]. Ecolog-ical Engineering， 2008， 34： 311-323.

[2]?遼寧遼河口國家級自然保護區.遼寧雙臺河口國家級自然保護區正式更名為遼寧遼河口國家級自然保護區[EB/OL].[2019-05-08].http：//www.shidi.org/sf_F29C061E763546C1BBF8DD4F3F015E6F_151_lnsthkgjjzrbhqglj.html.

[3]?常雄凱，劉淼，李春林，等.遼寧沿海土地利用變化的圖譜特征[J].生態學雜志，2015，34（12）：3459-3465.

[4]?申丹，常禹，胡遠滿，等.遼寧沿海地區景觀格局變化及驅動力分析[J].江西農業大學學報，2015，37（5）：874-884.

[5]?楊國強.昌邑國家海洋生態特別保護區檉柳地上生物量與地上碳儲量遙感估算研究[D].呼和浩特：內蒙古師范大學，2017：8-10.

[6]?趙麗瓊.北京山區森林碳儲量遙感估測技術研究[D].北京：北京林業大學，2010：1-6.

[7]?葉垚.SAR影像在濕地植被地上生物量估測中的應用：以洪河保護區為例[D].北京：首都師范大學，2011：12-13.

[8]?蘭天慧.遼河口濕地生態退化評價研究[D].沈陽：沈陽農業大學，2018：22-26.

【責任編輯?呂艷梅】