氣候變化背景下我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)研究

陳福軍，郭 英，李 倩

(1.河北省唐山市曹妃甸區(qū)國(guó)土資源局，河北唐山 063200；2.中國(guó)科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心，河北石家莊 050021)

?

氣候變化背景下我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)研究

陳福軍1，2，郭 英2*，李 倩1

(1.河北省唐山市曹妃甸區(qū)國(guó)土資源局，河北唐山 063200；2.中國(guó)科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心，河北石家莊 050021)

[目的]使我國(guó)生態(tài)系統(tǒng)盡可能地適應(yīng)全球性的氣候變化。[方法]對(duì)陸地植被固碳增匯的區(qū)域差異管理進(jìn)行分區(qū)，運(yùn)用歸一化植被指數(shù)(NDVI)對(duì)氣候變化的響應(yīng)關(guān)系，對(duì)我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)進(jìn)行研究。[結(jié)果]將我國(guó)陸地生態(tài)系統(tǒng)劃分為3個(gè)主動(dòng)適應(yīng)氣候變化固碳增匯區(qū)域。Ⅰ氣候暖濕化驅(qū)動(dòng)型區(qū)，包括華北、華南及西北等大部分區(qū)域；Ⅱ氣候暖干化驅(qū)動(dòng)型區(qū)，主要分布于東北大小興安嶺、長(zhǎng)白山等亞寒帶地區(qū)、青藏高原高寒區(qū)域大部、天山、祁連山及橫斷山脈等高海拔地區(qū)；Ⅲ氣候冷濕化驅(qū)動(dòng)型區(qū)，主要分布于內(nèi)蒙古東部及北部沙漠化嚴(yán)重地區(qū)、長(zhǎng)江流域下游水稻主產(chǎn)區(qū)及西藏林芝熱帶半濕潤(rùn)地區(qū)。[結(jié)論]以縣級(jí)行政界線(xiàn)為基本劃分單元的氣候變化背景下我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)，可為生態(tài)系統(tǒng)區(qū)域差異管理提供參考。

氣候變化；固碳增匯；NDVI；氣候因子；分區(qū)

2013年9月聯(lián)合國(guó)政府間氣候變化專(zhuān)門(mén)委員會(huì)(IPCC)發(fā)布的第五次全球氣候變化評(píng)估報(bào)告顯示，2003—2012年全球海陸表面平均氣溫比1850—1900年上升了0.78 ℃[1]，氣候變化已成為世界各國(guó)高度關(guān)注的全球性問(wèn)題。根據(jù)不同區(qū)域的氣候變化類(lèi)型，以及與之對(duì)應(yīng)的自然界響應(yīng)關(guān)系來(lái)制訂具體的適應(yīng)性措施，是人類(lèi)社會(huì)在全球氣候變化大環(huán)境下尋求可持續(xù)發(fā)展的重要途徑。歸一化植被指數(shù)(NDVI)是監(jiān)測(cè)大尺度生態(tài)系統(tǒng)植被的常用指標(biāo)之一，與植被吸收的光合有效輻射、生物固碳量呈線(xiàn)性關(guān)系[2]。樸世龍等[3]利用AVHRR/NDVI數(shù)據(jù)分析了1982—1999年我國(guó)植被覆蓋動(dòng)態(tài)變化；宋怡等[4]于2008年對(duì)我國(guó)寒區(qū)、旱區(qū)NDVI與氣候因子變化間的相互關(guān)系進(jìn)行了研究，結(jié)果表明我國(guó)寒區(qū)、旱區(qū)植被變化與降水、氣溫呈正相關(guān)關(guān)系。然而，結(jié)合氣候變化類(lèi)型對(duì)我國(guó)全域陸地生態(tài)系統(tǒng)固碳增匯區(qū)域劃分卻鮮見(jiàn)報(bào)道。筆者利用1981—2010年我國(guó)氣溫和降水的變化趨勢(shì)及生態(tài)系統(tǒng)植被與之對(duì)應(yīng)的響應(yīng)關(guān)系，定量識(shí)別植被與氣候的相互作用機(jī)制，結(jié)合我國(guó)地形特點(diǎn)，以縣級(jí)行政單元區(qū)劃為基本劃分單元，研究氣候變化背景下我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)，旨在為生態(tài)系統(tǒng)區(qū)域差異管理提供依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)來(lái)源 采用1981～2006年8 km分辨率的逐旬AVHRR/NDVI數(shù)據(jù)[5]，使用最大值合成法(MVC)求出逐月最大NDVI。2007～2010年遙感數(shù)據(jù)采用MODIS/NDVI產(chǎn)品[6]，AVHRR/NDVI數(shù)據(jù)持續(xù)至2006年，與MODIS衛(wèi)星數(shù)據(jù)在時(shí)間上有6個(gè)重合年份，運(yùn)用2001—2006年AVHRR/NDVI數(shù)據(jù)和MODIS/NDVI數(shù)據(jù)，從2種數(shù)據(jù)分別選取相同時(shí)間和相同地點(diǎn)的NDVI數(shù)值對(duì)多組，應(yīng)用線(xiàn)性回歸方法對(duì)該數(shù)據(jù)組進(jìn)行回歸，利用回歸方程對(duì)2007—2010年MODIS/NDVI數(shù)據(jù)進(jìn)行校正。

