陳福軍 張玉翠 李倩

摘要 [目的]分析我國不同區(qū)域氣候變化類型對當(dāng)?shù)厣鷳B(tài)系統(tǒng)植被長勢的長期影響。[方法]利用1981—2010年中國陸地生態(tài)系統(tǒng)植被指數(shù)NDVI與氣溫及降水之間的響應(yīng)關(guān)系，定量識別植被與氣候的相互作用機(jī)制，結(jié)合中國氣候變化區(qū)劃，完成氣候變化背景下的中國陸地植被覆蓋度預(yù)測。[結(jié)果]我國東北大小興安嶺、長白山、云貴高原等地區(qū)植被更適應(yīng)當(dāng)?shù)貧夂蚺苫厔?，西北地區(qū)大部、東南地區(qū)（長江下游除外）植被更適應(yīng)當(dāng)?shù)貧夂蚺瘽窕厔荩瑸闅夂蜃兓谥脖簧L區(qū)；我國內(nèi)蒙古東部及北部沙漠化嚴(yán)重地區(qū)植被不適應(yīng)當(dāng)?shù)貧夂蚺苫厔?，為氣候變化不利于植被生長區(qū)。我國其他大部分區(qū)域植被長勢與氣候變化無顯著響應(yīng)關(guān)系。[結(jié)果]該研究成果可為我國不同區(qū)域生態(tài)系統(tǒng)的差異化管理提供參考。

關(guān)鍵詞 氣候變化；陸地植被；NDVI；覆蓋度；適應(yīng)性

中圖分類號 S162；Q948文獻(xiàn)標(biāo)識碼 A文章編號 0517-6611（2018）01-0148-04

Abstract [Objective] The research aimed to quantitatively analyze the long-term effects for different types of regional climate change to the local ecosystem vegetation growth in China.[Method]Using the response relationship between Chinas terrestrial ecosystems vegetation index NDVI and the temperature，precipitation from 1981 to 2010，the interaction mechanism between vegetation and climate was quantitatively identified.Combined with Chinas climate change zoning，Chinas terrestrial vegetation coverage forecast under the background of climate change was completed.[Result]The vegetation of Xinganling Mountains，Changbai Mountain，Yunnan-guizhou Plateau regions was more adapt to the local warm and dry climatic change type，the vegetation of most northwest region，southeast region of China （except the lower reaches of Yangtze River） was more adapt to the local wet and warm climate change type，which were areas that climate change was favorable for vegetation growth.The vegetation in the serious desertification areas of eastern and north Inner Mongolia was not adapt to the local warm and dry climatic change type，which were areas that climate change was not beneficial to vegetation growth.At the other areas of China，there was no significant response relationship between vegetation growth and climate change.[Conclusion]The research results can provide a reference for the different regional ecological systems differentiation management.

Key words Climate change； Terrestrial vegetation； NDVI； Coverage；Adaptability

2003—2012年全球海陸表面平均氣溫比1850—1900年上升了0.8 ℃，有模型預(yù)測到21世紀(jì)末，全球地表氣溫度將升高2.0～4.8 ℃，不斷變化的氣候可導(dǎo)致前所未有的極端天氣和氣候事件更加頻繁的發(fā)生[1]。陸地以及陸地植被是人類賴以生存的環(huán)境主體，受到氣候變化的沖擊，植被覆蓋狀況已呈現(xiàn)了區(qū)域差異的變化[2]，在這種情況下，預(yù)測我國不同地區(qū)陸地植被的年際生長變化趨勢，對人們更好地適應(yīng)和應(yīng)對氣候變化尤為重要。丁勇等[3]認(rèn)為我國溫帶草原區(qū)呈現(xiàn)總體氣溫明顯升高，降水區(qū)域差異顯著的變化趨勢，受此影響，植物物候、覆蓋度等發(fā)生了明顯變化。趙茂盛等[4]應(yīng)用改進(jìn)的MAPSS模擬顯示，未來氣候變化可能導(dǎo)致我國東部森林植被帶的北移，西部的沙漠和草原可能略有退縮，相應(yīng)被草原和灌叢取代。牛建明[5]對內(nèi)蒙古草原的研究認(rèn)為，草原面積將顯著減少，南部界限大幅度北移，草原生產(chǎn)力明顯下降等。以上研究大多從物候?qū)W或氣象學(xué)的角度展開，筆者將從地理學(xué)分區(qū)的角度，利用中國1981—2010年氣溫和降水的變化趨勢以及生態(tài)系統(tǒng)植被與之對應(yīng)的響應(yīng)關(guān)系，定量識別植被與氣候的相互作用機(jī)制，結(jié)合中國氣候變化區(qū)劃，完成氣候變化背景下的中國陸地植被覆蓋度預(yù)測，為生態(tài)系統(tǒng)區(qū)域差異管理提供依據(jù)。

