郭金停，胡遠(yuǎn)滿，布仁倉(cāng)

（1.內(nèi)蒙古師范大學(xué) 生命科學(xué)與技術(shù)學(xué)院，內(nèi)蒙古 呼和浩特010022；2.中國(guó)科學(xué)院 沈陽(yáng)應(yīng)用生態(tài)研究所，遼寧 沈陽(yáng)110016）

多年凍土是冰凍圈的重要組成部分，近年來隨著全球氣候變暖，發(fā)生了不同程度的退化［1－3］。Taylor 等［4］發(fā)現(xiàn)北美北部多年凍土區(qū)的地表溫度隨著不斷升高的氣溫和降雪覆蓋狀況的改變呈增加的趨勢(shì)，導(dǎo)致北部廣泛分布的多年凍土不斷發(fā)生退化。顧鐘煒和周幼吾［5］對(duì)大興安嶺北坡阿爾木地區(qū)的研究發(fā)現(xiàn)，該區(qū)多年凍土退化趨勢(shì)明顯。由氣候變化引起的多年凍土退化過程對(duì)寒區(qū)植被的生長(zhǎng)、寒區(qū)水文特征以及整個(gè)寒區(qū)陸地生態(tài)系統(tǒng)都有著非常重要的影響。多年凍土退化導(dǎo)致凍土區(qū)水文模式、土壤中營(yíng)養(yǎng)物質(zhì)含量、土壤微生物活性、土壤質(zhì)地和組成發(fā)生改變，進(jìn)而影響地上植被［6］。郭正剛等［7］的研究表明隨著青藏高原北部地區(qū)多年凍土退化，植被覆蓋度降低，群落以及物種多樣性呈先增加后降低的趨勢(shì)。梁四海等［8］也發(fā)現(xiàn)黃河源區(qū)多年凍土發(fā)生退化后會(huì)使土壤含水量減少，造成植被物種更替。王根旭等［9］研究表明，隨著青藏高原多年凍土活動(dòng)層厚度的增加，高寒草甸草地的植被覆蓋度和生物量均顯著降低，同時(shí)高寒草甸與高寒沼澤草甸生態(tài)系統(tǒng)發(fā)生退化。譚俊和李秀華［10］預(yù)測(cè)當(dāng)年平均氣溫增加4 ℃，降水量增加10% 時(shí)，東北東部各森林地帶將向北移3～5 個(gè)緯度，屆時(shí)大興安嶺森林將有可能完全北移出境，取而代之的是以中溫性的針闊混交林以及草原為主的群落。孫廣友［11］認(rèn)為多年凍土的退化或者完全消失破壞了東北地區(qū)沼澤發(fā)生的物質(zhì)基礎(chǔ)，從而導(dǎo)致沼澤退化或者消失。毛德華等［12］對(duì)東北凍土區(qū)植被NDVI 對(duì)氣候變化的響應(yīng)進(jìn)行了研究，雖然探討了該區(qū)植被NDVI 與氣候因子的關(guān)系，但未對(duì)多年凍土退化對(duì)植被的影響進(jìn)行分析。

目前有關(guān)凍土退化的研究大都集中在分析多年凍土的分布以及退化成因等方面［13－16］，很少涉及區(qū)域尺度上關(guān)于多年凍土退化對(duì)植被NDVI 的影響研究，因此在氣候變化背景下，探討多年凍土退化對(duì)植被NDVI 的影響具有重要意義。

1 研究地區(qū)與研究方法

1.1 研究區(qū)概況

