結(jié)合地理國(guó)情普查成果的苕溪流域生態(tài)系統(tǒng)固碳釋氧價(jià)值估算

王 瑜，金姍姍，馮存均

(1. 浙江省地理信息中心，浙江 杭州 310012； 2. 浙江省地理國(guó)情監(jiān)測(cè)中心，浙江 杭州 310012； 3. 地理國(guó)情監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室，浙江 杭州 310012)

植被通過(guò)光合作用吸收CO2并釋放O2，調(diào)節(jié)大氣成分，在應(yīng)對(duì)氣候變化和大氣-陸地碳循環(huán)中發(fā)揮了重要作用，是生態(tài)系統(tǒng)重要的生態(tài)服務(wù)功能之一。多年來(lái)，已有不少學(xué)者對(duì)不同尺度的森林、草地等生態(tài)系統(tǒng)的服務(wù)價(jià)值進(jìn)行了研究。王兵等基于全國(guó)第7次森林資源清查數(shù)據(jù)及中國(guó)森林生態(tài)系統(tǒng)定位研究網(wǎng)絡(luò)長(zhǎng)期定位觀測(cè)數(shù)據(jù)，對(duì)全國(guó)各省區(qū)森林生態(tài)系統(tǒng)生態(tài)服務(wù)功能進(jìn)行了評(píng)價(jià)[1]。謝高地等以Costanza等對(duì)全球生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估的部分成果為參考，綜合對(duì)我國(guó)專(zhuān)業(yè)人士進(jìn)行的生態(tài)問(wèn)卷調(diào)查結(jié)果，建立了中國(guó)陸地生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值表，并在此基礎(chǔ)上評(píng)估了青藏高原生態(tài)資產(chǎn)的價(jià)值[2]。趙同謙等基于服務(wù)功能機(jī)制對(duì)中國(guó)草地的生態(tài)經(jīng)濟(jì)價(jià)值進(jìn)行了評(píng)價(jià)[3]。劉東等對(duì)浙江省森林生態(tài)服務(wù)功能價(jià)值進(jìn)行了估算，并對(duì)其逐月變化規(guī)律進(jìn)行了分析[4]。謝高地等在對(duì)青藏高原天然草地生態(tài)系統(tǒng)服務(wù)價(jià)值根據(jù)其生物量訂正的基礎(chǔ)上，根據(jù)高寒草地的地域分異特征對(duì)其生態(tài)價(jià)值進(jìn)行了評(píng)價(jià)[5]。也有學(xué)者專(zhuān)門(mén)就植被固碳釋氧價(jià)值進(jìn)行了探討。劉敏超等估算了三江源區(qū)植被固定CO2釋放O2的價(jià)值[6]。ZHANG Tao等評(píng)估了青海湖流域植被固碳釋氧價(jià)值，并分析了其1987—2010年的動(dòng)態(tài)變化規(guī)律[7]。固碳釋氧價(jià)值的估算除根據(jù)生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值估算[2，5]外，主要是先估算植被凈初級(jí)生產(chǎn)力(net primary productivity，NPP)，然后根據(jù)光合作用關(guān)系計(jì)算固碳釋氧價(jià)值[1，4，6-7]。植被凈初級(jí)生產(chǎn)力是固碳釋氧價(jià)值估算的關(guān)鍵步驟之一。根據(jù)前人的研究，CASA模型是一種廣泛用于植被凈初級(jí)生產(chǎn)力的估算模型[8-12]，具有參數(shù)較易獲取、利于動(dòng)態(tài)監(jiān)測(cè)等優(yōu)點(diǎn)。目前植被凈初級(jí)生產(chǎn)力估算模型使用的土地利用/覆蓋或植被分類(lèi)數(shù)據(jù)大多尺度較大，精度較低，制約了模型估算精度。

2013年7月，浙江省按照國(guó)務(wù)院的統(tǒng)一部署，開(kāi)始了第一次地理國(guó)情普查，歷時(shí)3年，圓滿完成了全省普查任務(wù)。本次普查采用優(yōu)于0.5 m分辨率的遙感影像作為主要數(shù)據(jù)源，搜集整合多行業(yè)專(zhuān)題數(shù)據(jù)，通過(guò)“室內(nèi)分析判讀、野外實(shí)地調(diào)查”等方法，獲取了多要素、全覆蓋、無(wú)縫隙的地理國(guó)情數(shù)據(jù)。地理國(guó)情普查成果具有精度高、時(shí)相新等特點(diǎn)，能夠?yàn)檎憬〉纳鷳B(tài)環(huán)境保護(hù)提供重要的數(shù)據(jù)支撐。