1 資料與方法

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

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

氣象數(shù)據(jù)基于1981—2010年中國境內(nèi)及周邊637個站點的氣象資料[8]，應(yīng)用Arcmap軟件克里金插值方法進(jìn)行月均溫、月降水量的空間插值，得到分辨率8 km氣象要素的空間柵格數(shù)據(jù)。地形數(shù)據(jù)來源于美國地質(zhì)勘探局（USGS）提供的Global 30 Arc-Second Elevation （GTOPO30）高程數(shù)據(jù)集，數(shù)據(jù)集的空間分辨率重采樣為8 km 分辨率。中國氣候變化區(qū)劃數(shù)據(jù)引自史培軍等[9]利用1961—2010年氣溫和降水量的變化趨勢值、波動特征值，結(jié)合中國地形特點，以縣級行政區(qū)劃為單元，完成的1961—2010年中國氣候變化區(qū)劃數(shù)據(jù)（圖1）。

1.2 計算方法

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

r=ni=1（xi-）（yi-）

ni=1（xi-）2ni=1（yi-）2（1）

其中，yi表示氣候因子如年降水量、年均溫第i年的某個柵格數(shù)值，代表氣候因子多年均值，xi表示在第i年 NDVI的某個柵格數(shù)值， 表示NDVI多年平均值。當(dāng)相關(guān)系數(shù)通過顯著性水平（ P <0.05）時，表示NDVI與氣候因子變化呈顯著線性相關(guān)關(guān)系。

1.3 預(yù)測方法 首先，對陸地植被長勢驅(qū)動類型進(jìn)行區(qū)域劃分，采用遙感數(shù)據(jù)提取感興趣區(qū)域（ROI）方法，對我國生態(tài)系統(tǒng)植被NDVI與氣候因子響應(yīng)顯著區(qū)域進(jìn)行提取，若某區(qū)域植被NDVI與過去30年均溫變化呈顯著正相關(guān)、與年降水變化呈顯著正相關(guān)，則該區(qū)域為氣候暖濕化驅(qū)動型區(qū)；若某區(qū)域植被NDVI與過去30年均溫變化呈顯著正相關(guān)、與年降水變化呈顯著負(fù)相關(guān)，則該區(qū)域為氣候暖干化驅(qū)動型區(qū)；若某區(qū)域植被NDVI與過去30年均溫變化呈顯著負(fù)相關(guān)、與年降水變化呈顯著正相關(guān)，則該區(qū)域為氣候冷濕化驅(qū)動型區(qū)[10]。

其次，將陸地植被長勢驅(qū)動區(qū)域劃分與中國氣候變化區(qū)劃數(shù)據(jù)進(jìn)行GIS空間疊加，提取我國陸地范圍內(nèi)：①植被適應(yīng)氣候變化的區(qū)域，即植被覆蓋度上升區(qū)域；②植被不適應(yīng)氣候變化的區(qū)域，即植被覆蓋度下降區(qū)域。

1.4 區(qū)域劃分原則

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

2 中國陸地植被覆蓋變化趨勢預(yù)測

2.1 植被NDVI對氣候因子響應(yīng)的區(qū)域分異

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

2.2 中國陸地植被長勢驅(qū)動區(qū)域劃分

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

2.3 中國陸地植被覆蓋度變化趨勢預(yù)測