我國(guó)東北多年凍土區(qū)（46°30'～53°30' N，115°52'～135°09' E）處于歐亞大陸多年凍土帶的南緣，同時(shí)也是我國(guó)第二大多年凍土分布區(qū)，地跨內(nèi)蒙古自治區(qū)和黑龍江省的最北部。 依據(jù)國(guó)內(nèi)學(xué)者郭東信、金會(huì)軍等［17－18］的研究成果，經(jīng)過掃描數(shù)字化以及屬性添加，確定本研究區(qū)多年凍土南界以及不同多年凍土類型區(qū)范圍。研究區(qū)內(nèi)包括郭東信等［17］在20 世紀(jì)70 年代所考察出的多年凍土南界以及金會(huì)軍等［18］在21 世紀(jì)初期推斷出的南界范圍。按照影響多年凍土空間分布的地理以及氣候條件，將東北多年凍土區(qū)劃分為連續(xù)多年凍土區(qū)（6.2×104km2）、不連續(xù)多年凍土區(qū)（6.6×104km2）、稀疏島狀多年凍土區(qū)（1.3×105km2）以及多年凍土完全退化區(qū)（1.3×105km2）。該區(qū)屬于寒溫帶大陸性氣候，冬季漫長(zhǎng)干冷，夏季短暫濕熱。研究區(qū)的年平均氣溫為－5 ℃～2 ℃，年降水量為260～600 mm［19］。主要植被類型為森林，且大都分布在大小興安嶺地區(qū)，樹種以興安落葉松（Larix gmelinii Kuzen.）、白樺（Betula platyphylla Suk.）、蒙古櫟（Quercus mongolica Fisch.ex Ledeb.）等混交樹種為主。

1.2 數(shù)據(jù)來源與處理

獲取LTDR（land term data record）NDVI 的AVH13C1 以及MODIS NDVI 的MOD13C2 兩種NDVI產(chǎn)品數(shù)據(jù)集進(jìn)行研究，其中AVH13C1 產(chǎn)品的時(shí)間分辨率為1 d，空間分辨率為0.05°，本文使用該數(shù)據(jù)集的時(shí) 間 為1982－1999 年，該 產(chǎn) 品 獲 取 網(wǎng) 址http：//ltdr.nascom.nasa.gov/cgi－bin/ltdr/ltdrPage.cgi。MOD13C2 產(chǎn)品屬于半月合成數(shù)據(jù)，空間分辨率0.05°，使用該數(shù)據(jù)集時(shí)間為2000－2014 年，該產(chǎn)品獲取網(wǎng)址https：//ladsweb.nascom.nasa.gov/。上述兩種NDVI 產(chǎn)品均已經(jīng)過遙感圖像預(yù)處理，如幾何校正、輻射校正、大氣校正等［20］，通過采用最大合成法（MVC）［21］來減少云、大氣、太陽(yáng)高度角等影響，且兩個(gè)NDVI 產(chǎn)品數(shù)據(jù)具有良好的一致性，可以用于后續(xù)研究［22］。受東北多年凍土區(qū)冬季降雪的影響，本研究分別選擇植被春季（4 月－5 月）、夏季（6 月－8 月）、秋季（9 月－10 月）平均NDVI 進(jìn)行研究，且NDVI 值小于0.05 的像元被定義為非植被區(qū)，不參與本文的結(jié)果分析。東北多年凍土區(qū)范圍內(nèi)的35 個(gè)氣象站點(diǎn)1982－2014 年的月平均氣溫、月降水量以及月平均地溫?cái)?shù)據(jù)來源于中國(guó)氣象數(shù)據(jù)網(wǎng)（http：//cdc.cma.gov.cn/）。

1.3 研究方法

利用“空間替代時(shí)間”方法，從連續(xù)多年凍土區(qū)－不連續(xù)多年凍土區(qū)－稀疏島狀多年凍土區(qū)－多年凍土完全退化區(qū)視為多年凍土區(qū)退化過程，由此分析多年凍土退化對(duì)植被產(chǎn)生的影響。

2 結(jié)果與分析

2.1 多年凍土退化對(duì)春季NDVI 影響

2.1.1 不同類型多年凍土區(qū)春季NDVI 變化 四種類型區(qū)多年凍土植被春季NDVI 均呈增加趨勢(shì)，除多年凍土完全退化區(qū)外均達(dá)到了顯著增加水平（圖1）。不同多年的凍土類型區(qū)春季NDVI 增加的幅度不同，連續(xù)多年凍土區(qū)植被春季NDVI 增加的幅度最大，為0.002 a－1，而多年凍土完全退化區(qū)增加幅度最小，為0.001 a－1。不連續(xù)多年凍土區(qū)植被NDVI 增加的幅度與連續(xù)多年凍土區(qū)相近，為0.002 a－1，稀疏島狀多年凍土區(qū)植被春季NDVI 增加幅度為0.001 a－1。