本文利用地理國(guó)情普查成果中的地表覆蓋數(shù)據(jù)，基于CASA模型估算浙江省苕溪流域植被凈初級(jí)生產(chǎn)力，進(jìn)而評(píng)估苕溪流域生態(tài)系統(tǒng)的固碳釋氧價(jià)值。

1 研究區(qū)概況

苕溪為浙江省八大水系之一，位于浙江省西北部，是太湖流域重要支流。流域地跨杭州、湖州兩市，流經(jīng)臨安、余杭、湖州城區(qū)、德清、安吉、長(zhǎng)興等縣(市、區(qū))。水系有東、西苕溪兩大源流，流域面積達(dá)4579 km2。境內(nèi)主要山脈為天目山，地勢(shì)自西南向東及東北逐步遞減和傾斜。全流域山地丘陵面積占88%，平原占12%左右。流域地處中熱帶季風(fēng)區(qū)北緣和北亞熱帶季風(fēng)區(qū)南緣，氣候溫和濕潤(rùn)，水熱同步，雨量充沛，四季分明。多年平均降雨量1460 mm，年均氣溫15.5℃～15.8℃。區(qū)內(nèi)生態(tài)環(huán)境和植被良好，如圖1所示。

圖1 苕溪流域地形暈渲圖

2 研究方法

2.1 CASA模型

CASA模型屬于參數(shù)模型，即光能利用率模型。CASA模型中NPP主要由植被所吸收的光合有效輻射(APAR)和光能利用率(ε)兩個(gè)變量來(lái)確定，其估算公式如下

NPP(x,t)=APAR(x,t)×ε(x,t)

(1)

式中，APAR(x,t)為t月在像元x處植被吸收的光合有效輻射(MJ·m-2)；ε(x,t)為t月在像元x處的光能轉(zhuǎn)化率(gC·MJ-1)。

2.1.1 APAR

太陽(yáng)輻射能是植物進(jìn)行光合作用的能量來(lái)源，其中能被植物有效吸收并進(jìn)行光合作用的部分叫光合有效輻射(PAR)，其波長(zhǎng)范圍約在0.38～0.78 μm。光合有效輻射一般約占太陽(yáng)總輻射的50%左右。其中，通過(guò)光合作用，并最終能固定為植物有機(jī)質(zhì)的這部分光能叫APAR(absorbed photosynthetic ally active radiation)，太陽(yáng)總輻射的大小和植物本身的生理生態(tài)特征共同決定了APAR的大小。

APAR(x,t)=SOL(x,t)×FPAR(x，t)×r

(2)

式中，SOL(x,t)為t月在像元x處的太陽(yáng)總輻射量(MJ·m-2)；FPAR(x,t)為t月在像元x處植被層對(duì)入射的光合有效輻射的吸收比例；r為植被所能利用的太陽(yáng)輻射占太陽(yáng)總輻射的比例，在此取值0.5。

FPARmin

(3)

(4)

式中，F(xiàn)PARmax=0.950，F(xiàn)PARmin=0.001，且FPARmax和FPARmin不隨植被類(lèi)型變化而變化；SR(x,t)為比值植被指數(shù)；SRmin為最小比值植被指數(shù)；SRmax為最大比值植被指數(shù)；NDVI(x,t)為歸一化植被指數(shù)。

2.1.2 光能利用率(ε)

光能利用率是指植被把所吸收的入射光合有效輻射(APAR)轉(zhuǎn)化為有機(jī)碳的比率(gC/MJ)。光能利用率因不同植被類(lèi)型及不同的生態(tài)環(huán)境而差異很大，包括氣溫、水分、土壤、營(yíng)養(yǎng)、疾病、個(gè)體發(fā)育、基因差異和植被維持與生長(zhǎng)的不同能量分配等影響因子。光利用率主要受溫度和水分的影響，其計(jì)算如下

ε(x,t)=Tε1(x,t)×Tε2(x,t)×Wε(x,t)×εmax

(5)

式中，Tε1(x,t)和Tε2(x,t)反映溫度的影響；Wε(x,t)反映水分條件的影響；εmax為理想條件下的最大光能轉(zhuǎn)化率。根據(jù)朱文泉等[13]，不同植被類(lèi)型最大光能利用率見(jiàn)表1。

