西南喀斯特?cái)嘞菖璧刂脖晃锖騽討B(tài)變化及其與氣候因子的響應(yīng)

李 丹, 吳秀芹, 張靖宙, 于 洋

(1.北京林業(yè)大學(xué) 水土保持學(xué)院, 北京 100083; 2.北京林業(yè)大學(xué)水土保持國家林業(yè)局重點(diǎn)實(shí)驗(yàn)室, 北京 100083; 3.中國水利水電科學(xué)研究院 泥沙研究所, 北京 100048)

植被物候反映了植物在生命周期中不斷適應(yīng)季節(jié)變化形成與環(huán)境條件密切相關(guān)的生長發(fā)育節(jié)律[1-2]。植被是自然環(huán)境的一面鏡子，環(huán)境中水、土、氣、生等組成成分的任何變化都會被覆蓋在地表的植物群落以不同的形式表現(xiàn)出來，植被的動態(tài)變化蘊(yùn)藏著氣候和人類活動的變遷。氣候的變化對生態(tài)系統(tǒng)的構(gòu)造與性能產(chǎn)生顯著影響，當(dāng)代以增溫為主要特征的全球氣候變化問題，對地球環(huán)境產(chǎn)生了前所未有的影響，成為當(dāng)代人類面臨的重大難題之一[3]。目前國內(nèi)外針對植被物候觀測的方法有多種，不同植被物候觀測技術(shù)具有各自適宜的時間和空間尺度，不同于以野外觀測為基礎(chǔ)的傳統(tǒng)目視觀察法，近20 a來，遙感技術(shù)的廣泛應(yīng)用為物候研究提供了新的角度和機(jī)遇。不同觀測技術(shù)或不同取方法使獲取的植被物候指標(biāo)之間存在顯著的差異性，在植被秋季物候指標(biāo)上尤其明顯。但基于地面監(jiān)測和遙感觀測的研究一致顯示，全球尺度下氣溫的持續(xù)升高已經(jīng)導(dǎo)致了許多地區(qū)植被春季物候始期提前、秋季物候末期推遲[4-6]。

西南喀斯特?cái)嘞菖璧氐靥幍途暩咴匦螐?fù)雜多樣，生物資源豐富，生態(tài)環(huán)境對水熱變化相當(dāng)敏感?；谙嚓P(guān)文獻(xiàn)的統(tǒng)計(jì)，該地區(qū)近50 a的氣候變化與全球氣候變暖趨勢相符，年平均氣溫呈上升趨勢，升溫速率略低于全國平均水平，略大于全球氣候變化平均水平[7]。但針對西南喀斯特?cái)嘞菖璧匚锖蚍矫娴难芯恐饕匀斯び^測法集中在小尺度、特定物種、短時序的研究[8-10]，缺少對區(qū)域尺度上宏觀全面的認(rèn)識。因此，本研究基于2001—2016年的MOD13Q1產(chǎn)品和氣象數(shù)據(jù)，利用Matlab對長時間序列NDVI進(jìn)行重構(gòu)并提取物候參數(shù)，分析典型物候參數(shù)(SOS，EOS，LOS)的時空變化特征，揭示植被物候時空動態(tài)變化與氣候因子的響應(yīng)機(jī)制，重點(diǎn)圍繞兩個問題展開：(1) 喀斯特?cái)嘞菖璧刂脖晃锖蜃兓绾巍?2) 氣候因子對該地區(qū)植被物候產(chǎn)生了怎樣的影響。研究結(jié)果為分析喀斯特?cái)嘞菖璧刂脖晃锖虻臅r空變化動態(tài)特征及其與氣溫、降水的響應(yīng)關(guān)系提供科學(xué)依據(jù)，同時為該地區(qū)根據(jù)物候變化開展植被恢復(fù)提供理論支撐。

1 研究區(qū)概況

西南喀斯特?cái)嘞菖璧?100°03′—105°10′E，22°39′—29°10′N)(圖1)地處云貴高原，南面中越邊境，西臨橫斷山脈，北接青藏高原，是青藏高原隆升產(chǎn)生的斷裂活動所引起的斷塊差異沉陷及溶蝕作用共同形成的山間盆地[11]，是中國西南地區(qū)重要的生態(tài)屏障和生態(tài)脆弱區(qū)。本研究區(qū)橫跨云南、四川和貴州3省共43個區(qū)縣，總面積達(dá)11.04萬km2。研究區(qū)屬亞熱帶季風(fēng)氣候，全年分干濕兩季，雨熱同期，日照充足，部分高海拔山區(qū)的立體氣候特征明顯。其中，有滇東地面定株實(shí)測研究表明，該區(qū)域各物候期的平均溫度和≥5℃積溫對主要鄉(xiāng)土物種的物候表現(xiàn)較為密切，植被多于2月氣溫達(dá)到10℃時開始萌動，至3月氣溫增高至15℃以上開始展葉，在5月溫度超過20℃時，絕大部分開始開花，果實(shí)成熟期分散在4—12月[8]。

