浙江省近20年植被變化趨勢及驅(qū)動因子分析

賀忠華，張育慧*，何月，張小偉，蔡菊珍，雷莉萍

1.浙江省氣候中心，浙江 杭州 310017；2.中國科學院空天信息創(chuàng)新研究院/數(shù)字地球重點實驗室，北京 100094

植被是陸地生態(tài)系統(tǒng)的重要組成部分，不僅參與生物圈碳水循環(huán)等調(diào)節(jié)區(qū)域氣候，也能反映區(qū)域自然和人類活動的變化（劉啟興等，2019；趙維清等，2019）。遙感觀測是當前植被變化監(jiān)測的重要手段。研究表明，自1982年以來，全球和區(qū)域植被活動發(fā)生了一系列變化，包括21世紀前后的降低和增長趨勢等（Park et al.，2010；Yuan et al.，2018）。大量研究將植被變化與氣候因子進行關聯(lián)分析，探討其受氣候因子及其變化的影響，特別是中國南方的溫度變化和北方的降水限制等（Liu et al.，2015；Piao et al.，2011；Zhang et al.，2011）。

植被變化的氣候響應存在顯著的區(qū)域差異。趙鴻雁等（2019）研究發(fā)現(xiàn)甘肅中東部地區(qū)植被變化與溫度存在顯著相關，降水次之、太陽輻射最小。李思慧（2019）發(fā)現(xiàn)氣候變暖和植被生長季變長是影響通遼2003年以來植被增長的主要因素。賈路等（2019）研究發(fā)現(xiàn)降水增加對西安植被覆蓋具有顯著促進作用。另外，熊巧利等（2019）研究發(fā)現(xiàn)植被覆蓋和生長對氣候變化的響應受到地形因素的影響。而在經(jīng)濟快速發(fā)展區(qū)域，植被年際變化與氣象因子則可能沒有顯著關系（何全軍，2019）。因此，局部地區(qū)的植被生長對氣候因子的響應存在顯著差異，理解區(qū)域植被變化及其驅(qū)動因子的敏感性差異有利于區(qū)域生態(tài)保護和治理。

除氣候影響外，人為因素也是影響區(qū)域植被變化的重要因子。劉啟興等（2019）通過歸一化植被指數(shù)（Normalized Difference Vegetation Index：NDVI）與氣象因子的相關和復回歸分析發(fā)現(xiàn)，氣象因素和人類活動在一定程度上共同促進黃河源區(qū)植被2000年以來的生長。趙維清等（2019）研究表明，城市擴張是影響湖北省植被變化的主要原因。王瑾等（2019）研究表明人類的環(huán)境保護意識提升對內(nèi)蒙草地有一定的積極作用。董晨煒等（2017）發(fā)現(xiàn)2007年以來杭州灣地區(qū)城市擴張對建成區(qū)植被產(chǎn)生了不利影響。所以，探索植被變化的氣候響應同時，分析人為因素對植被趨勢變化的影響有利于分析城市擴張的生態(tài)影響和區(qū)域的可持續(xù)發(fā)展。

研究顯示，二十一世紀以來全球出現(xiàn)植被綠化，且中國對其貢獻了近25%（Chen et al.，2019；Zhu et al.，2016）。中國的植被綠化主要是由于東部地區(qū)森林（Ahrends et al.，2017）和農(nóng)田（Zuo et al.，2018）植被的持續(xù)增長。浙江是位于中國東南部的省份，以丘陵地帶為主，地處亞熱帶季風氣候。該區(qū)域土地覆蓋以森林和農(nóng)業(yè)區(qū)為主，而浙北地區(qū)是長江三角洲經(jīng)濟區(qū)的重要組成部分。因此，浙江植被變化及其驅(qū)動因子解析對區(qū)域生態(tài)文明建設具有重要意義。高大偉等（2010）和何月等（2012）研究了浙江省1982年至21世紀初植被覆蓋的趨勢變化及其氣候響應，發(fā)現(xiàn)區(qū)域的溫度抑制和濕度促進作用。董晨煒等（2017）利用三期（2000、2007、2013年）夜光數(shù)據(jù)與植被指數(shù)的閾值關聯(lián)分析，探討了杭州灣城市擴張對植被的不利影響。由于植被變化趨勢受到不同主導因素影響，其在遠離人類活動的森林區(qū)域、人類參與種植的農(nóng)業(yè)區(qū)域以及以人類活動為主的城市及周邊區(qū)域存在明顯差異。然而，浙江區(qū)域2000年以來的植被趨勢變化的空間分布及其主導驅(qū)動歸因研究還未開展。

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