Tε1(x,t)的計(jì)算公式為

Tε1(x,t)=0.8+0.02Topt(x)-0.000 5[Topt(x)]2

(6)

式中，Topt(x)為某一區(qū)一年內(nèi)NDVI值達(dá)到最高時(shí)月份的平均氣溫，認(rèn)為此溫度為植被生長(zhǎng)的最適溫度。當(dāng)某一月平均溫度≤-10℃時(shí)，Tε1取0，認(rèn)為光合生產(chǎn)力為零。Tε2(x,t)的計(jì)算公式為

(7)

表1 不同植被類(lèi)型最大光能利用率 gC/MJ

式中，若月平均氣溫T(x,t)比最適溫Topt(x)高10℃或低13℃時(shí)，該月Tε2(x,t)值為最適溫度Topt(x)時(shí)Tε2(x,t)的一半。Wε(x,t)的計(jì)算公式為

Wε(x,t)=0.5+0.5×EET(x,t)/PET(x,t)

(8)

式中，EET(x,t)為t月在像元x處的區(qū)域估計(jì)蒸散量(mm)；PET(x,t)為t月在像元x處的區(qū)域可能蒸散量(mm)。

當(dāng)該年降水量PPT(x,t)≥PET(x,t)時(shí)，EET(x,t)=PET(x,t)，即Wε(x,t)=1。

當(dāng)該年降水量PPT(x,t)

(9)

(10)

式中，PPT(x,t)為像元x處的年降水量(mm)；Rn(x,t)為像元x處的凈輻射量(mm·mon-1)；Ep(x,t)為局地可能蒸散量(mm)。計(jì)算模型如下

Rn(x,t)=[Ep(x,t)×PPT(x,t)]0.5×{0.369+

0.598×[Ep(x,t)/PPT(x,t)]0.5}

(11)

Ep(x,t)=16×[10×T(x,t)/I]a

(12)

a=(0.675I3-77.1I2+17 920I+492 390)×10-6

(13)

(14)

2.2 固碳釋氧價(jià)值計(jì)算方法

根據(jù)《森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范》(LY/T 1721—2008)，植被固碳量計(jì)算公式為

G固碳=1.63R碳·A·B年

(15)

式中，G固碳為植被年固碳量，單位為t·a-1；R碳為CO2中碳的含量，為27.27%；B年為植被凈初級(jí)生產(chǎn)力，單位為t·hm-2·a；A為面積，單位為hm2。

植被釋氧量計(jì)算公式為

G氧氣=1.19A·B年

(16)

式中，G氧氣為植被年釋氧量，單位為t·a-1；B年為植被凈初級(jí)生產(chǎn)力，單位為t·hm-2·a；A為面積，單位為hm2。

固碳釋氧價(jià)值計(jì)算，采用直接市場(chǎng)法、影子市場(chǎng)法，依據(jù)《森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范》中瑞典碳稅約為1200元/噸，我國(guó)國(guó)家衛(wèi)生部網(wǎng)站中2007年春季氧氣平均價(jià)格為1000元/噸。

3 數(shù)據(jù)及處理

本文中使用到的數(shù)據(jù)包括浙江省第一次地理國(guó)情普查成果地表覆蓋數(shù)據(jù)、氣象站點(diǎn)數(shù)據(jù)、MOD13Q1等。

地表覆蓋數(shù)據(jù)具有全覆蓋、高精度的特征，包括耕地、園地、林地、草地、房屋建筑(區(qū))、道路、構(gòu)筑物、人工堆掘地、荒漠與裸露地表和水域十大類(lèi)，各大類(lèi)又根據(jù)情況相應(yīng)細(xì)分為二級(jí)類(lèi)或三級(jí)類(lèi)。地表覆蓋類(lèi)型主要影響光能利用率。將地表覆蓋類(lèi)型與表1中指標(biāo)類(lèi)型進(jìn)行對(duì)應(yīng)，得到苕溪流域不同地表覆蓋類(lèi)型的最大光能利用率。植被光能利用率柵格化為250 m分辨率的柵格。

氣象站點(diǎn)數(shù)據(jù)包括12個(gè)月的月均溫、年降水量、年蒸發(fā)量。使用的氣象站點(diǎn)有長(zhǎng)興、湖州、安吉、德清、臨安、富陽(yáng)和杭州站共7個(gè)站點(diǎn)，采用IDW插值法進(jìn)行插值，柵格分辨率為250 m。