氣象數(shù)據(jù)基于1981～2010年我國(guó)境內(nèi)及周邊637個(gè)站點(diǎn)的氣象資料[7]，運(yùn)用Arcmap軟件克里金插值方法進(jìn)行月均溫、月降水量的空間插值，得到分辨率8 km氣象要素的空間柵格數(shù)據(jù)。

地形數(shù)據(jù)來(lái)源于美國(guó)地質(zhì)勘探局(USGS)提供的Global 30 Arc-Second Elevation(GTOPO30)高程數(shù)據(jù)集，數(shù)據(jù)集的空間分辨率重采樣為8 km 分辨率。

1.2 計(jì)算方法 為計(jì)算分析NDVI對(duì)氣溫和降水的響應(yīng)關(guān)系，逐像元計(jì)算我國(guó)陸地生態(tài)系統(tǒng)近30年間NDVI年累積值與標(biāo)準(zhǔn)化的年降水、年均溫柵格數(shù)據(jù)的相關(guān)系數(shù)，并對(duì)計(jì)算結(jié)果進(jìn)行顯著性檢驗(yàn)，當(dāng)相關(guān)系數(shù)通過(guò)顯著性水平(P<0.05)時(shí)，表示NDVI與氣候因子(氣溫或降水)變化呈顯著線(xiàn)性相關(guān)關(guān)系。計(jì)算公式為：

按照以上計(jì)算結(jié)果，若某區(qū)域植被NDVI與年均溫變化呈顯著正相關(guān)，與年降水變化呈顯著正相關(guān)，則該區(qū)域?yàn)棰?氣候暖濕化驅(qū)動(dòng)型區(qū)；若某區(qū)域植被NDVI與年均溫變化呈顯著正相關(guān)，與年降水變化呈顯著負(fù)相關(guān)，則該區(qū)域?yàn)棰?氣候暖干化驅(qū)動(dòng)型區(qū)；若某區(qū)域植被NDVI與年均溫變化呈顯著負(fù)相關(guān)，與年降水變化呈顯著正相關(guān)，則該區(qū)域?yàn)棰髿夂蚶錆窕?qū)動(dòng)型區(qū)。