1.1 研究區(qū)概況

圖1 浙江省2015年土地覆蓋（a）和數(shù)字高程及區(qū)縣代表氣象站點分布（b）Fig.1 Land cover in 2015 (a) and digital elevation and district and county representative meteorological station distribution (b) over Zhejiang Province

本研究以浙江省區(qū)域為例，分析氣候變化和城市發(fā)展對區(qū)域生態(tài)植被變化的影響。浙江省土地覆蓋和數(shù)字高程及區(qū)縣代表氣象站點分布見圖1，其地表覆蓋以森林、農(nóng)田和城市為主（圖1a）。森林區(qū)域主要分布在浙南和浙西等較高海拔地區(qū)，農(nóng)田區(qū)域主要分布在浙中洼地和浙北平原地區(qū)（圖1b），而城市區(qū)域主要分布在杭州灣、寧波、臺州和溫州等沿海地區(qū)。其中，土地覆蓋數(shù)據(jù)來源于歐空局土地覆蓋產(chǎn)品（http://maps.elie.ucl.ac.be/CCI/viewer/index.php），數(shù)字高程數(shù)據(jù)來源于美國國家航空航天局提供的航天飛機雷達地形任務（Shuttle Radar Topography Mission: SRTM）數(shù)字高程數(shù)據(jù)集。浙江屬于亞熱帶濕潤季風氣候，年平均氣溫15—18 ℃，年平均降水量980—2 000 mm，年平均日照時數(shù)1 710—2 100 h。

1.2 數(shù)據(jù)和預處理

本研究主要使用數(shù)據(jù)包括表征植被活動的中分辨率成像光譜儀（Moderate-resolution Imaging Spectroradiometer：MODIS）觀測增強型植被指數(shù)（MOD13A2，Enhanced Vegetation Index：EVI）、總初級生產(chǎn)力（MOD17A2H，Gross Primary Production：GPP）等植被生態(tài)參量數(shù)據(jù)；表征氣候變化的區(qū)縣代表氣象站觀測氣溫、降水和日照等氣象和氣候因子；以及表征人類活動特征的美國國防軍事氣象衛(wèi)星（Defense Meteorological Satellite Program-Operational Linescan System：DMSP/OLS）和Soumi-NPP對地觀測衛(wèi)星（Suomi National Polar-Orbiting Partnership/Visible Infrared Imaging Radiometer Suite：Soumi-NPP/VIIRS）觀測夜間燈光數(shù)據(jù)等。

1.2.1 植被生態(tài)參量

由于研究區(qū)域植被覆蓋茂密，為減少植被指數(shù)飽和對研究結(jié)果的影響（Huete et al.，2002），本研究采用EVI數(shù)據(jù)進行植被趨勢變化研究。該數(shù)據(jù)來源于美國國家航空航天局2000—2019年Terra-MODIS觀測的16天合成產(chǎn)品，空間分辨率為1 km（Huete et al.，2011）。

為探討區(qū)域植被生長相關的植被光合作用變化，本研究采用EVI對應時期（2000—2019年）的Terra-MODIS觀測結(jié)果反演得到GPP數(shù)據(jù)，其時空分辨率為8 d和500 m（Running et al.，2015）。該產(chǎn)品為美國國家航空航天局（National Aeronautics and Space Administration：NASA）發(fā)布的官方產(chǎn)品，已經(jīng)過多次的地面資料驗證（Turner et al.，2006；Liu et al.，2015）。

文中使用的MOD11A2-EVI和MOD17A2HGPP產(chǎn)品數(shù)據(jù)均通過谷歌地球引擎（Google earth engine）獲取（Gorelick et al.，2017）。數(shù)據(jù)預處理包括：（1）云掩膜處理，獲取高質(zhì)量觀測數(shù)據(jù)結(jié)果；（2）通過20年觀測數(shù)據(jù)進行年內(nèi)不同時期的中值求解，表征植被年內(nèi)常態(tài)變化；（3）在區(qū)域或不同土地覆蓋類型內(nèi)的參數(shù)均值求解中，原始影像數(shù)據(jù)空白值由對應時期中值數(shù)據(jù)填補；（4）GPP數(shù)據(jù)采用雙線性內(nèi)插進行空間重采樣至1 km，累積計算至16天的時間分辨率。

