易 揚,胡昕利,史明昌,康宏樟,王 彬,張 辰,劉春江,*

1 上海市園林科學規(guī)劃研究院,上海 200232 2 上海交通大學農(nóng)業(yè)與生物學院,上海 200240 3 國家林業(yè)和草原局上海城市森林生態(tài)系統(tǒng)國家定位觀測研究站,上海 200240 4 北京林業(yè)大學水土保持學院,北京 100083 5 上海礎鼎環(huán)境科技有限公司,上海 200433

植被可定義為陸地地表上所覆蓋的所有植物群落的總稱,是陸地生態(tài)系統(tǒng)中的主體,其生長受到自然和人為的共同影響,對環(huán)境變化和氣候變化十分敏感[1-2]。植被長期處在動態(tài)演變之中,并在生物圈的物質循環(huán)、能量流動、生態(tài)系統(tǒng)穩(wěn)定性上發(fā)揮著重要的作用,氣候變化是植被活動的重要影像因素,研究植被的動態(tài)變化和氣候因子的關系,對區(qū)域生態(tài)平衡維護有重要的現(xiàn)實意義[3-4]。歸一化植被指數(shù)(Normalized Difference Vegetation Index,NDVI)是反應和監(jiān)測植被生長狀況和動態(tài)變化的重要指示因子,對地表植被的動態(tài)變化非常敏感(包括植被覆蓋度、生物量和葉面積指數(shù)等),在時間序列的植被監(jiān)測方面有著廣泛的運用,已成為全球變化研究的重要方向之一[5- 7]。

國內(nèi)外學者利用NDVI數(shù)據(jù),從不同的時間和空間尺度上對各個區(qū)域的植被動態(tài)變化規(guī)律進行了研究,探討了各個不同研究區(qū)的NDVI變化特征、演變趨勢和對氣候變化的響應機制[8- 10]。陶帥等通過分析2000—2015年長江上游植被變化特征,指出該區(qū)NDVI未來可能出現(xiàn)波動的區(qū)域占一半以上,海拔和年均溫是影響NDVI變化的主要因素[11]。袁喆等探討了2000—2015年長江流域NDVI與氣象因子的關系,結果表明NDVI呈逐年增加的趨勢,長江流域水量豐沛,熱量條件是影響植被生長的限制性因子[12]。Mao等和Stenhman等研究表明全球的植被變化受人為活動影響較大,特別是一些大型的植樹造林活動[13-14]。John和Ichii等分析全球范圍內(nèi)NDVI與氣候變量之間的關系,結果顯示區(qū)域北部中高緯度地區(qū),植被NDVI與氣溫呈顯著正相關,而在北部和南部半干旱地區(qū),植被NDVI與降水呈顯著負相關[15]。Fensholt等利用NDVI數(shù)據(jù)分析了半干旱地區(qū)植被的季節(jié)變化趨勢和氣候驅動因素,表明全球平均半干旱地區(qū)的植被覆蓋呈增加趨勢,降水是該地區(qū)植物生產(chǎn)的主要限制因素[16]。王桂鋼等通過分析天山地區(qū)氣候參數(shù)對草地季節(jié)變化影響的滯后性特征,指出當旬降水和前一個月降水累積量對于草地長勢的影響最大[17]。

圖1 研究區(qū)Fig.1 Study area

不同區(qū)域尺度上的NDVI響應特征存在差別,已有研究對長江流域植被變化進行了相關分析,然而針對長江中游地區(qū)植被演變趨勢和響應機制還未做充分探討。本研究選取長江中游地區(qū)(湖北省、湖南省和江西省),基于1999—2015年MODIS NDVI數(shù)據(jù)和氣象資料,分析長江中游地區(qū)NDVI的變化特征和演變趨勢,辨析了氣象因子與NDVI在年際尺度和月際尺度上的關系,從而全面揭示了長江中游地區(qū)植被動態(tài)變化規(guī)律,探討了植被對氣候變化的響應機制,以期為長江流域應對氣候變化影響、區(qū)域生態(tài)環(huán)境建設與保護措施的制定,提供科學參考依據(jù)。