1.3 分區(qū)原則及方法 區(qū)域劃分主要遵循以下原則：①行政單元完整原則。保持縣級(jí)行政單元界線(xiàn)的完整，在分區(qū)時(shí)沿縣級(jí)行政單元的邊界線(xiàn)劃定分區(qū)界線(xiàn)。②主導(dǎo)因素原則。在分區(qū)時(shí)以植被對(duì)氣候變化響應(yīng)顯著為主導(dǎo)因素，若某個(gè)行政單元中有幾種顯著響應(yīng)類(lèi)型，以面積占優(yōu)者為準(zhǔn)進(jìn)行劃定。③空間分布連續(xù)性原則。在分區(qū)過(guò)程中要根據(jù)分區(qū)空間范圍的大小進(jìn)行適當(dāng)?shù)娜∩幔员３址謪^(qū)結(jié)果的完整性。④大尺度地形單元一致性原則。我國(guó)緯度跨度較大、地貌階梯格局復(fù)雜， 使得不同地域植被類(lèi)型響應(yīng)氣候變化的規(guī)律受大尺度地形單元的影響，因此在分區(qū)中將大尺度地形相對(duì)一致的縣級(jí)單元?jiǎng)澐衷谝粋€(gè)氣候變化區(qū)內(nèi)。

根據(jù)上述分區(qū)原則，采用遙感數(shù)據(jù)提取感興趣區(qū)域(ROI)方法，對(duì)我國(guó)生態(tài)系統(tǒng)植被NDVI與氣候因子顯著響應(yīng)區(qū)域進(jìn)行提取，再采用GIS空間疊置方法，將上述提取的ROI區(qū)域與我國(guó)縣級(jí)行政單元矢量數(shù)據(jù)以及大范圍地形(高程與地貌單元)進(jìn)行疊加，綜合分析確定我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)。

2 結(jié)果與分析

2.1 植被NDVI對(duì)氣候因子響應(yīng)的區(qū)域分異 從圖1可見(jiàn)，年均氣溫升高對(duì)植被(NDVI)生長(zhǎng)不利的地區(qū)包括內(nèi)蒙古東部及北部沙漠化嚴(yán)重地區(qū)(呼倫貝爾沙地、科爾沁沙地、渾沙達(dá)克沙地和巴丹吉林沙漠等)、長(zhǎng)江流域下游水稻主產(chǎn)區(qū)及西藏林芝熱帶半濕潤(rùn)地區(qū)，其余大部分區(qū)域年均氣溫升高對(duì)生態(tài)系統(tǒng)植被生長(zhǎng)有利。年降水量升高對(duì)植被(NDVI)生長(zhǎng)不利的地區(qū)包括我國(guó)東北大小興安嶺、長(zhǎng)白山等亞寒帶地區(qū)、青藏高原高寒區(qū)域大部、天山、祁連山及橫斷山脈等高海拔地區(qū)，其余我國(guó)大部分區(qū)域年降水量升高對(duì)生態(tài)系統(tǒng)植被生長(zhǎng)有利。

注：a為植被NDVI對(duì)年均氣溫的響應(yīng)；b為植被NDVI對(duì)年降水量的響應(yīng)。Note：a response of vegetation NDVI to annual average temperature；b response of vegetation NDVI to annual precipitation.圖1 植被NDVI 對(duì)氣候因子響應(yīng)的區(qū)域分異Fig.1 Response of vegetation NDVI to climatic factors in different regions

圖2 1981—2010年我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)Fig.2 Carbon sequestration increase zoning of China terrestrial ecosystem during 1981-2010