1.2.2 氣象和氣候因子

本研究使用的同期氣象數(shù)據(jù)來源于浙江省氣象局，選取省內(nèi)75個區(qū)縣代表氣象站觀測2000—2019年逐日氣溫、降水和日照時數(shù)結(jié)果。氣象站的空間分布參見圖1b。為探討氣候因子與植被變化的對應關系，本研究將氣象要素的日均值數(shù)據(jù)重采樣到16天時間分辨率。

1.2.3 夜間燈光數(shù)據(jù)

為分析生態(tài)植被變化與人類活動的關系，本研究采用DMSP/OLS和Soumi-NPP/VIIRS觀測夜間燈光數(shù)據(jù)表征人類活動變化（Elvidge et al.，2010）。DMSP/OLS夜光指數(shù)為年數(shù)據(jù)，空間分辨率為30′′（弧秒，浙江地區(qū)約為0.8—1 km）。該數(shù)據(jù)來源于不同傳感器：F14（2000—2003）、F15（2000—2008）、F16（2004—2009）和F18（2010—2013）。由于不同傳感器和成像條件的差異，影響數(shù)據(jù)的連續(xù)性和可比性。本研究采用DMSP/OLS夜光數(shù)據(jù)校正方法（鄒進貴等，2014；Liu et al.，2012），獲取得到2000—2013年的DMSP/OLS夜光數(shù)據(jù)。該處理過程包括：（1）穩(wěn)定像元篩選；（2）參考年份數(shù)據(jù)及夜光穩(wěn)定區(qū)域選取；（3）不同傳感器觀測數(shù)據(jù)融合和參考校正；（4）夜光數(shù)據(jù)的時序校正等步驟。

相比DMSP/OLS，Soumi-NPP/VIIRS夜光數(shù)據(jù)在時空分辨率和觀測信號敏感性等方面有了較大提升。其觀測時間為2012年4月至今，數(shù)據(jù)時空分辨率為月均值和15′′（弧秒，浙江地區(qū)約為0.4—0.5 km）。由于Soumi-NPP/VIRR和DMSP/OLS的夜光數(shù)據(jù)差異，本文采用楊任飛（2018）引用的數(shù)據(jù)整合方法，得到2000—2019年遙感觀測夜間燈光指數(shù)產(chǎn)品。

1.3 研究方法

本研究采用的方法主要包括：（1）一元線性回歸和Mann-Kendall趨勢檢驗；（2）年際變化提取和數(shù)據(jù)歸一化；（3）相關分析與多元線性回歸。通過一元線性回歸和趨勢檢驗獲取生態(tài)植被變化趨勢及其顯著性。利用年際變化提取和數(shù)據(jù)歸一化得到植被和氣候因子的年際變化特征，并使其具有可比性。使用相關分析與多元線性回歸方法進行生態(tài)植被變化的自然和人為因素的驅(qū)動分析。

1.3.1 一元線性回歸和Mann-Kendall趨勢檢驗

本研究采用一元線性回歸計算EVI數(shù)據(jù)的趨勢變化特征，使用Mann-Kendall對計算趨勢進行顯著性檢驗。該方法對EVI進行逐像元應用，獲取EVI數(shù)據(jù)的時空變化特征。研究中使用EVI數(shù)據(jù)2000—2019年的16 d時間分辨率的數(shù)據(jù)進行16 d趨勢計算，進而轉(zhuǎn)換成年趨勢。趨勢計算見公式（1）、（2）。

式中，Slope(16d)為2000—2019年間時間分辨率為16 d的EVI數(shù)據(jù)的變化趨勢；n為數(shù)據(jù)總量；i為數(shù)據(jù)序列號；EVIi表示序列號為i的像元EVI數(shù)據(jù)；Slopea為年趨勢（MODIS-EVI數(shù)據(jù)每年設置為23個16 d的時間單位）。

