馬彩虹, 任志遠(yuǎn)

(1.陜西師范大學(xué) 旅游與環(huán)境學(xué)院, 陜西 西安 710062; 2.陜西理工學(xué)院 秦嶺與蜀道地理研究所， 陜西 漢中 723001)

?

陜西黃土臺塬區(qū)植被-土壤系統(tǒng)碳收支估算

馬彩虹1,2, 任志遠(yuǎn)1

(1.陜西師范大學(xué) 旅游與環(huán)境學(xué)院, 陜西 西安 710062; 2.陜西理工學(xué)院 秦嶺與蜀道地理研究所， 陜西 漢中 723001)

摘要：[目的] 為土地利用的低碳戰(zhàn)略提供決策依據(jù)。 [方法] 基于GIS和RS分析與處理NDVI數(shù)據(jù)和Landsat遙感影像數(shù)據(jù)源,采用CASA光能利用模型對陜西黃土臺塬區(qū)植被—土壤系統(tǒng)的碳收支進(jìn)行估算。 [結(jié)果] 1990—2010年陜西黃土臺塬區(qū)植被—土壤系統(tǒng)碳源效應(yīng)明顯。其中,咸陽臺塬區(qū)、寶雞臺塬區(qū)和西安臺塬區(qū)有所減弱,但渭南臺塬區(qū)和銅川臺塬區(qū)有所增強(qiáng)。研究時段內(nèi),陜西黃土臺塬區(qū)的高碳源區(qū)的面積由38.22%縮減為21.13%,中碳源區(qū)由57.40%擴(kuò)大為67.71%,低碳塬區(qū)由3.43%擴(kuò)大為7.76%;碳匯區(qū)由0.96%擴(kuò)大為3.40%。 [結(jié)論] 1990—2010年該區(qū)建設(shè)用地占用耕地、林地、草地,土地利用結(jié)構(gòu)變化不利于碳固定,提高復(fù)種指數(shù)有利于降低植被—土壤系統(tǒng)的碳排放。

關(guān)鍵詞：CASA; 碳源; 碳匯; 陜西黃土臺塬區(qū); 植被—土壤系統(tǒng)

工業(yè)革命以來，以二氧化碳為代表的溫室氣體排放已成為全球氣候變暖的主要原因[1]。作為世界上碳排放大國之一，中國政府承諾到2020年單位國內(nèi)生產(chǎn)總值二氧化碳排放比2005年下降40%～45%[2]。如何有效減少碳排放、增加生態(tài)系統(tǒng)的碳匯能力是當(dāng)前研究的熱點(diǎn)問題[3]。土地利用與覆被變化直接影響陸地生態(tài)系統(tǒng)的碳儲量和碳通量過程[4-5]。20世紀(jì)80年代中后期以來，學(xué)界對不同區(qū)域不同類型的生態(tài)系統(tǒng)碳收支作了大量研究，并取得了可喜的成果[6]，但評估結(jié)果仍然存在較大差異。如葛全勝[7]等的研究發(fā)現(xiàn),過去300年中土地利用變化導(dǎo)致中國陸地生態(tài)系統(tǒng)的碳排放達(dá)4.50～9.54 Pg;方精云[8]等的研究卻表明中國土地利用活動引起陸地植被—土壤系統(tǒng)的碳吸收大約為0.45 Pg；張秀梅[9]等對江蘇省土地利用的碳排放研究發(fā)現(xiàn)，土地利用和碳排放之間存在復(fù)雜的耦合關(guān)系，土地利用的碳源/匯效應(yīng)存在較大的空間差異性和不確定性[10]，需要對不同地方的不同生態(tài)系統(tǒng)碳收支情況進(jìn)行研究，以探究其內(nèi)在的規(guī)律和特征。陜西黃土高原臺塬區(qū)是西北黃土高原臺塬區(qū)的主體，既處于關(guān)中平原生態(tài)保護(hù)的屏障地帶，又是糧食生產(chǎn)的主體功能區(qū)，植被—土壤系統(tǒng)的碳收支對區(qū)域碳排放有較大影響。本文對1990—2010年陜西黃土臺塬區(qū)植被—土壤系統(tǒng)的碳收支情況進(jìn)行估算，以期為土地利用的低碳戰(zhàn)略提供一定的參考依據(jù)。

