氣候變化背景下我國陸地生態系統固碳增匯分區研究

陳福軍，郭 英，李 倩

(1.河北省唐山市曹妃甸區國土資源局，河北唐山 063200；2.中國科學院遺傳與發育生物學研究所農業資源研究中心，河北石家莊 050021)

?

氣候變化背景下我國陸地生態系統固碳增匯分區研究

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

(1.河北省唐山市曹妃甸區國土資源局，河北唐山 063200；2.中國科學院遺傳與發育生物學研究所農業資源研究中心，河北石家莊 050021)

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

氣候變化；固碳增匯；NDVI；氣候因子；分區

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

1 材料與方法

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

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

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

1.2 計算方法 為計算分析NDVI對氣溫和降水的響應關系，逐像元計算我國陸地生態系統近30年間NDVI年累積值與標準化的年降水、年均溫柵格數據的相關系數，并對計算結果進行顯著性檢驗，當相關系數通過顯著性水平(P<0.05)時，表示NDVI與氣候因子(氣溫或降水)變化呈顯著線性相關關系。計算公式為：

按照以上計算結果，若某區域植被NDVI與年均溫變化呈顯著正相關，與年降水變化呈顯著正相關，則該區域為Ⅰ 氣候暖濕化驅動型區；若某區域植被NDVI與年均溫變化呈顯著正相關，與年降水變化呈顯著負相關，則該區域為Ⅱ 氣候暖干化驅動型區；若某區域植被NDVI與年均溫變化呈顯著負相關，與年降水變化呈顯著正相關，則該區域為Ⅲ氣候冷濕化驅動型區。

1.3 分區原則及方法 區域劃分主要遵循以下原則：①行政單元完整原則。保持縣級行政單元界線的完整，在分區時沿縣級行政單元的邊界線劃定分區界線。②主導因素原則。在分區時以植被對氣候變化響應顯著為主導因素，若某個行政單元中有幾種顯著響應類型，以面積占優者為準進行劃定。③空間分布連續性原則。在分區過程中要根據分區空間范圍的大小進行適當的取舍，以保持分區結果的完整性。④大尺度地形單元一致性原則。我國緯度跨度較大、地貌階梯格局復雜， 使得不同地域植被類型響應氣候變化的規律受大尺度地形單元的影響，因此在分區中將大尺度地形相對一致的縣級單元劃分在一個氣候變化區內。

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

2 結果與分析

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

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

圖2 1981—2010年我國陸地生態系統固碳增匯分區Fig.2 Carbon sequestration increase zoning of China terrestrial ecosystem during 1981-2010

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

3 結論與討論

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

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

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

(4)以上我國陸地生態系統固碳增匯分區是依托1981—2010年氣候變化背景下植被長勢情況計算得到，為了在該研究成果的基礎上主動適應氣候變化，提高生態系統植被固碳增匯能力，今后需要結合我國氣候變化的基本區域分異，優化生態系統的人為管理。

[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] 趙英時.遙感應用分析原理與方法[M].北京：科學出版社，2004.

[3] 樸世龍，方精云.最近18年來中國植被覆蓋的動態變化[J].第四紀研究，2001，21(4)：294—302.

[4] 宋怡，馬明國.基于GIMMS AVHRR/NDVI數據的中國寒旱區植被動態及其與氣候因子的關系[J].遙感學報，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] 任國玉，張愛英，初子瑩，等.我國地面氣溫參考站點遴選的依據、原則和方法[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- )，男，河北唐山人，工程師， 碩士，從事氣候變化與生態系統碳循環研究。*通訊作者，助理研究員，博士，從事生態水文模型與遙感研究。

2016-08-29

S 181；P 966.7

A

0517-6611(2016)31-0064-03