MOD13Q1來(lái)源于NASA(http:∥ladsweb.nascom.nasa.gov/data/search.html)，包括2015年24期16 d合成植被指數(shù)影像，分辨率為250 m，影像經(jīng)過(guò)坐標(biāo)轉(zhuǎn)換、鑲嵌、波段提取、裁剪等處理，形成苕溪流域2015年逐月NDVI。

太陽(yáng)總輻射通過(guò)ArcGIS軟件的Solar Radiation模塊進(jìn)行模擬[14]，該模塊輸入DEM數(shù)據(jù)，能夠根據(jù)研究區(qū)緯度、月份和地形去模擬研究區(qū)太陽(yáng)總輻射。柵格分辨率為250 m。

4 結(jié)果分析

經(jīng)CASA模型計(jì)算的苕溪流域凈初級(jí)生產(chǎn)力分布如圖2所示，苕溪流域平均NPP約為576.83 gC/m2。低值區(qū)多集中在城鎮(zhèn)、水域和平原水網(wǎng)區(qū)，這些地區(qū)植被覆蓋較少，NPP大多低于100 gC/m2。高值區(qū)多為山地，植被茂密，NPP集中在400～850 gC/m2及以上。境內(nèi)天目山和莫干山地區(qū)為NPP高值區(qū)。

圖2 苕溪流域NPP分布

根據(jù)NPP計(jì)算，苕溪流域全年植被固碳116.72萬(wàn)噸，按《森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范》中瑞典碳稅約1200元/噸算，苕溪流域全年植被固碳價(jià)值約14.01億元；苕溪流域全年植被釋放氧氣312.48萬(wàn)噸，按2007年春季氧氣平均價(jià)格1000元/噸計(jì)算，苕溪流域全年植被釋放氧氣價(jià)值31.25億元。苕溪流域全年固碳釋氧總價(jià)值45.26億元。

全年中5—9月水熱條件優(yōu)越，植被茂密，光合作用旺盛，固碳釋氧總價(jià)值與NPP同步，約占全年固碳釋氧總價(jià)值的79.90%，如圖3所示。

圖3 苕溪流域不同月份固碳釋氧總價(jià)值

各種地表覆蓋類(lèi)型中NPP從大到小依次為林地、園地、耕地、草地、其他和水域，其中其他部分包括房屋建筑(區(qū))、道路、構(gòu)筑物、人工堆掘地和荒漠與裸露地表。林地平均NPP達(dá)到736.04 gC/m2。水域是NPP低值區(qū)，平均NPP為194.39 gC/m2，主要來(lái)源于水生植物及由于混合像元效應(yīng)產(chǎn)生的影響。林地固碳釋氧總價(jià)值為35.10億元，占苕溪流域固碳釋氧總價(jià)值的77.57%。如圖4、圖5所示。

圖4 不同地表覆蓋類(lèi)型平均NPP統(tǒng)計(jì)

圖5 不同地表覆蓋類(lèi)型固碳釋氧價(jià)值占比

5 結(jié) 語(yǔ)

本文利用浙江省第一次地理國(guó)情普查成果，將GIS和RS技術(shù)運(yùn)用到NPP和固碳釋氧價(jià)值估算中，得出浙江省苕溪流域生態(tài)系統(tǒng)NPP與固碳釋氧價(jià)值的分布。本文所利用的地理國(guó)情普查成果地表覆蓋相對(duì)于傳統(tǒng)的土地利用或植被分類(lèi)數(shù)據(jù)，具有精度高、時(shí)相新的特點(diǎn)，更能反映地表植被的真實(shí)情況，提高NPP估算精度。

由于數(shù)據(jù)源的限制，本文所用的氣象數(shù)據(jù)為插值數(shù)據(jù)，站點(diǎn)較少，且為常年數(shù)據(jù)，非2015年當(dāng)年測(cè)量值；因太陽(yáng)輻射數(shù)據(jù)缺乏而采用ArcGIS軟件模擬，這些都會(huì)增加NPP和固碳釋氧價(jià)值的估算誤差。由于數(shù)據(jù)的限制，本文僅討論了植被光合作用所固定的CO2和釋放O2的價(jià)值，而土壤和凋落物呼吸作用還有待進(jìn)一步研究。