1數(shù)據(jù)與方法

1.1　研究區(qū)概況

研究區(qū)陜西黃土臺塬區(qū)位于34°08′—35°52′N，106°20′—110°36′E，是介于生態(tài)脆弱區(qū)和糧食主產(chǎn)區(qū)之間的過渡地帶，呈帶狀展布于渭河谷地南北兩側(cè)，總面積1 263 773 hm2。按行政隸屬關(guān)系可分為寶雞臺塬區(qū)、咸陽臺塬區(qū)、西安臺塬區(qū)、銅川臺塬區(qū)和渭南臺塬區(qū)[11]。處于季風(fēng)向干旱半干旱氣候過渡帶，多年均溫12.46 ℃，年均降水量456.3～735.9 mm，主要集中于夏季。土壤以肥力較高的土分布最廣，其次是黃墡土。水資源短缺，降水是地表水、地下水以及土壤水的主要補(bǔ)給來源；地下水埋深達(dá)50 m以下，開發(fā)利用較困難。人口分布較多，產(chǎn)業(yè)相對比較密集。根據(jù)第六次人口普查數(shù)據(jù)資料，截至2010年，陜西黃土臺塬區(qū)的人口為1.01×107人，人口密度約為845.84人/km2。近年來，經(jīng)濟(jì)發(fā)展較快，城市化水平也有較大提高。

1.2　研究方法

基于碳循環(huán)理論，生物生產(chǎn)力是生態(tài)系統(tǒng)碳循環(huán)過程以及各種碳吸收和排放分量平衡的結(jié)果，因此，植被—土壤系統(tǒng)的碳收支估算主要結(jié)合生物生產(chǎn)力進(jìn)行。早期的研究主要關(guān)注植被的凈初級生產(chǎn)力NPP[12]，即綠色植物在單位時間、單位面積上所積累的有機(jī)物數(shù)量。后來研究發(fā)現(xiàn)NPP只反映植物固定和轉(zhuǎn)化光合產(chǎn)物的效率，不能反映生態(tài)系統(tǒng)的碳源/匯功能。CASA模型是一個充分考慮環(huán)境條件和植被本身特征的光能利用率模型，適合于估算區(qū)域尺度上的植被NPP[13-15]。NPP減去土壤異養(yǎng)呼吸Rh所消耗的光合作用產(chǎn)物之后，即可獲得凈生態(tài)系統(tǒng)生產(chǎn)力NEP，可以直接定性定量地描述植被—土壤系統(tǒng)的碳收支情況。若NEP>0為碳匯，反之則為碳源。CASA模型系列公式如下：

NEP=NPP-Rh

(1)

式中：NEP——植被凈生態(tài)系統(tǒng)生產(chǎn)力； NPP——植被凈初級生產(chǎn)力；Rh——土壤微生物呼吸量。本文中Rh的取值來自于前人研究結(jié)果，研究區(qū)Rh為0.14～0.19，本文取其均值0.17。下同。

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

(2)

式中：APAR——吸收光合有效輻射(MJ/m2)。下同。

(4)

式中：FPAR(x,t)NDVI——由NDVI計算所得的植被層對入射光合有效輻射的吸收比例； FPAR(x,t)SR——比值植被指數(shù)(SR)計算所得的植被層對入射光合有效輻射的吸收比例；FPARmax及FPARmin的值為常數(shù)，其值與植被類型無關(guān)，F(xiàn)PARmax=0.95，F(xiàn)PARmin=0.001；NDVIi,max以及NDVIi,min——第i種植被類型的NDVI最大值和最小值。

(5)

式中：SR(x,t)——t月份在像元x處的太陽總輻射(MJ/m2)，常數(shù)0.48表示植被所能利用的波長為0.38～0.71 μm的太陽有效輻射占太陽總輻射的比重大小。SRi,max，SRi,min——第i種植被類型的比值植被指數(shù)SR的最大值及最小值。NDVI(x,t)——t月份像元x處的歸一化植被指數(shù)和比值植被指數(shù)。下同。