圖1 研究區(qū)位置

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

2.1 研究數(shù)據(jù)

歸一化植被指數(shù)(NDVI)采用MODIS/Terra Vegetation Indices 16-Day L3 Global 250 m SIN Grid產(chǎn)品(簡稱MOD13Q1)，數(shù)據(jù)源從2001年1月1日—2016年12月19日，每16 d提供一景，共690景影像。為消除NDVI時間序列的云影噪聲影響，采用TIMESAT軟件的S-G濾波法對數(shù)據(jù)進(jìn)行降噪處理[12]。氣象產(chǎn)品是基于中國氣象數(shù)據(jù)共享網(wǎng)站(http:∥data.cma.cn)得到的2001—2016年西南喀斯特?cái)嘞菖璧貐^(qū)域及周邊的43個氣象站點(diǎn)氣溫和降水?dāng)?shù)據(jù)，通過克里金插值(Kriging)對年降水和年均氣溫?cái)?shù)據(jù)進(jìn)行插值為空間精度250 m的柵格數(shù)據(jù)集。

表1 數(shù)據(jù)概況

2.2 研究方法

(1) 植被物候指標(biāo)提取方法。衛(wèi)星遙感監(jiān)測數(shù)據(jù)提取物候指標(biāo)的常用方法有NDVI閾值法[13]、Gu法[14]、移動平均法[15]、最大斜率法[16]。不同提取方法各有利弊，一定程度上體現(xiàn)出物候研究的不確定性，基于西南喀斯特區(qū)獨(dú)特的地理特性，本研究選用NDVI曲線反演——最大斜率法進(jìn)行物候參數(shù)提取，該方法是基于時間序列的NDVI曲線變化特征確定植被物候期[17-18]，如植被開始迅速生長或植被葉片脫落的時間點(diǎn)對NDVI曲線迅速增長和減小的時間點(diǎn)[19]。計(jì)算公式如下：

NDVIratio(t)=[NDVI(t+1)-NDVI(t)/NDVI(t)]

(1)

式中：t為具體日期，當(dāng)NDVIratio達(dá)到最大時對應(yīng)的NDVI(t)即為生長季始期，當(dāng)NDVIratio達(dá)到最小時對應(yīng)的NDVI(t)即為生長季末期。

NDVI(t)=a0+a1x1+a2x2+…+anxn

(2)

式中：NDVI(t)是x的函數(shù)；x為每年的第多少天，采用最小二乘法來計(jì)算天數(shù)。

(2) 趨勢分析方法。采用一元線性回歸方法建立氣溫、降水(y)與時間序列(x)；物候參數(shù)(y)與時間序列(x)的一元線性回歸方程，方程的斜率b用來表示方程的變化趨勢。

y=a+bt+ε

(3)

式中：y為生長季始期或氣象因子變量；t為年時間序列；a和b為系數(shù)，是擬合的殘差。

(4)

式中：n為累積年數(shù)，本研究為16；變量i為年序號；M為第i年的物候參數(shù)/氣象數(shù)據(jù)值；slope為一元線性回歸方程的斜率，某像元的斜率即為該像元多年的生長季始期氣象數(shù)據(jù)均值變化的總趨勢。

(3) 相關(guān)性分析方法。采用空間相關(guān)分析方法研究兩個要素之間的相關(guān)關(guān)系，相關(guān)系數(shù)取值為[-1，1]。

(5)

式中：n為研究時段年數(shù)；x，y為相關(guān)分析的兩個變量；xi，yi分別為它們的樣本值。

3 結(jié)果與分析

3.1 喀斯特?cái)嘞菖璧刂脖晃锖蜃兓厔菁翱臻g分異

2001—2016年，研究區(qū)內(nèi)3項(xiàng)物候指標(biāo)均呈推遲或延長趨勢，變化率依次為EOS(10 d/10 a)>SOS(7.7 d/10 a)>LOS(2.3 d/10 a)(圖2)。

圖2 物候變化趨勢