Mann-Kendall是一種非參數(shù)統(tǒng)計方法，被廣泛應用于氣象、水文等數(shù)據(jù)的趨勢判別，它對序列數(shù)據(jù)沒有強制的分布要求，并不受少量異常數(shù)據(jù)的干擾。本研究對EVI逐像元時間序列進行Mann-Kendall趨勢檢驗，判別不同像元趨勢的顯著性差異。趨勢計算見公式（3）—（5）。

式中，ZMK為單個序列的Mann-Kendal標準化檢驗統(tǒng)計量；S統(tǒng)計量是由前后影像的EVI數(shù)據(jù)計算得到。var(S)為S統(tǒng)計量的方差，求解公式可簡化為var(S)=n(n-1)(2n+5)/18。sgn(θ)提供邏輯判別函數(shù)。本文分別進行不同置信水平（0.01、0.05和0.1）上的EVI趨勢檢驗，即ZMK＞2.58，ZMK＞1.96和ZMK＞1.65。

1.3.2 年際變化提取和數(shù)據(jù)歸一化

為了獲取植被變化和氣候因子的年際變化特征，本研究采用多年數(shù)據(jù)不同時期的中值表征植被和氣候的年內(nèi)常態(tài)變化特征，用于年內(nèi)變化的特征去除。年際變化特征提取見公式（6）、（7）

式中，EVIChangei,j表示EVI年際變化時間序列；EVIi,j表示第i年、第j期的EVI數(shù)據(jù)。EVImedianj表示2000—2019年第j期的中值數(shù)據(jù)。

由于植被指數(shù)和多種氣候因子表征內(nèi)容不同，且在不同區(qū)域像元本底和變化值存在較大差異，故在數(shù)據(jù)間的對比分析中，需要進行數(shù)據(jù)的歸一化處理。本研究采用Z分數(shù)進行數(shù)據(jù)的歸一化處理，公式如下：

式中，EVIChangeZi,j表示EVI變化的歸一化序列。mean(EVIChangei,j)和std(EVIChangei,j)分別表示EVI變化時序數(shù)據(jù)的均值和標準差。

年際變化的累積變化量則為多年數(shù)據(jù)與初始年份（2000）對應時期數(shù)據(jù)進行差值處理，再進行累積計算，獲取物理量的近20年的累積變化量。

1.3.3 相關分析與多元線性回歸

本研究通過對植被變化和氣候因子、夜間燈光數(shù)據(jù)等進行相關分析，分別從自然氣候因子和人類活動兩方面進行植被變化驅(qū)動因子分析。應用皮爾斯相關系數(shù)和t檢驗進行相關分析。皮爾斯相關系數(shù)計算公式：

式中，ra,b為序列a和b的相關系數(shù)；ai為序列第i個數(shù)據(jù)的EVI值；bi為氣候因子或夜間燈光指數(shù)的第i個數(shù)據(jù)值；和分別為EVI序列和氣候因子或夜間燈光序列的均值。

本研究利用EVI及與其具有顯著相關關系（P＜0.1）的變量進行多元線性回歸，獲取標準化系數(shù)（Beta），分析多變量對EVI年際變化的貢獻強弱及其分布。其中，EVI數(shù)據(jù)為植被生長季期間的年均值，夜間燈光指數(shù)為年值，氣候和氣象因子為像元植被生長季期間的鄰近站點觀測均值。像元的鄰近站點采用距離小于50 km的多個站點或者大于50 km的最近站點。

2 結(jié)果與分析

2.1 浙江省區(qū)域生態(tài)植被時空變化

圖2 浙江省2000—2019年EVI年趨勢變化（a）、及GPP累積變化（b）Fig.2 Yearly trends of EVI (a) and cumulative change of GPP (b) from 2000 to 2019 over Zhejiang Province

圖3 浙江省2000—2019年EVI趨勢顯著性Fig.3 Significance of EVI trends from 2000 to 2019 over Zhejiang Province