2.2 我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)及分布特征 根據(jù)陸地生態(tài)系統(tǒng)植被長(zhǎng)勢(shì)情況(以NDVI為標(biāo)識(shí)因子)對(duì)氣候變化的響應(yīng)關(guān)系，將我國(guó)陸地劃分為3個(gè)主動(dòng)適應(yīng)氣候變化固碳增匯區(qū)域：Ⅰ 氣候暖濕化驅(qū)動(dòng)型區(qū)、Ⅱ 氣候暖干化驅(qū)動(dòng)型區(qū)、Ⅲ氣候冷濕化驅(qū)動(dòng)型區(qū)(圖2)。其中，Ⅱ 氣候暖干化驅(qū)動(dòng)型區(qū)主要分布于我國(guó)東北大小興安嶺、長(zhǎng)白山等亞寒帶地區(qū)、青藏高原高寒區(qū)域大部、天山、祁連山及橫斷山脈等高海拔地區(qū)，該區(qū)域若年均氣溫升高、年降水量減少對(duì)陸地生態(tài)系統(tǒng)植被生長(zhǎng)有利，相反則對(duì)植被生長(zhǎng)不利；Ⅲ 氣候冷濕化驅(qū)動(dòng)型區(qū)主要分布于內(nèi)蒙古東部及北部沙漠化嚴(yán)重地區(qū)、長(zhǎng)江流域下游水稻主產(chǎn)區(qū)及西藏林芝熱帶半濕潤(rùn)地區(qū)，該區(qū)域若年均氣溫降低、年降水量增多對(duì)陸地生態(tài)系統(tǒng)植被生長(zhǎng)有利，相反則對(duì)植被生長(zhǎng)不利；其余大部分區(qū)域均屬于Ⅰ 氣候暖濕化驅(qū)動(dòng)型區(qū)，該區(qū)域年均氣溫升高，同時(shí)年降水量增多對(duì)陸地生態(tài)系統(tǒng)植被生長(zhǎng)有利，相反則對(duì)植被生長(zhǎng)不利。

3 結(jié)論與討論

(1)該研究依托的空間數(shù)據(jù)有植被NDVI數(shù)據(jù)、年均氣溫?cái)?shù)據(jù)和年降水量數(shù)據(jù)，其中NDVI數(shù)據(jù)來(lái)自2種不同渠道，運(yùn)用線(xiàn)性回歸方法對(duì)2007～2010年MODIS/NDVI數(shù)據(jù)進(jìn)行校正，可以避免數(shù)據(jù)在時(shí)間尺度上的突變。氣象數(shù)據(jù)依據(jù)我國(guó)境內(nèi)及周邊637個(gè)站點(diǎn)的氣象資料進(jìn)行空間插值取得，因我國(guó)氣象站點(diǎn)東西部分布不均、坐落地形或與城市距離不同、插值方法不同等，空間插值結(jié)果可能存在誤差，今后應(yīng)選定一批具有代表性、分布均勻的氣象站作為參考站，對(duì)插值結(jié)果進(jìn)行驗(yàn)證[8]。

(2)固碳增匯分區(qū)選取縣級(jí)行政區(qū)為基本單元，目的是盡量詳實(shí)細(xì)致地對(duì)不同類(lèi)型的區(qū)域進(jìn)行分類(lèi)，同時(shí)也考慮到目前縣級(jí)區(qū)域是我國(guó)行政管理的基本單元。而實(shí)際分類(lèi)過(guò)程中，我國(guó)西部省份有些縣域面積較大，或某些縣域包括了幾種不同植被NDVI對(duì)氣候變化的響應(yīng)類(lèi)型。這種情況只能按照主導(dǎo)因素原則將面積較大的類(lèi)型作為該縣域固碳增匯的基礎(chǔ)類(lèi)型進(jìn)行劃分。

(3)該研究應(yīng)用植被NDVI對(duì)氣候變化的響應(yīng)關(guān)系，將我國(guó)陸地劃分為3個(gè)主動(dòng)適應(yīng)氣候變化固碳增匯區(qū)域：Ⅰ氣候暖濕化驅(qū)動(dòng)型區(qū)，包括除Ⅱ、Ⅲ類(lèi)型以外我國(guó)大部分區(qū)域；Ⅱ氣候暖干化驅(qū)動(dòng)型區(qū)，主要分布于我國(guó)東北大小興安嶺、長(zhǎng)白山等亞寒帶地區(qū)、青藏高原高寒區(qū)域大部、天山、祁連山及橫斷山脈等高海拔地區(qū)；Ⅲ氣候冷濕化驅(qū)動(dòng)型區(qū)，主要分布于內(nèi)蒙古東部及北部沙漠化嚴(yán)重地區(qū)、長(zhǎng)江流域下游水稻主產(chǎn)區(qū)及西藏林芝熱帶半濕潤(rùn)地區(qū)。