研究區(qū)多年SOS均值為第192.41天，空間上SOS表現(xiàn)出東早西晚的趨勢(圖3)。約有45.04%的地區(qū)生長季開始于第170～190天(6月中旬—7月上旬)，如木里河谷地帶、鹽源盆地、川滇邊界、云南省東北部、昆明玉溪以東的南北向谷地和丘北南部等。24.91%的地區(qū)SOS晚于第200天(7月下旬)，大多分布于青藏高原東部邊緣的高海拔山區(qū)。研究區(qū)東南部、南盤江河谷地區(qū)也明顯表現(xiàn)出SOS較早的現(xiàn)象。研究區(qū)多年EOS均值為第339.31天，空間上EOS表現(xiàn)出東早西晚的趨勢。各地EOS時間分異跨度較長，近100 d。最早的如鹽源盆地、安寧河谷、瀘西小江盆地、陸良南盤江水系，EOS多發(fā)生在300 d(10月末)之前；最遲的如木里高海拔區(qū)、哀牢山東部等地，EOS多晚于400 d(次年2月初)。56.28%的地區(qū)EOS發(fā)生在第300～340天(10月末—12月上旬)，植被停止發(fā)育，生長季結(jié)束。海拔與物候期體現(xiàn)出較明顯的相關(guān)性，隨著海拔升高LOS會逐漸縮短，并且SOS會呈現(xiàn)明顯推遲的趨勢。這主要是高海拔地區(qū)分布的山地植被類型對季相變化的響應(yīng)不強(qiáng)烈的表現(xiàn)。

圖3 16 a間SOS，EOS，LOS空間分布

研究區(qū)多年LOS平均為146.91 d，范圍在140～160 d的地區(qū)分布較多，占總面積的37.22%。人類活動較少、氣候條件優(yōu)越的南盤江河谷、紅河河谷、盤龍江河谷地區(qū)LOS長于160 d。研究區(qū)內(nèi)高海拔地區(qū)、部分城鎮(zhèn)開發(fā)區(qū)(安寧河谷、蒙自市西北部等)、農(nóng)事活動區(qū)、石漠化發(fā)育區(qū)的LOS多短于120 d。

在空間上，大部分地區(qū)SOS發(fā)生了不同程度的推遲現(xiàn)象(圖4)，但相反，安寧河谷等部分地區(qū)SOS提前較為明顯。除木里縣低海拔谷地、蒙自市西北部等地區(qū)EOS發(fā)生提前，其他地區(qū)EOS均呈現(xiàn)出推遲現(xiàn)象。統(tǒng)計(jì)得出LOS發(fā)生延長的區(qū)域面積大于縮短的區(qū)域，延長的區(qū)域主要集中分布在鹽源盆地、安寧河谷及研究區(qū)東北部等地，縮短主要分布在木里山谷低地、攀枝花東部、川滇邊界周圍地區(qū)等地。

圖4 物候參數(shù)年變化的空間分布差異

3.2 2001-2016年喀斯特?cái)嘞菖璧貧鉁睾徒邓兓厔?/h3>

研究區(qū)多年年均氣溫為15.06℃，在過去16 a中氣溫以0.33℃/10 a的速度增加。16 a各年總降水量947.04 mm，降水量波動較大，總體呈下降趨勢，年降水量差異明顯，2011年當(dāng)?shù)睾禐?zāi)是16 a來降水量最少的年份。2001—2016年，在春季、夏季和秋季，研究區(qū)的年均氣溫均呈現(xiàn)上升趨勢，變化幅度依次為0.66，0.44，0.36℃/10 a。冬季氣溫呈下降趨勢，降幅為-0.43℃/10 a。春、夏兩季降水量總體呈下降趨勢，春季變化率為-19.12 mm/10 a，夏季變化率為-18.21 mm/10 a。秋冬兩季降水量呈現(xiàn)上升趨勢，秋季變化率14.63 mm/10 a，冬季變化率為2.85 mm/10 a(圖5)。

3.3 喀斯特?cái)嘞菖璧刂脖晃锖蚺c氣候因子相關(guān)性分析

研究分析氣象要素與關(guān)鍵物候期參數(shù)的相關(guān)性有助于進(jìn)一步認(rèn)識喀斯特?cái)嘞菖璧刂脖粚夂蜃兓捻憫?yīng)機(jī)制。結(jié)合氣候和物候動態(tài)變化趨勢的研究結(jié)果，可以發(fā)現(xiàn)在16 a間，喀斯特?cái)嘞菖璧貧鉁睾徒邓淖兓瘜ιL季始期和生長季末期的推遲都產(chǎn)生了較大影響(表2)。