(6)

ε(x,t)=f1(x,t)×f2(x,t)×w(x,t)×εmax

(7)

式中：ε(x,t)——像元x在t月份的實(shí)際光能利用率。FPAR值(植被層對入射光合有效輻射的吸收比例)；f1(x,t)——低溫脅迫指數(shù)；f2(x,t)——高溫脅迫指數(shù)；w(x,t)——水分脅迫系數(shù)；εmax——理想條件下的最大光能利用率(g C/MJ)，取值根據(jù)朱文泉模擬值[16]。下同。

(8)

式中：Topt(x)——研究區(qū)域內(nèi)一年之中NDVI值為最高值時的當(dāng)月平均氣溫(℃)。

(9)

式中：T(x,t)——t月的平均氣溫(℃)。下同。

w(x,t)=0.5+0.5×E(x,t)/Ep(x,t)

(10)

式中：EP(x,t)——研究區(qū)的潛在蒸散量(mm/月)；E(x,t)——研究區(qū)的實(shí)際蒸散量(mm/月)，潛在蒸散量及實(shí)際蒸散量的值也取朱文泉模擬值[15]。

RH=0.22×〔Exp(0.0913T)〕+ln〔0.3145R+1)〕×30×46.5%

(11)

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

黃土臺塬區(qū)的范圍較小，所需氣象數(shù)據(jù)主要來源于該范圍內(nèi)的氣象站點(diǎn)如銅川站點(diǎn)(53 947)、耀縣站點(diǎn)(57 037)、鳳翔站點(diǎn)(57 025)、寶雞站點(diǎn)(57 016)、武功站點(diǎn)(57 034)、西安站點(diǎn)(57 036)、華縣站點(diǎn)(57 046)。所需1990年數(shù)據(jù)通過遙感影像提取，2010年NDVI數(shù)據(jù)通過SPOT 衛(wèi)星NDVI數(shù)據(jù)獲取。二者精度有一定差別，但不影響總體NPP的計算精度的比較效果。土地利用變化數(shù)據(jù)主要來源于Landsat陸地資源衛(wèi)星影像，基于ERDAS9.2進(jìn)行遙感影像的校正、投影轉(zhuǎn)換及解譯；基于ArcGIS10.0支持進(jìn)行數(shù)據(jù)處理、統(tǒng)計分析以及制圖輸出。

2結(jié)果與分析

2.1　植被-土壤系統(tǒng)的碳源/匯效應(yīng)估算

1990年，臺塬區(qū)的NEP位于碳含量-99.9~75 g/(m2·a)和-124.9~100 g/(m2·a)之間的范圍較大，其中，前一個區(qū)間的值主要分布在陜西黃土臺塬區(qū)中部和北部，后一個區(qū)間的值的范圍主要出現(xiàn)在渭南臺塬區(qū)的西南部和咸陽臺塬區(qū)大部分范圍，反映出臺塬區(qū)的植被—土壤系統(tǒng)的碳源效應(yīng)明顯。2010年，臺塬區(qū)NEP的值位于碳含量-124.9~100 g/(m2·a)范圍的面積有所縮小，而 -49.9~725 g/(m2·a)之間區(qū)域有所增大。其中，渭南臺塬區(qū)植被—土壤系統(tǒng)的碳源效應(yīng)有所減弱，咸陽臺塬區(qū)的的碳源效應(yīng)有所增強(qiáng)(附圖1)。分析可知，研究時段內(nèi)陜西黃土臺塬區(qū)植被—土壤系統(tǒng)在大部分區(qū)域呈現(xiàn)出碳源效應(yīng)，只有東北部邊緣區(qū)很小的范圍表現(xiàn)出碳匯效應(yīng)。

2.2　植被-土壤系統(tǒng)的碳源/匯效應(yīng)分區(qū)評價

根據(jù)陜西黃土臺塬區(qū)NEP出現(xiàn)區(qū)間，采用等間距法劃分為4個等級，即高碳源區(qū)、中碳源區(qū)、低碳源區(qū)和碳匯區(qū)(表1)。