(4)以上我國(guó)陸地生態(tài)系統(tǒng)固碳增匯分區(qū)是依托1981—2010年氣候變化背景下植被長(zhǎng)勢(shì)情況計(jì)算得到，為了在該研究成果的基礎(chǔ)上主動(dòng)適應(yīng)氣候變化，提高生態(tài)系統(tǒng)植被固碳增匯能力，今后需要結(jié)合我國(guó)氣候變化的基本區(qū)域分異，優(yōu)化生態(tài)系統(tǒng)的人為管理。

[1] IPCC.Climate change 2013：The physical science basis[M]//Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change.cambridge：Cambridge University Press，2013.

[2] 趙英時(shí).遙感應(yīng)用分析原理與方法[M].北京：科學(xué)出版社，2004.

[3] 樸世龍，方精云.最近18年來(lái)中國(guó)植被覆蓋的動(dòng)態(tài)變化[J].第四紀(jì)研究，2001，21(4)：294—302.

[4] 宋怡，馬明國(guó).基于GIMMS AVHRR/NDVI數(shù)據(jù)的中國(guó)寒旱區(qū)植被動(dòng)態(tài)及其與氣候因子的關(guān)系[J].遙感學(xué)報(bào)，2008，12(3)：500-505.

[5] TUCKER C J，PINZON J E，BROWN M E.Global inventory modeling and mapping studies，NA94apr15b.n11-VIg，2.0[M].Maryland：Global Land Cover Facility University of Maryland，College Park，2004.

[6] Land Processes Distributed Active Archive Center.ASTER Level 1B Registered Radiance at the Sensor.Version 3.NASA EOSDIS Land Processes DAAC，USGS Earth Resources Observation and Science(EROS)Center，Sioux Falls，South Dakota(https：//lpdaac.usgs.gov)，accessed January 1，2015，at http：//dx.doi.org/10.5067/ASTER/AST_L1B.003.

[7] MENNE M J，DURRE I，VOSE R S，et al.An overview of the global historical climatology network-daily database[J]. Journal of atmospheric and oceanic technology，2012，29(7)：897-910.

[8] 任國(guó)玉，張愛(ài)英，初子瑩，等.我國(guó)地面氣溫參考站點(diǎn)遴選的依據(jù)、原則和方法[J].氣象科技，2010，38(1)：78-85.

Study on Carbon Sequestration Increase Zoning of China Terrestrial Ecosystem under Climate Change

CHEN Fu-jun1,2, GUO Ying2*, LI Qian1

(1. Tangshan Caofeidian Land Resources Bureau, Tangshan, Hebei 063200; 2. Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang, Hebei 050021)

[Objective] The aim was to provide theoretical basis for the ecological system in China to adapt to global climate change as much as possible. [Method] The response of vegetation NDVI to climate change was applied to study on carbon sequestration increase zoning of China terrestrial ecosystem. [Result] The terrestrial ecosystem in China was divided into three regions: Ⅰ warm and humid climate driven type region, including most areas of China, as North China, South China and the northwest of China, Ⅱ warm and dry climate driven type region, mainly distributed in Daxing’an and Xiaoxing’an Mountains, Changbai Mountain in the subarctic region of Northeast China, high altitude areas, as Qinghai Tibet Plateau alpine region, Tianshan Mountains, Qilian Mountains and Hengduan Mountains, Ⅲ cold and humid climate driven region, mainly distributed in the eastern part and the northern desert region of Inner Mongolia, the main produce areas of rice in the lower reaches of the Yangtze River and the tropical semi humid regions of Linzhi, Tibet. [Conclusion] Carbon sequestration increase zoning of China terrestrial ecosystem under climate change with county administrative boundaries as the basic unit can provide reference for regional difference management of ecosystem.

Climate change; Carbon sequestration; NDVI; Climatic factor; Zoning

陳福軍(1983- )，男，河北唐山人，工程師， 碩士，從事氣候變化與生態(tài)系統(tǒng)碳循環(huán)研究。*通訊作者，助理研究員，博士，從事生態(tài)水文模型與遙感研究。

2016-08-29

S 181；P 966.7

A

0517-6611(2016)31-0064-03