浙江省2000—2019年EVI趨勢及近20年GPP累積變化量見圖2，EVI趨勢顯著性見圖3。參考圖1中土地類型分布圖，可以發(fā)現(xiàn)浙江省內(nèi)EVI（圖2a）存在顯著趨勢（P＜0.05）變化的植被主要覆蓋于森林及部分天養(yǎng)農(nóng)田區(qū)，且為顯著增長。而在城市區(qū)域的EVI顯示出降低趨勢，特別是杭州灣、寧波、臺州和溫州城市及周邊地區(qū)。圖2b結(jié)果顯示，GPP累積量在EVI顯著增長和減少趨勢區(qū)域分別出現(xiàn)明顯增長（＞2.0 kg·m-2）和減少（＜-1.0 kg·m-2）。總體結(jié)果顯示，浙江地區(qū)近20年出現(xiàn)植被的顯著增長，并引起植被光合作用碳吸收的增強，而浙江城市及周邊地區(qū)的植被生態(tài)出現(xiàn)明顯退化。

圖4 浙江省2000—2019年EVI年均值（a）、及GPP年均值（b）Fig.4 Yearly mean values of EVI (a) and GPP (b) from 2000 to 2019 over Zhejiang Province

圖5 浙江省2015年典型土地類型劃分Fig.5 Typical land cover in 2015 over Zhejiang Province

表1 浙江省典型土地類型上EVI和GPP的統(tǒng)計Table 1 EVI and GPP statistics on typical land cover over Zhejiang Province

如圖4所示，浙江EVI和GPP多年均值的高值主要分布在浙南和浙西的高海拔森林區(qū)域（EVI＞0.4；GPP＞1.8 kg·m-2），低值則分布在浙東沿海的平原城市及周邊區(qū)域（EVI＜0.1；GPP＜0.6 kg·m-2）。另外，本文對浙江地區(qū)植被指數(shù)和總初級生產(chǎn)力的均值和變化按典型土地類型劃分（圖5）進行統(tǒng)計見表1。結(jié)果顯示，全省呈現(xiàn)植被增長趨勢以及近20年的植被光合作用碳吸收的增加，植被指數(shù)增長和光合作用碳吸收最高值在森林區(qū)[EVI趨勢：(1.78±0.94) 10-3·a-1，GPP變化量：(2.47±2.12) kg·m-2]，最低值出現(xiàn)在城區(qū)。全省EVI增加趨勢達年均值的0.40%，2000—2019的EVI總變化趨近年均值的7.66%。而近20年所有陸地像元GPP與2000年數(shù)據(jù)的相對變化統(tǒng)計總增加量約為(1.92±2.34) kg·m-2，為多年均值 (1.30±0.41) kg·m-2的1.48倍。另外，圖6顯示了典型土地類型植被指數(shù)EVI均值的時序變化特征。其中，林區(qū)EVI均值最高且存在顯著增長趨勢（0.60×10-3a-1），2016年夏季達到峰值；而城區(qū)EVI均值呈現(xiàn)顯著減少趨勢（-1.11×10-3a-1）。

2.2 植被變化與氣候因子的關系

EVI與浙江省區(qū)縣代表氣象站觀測氣溫、降水和日照時數(shù)的省內(nèi)均值時序分布見圖7。省內(nèi)EVI均值為主要植被覆蓋區(qū)（EVI＞0.15）的16天合成數(shù)據(jù)，氣候因子為75個站點觀測日均產(chǎn)品的16天均值合成數(shù)據(jù)。圖7a結(jié)果顯示，EVI與氣溫具有良好的一致性分布，主要表現(xiàn)為季節(jié)變化和趨勢增長的相似性，其時序相關系數(shù)達0.94（P＜0.01）。而EVI與日均日照和日均降水量時長存在較弱的相關關系（r=0.55，P＜0.01；r=0.27，P＜0.01）。植被生長季（4—9月）期間的EVI和氣候因子年際變化（Z分數(shù)年均值）（圖7b）結(jié)果顯示，EVI年際變化與氣溫變化存在顯著正相關關系（r=0.54，P＜ 0.05），而與降水和日照時數(shù)變化無顯著相關。EVI與氣溫、降水和日照時數(shù)的多元線性回歸分析結(jié)果顯示，氣溫為顯著影響因子，而降水和日照時數(shù)為排除變量，回歸模型的Beta為0.54。以上結(jié)果表明，溫度變化是浙江省區(qū)域近20年植被年際變化的主導氣候因子。

2.3 植被變化與夜光指數(shù)的關系