從分布區(qū)域上來看，1990年高碳塬區(qū)主要分布在咸陽臺塬區(qū)、西安臺塬區(qū)和渭南臺塬區(qū)的中部，中碳源區(qū)主要分布在渭南臺塬區(qū)北部、東部和南部、咸陽臺塬區(qū)東部和西部、以及寶雞臺塬區(qū)北部區(qū)域，低碳源區(qū)和碳匯區(qū)主要分布在渭南臺塬區(qū)北部邊緣地帶。2010年，高碳源區(qū)范圍縮減了很多，主要分布在寶雞臺塬區(qū)西部和東部、咸陽臺塬區(qū)西部和東部、西安臺塬區(qū)、渭南臺塬區(qū)中部也有零星分布區(qū)；中碳源區(qū)主要分布于咸陽臺塬區(qū)中部和東部、銅川臺塬區(qū)、渭南臺塬區(qū)南部、西部和中部大部分區(qū)域；低碳塬區(qū)主要分布為渭南臺塬區(qū)，碳匯區(qū)面積擴(kuò)展明顯，主要分布于渭南臺塬區(qū)北部邊緣區(qū)(附圖2)。

如表2所示，2010年高碳源區(qū)大幅降低，中碳源區(qū)是陜西黃土臺塬區(qū)面積最大的一類，低碳塬區(qū)比例較小?？梢钥闯?，2010年土地利用的碳源效應(yīng)有所降低，表明1990—2010年陜西黃土臺塬區(qū)的植被覆蓋狀況改善程度較大。

表1　陜西黃土臺塬區(qū)碳/塬匯等級劃分標(biāo)準(zhǔn)　g/m2

表2　陜西黃土臺塬區(qū)碳源/匯分區(qū)面積及比例

2.3　植被-土壤系統(tǒng)碳收支變化原因分析

闊葉林、針葉林、灌叢、草原、草甸、荒漠以及高山稀疏植被等是陜西黃土臺塬區(qū)主要的植被類型。陜西黃土臺塬區(qū)農(nóng)業(yè)歷史悠久，天然植被很少，主要以人工植被為主。從土地利用變化的基本情況來看，研究區(qū)以耕地為主，其次為草地，林地比例很低(表3)。測算的NEP密度值，發(fā)現(xiàn)大部分區(qū)域表現(xiàn)為NEP<0，植被—土壤系統(tǒng)表現(xiàn)為碳排放效應(yīng)，只有臺塬區(qū)東北部的林地集中分布區(qū)表現(xiàn)為微弱的碳匯效應(yīng)。2010年，陜西黃土臺塬區(qū)耕地、林地、草地等貢獻(xiàn)NPP的地類均有所下降，而建設(shè)用地凈增加了26 538.3 hm2。因此，從土地利用情況來看，陜西黃土臺塬區(qū)的土堆利用變化趨勢不利于碳固定。

從NDVI的變化情況來看，1990年陜西黃土臺塬區(qū)的東北部、東南部及西南部的植被覆蓋狀況較好，即渭南市及寶雞市的植被覆蓋狀況較好，西安市較差。2010年，寶雞臺塬區(qū)、中部的咸陽臺塬區(qū)植被覆蓋比1990年明顯改善，渭南臺塬區(qū)和銅川臺塬區(qū)植被覆蓋有所下降。

由于本文研究范圍不是很大，因此降水和氣溫變化的空間差異性對植被—土壤系統(tǒng)的碳收支影響不是很明顯。因此，這種以農(nóng)田生態(tài)系統(tǒng)為主導(dǎo)的植被覆蓋情況下，研究區(qū)的NEP應(yīng)該主要受復(fù)種指數(shù)的影響。已有研究表明，研究時段內(nèi)陜西黃土臺塬區(qū)的復(fù)種指數(shù)總體表現(xiàn)為升高態(tài)勢[17]?？梢姡瑥?fù)種指數(shù)的提高對于提高農(nóng)田生態(tài)系統(tǒng)的碳匯功能有重要貢獻(xiàn)。

表3　1990-2010年陜西黃土臺塬區(qū)土地利用變化統(tǒng)計

3結(jié)論與討論