圖1 1982－2014 年?yáng)|北不同類型多年凍土區(qū)春季平均NDVI 年際變化趨勢(shì)Fig.1 Interannual variation of spatial average spring average NDVI values of variable types of permafrost zones in northeastern China during 1982－2014

NDVI 增加趨勢(shì)最大值（≥0.004）所占的像元比例依次為連續(xù)多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞稀疏島狀多年凍土區(qū)＞多年凍土完全退化區(qū)；NDVI 增加趨勢(shì)（0.002～＜0.004）所占的像元比例依次為連續(xù)多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞稀疏島狀多年凍土區(qū)＞多年凍土完全退化區(qū)；NDVI 增加趨勢(shì)（0～＜0.002）所占的像元比例依次為稀疏島狀多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞連續(xù)多年凍土區(qū)＞多年凍土完全退化區(qū)；而NDVI 減少（＜0）所占的像元比例依次為多年凍土完全退化區(qū)＞稀疏島狀多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞連續(xù)多年凍土區(qū)（表1）。

表1 1982－2014 年?yáng)|北不同多年凍土分區(qū)春季NDVI 變化趨勢(shì)的面積比例Tab.1 Proportions of pixels for the variations trend of spring NDVI in the different permafrost zones of northeastern China during 1982－2014 %

2.1.2 不同類型多年凍土區(qū)春季NDVI 與氣候因子相關(guān)性 不同類型多年凍土區(qū)春季NDVI 與春季平均氣溫呈顯著正相關(guān)關(guān)系（見表2），其中稀疏島狀多年凍土區(qū)以及多年凍土完全退化區(qū)達(dá)到了極顯著水平（P＜0.01），連續(xù)多年凍土區(qū)與島狀多年凍土區(qū)植被春季NDVI 與春季降水呈負(fù)相關(guān)關(guān)系，而稀疏島狀多年凍土區(qū)與多年凍土完全退化區(qū)植被NDVI 與春季降水呈正相關(guān)關(guān)系。NDVI 與氣溫相關(guān)性系數(shù)依次為多年凍土完全退化區(qū)＞稀疏島狀多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞連續(xù)多年凍土區(qū)。

表2 1982—2014 年?yáng)|北不同類型多年凍土區(qū)植被春季平均NDVI 與氣候因子的相關(guān)系數(shù)Tab.2 Correlation coefficients between mean spring NDVI and climatic factors during growing season in different permafrost zones of northeastern China during 1982－2014

2.2 多年凍土退化對(duì)夏季NDVI 影響

2.2.1 不同類型多年凍土區(qū)夏季NDVI 變化 四種類型區(qū)植被NDVI 均呈極顯著增加趨勢(shì)（P＜0.01），然而增加的幅度不盡相同（圖2）。連續(xù)多年凍土區(qū)植被夏季NDVI 增加的幅度最大，為0.005 a－1，而多年凍土完全退化區(qū)增加幅度最小，為0.002 a－1。不連續(xù)多年凍土區(qū)為0.004 a－1，稀疏島狀多年凍土區(qū)植被夏季NDVI 增加趨勢(shì)為0.003 a－1。

圖2 1982－2014 年?yáng)|北不同類型多年凍土區(qū)夏季平均NDVI 年際變化趨勢(shì)Fig.2 Interannual variation of spatial average summer average NDVI values of variable types of permafrost zones in northeastern China during 1982－2014

空間像元尺度對(duì)不同多年凍土區(qū)夏季NDVI 變化幅度所占的像元比例見表3，NDVI 增加趨勢(shì)最大值（＞0.004）所占的像元比例依次為連續(xù)多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞稀疏島狀多年凍土區(qū)＞多年凍土完全退化區(qū)；NDVI 增加趨勢(shì)（0.002～＜0.004）所占的像元比例依次為稀疏島狀多年凍土區(qū)＞多年凍土完全退化區(qū)＞不連續(xù)多年凍土區(qū)＞連續(xù)多年凍土區(qū)；NDVI 增加趨勢(shì)（0～＜0.002）所占的像元比例依次為多年凍土完全退化區(qū)＞稀疏島狀多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞連續(xù)多年凍土區(qū)；NDVI 減少（＜0）所占的像元比例依次為多年凍土完全退化區(qū)＞稀疏島狀多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞連續(xù)多年凍土區(qū)。