圖8所示為浙江省近20年夜光指數(shù)均值和變化量。圖8a顯示出浙江省內(nèi)夜光高亮值主要分布在杭州灣周邊城市群，主要包括杭州、寧波等市區(qū)，金華、臺州和溫州城區(qū)等，亮度值向城區(qū)周邊降低。圖8b顯示2000—2019年間，各城市周邊均出現(xiàn)一定的夜光增強，表征了區(qū)域城市擴張范圍。

圖6 浙江省典型土地覆蓋區(qū)域EVI均值的時序變化Fig.6 Temporal variation of EVI in typical land cover over Zhejiang Province

圖7 浙江省EVI與氣候因子16天均值（a）和年均值（b）的時序變化Fig.7 Temporal variation of mean values of EVI and climate factors in 16 days (a) and yearly mean (b) over Zhejiang Province

圖8 浙江省夜間燈光指數(shù)近20年均值（a）和變化量（b）Fig.8 The mean value (a) and changes (b) of night light index in recent 20 years over Zhejiang Province

本研究對城市夜光多年均值進行數(shù)據(jù)篩選，評估夜光與植被指數(shù)的相關關系。通過對夜間燈光亮度值進行分層，分別分為0—10、10—20、…、50—63共6層，從每層數(shù)據(jù)中隨機篩選100個采樣點，獲取夜光多年均值及對應EVI多年均值數(shù)據(jù)，其散點圖如圖9a所示。結(jié)果顯示，夜間燈光高低值分別對應EVI數(shù)據(jù)低高值，呈現(xiàn)顯著負相關關系（r=0.62，P＜0.01）。為探索城市夜光與植被指數(shù)變化的區(qū)域差異性，本研究利用夜光數(shù)據(jù)年均值與EVI生長季（4—9月）年均值進行逐像元相關分析，相關結(jié)果（P＜0.1）見圖9b。結(jié)果表明，杭州灣周邊城區(qū)，特別是杭州市東北和東南區(qū)域、湖州東部、寧波市周邊城區(qū)等地，植被指數(shù)年際變化與夜光指數(shù)年均值呈現(xiàn)顯著負相關關系（r= -0.60±0.13）。而在部分遠離城市區(qū)域，夜光指數(shù)與植被指數(shù)呈現(xiàn)正向相關關系（r=0.56±0.13）。因此，城市發(fā)展對周邊植被具有一定的抑制作用，主要表現(xiàn)為城市擴張引起的土地利用變化等。而在遠離城市區(qū)域，夜間燈光對本地植被生長存在一定促進作用，相關研究有待進一步探索。

2.4 植被變化的人為和氣候貢獻率分布

利用植被生長季期間EVI的多年變化與夜光指數(shù)、鄰近站點觀測氣象和氣候因子進行多元回歸，計算得到的標準化參數(shù)空間分布見圖10。結(jié)果顯示，表征人類活動的夜光指數(shù)對EVI的年際變化影響廣泛。浙東、浙北地區(qū)夜光指數(shù)對城市周邊區(qū)域植被生長起到顯著抑制作用（Beta為-0.53±0.22），而在浙中、浙西和浙南地區(qū)，夜光指數(shù)具有一定的促進作用（Beta為0.46±0.21）。各種氣候因子的影響則分布在不同的局部地區(qū)。站點觀測氣溫主要對浙西南植被顯著增長區(qū)（圖2a和圖3）的植被變化具有明顯的促進作用（Beta為0.52±0.22）。降水對浙西部分地區(qū)植被具有明顯的促進作用，而對浙西北和浙東南地區(qū)植被具有較弱的抑制作用。另外，日照時數(shù)對浙北部分植被區(qū)域具有一定的促進作用。

圖9 EVI與夜間燈光指數(shù)的近20年均值的散點圖（a）及像元年均值時序相關的空間分布（b）Fig.9 The scatter diagram of average EVI and night light index for 20 year (a) and spatial distribution of time series correlation coefficients between EVI and night light index (b)

圖10 植被生長季期間EVI年際變化的多變量貢獻分布Fig.10 Distribution of the multivariate contribution to EVI’s inter-annual change during the growing period

3 討論