將以上陸地植被長勢驅(qū)動區(qū)域劃分?jǐn)?shù)據(jù)與中國氣候變化區(qū)劃數(shù)據(jù)進(jìn)行GIS空間疊加（圖4），提取我國陸地范圍內(nèi)2種植被覆蓋變化明顯區(qū)域：①植被適應(yīng)氣候變化的區(qū)域，即植被覆蓋度將上升區(qū)域，我國東北大小興安嶺、長白山、云貴等地區(qū)為氣候暖干化驅(qū)動區(qū)，疊加氣候暖干趨勢帶；西北地區(qū)大部、東南地區(qū)（長江下游除外）為氣候暖濕化驅(qū)動區(qū)，疊加氣候暖濕趨勢帶，氣候變化趨勢對陸地植被生長有利。②植被不適應(yīng)氣候變化的區(qū)域，即植被覆蓋度下降區(qū)域。內(nèi)蒙古東部及北部沙漠化嚴(yán)重地區(qū)為氣候冷濕化驅(qū)動區(qū)，疊加氣候暖干化趨勢帶，氣候變化對陸地植被生長不利。圖4中空白區(qū)域為氣候變化背景下植被覆蓋度預(yù)測變化不顯著區(qū)域。

3 結(jié)論

該研究首先應(yīng)用植被NDVI對氣候變化的響應(yīng)關(guān)系，將中國陸地劃分為3個陸地植被長勢驅(qū)動類型區(qū)域，其次將這3個陸地植被長勢驅(qū)動區(qū)域與中國氣候變化區(qū)劃數(shù)據(jù)進(jìn)行GIS空間疊加，提取我國陸地范圍內(nèi)2種植被覆蓋變化明顯區(qū)域。基于以上分析，預(yù)測我國陸地植被覆蓋度將上升的地區(qū)包括東北大小興安嶺、長白山、云貴等地區(qū)，以及西北地區(qū)大部、東南地區(qū)（長江下游除外），氣候變化趨勢對陸地植被生長有利；預(yù)測我國陸地植被覆蓋度將下降的地區(qū)包括內(nèi)蒙古東部及北部沙漠化嚴(yán)重地區(qū)，氣候變化對陸地植被生長不利。

以上分析并未對混合情況進(jìn)行深入探討，如陸地植被氣候暖干化驅(qū)動區(qū)域疊加氣候暖濕趨勢帶，或陸地植被暖濕化驅(qū)動區(qū)域疊加氣候暖干趨勢帶等，該研究對這種區(qū)域定性分析為植被覆蓋度預(yù)測變化不顯著區(qū)域，需進(jìn)一步詳細(xì)分析。

參考文獻(xiàn)

[1] IPCC.Climate change 2013：the physical science basis.Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[M].Cambridge：Cambridge University Press，2013.

[2] 于貴瑞.人類活動與生態(tài)系統(tǒng)變化的前沿科學(xué)問題[M].北京：高等教育出版社，2009：16-21.

[3] 丁勇，侯向陽，UBUGUNOV L，等.溫帶草原區(qū)氣候變化及其對植被影響的研究進(jìn)展[J].中國農(nóng)學(xué)通報，2012，28（17）：310-316.

[4] 趙茂盛，NEILSON R P，延曉冬，等.氣候變化對中國植被可能影響的模擬[J].地理學(xué)報，2002，57（1）：28-38.

[5] 牛建明.氣候變化對內(nèi)蒙古草原分布和生產(chǎn)力影響的預(yù)測研究[J].草地學(xué)報，2001，9（4）：276-282.

[6] TUCKER C J，PINZON J E，BROWN M E.Global Inventory Modeling and Mapping Studies[M].Mary land：Global Land Cover Facility，University of Maryland，2004.

[7] Land Processes Distributed Active Archive Center （LP DAAC）.ASTER Level 1B registered radiance at the Sensor[M].Version 3.Sioux Falls，South Dakota：NASA EOSDIS Land Processes DAAC，USGS Earth Resources Observation and Science （EROS） Center，2001.

[8] 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：897-910.

[9] 史培軍，孫劭，汪明，等.中國氣候變化區(qū)劃（1961～2010年）[J].中國科學(xué)（地球科學(xué)），2014，44（10）：2294–2306.

[10] 陳福軍，沈彥俊，胡喬利，等.海河流域NDVI對氣候變化的響應(yīng)研究[J].遙感學(xué)報，2011，15（2）：401-414.