表3 1982－2014 年?yáng)|北不同多年凍土分區(qū)夏季NDVI 變化趨勢(shì)的面積比例Tab.3 Proportions of pixels for the variations trend of summer NDVI in the different permafrost zones of northeastern China during 1982－2014 %

2.2.2 不同類型多年凍土區(qū)夏季NDVI 與氣候因子相關(guān)性 不同類型多年凍土區(qū)夏季NDVI 與平均氣溫呈顯著正相關(guān)關(guān)系（表4），其中不連續(xù)多年凍土區(qū)、稀疏島狀多年凍土區(qū)以及多年凍土完全退化區(qū)達(dá)到了極顯著水平（P＜0.01）。除多年凍土完全退化區(qū)外，夏季NDVI 與降水呈負(fù)相關(guān)關(guān)系，其中不連續(xù)多年凍土區(qū)達(dá)到了顯著性水平（P＜0.05）。NDVI 與氣溫相關(guān)性系數(shù)依次為不連續(xù)多年凍土區(qū)＞稀疏島狀多年凍土區(qū)＞多年凍土完全退化區(qū)＞連續(xù)多年凍土區(qū)。多年凍土由連續(xù)多年凍土區(qū)退化成不連續(xù)多年凍土區(qū)，NDVI 與氣溫之間的相關(guān)性系數(shù)增加，即氣溫升高可以促進(jìn)夏季NDVI 的增加，但隨著多年凍土區(qū)由不連續(xù)多年凍土區(qū)退化成稀疏島狀多年凍土區(qū)，甚至退化為多年凍土完全退化區(qū)，植被NDVI 與氣溫之間的相關(guān)性系數(shù)均逐漸降低。

表4 1982－2014 年?yáng)|北不同類型多年凍土區(qū)植被夏季平均NDVI 與氣候因子的相關(guān)系數(shù)Tab.4 Correlation coefficients between mean summer NDVI and climatic factors during growing season in different permafrost zones of northeastern China during 1982－2014

2.3 多年凍土退化對(duì)秋季NDVI 影響

2.3.1 不同類型多年凍土區(qū)秋季NDVI 變化 四種類型區(qū)植被秋季NDVI 均呈增加趨勢(shì)（圖3），其中稀疏島狀多年凍土區(qū)與多年凍土完全退化區(qū)未達(dá)到顯著性水平。不同多年凍土類型區(qū)秋季NDVI 增加的幅度不盡相同，連續(xù)多年凍土區(qū)植被秋季NDVI 增加的幅度最大，為0.002 a－1，而多年凍土完全退化區(qū)增加幅度最小，為0.0002 a－1。不連續(xù)多年凍土區(qū)植被NDVI 增加的幅度與連續(xù)多年凍土區(qū)相近，為0.002 a－1，稀疏島狀多年凍土區(qū)植被秋季NDVI 增加趨勢(shì)為0.001 a－1。

圖3 1982－2014 年?yáng)|北不同類型多年凍土區(qū)秋季平均NDVI 年際變化趨勢(shì)Fig.3 Interannual variation of spatial average autumn average NDVI values of variable types of permafrost zones in northeastern China during 1982－2014

空間像元尺度不同對(duì)多年凍土區(qū)秋季NDVI 變化幅度所占的像元比例見表5，NDVI 增加趨勢(shì)最大值（≥0.004）所占的像元比例依次為連續(xù)多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞稀疏島狀多年凍土區(qū)＞多年凍土完全退化區(qū)；NDVI 增加趨勢(shì)（0.002～＜0.004）所占的像元比例依次為連續(xù)多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞稀疏島狀多年凍土區(qū)＞多年凍土完全退化區(qū)；NDVI 增加趨勢(shì)（0～＜0.002）所占的像元比例依次為多年凍土完全退化區(qū)＞不連續(xù)多年凍土區(qū)＞連續(xù)多年凍土區(qū)＞稀疏島狀多年凍土區(qū)；而NDVI 減少（＜0）所占的像元比例依次為多年凍土完全退化區(qū)＞稀疏島狀多年凍土區(qū)＞不連續(xù)多年凍土區(qū)＞連續(xù)多年凍土區(qū)。