1 研究區(qū)概況

長江中游區(qū)域包含湖北省、湖南省和江西省(24°25′N—33°16′N,108°24′E—118°23′E),總面積為56.46萬km2(圖1)。區(qū)域氣候類型屬于亞熱帶季風氣候,降水量達到1000—1600 mm/a,年平均氣溫在16—18℃[18]。長江中游區(qū)域地形主要以山地丘陵與平原為主,海拔最高處位于鄂西山地一帶,最低處位于中部平原地區(qū),其中山地丘陵占50%以上,地勢西高東低,平均海拔約為1497 m[19]。2018年地區(qū)總人口為17438萬人,占全國總人口12.7%,生產(chǎn)總值為9.8萬億元,城市化水平達50%以上,人類經(jīng)濟活動頻繁,對區(qū)域擾動較大[20]。

長江中游地區(qū)林木資源豐富,2015年湖北省、湖南省和江西省的森林覆蓋率分別為39.61%、59.82%和64.00%[21]。主要植被資源有水杉(Metasequoiaglyptostroboides)、冬青(IlexchinensisSims)、銀杉(Cathayaargyrophylla)、珙桐(DavidiainvolucrataBaill.)、銀杏(GinkgobilobaL.)、杜仲(Eucommiaulmoides)、山茶(Camelliajaponica)、南方鐵杉(Tsugachinensisvar.tchekiangensis)竹柏(Podocarpusnagi)、江西槭(Acerkiangsiense)等。

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

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