浙江省近20年植被年際變化與二十世紀末存在顯著差異，整體趨勢由顯著減少（何月等，2012；高大偉等，2010）轉(zhuǎn)變?yōu)轱@著增長，這與全球植被整體變化趨勢（Yuan et al.，2018）具有良好的一致性。浙江植被指數(shù)近20年總增長趨近年均值的7.66%，接近于中國東部2003—2016年間的植被指數(shù)增長比例（He et al.，2020）；其植被總生產(chǎn)力近20年總增量為 (1.92±2.34) kg·m-2，即年增長趨勢約為 (9.14±11.14) g·m-2·a-1，與2000—2016年的浙江地區(qū)GPP增長趨勢分布一致（Ma et al.，2019）。研究結(jié)果表明，溫度是促進近20年區(qū)域植被生長的主要氣候因子，這與20世紀末的影響因子存在顯著差異。降水量和干濕程度的大小對二十世紀末的浙江植被變化具有正向作用（何月等，2012），而熱量則對植被生長具有明顯的抑制作用（高大偉等，2010）。這與該地區(qū)近20年降水量逐年增加有關：水汽充沛使得其不再成為區(qū)域植被年際變化的主導因子。浙江地區(qū)被發(fā)現(xiàn)是中國東部植被綠化的重要區(qū)域，它與相關的農(nóng)業(yè)發(fā)展和水土保持具有重要關系（Chen et al.，2019），這也印證了遠離城市區(qū)域植被變化的人類活動的促進作用（夜光指數(shù)正相關于EVI）。浙江省內(nèi)大部分城市的邊緣區(qū)域植被呈現(xiàn)顯著的減少趨勢，且與夜光指數(shù)具有顯著負相關關系，表明城市的發(fā)展和擴張對于周邊植被會產(chǎn)生不利影響（劉沁萍等，2014），這與杭州灣城市擴張對植被的影響結(jié)果研究具有一致性（董晨煒等，2017）。該研究結(jié)果顯示，氣候條件及其變化是影響浙江整體植被年際變化的重要因子，而人類活動對省內(nèi)植被的影響存在區(qū)域差異。

本研究采用了多元線性回歸的標準化系數(shù)探討多種因子對植被年際變化的貢獻率，半定量地描述了各變量的貢獻強度及其空間分布。但是，當前變量對浙江植被年際變化的代表性有待進一步探討，各變量之間的自相關及其正交化也需要進一步完善。當前的貢獻強度可用于變量間的相對比較，但難以定量刻畫變量對植被變化的絕對貢獻。另外，本研究使用的遙感數(shù)據(jù)空間分辨率較粗（1 km），難以滿足市縣等局部地區(qū)生態(tài)植被變化的評估需求。當前的夜光指數(shù)和氣候因子對植被變化的整體貢獻不高，更多影響因素的共同評估有利于區(qū)域植被變化的氣候和人為貢獻率評價。

4 結(jié)論

本研究使用多源衛(wèi)星遙感和地面觀測生態(tài)植被參量、氣象和氣候因子、人類活動相關數(shù)據(jù)等分析浙江區(qū)域生態(tài)植被變化及其驅(qū)動歸因，通過線性回歸和趨勢分析等探討了浙江區(qū)域生態(tài)植被的顯著增長。結(jié)果發(fā)現(xiàn)，2000—2019浙江省內(nèi)EVI增加比例趨近多年均值的7.66%，而近20年GPP的增加量約為多年均值的1.48倍。其中主要增長區(qū)域為浙南、浙中和浙西區(qū)域的森林覆蓋區(qū)，而浙北和浙東區(qū)域的城市區(qū)域為明顯減少區(qū)域。通過年際變化提取和相關分析等發(fā)現(xiàn)氣溫升高對省內(nèi)植被增長趨勢具有顯著貢獻（r=0.54，P＜0.01）。浙江城市周邊區(qū)域植被趨勢變化受到人類活動的顯著抑制影響（r= -0.60±0.13），而在遠離城市區(qū)域夜光數(shù)據(jù)呈現(xiàn)正向驅(qū)動（r=0.56±0.13）。植被年際變化的多變量貢獻結(jié)果顯示，表征人類活動的夜光數(shù)據(jù)在全省范圍具有重要影響，而氣候因素則會影響局部地區(qū)。因此，氣候條件和人類活動分別在不同區(qū)域影響植被的年際趨勢變化，其中人類活動影響范圍更廣。