表5 1982－2014 年?yáng)|北不同多年凍土分區(qū)秋季NDVI 變化趨勢(shì)的面積比例Tab.5 Proportions of pixels for the variations trend of autumn NDVI in the different permafrost zones of northeastern China during 1982－2014 %

2.3.2 不同類型多年凍土區(qū)秋季NDVI 與氣候因子相關(guān)性 不同類型多年凍土區(qū)秋季NDVI 與平均氣溫相關(guān)性不顯著（表6），連續(xù)多年凍土區(qū)與不連續(xù)多年凍土區(qū)植被秋季NDVI 與降水量表現(xiàn)出顯著負(fù)相關(guān)關(guān)系，相關(guān)系數(shù)均為－0.426。對(duì)于多年凍土完全退化區(qū)而言，秋季NDVI 與氣溫呈負(fù)相關(guān)關(guān)系，為－0.107。

表6 1982—2014 年?yáng)|北不同類型凍土區(qū)植被秋季平均NDVI 與氣候因子的相關(guān)系數(shù)Tab.6 Correlation coefficients between mean autumn NDVI and climatic factors during growing season in different permafrost zones of northeastern China during 1982－2014

3 結(jié)論與討論

本研究利用景觀生態(tài)學(xué)中“空間代時(shí)間”方法，將連續(xù)多年凍土－不連續(xù)多年凍土－稀疏島狀多年凍土－多年凍土完全退化區(qū)視為多年凍土的退化過程。隨著多年凍土不斷退化，春季NDVI 受地區(qū)氣溫的影響其相關(guān)系數(shù)逐漸增高，這表明多年凍土退化過程有助于春季植被的生長(zhǎng)；而對(duì)于夏季NDVI，多年凍土由連續(xù)多年凍土區(qū)退化成不連續(xù)多年凍土區(qū)，NDVI 與氣溫之間的相關(guān)性系數(shù)增加，即氣溫升高可以促進(jìn)夏季NDVI 的增加，但隨著多年凍土區(qū)由不連續(xù)多年凍土區(qū)退化成稀疏島狀多年凍土區(qū)，直至多年凍土完全退化區(qū)，植被NDVI 與氣溫之間的相關(guān)性系數(shù)逐漸降低，植被受氣溫影響的敏感程度下降。這在一定程度上表明，短期來看，多年凍土區(qū)退化可以促進(jìn)夏季植被生長(zhǎng)，但長(zhǎng)期來看，多年凍土退化甚至消失可能會(huì)阻礙夏季植被生長(zhǎng)。對(duì)于秋季植被NDVI 而言，在多年凍土完全退化區(qū)，秋季NDVI 與氣溫呈負(fù)相關(guān)關(guān)系，在一定程度上表明隨著氣溫升高，多年凍土區(qū)退化會(huì)阻礙秋季植被生長(zhǎng)。

多年凍土退化在植被生長(zhǎng)過程中起到積極作用，這與一些研究結(jié)果一致［19，23］。多年凍土區(qū)隨著地表溫度的逐漸升高，多年凍土?xí)l(fā)生融化，導(dǎo)致凍土區(qū)活動(dòng)層厚度增加，進(jìn)而可以為多年凍土區(qū)的植被生長(zhǎng)提供更多的水分和營(yíng)養(yǎng)物質(zhì)［24］。多年凍土退化屬于長(zhǎng)期的動(dòng)態(tài)過程，而且也存在一定滯后效應(yīng)，多年凍土退化的初期階段，升高的溫度減弱了該區(qū)低溫環(huán)境對(duì)植被生長(zhǎng)的消極作用，增加了土壤水分，促進(jìn)植被生長(zhǎng)。隨著多年凍土進(jìn)一步退化，植被NDVI 與氣溫的正相關(guān)關(guān)系逐漸減弱，可能是長(zhǎng)期退化過程中，土壤水分供應(yīng)不斷減少。因此，短期來看，多年凍土發(fā)育或者維持相對(duì)較好的區(qū)域，隨著溫度的增升高，多年凍土發(fā)生退化，可以促進(jìn)植被生長(zhǎng)，增加植被覆蓋，但長(zhǎng)期的持續(xù)退化甚至消失會(huì)阻礙植被生長(zhǎng)。