由表2可知，春季氣候因子與SOS具有較強(qiáng)的相關(guān)性，春季降水對植被生長季提前的影響大于溫度對其推遲的影響，說明研究區(qū)植被生長前期對水分的需求是促使植被生長開始的主要因素，16 a來研究區(qū)春季溫度升高加劇了植被的蒸騰作用，結(jié)合春季降水減少進(jìn)一步導(dǎo)致了植被SOS的推遲。全年來說溫度對SOS的影響遠(yuǎn)大于降水對SOS的影響，SOS的推遲與當(dāng)年溫度呈顯著相關(guān)(p<0.05)。EOS受溫度和降水的影響均不顯著，總體來說，降水對EOS的影響大于溫度對其的影響，且夏季降水減少會導(dǎo)致生長季推遲結(jié)束，但研究區(qū)EOS常開始發(fā)生于10月末之后，表明降水對土壤水分貯存及生長季末期均存在一定的時滯作用[20-21]。綜合來看，SOS對氣候因子的敏感性高于EOS。

圖5 2001-2016年溫度、降水變化

表2 植被物候與氣候因子的相關(guān)分析

4 討 論

(1) 遙感技術(shù)的發(fā)展為植被物候研究提供了新手段，但依據(jù)遙感影像計(jì)算所得物候指標(biāo)與人工目視實(shí)測數(shù)據(jù)存在一定差異[22]，且遙感多數(shù)據(jù)源的不同物候提取方法各有利弊，一定程度上體現(xiàn)了植被物候研究的不確定性。本研究區(qū)為亞熱帶地區(qū)，常綠闊葉林是常見的植被類型，終年常綠和季相變化不明顯是該植被類型的特征，也是該區(qū)域利用遙感手段研究其物候變化的疑難點(diǎn)之一。本研究選取NDVI最大斜率法提取物候指標(biāo)，最大化彰顯了該區(qū)域植被生長變化趨勢，宏觀尺度上彌補(bǔ)了對西南喀斯特?cái)嘞菖璧刂脖晃锖蛑笜?biāo)變化的規(guī)律性認(rèn)知，但如何決策選取合適的方法提取物候參數(shù)，還需要通過研究者們深入挖掘植被物候自身特性，并盡可能考慮氣候變化、地貌特征、人類活動等多方面的影響因子，做出更精確的判斷。

(2) 本研究結(jié)果所得研究區(qū)西北部高海拔地區(qū)植被生長季持續(xù)到了次年2月，高海拔地區(qū)植被生長代謝限制因子之間的相互關(guān)系復(fù)雜，關(guān)于植被生長受限機(jī)理存在眾多爭議，其中生長受限假說認(rèn)為高海拔地區(qū)增溫對植被呼吸的促進(jìn)遠(yuǎn)大于光合的促進(jìn)[23]，同時很多學(xué)者實(shí)測研究表明：在植物的生長季期間,植物葉片的凈光合速率并沒有受到生長季低溫環(huán)境的影響而顯著下降，這可能是研究區(qū)西北部高海拔地區(qū)植被生長季結(jié)束較晚的原因之一[24-25]。

喀斯特?cái)嘞菖璧匾蚱洫?dú)特且脆弱的生態(tài)環(huán)境，成為我國長期開展生態(tài)修復(fù)的重點(diǎn)地區(qū)，研究植被物候動態(tài)變化及其與氣候因子的響應(yīng)，對理解喀斯特?cái)嘞菖璧丨h(huán)境下的植被物候格局、探明植被物候變化的驅(qū)動機(jī)制有著重要意義。植被物候作為自然生態(tài)環(huán)境變化的重要指示因子，將其納入生態(tài)修復(fù)成效的評估指標(biāo)將會促進(jìn)評價(jià)體系進(jìn)一步完善。

5 結(jié) 論

(1) 西南喀斯特?cái)嘞菖璧刂脖晃锖騍OS和EOS均呈現(xiàn)東早西晚的空間分異特征，LOS分異特征受海拔和人類活動的影響較大，三類物候指標(biāo)均呈現(xiàn)不同程度的推遲或延長趨勢。(2) 喀斯特?cái)嘞菖璧厣L季始期對氣候的敏感性高于生長季末期，且溫度是生長季始期發(fā)生推遲的主導(dǎo)自然因子。區(qū)分季節(jié)維度分析，春季氣候因子對生長季始期的推遲影響最為顯著，夏季氣候因子對生長季末期的推遲影響尤為關(guān)鍵。研究揭示了氣候變化背景下西南喀斯特?cái)嘞菖璧厮疅釛l件對植被物候變化特征的影響趨勢，有助于探明氣候因子對區(qū)域植被變化的驅(qū)動機(jī)制，了解區(qū)域自然生態(tài)系統(tǒng)的結(jié)構(gòu)和功能特性。