本研究所用的NDVI數(shù)據(jù)產(chǎn)品為1999—2015年MODIS NDVI數(shù)據(jù),來源于美國國家航空航天局(NASA)的EOS/MODIS植被指數(shù)產(chǎn)品數(shù)據(jù)MOD13A2(https：//wist.echo.nasa.gov/ap)。數(shù)據(jù)空間分辨率為1 km×1 km,時間分辨率為16d,獲取的數(shù)據(jù)利用Envi 5.3對原始數(shù)據(jù)進行預處理。氣象數(shù)據(jù)來源于國家氣象科學數(shù)據(jù)中心(http://data.cma.cn),研究區(qū)附近的57個基準氣象臺站的逐月數(shù)據(jù)(月平均氣溫、月降水及月輻射等),時間序列為1999—2015年。

2.2 研究方法

2.2.1Theil-Sen Median趨勢分析和Mann-Kendall檢驗

Theil-Sen Median趨勢分析計算時間序列數(shù)據(jù)的所有n×(n-1)/2成對組合之間的斜率的中位數(shù)。其計算公式為[22]：

(1)

式中,?是植被變化趨勢,NDVIj是第j年的NDVI值,NDVIi是第i年的NDVI值。如果?>0,則NDVI有增加的趨勢,表明在時間段內(nèi)植被一直在改善或恢復。如果?<0,則NDVI呈下降趨勢,表明在時間段內(nèi)植被呈現(xiàn)退化的趨勢。

Mann-Kenddall檢驗可以判斷趨勢顯著性,其樣本無特定的分布要求,其結果也更為準確,不受異常值的干擾[23]。計算公式如下：

設定NDVIj,j=1999,2000,…,2015,定義Z統(tǒng)計量為：

(2)

(3)

(4)

(5)

式中,NDVIi和NDVIj表示像元i年和j年的NDVI值；n是時間序列的長度；sgn是一個符號函數(shù)。統(tǒng)計量Z的值在(-∞,+∞)范圍內(nèi)。在給定顯著性水平α下,當|Z|>Z1-α/2時,表示時間序列在α水平上存在顯著的變化。在本研究中,取α=0.05,判斷在0.05置信水平上即|Z|>1.96時,1999—2015年區(qū)域NDVI變化趨勢的顯著性。

2.2.2變異系數(shù)分析

變異系數(shù)通過計算標準差與平均值的比值而來,又稱為離散系數(shù),是量度各觀測值與單位均值離散程度的數(shù)學指標[24]。

(6)

2.2.3未來變化趨勢分析

基于重標極差法(R/S)的Hurst指數(shù)是定量描述時間序列長期依賴性的有效方法[25],其計算公式是：

定義時間序列NDVI(t),t=1,2,…,n.,對于任意正整數(shù)τ,定義其均值序列：

(7)

計算累積離差：

(8)

計算極差：

(9)

計算標準差：

(10)

計算Hurst指數(shù)：

(11)

式中,通過使用最小二乘法擬合等式log(R/S)n=a+H×log (n)獲得H值,其中H為赫斯特指數(shù)。根據(jù)Hurst[25]和Mandelbrot[26]的觀點,Hurst指數(shù)值的范圍為(0,1),有3種類型：當H=0.5時,表明NDVI時間序列是隨機的序列,沒有可持續(xù)性,表明了此變化趨勢未來的時間序列與研究的時間序列無關；當H>0.5時,反映NDVI時間序列具有可持續(xù)性,表明未來趨勢變化與當前NDVI時間序列趨勢相同；當H<0.5時,此類未來變化趨勢有可能保持基本穩(wěn)定或者發(fā)生逆轉,本研究中在空間變化上將其歸為變化趨勢不確定[26]。

2.2.4氣象因子數(shù)據(jù)插值

利用1999—2015年研究區(qū)及周邊57個國家級氣象站逐日記錄的氣象數(shù)據(jù)(氣溫、降水、相對濕度、日照時數(shù))進行最大值合成計算出月最大值,基于月最大值求取年均值,形成站點逐年數(shù)據(jù)集,利用ArcGIS軟件的克里金插值法得出氣象因子空間分布。

2.2.5偏相關分析

地球系統(tǒng)是由多要素所構成的復雜系統(tǒng),系統(tǒng)中任何一個要素的變化都會引起其它要素也產(chǎn)生變化,偏相關分析便可有效解決上述問題,單獨研究兩要素之間的相關性時,偏相關將其它要素的影響視為常數(shù)[27]。計算公式為

(12)

(13)

式中,rxy、rxz、ryz分別為要素x與y、x與z、y與z之間的相關系數(shù)；rxy·z為將要素z固定后要素x與y的偏相關系數(shù),通常采用t檢驗法對偏相關系數(shù)進行顯著性檢驗。

圖2 長江中游地區(qū)NDVI平均值年際變化 Fig.2 Interannual variation of average NDVI in the middle reaches of the Yangtze RiverNDVI：表示歸一化植被指數(shù) Normalized difference vegetation index

3 結果

3.1 NDVI的時間變化及空間分布特征

3.1.1NDVI的時間變化特征

1999—2015年NDVI數(shù)據(jù)的年平均值變化成體呈現(xiàn)增加的趨勢,最低值為2001年的0.72,最高值為2015年的0.80,整個研究時段內(nèi)的平均值為0.76。研究區(qū)植被覆蓋向好的趨勢發(fā)展,可大致分為3個階段：第一階段為緩慢下降期(1999—2001年)、第二階段為波動上升期(2002—2011年)、第三階段為平穩(wěn)上升期(2011—2015年)。區(qū)域NDVI在2000—2005年、2006—2010年、2011—2015年的增長率分別為4.86%、2.23%、4.81%。在2000—2005年與2011—2015年時間段內(nèi),研究區(qū)植被得到了顯著的改善(圖2)。

3.1.2NDVI的空間分布特征

利用1999—2015年的年平均NDVI數(shù)據(jù),計算16年間平均值得到長江中游地區(qū)NDVI空間分布圖(圖3)。研究區(qū)植被覆蓋格局整體呈現(xiàn)“西高東低、北高南低”的分布特征,西部地區(qū)海拔較高,以山區(qū)為主,天然林、灌木和草地等分布廣泛,植被生長良好,NDVI值較高。北部地區(qū)以平原為主,是中國重要的商品糧基地,NDVI值同樣較高。NDVI低值區(qū)主要分布在三生省會城市(湖南省長沙市、湖北省武漢市和江西省南昌市)及其附近城市地區(qū),此部分區(qū)域多數(shù)經(jīng)濟發(fā)展迅速,人工擾動較大。16年間,NDVI值在0.1—0.5段占比1.40%,NDVI值在0.5—0.7段的區(qū)域分別占11.45%,NDVI值從0.7—0.8的區(qū)域占55.06%,NDVI值大于0.8的區(qū)域占32.09%。

圖3 1999—2015年長江中游地區(qū)平均NDVI空間分布及像元特征Fig.3 Spatial distribution of average NDVI and pixel distribution in the middle reaches of the Yangtze River from 1999 to 2015

3.2 NDVI趨勢變化特征

3.2.1NDVI動態(tài)趨勢分析

長江中游地區(qū)1999—2015年NDVI空間格局變化特征表現(xiàn)為,植被明顯改善區(qū)域占最大比例(67.39%),主要分布在研究區(qū)西部及東南部；穩(wěn)定不變及輕微改善區(qū)域占區(qū)域面積的四分之一(25.79%),主要以江漢平原、洞庭湖平原、鄱陽湖平原為中心分布；退化區(qū)域占區(qū)域面積的6.82%,零星散布在各大城市的中心區(qū)位置,并以中心區(qū)向外擴散分布；嚴重退化位置以三省省會武漢市、長沙市和南昌市為中心向外輻射,表明城市規(guī)模的擴大對NDVI帶來顯著影響(表1和圖4)。

表1 長江中游地區(qū)NDVI變化趨勢

3.2.2NDVI動態(tài)穩(wěn)定性分析

長江中游區(qū)域1999—2015年NDVI穩(wěn)定性整體表現(xiàn)為“高低波動并存,低波動占比高,區(qū)域性明顯”。空間格局變化特征表現(xiàn)為：低波動區(qū)(0.014

圖4 1999—2015年長江中游地區(qū)NDVI變化趨勢 Fig.4 Trends of NDVI in the Middle Reaches of the Yangtze River from 1999 to 2015

圖5 1999—2015年長江中游地區(qū)植被年均NDVI變異程度 Fig.5 Annual average NDVI variation of vegetation in the middle reaches of the Yangtze River from 1999 to 2015

表2 長江中游地區(qū)NDVI變異系數(shù)分類表

3.2.3NDVI動態(tài)持續(xù)性分析

長江中游地區(qū)1999年至2015年間,植被Hurst小于0.5的區(qū)域占總面積的60.54%,大于0.5的區(qū)域占39.46%。將NDVI變化趨勢與Hurst指數(shù)結果耦合疊加,得到研究區(qū)NDVI變化趨勢與持續(xù)性的耦合空間圖(圖6),將結果分為6類(表3)：持續(xù)性嚴重退化和持續(xù)性輕微退化,其比例為0.79%和2.47%,主要以大斑塊分布在長沙、南昌、武漢省會發(fā)達城市及其周邊；持續(xù)性穩(wěn)定不變,其比例達1.42%；持續(xù)性輕微改善和持續(xù)性明顯改善,其比例為9.78%和25%,主要分布在低山丘陵區(qū)；未來變化趨勢不確定,此類占比高達60.54%,均勻分布在各個省市。

3.3 NDVI與氣候變化的關系

3.3.1NDVI與氣候因子的時間相關性分析

對長江中游地區(qū)1999—2015年NDVI與各氣象因子在年際變化尺度上,分別做了相關、偏相關分析(表4)。結果表明,區(qū)域整體NDVI與氣溫、降水呈正相關,偏相關系數(shù)分別為0.488、0.277,與相對濕度、日照時數(shù)呈負相關,偏相關系數(shù)分別為-0.646、-0.532。NDVI與氣溫的相關性大于與降水的相關性。NDVI與相對濕度的相關性大于日照時數(shù)的相關性,NDVI與各氣候因子的相關系數(shù)及偏相關系數(shù)均未達到顯著水平(P>0.05)。表明在年際變化尺度上,長江中游地區(qū)整體區(qū)域NDVI對同期氣候變化的響應不是特別強烈。

表3 長江中游地區(qū)NDVI變化趨勢及Hurst指數(shù)分類

表4 1999—2015年長江中游地區(qū)NDVI與氣候因子的相關系數(shù)

圖6 長江中游地區(qū)NDVI未來變化趨勢 Fig.6 Future trends of NDVI in the middle reaches of the Yangtze River

分析每月氣候因子與NDVI之間的相關性,可以更好地揭示氣候因子是如何在月際尺度上影響NDVI變化的。對1999—2015年12個月各月份的NDVI年際變化與當月及前1—3月的氣溫、降水、相對濕度、日照時數(shù)進行雙變量相關性分析,得出對應的相關系數(shù)表(表5)。結果顯示,1—4月份NDVI生長受降水和日照影響較為顯著,并且具有明顯的滯后性；6月NDVI與當月相對濕度具有顯著的相關性；7月NDVI受5月降水的滯后性影響；8月與當月的氣溫、相對濕度、日照時數(shù)均存在顯著的相關性,受氣候因子的影響最大；9月受當月降水和相對濕度的共同影響；10月受9月的降水、相對濕度和日照的影響；11月與各個氣象因子均無顯著相關關系；12月受當月相對濕度與日照時數(shù)的影響。降水、相對濕度和日照時數(shù)是影響長江中游區(qū)域NDVI季節(jié)變化的主要氣候因子,降水和日照時數(shù)的時滯性較為明顯。

3.3.2NDVI與氣候因子的空間相關性分析

從空間相關性來看,研究區(qū)絕大部分地區(qū)NDVI年際變化與氣候因子的相關性并不顯著,NDVI與氣溫、降水呈正相關的面積比例大于負相關的面積比例(表6和圖7)。NDVI與氣溫和降水呈正相關的面積比例占71.82%和61.56%,其中,在P<0.01(氣溫：1.85%；降水：0.63%)和P<0.05(氣溫：5.39%；降水：2.82%)水平上達到顯著。NDVI與相對濕度和日照時數(shù)呈負相關的面積比例占66.38%和60.51%,其中,在P<0.01(相對濕度：1.90%；日照：0.32%)和P<0.05(相對濕度：4.66%；日照：2.28%)水平上達到顯著?？傮w來看,長江中游地區(qū)NDVI與各氣候因子達到顯著的相關性比例分別為(圖7和圖8)：7.53%(氣溫)、7.31%(相對濕度)、4.12%(降水)、3.76%(日照時數(shù))。表明在研究區(qū)大面積范圍內(nèi)氣溫升高、降水增多有利于植被的動態(tài)生長。

表5 長江中游地區(qū)當月NDVI與前0—3月各氣候因子的相關系數(shù)

表6 長江中游地區(qū)NDVI年際變化與氣候因子的相關性面積分布/%

圖7 1999—2015年長江中游地區(qū)NDVI與氣候因子空間相關性分類Fig.7 Spatial correlation of NDVI and climatic factors in the middle reaches of the Yangtze River from 1999 to 2015

圖8 1999—2015年長江中游地區(qū)NDVI與氣候因子空間相關性系數(shù)分布Fig.8 Spatial correlation of NDVI and climatic factors in the middle reaches of the Yangtze River from 1999 to 2015

4 討論

1999—2015年,長江中游地區(qū)NDVI呈總體上升趨勢,這與湖北省[28]、湖南省[29]、江西省[30]、長江流域[12]等區(qū)域的研究果相一致。長三角、珠三角等沿海城市群的NDVI年均值多在0.3—0.5之間,與之相比,長江中游地區(qū)NDVI相對較高,多年均值為0.76。從NDVI的變化趨勢來看,NDVI穩(wěn)定不變及輕微改善區(qū)域占區(qū)域面積的四分之一(25.79%),主要以江漢平原、洞庭湖平原、鄱陽湖平原為中心展布,這三大平原是我國主要商品糧基地,符合農(nóng)作物輕微變化特征。退化區(qū)域占區(qū)域面積的6.82%,以三省省會武漢市、長沙市、南昌市為中心向外輻射,表明城市規(guī)模的擴大對植被覆蓋帶來顯著影響。

從NDVI的波動性來看,NDVI較高和高波動區(qū)占7.00%,在區(qū)域中部形成連接三省省會的帶狀分布(長沙市-岳陽市-咸寧市-武漢市-黃岡市-九江市-上饒市-南昌市),其余呈斑塊狀分布在各城市中心,對應輕微退化及嚴重退化區(qū)域。中度波動區(qū)占全區(qū)32.59%,環(huán)簇于較高波動區(qū)邊界并向外延伸,低波動區(qū)與較低波動區(qū)占全區(qū)面積的59.32%,植被變化較小,與NDVI穩(wěn)定不變及輕微改善相吻合。從NDVI的未來變化趨勢來看,未來變化趨勢不確定的區(qū)域占研究區(qū)一半以上,這與前人在長江流域的結果一致[11]。持續(xù)性明顯改善區(qū)域占研究區(qū)四分之一,其斑塊分布在各個省市內(nèi),主要分布在區(qū)域西部山地區(qū),十堰市、神農(nóng)架林區(qū)等地。持續(xù)性嚴重退化區(qū)域(0.79%)則主要以大斑塊分布在長沙、南昌、武漢省會發(fā)達城市,在襄樊市、宜昌市、常德市等人類活動頻繁的較發(fā)達城市呈小斑塊分布。

Nemani等和Liu等認為水熱氣候條件是影響陸地植被覆蓋空間格局的驅動因素[31-32]。本文研究表明,在年際變化上氣溫是影響研究區(qū)植被生長的主要氣候因素,這與崔林麗等[33]、李本綱等[34]和陳云浩等[35]對我國華東及周邊地區(qū)植被NDVI對氣候響應的研究結果基本一致,其原因可能是研究區(qū)降水較豐富,降水足夠滿足植被生長需求,而熱量的差異便是造成NDVI差異的主要驅動因素。本文研究表明,降水對該區(qū)植被生長的影響為滯后2—3個月,日照時數(shù)為滯后2個月,相對濕度為滯后3個月,表明植被NDVI變化也受土壤水分和養(yǎng)分有效性的影響。這與Bao等[36]從全球、區(qū)域等多尺度研究證實結果一致,植被覆蓋對氣候變化的反饋存在一定的滯后效應。從空間相關性來看,研究區(qū)與各氣候因子達到顯著相關的比例為：氣溫>相對濕度>降水>日照時數(shù)。

5 結論

(1)在年際尺度上,研究區(qū)近16年NDVI呈上升趨勢平均值為0.76,植被狀況較好。在空間尺度上,分布呈現(xiàn)“西高東低、北高南低“的分布特征。

(2)從變化趨勢來看,研究區(qū)NDVI改善區(qū)域占比一半以上,穩(wěn)定性表現(xiàn)為“高低波動并存,低波動占比高,區(qū)域性明顯“。區(qū)域NDVI未來變化趨勢不確定的區(qū)域占比60.54%,需要引起持續(xù)關注和重視。

(3)從植被與氣候因子的相關性來看,年際尺度上,研究區(qū)NDVI與各氣象因子關系均不顯著；月際尺度上,降水、相對濕度和日照時數(shù)是影響NDVI變化的主要因子,降水和日照時數(shù)有明顯的時滯性；空間相關性上,植被變化與氣溫、降水呈正相關的面積大于負相關,與相對濕度、日照時數(shù)呈負相關的面積大于正相關。