陜西黃土臺塬區(qū)植被—土壤系統(tǒng)總體表現(xiàn)為碳源效應(yīng)，只有東北部邊緣區(qū)很小的范圍表現(xiàn)為碳匯效應(yīng)。1990—2010年，咸陽臺塬區(qū)、寶雞臺塬區(qū)和西安臺塬區(qū)的碳源效應(yīng)有所減弱，渭南臺塬區(qū)和銅川臺塬區(qū)的碳源效應(yīng)有所增強(qiáng)。高碳源區(qū)由1990年的38.22%降低為2010年的21.13%，同期中碳源區(qū)面積最大且由57.40%增大為67.71%，低碳塬區(qū)由3.43%提高為7.76%，碳匯區(qū)由0.96%擴(kuò)大為3.40%。

植被—土壤系統(tǒng)的碳收支估算問題，對于指導(dǎo)區(qū)域調(diào)整土地利用方式具有重要作用。本研究中發(fā)現(xiàn)陜西黃土臺塬區(qū)土地利用變化呈現(xiàn)出建設(shè)用地占用耕地、林地、草地的土地利用基本態(tài)勢，這種態(tài)勢不利于碳固定。但復(fù)種指數(shù)的提高對植被—土壤系統(tǒng)碳固定具有明顯的正向作用。另外，植被—土壤系統(tǒng)的碳收支計算存在一定的不確定性，測算精度受土壤系統(tǒng)的RH的影響明顯。本研究中采用了前人研究的RH均值，在一定程度上影響了本文測評結(jié)果的精度。但研究區(qū)為半干旱區(qū)，RH效應(yīng)的空間差異較小，總體上不影響結(jié)果的時空分布特征，后續(xù)研究中擬考慮不同土壤類型的RH對植被—土壤系統(tǒng)的碳收支影響。

[參考文獻(xiàn)]

[1]Kaku K. Global warming and climate change of Asian countries including Japanese domestic greenhouse gas(GHG) reduction in the filed of poultry and swine industries[J]. Procedia Engineering, 2011(8):511-514.

[2]劉春蘭,蔡博峰,陳操操,等.中國碳減排目標(biāo)的地區(qū)分解方法研究述評[J].地理科學(xué),2013,33(9):1089-1096.

[3]Tick Hui Oh. Carbon capture and storage potential in coal-fired plant in Malaysia: A review[J]. Renewable and Sustainable Energy Reviews, 2010,14(9):2697-2709.

[4]方精云,郭兆迪,樸世龍,等.1981—2000年中國陸地植被碳匯的估算[J].中國科學(xué):D輯,2007,37(6):804-812.

[5]高志強(qiáng),劉紀(jì)遠(yuǎn),曹明奎,等.土地利用和氣候變化對農(nóng)牧過渡區(qū)生態(tài)系統(tǒng)生產(chǎn)力和碳循環(huán)的影響[J].中國科學(xué):D輯,2005,34(10):946-957.

[6]Li Wenhua. Progresses and perspectives of ecological research in China[J]. Journal of Resources and Ecology, 2010,1(1):3-14.

[7]葛全勝,戴君虎,何凡能,等.過去300年中國土地利用,土地覆被變化與碳循環(huán)研究[J].中國科學(xué):D輯,2008,38(2):197-210.

[8]Fang J, Chen A, Peng C, et al. Changes in forest biomass carbon storage in China between 1949 and 1998[J]. Science, 2001,292(5525):2320-2322.

[9]張秀梅,李升峰,黃賢金,等.江蘇省1996年至2007年碳排放效應(yīng)及時空格局分析[J].資源科學(xué),2010,32(4):768-775.

[10]于貴瑞,王秋鳳,朱先進(jìn).區(qū)域尺度陸地生態(tài)系統(tǒng)碳收支評估方法及其不確定性[J].地理科學(xué)進(jìn)展,2011,30(1):103-113.

[11]馬彩虹,任志遠(yuǎn),李小燕.黃土臺塬區(qū)土地利用轉(zhuǎn)移流及空間集聚特征分析[J].地理學(xué)報,2013,68(2):257-267.

[12]孫睿,朱啟疆.中國陸地植被凈第一性生產(chǎn)力及季節(jié)變化研究[J].地理學(xué)報,2000,55(1):36-45.

[13]董丹,倪健.利用CASA模型模擬西南喀斯特植被凈第一性生產(chǎn)力[J].生態(tài)學(xué)報,2011,31(7):1855-1866.

[14]湯潔,姜毅,李昭陽,等.基于CASA模型的吉林西部植被凈初級生產(chǎn)力及植被碳匯量估測[J].干旱區(qū)資源與環(huán)境,2013,27(4):1-7.

[15]毛德華,王宗明,羅玲,等.1982—2009年東北多年凍土區(qū)植被凈初級生產(chǎn)力動態(tài)及其對全球變化的響應(yīng)[J].應(yīng)用生態(tài)學(xué)報,2012,23(6):1511-1519.

[16]朱文泉,潘耀忠,張錦水.中國陸地植被凈初級生產(chǎn)力遙感估算[J].植物生態(tài)學(xué)報,2007,31(3):413-424.

[17]李晶,任志遠(yuǎn).基于SPOTNDVI的陜西省耕地復(fù)種指數(shù)時空變化[J].干旱區(qū)資源與環(huán)境,2011,25(10):86-91.

Evaluation of Carbon Source or Sink of Vegetation and Soil System in Shaanxi Loess Platform Region

MA Caihong1,2, REN Zhiyuan1

(1.CollegeofTourismandEnvironmentScience,ShaanxiNormalUniversity,Xi’an,Shaanxi710062,China; 2.QinlingandIntonesGeographyResearchInstitute,SchoolofHistoryandTourism,ShaanxiUniversityofTechnology；Hanzhong,Shaanxi723001,China)

Abstract：[Objective] To provide a decision basis for low-carbon strategy of land use and land cover change. [Methods] Based on RS and GIS, taking the data source of NDVI and Landsat, the carbon of vegetation and soil system in Shaanxi Loess platform region during 1990—2010 was calculated by the Carnegie-Ames-Stanford Approach(CASA) model. [Results] Large part in the research area was carbon source, however, the carbon sink area were small and mainly distributed in the northern part in Weinan loess platform. During the research period, the carbon source degree in Xianyang loess platform, Baoji loess platform and Xi’an loess platform decreased, whereas, of which in Weinan loess platform and Tongchuan platform increased in some degree. From 1990 to 2010, the high degree carbon source area decreased from 38.22% to 21.13%, middle degree carbon source area were the largest and improved from 57.40% to 67.71%, and the area of low degree carbon source increased from 3.43% to 7.76%.At the same period, the area of carbon sink enlarged from 0.96% to 3.40%. [Conclusion] Cultivated land, forest land and grass land were changed into built-up land were not benefit to the carbon sink, and the enhancing of cropping index benefit to the carbon sink of vegetation and soil system.

Keywords:Carnegie-Ames-Stanford Approach(CASA); carbon source; carbon sink; Shaanxi Loess platform region; vegetation and soil system

文獻(xiàn)標(biāo)識碼：B

文章編號：1000-288X(2015)03-0231-04

中圖分類號：X51;F301

通信作者：任志遠(yuǎn)(1953—)，男(漢族)，陜西省興平市人，教授，博導(dǎo)，主要從事國土資源開發(fā)與生態(tài)環(huán)境評價研究。E-mail: renzhy@snnu.edu.cn。

收稿日期：2014-03-27修回日期：2014-04-21

資助項目：國家自然科學(xué)基金“我國生態(tài)脆弱區(qū)能源開發(fā)生態(tài)效應(yīng)測評與調(diào)控研究”(41371523); 國家社會科學(xué)基金(14XKS019); 教育部人文社會科學(xué)重點(diǎn)研究基地項目(14JJD840004)

第一作者：馬彩虹(1974—),女(漢族),寧夏回族自治區(qū)西吉縣人,博士,副教授,主要從事資源環(huán)境評價與GIS/RS應(yīng)用、生態(tài)經(jīng)濟(jì)與區(qū)域可持續(xù)發(fā)展研究。E-mail:mchyanni@aliyun.com。