李 珍，王永波*，牟乃夏，劉文寶

(1.首都師范大學(xué)資源環(huán)境與旅游學(xué)院,中國 北京 100048；2.山東科技大學(xué)測繪與空間信息學(xué)院，中國 青島 266590)

陸地植被作為地球生態(tài)系統(tǒng)的重要組成部分，是人類生存和發(fā)展不可或缺的自然資源[1,2]。植被覆蓋度(FVC)可以反映單位面積內(nèi)植被在地面垂直投影面積占比[3]，是衡量生態(tài)環(huán)境變化的一個(gè)重要標(biāo)準(zhǔn)，也是國家進(jìn)行自然資源管理和宏觀調(diào)控的重要依據(jù)[4]。近年來，隨著全球變暖，生態(tài)環(huán)境惡化、資源短缺等全球性問題的陸續(xù)出現(xiàn)，了解地表覆被情況對(duì)于土地資源利用、生態(tài)環(huán)境恢復(fù)等具有重要意義[5-8]。另一方面，長時(shí)間尺度植被演化同樣是過去全球變化研究的重要手段，在探討和重建過去氣候環(huán)境變遷當(dāng)中發(fā)揮了重要作用[9,10]。因而，深入理解現(xiàn)代氣候變化過程對(duì)植被覆蓋度變化的影響，探討植被覆蓋度與氣候因子的相關(guān)性，對(duì)了解現(xiàn)代生態(tài)環(huán)境變化以及理解過去植被變化過程均有重要的科學(xué)意義[3,11]。

植被覆蓋度調(diào)查的傳統(tǒng)方法主要是實(shí)地測量，包括目估法、采樣法、儀器法和模型法等[12]，但是地面測量往往只能進(jìn)行小區(qū)域的植被覆蓋度監(jiān)測，且耗時(shí)較長、工作量大、受條件限制多，尤其是結(jié)果受人為因素影響較大，具有一定的局限性[13,14]。近年來，隨著遙感技術(shù)不斷發(fā)展與成熟，諸多學(xué)者利用遙感數(shù)據(jù)研究了不同空間尺度植被覆蓋度變化并對(duì)其與氣候因子的相關(guān)性進(jìn)行了探討[15-18]。比如，基于MODIS-NDVI(2000—2016)數(shù)據(jù)集，丁玥等研究發(fā)現(xiàn)和田地區(qū)近17年的植被覆蓋度呈上升趨勢，并且結(jié)合氣溫、降水?dāng)?shù)據(jù)獲取了不同時(shí)空尺度上NDVI與氣候因子的相關(guān)性[19]。趙舒怡等利用MODIS-NDVI數(shù)據(jù)和混合像元分解模型法發(fā)現(xiàn)2001—2013年華北地區(qū)植被覆蓋度呈上升趨勢且與氣候干旱程度存在明顯相關(guān)關(guān)系[20]?；贚andsat TM數(shù)據(jù)，李苗苗等利用像元二分模型法估算了2002年5月和7月北京密云水庫上游的植被覆蓋度，發(fā)現(xiàn)植被覆蓋度的變化一定程度上也受到地形條件差異的影響[21]。因而，利用歸一化植被指數(shù)(NDVI)研究植被覆蓋度變化的方法已較為成熟，尤其是基于Landsat數(shù)據(jù)的像元二分模型法削弱了大氣、土壤與植被類型等的影響，操作簡單且不受地區(qū)限制，得到了廣泛應(yīng)用[22,23]。

青島市地處山東半島東南部沿海，地形復(fù)雜，受到典型溫帶季風(fēng)氣候影響，植被類型豐富[24]。青島市近年來植被覆蓋度變化研究一方面可以為城市生態(tài)環(huán)境治理提供參考，另一方面，對(duì)探討我國東部季風(fēng)區(qū)植被對(duì)氣候變化的響應(yīng)過程也具有重要意義[18]。然而，該地區(qū)已有研究主要關(guān)注過去10年的植被變化。比如，宋鵬飛等利用MODIS數(shù)據(jù)獲取了山東省過去10年植被覆蓋度變化趨勢，發(fā)現(xiàn)其與氣候因子變化具有一定相關(guān)性，并且對(duì)溫度變化的敏感度高于降水[25]；同樣基于MODIS數(shù)據(jù)，趙鑫等對(duì)青島市2008—2017年間的植被覆蓋度動(dòng)態(tài)變化進(jìn)行了研究[24]。另一方面，這些研究中所采用的時(shí)間、空間分辨率較低，相對(duì)于MODIS數(shù)據(jù)而言，Landsat數(shù)據(jù)具有更高的空間分辨率(30 m)，可以獲取區(qū)域范圍內(nèi)更為精確的植被覆蓋度數(shù)據(jù)[21]。

因此，本文利用青島市1998—2019年的Landsat衛(wèi)星影像，依托ENVI軟件平臺(tái)，使用像元二分模型法分析青島市近21年來年際尺度的植被覆蓋度變化情況，并結(jié)合氣象站點(diǎn)的觀測氣象數(shù)據(jù)，探討植被覆蓋度與氣候因子的相關(guān)性，以期在全球氣候變化背景下對(duì)青島市生態(tài)環(huán)境的治理提供科學(xué)依據(jù)，為季風(fēng)區(qū)長時(shí)間尺度植被和氣候演化研究工作提供參考。

1 研究區(qū)概況

青島市地處山東半島東南部，位于東經(jīng)119°30′～121°00′、北緯35°35′～37°09′，東、南瀕臨黃海，東北與煙臺(tái)市毗鄰，西與濰坊市相連，西南與日照市接壤，總面積[24]為1.1 萬km2。青島為海濱丘陵城市，地勢東高西低，南北兩側(cè)隆起，中間低凹(圖1A)。山地約占青島市總面積的15.5%，丘陵占2.1%，平原占37.7%，洼地占21.7%[24]。青島地處北溫帶，屬溫帶季風(fēng)氣候，溫度適中，年平均氣溫12.7 ℃，全年8月份氣溫最高(平均氣溫25.3 ℃)，1月份溫度最低(平均氣溫-0.5 ℃)(圖1B)；受夏季風(fēng)影響雨量充沛，年平均降水量約為660 mm，且主要集中在7～8月(圖1C)。春季氣溫回升緩慢，夏季濕熱多雨，秋季降水少，蒸發(fā)強(qiáng)，冬季風(fēng)大溫度低，持續(xù)時(shí)間較長[26]。

圖1 研究區(qū)位置(A)及各月氣溫(B)、降水(C)分布圖

受地形、氣候等因素影響，青島地區(qū)植物種類豐富，為同緯度地區(qū)植被建群種較多的地區(qū)[27]。全市以華北植物區(qū)系為主，主要為溫帶落葉闊葉林和溫帶針葉闊葉林，灌叢和草原較少，中部平原地區(qū)以分布栽培植物為主，南部分布少量常綠闊葉喬灌木成分[27]。區(qū)域內(nèi)共有種子植物151科，768屬，1 469種，其中栽培種占40%。常見喬木主要為櫟、槭、楊、榆、柳、楓、刺槐、黃連木和泡桐等，裸子植物主要為松屬，灌木以荊條、酸棗等為主，草本植物分布最廣的為菅草、白羊草，等。此外，在東南沿海及海島上分布少量亞熱帶植物，如耐冬、楠木、胡枝子，等[28]。平原地區(qū)農(nóng)業(yè)栽培植物分布廣泛，以冬小麥、玉米、高粱，等為主，山麓地區(qū)分布部分栽培果樹。濱海地區(qū)分布有少量鹽生植被、砂生植被及沼生植被等植物[28]。

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

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

本研究所用Landsat數(shù)據(jù)來源于美國地質(zhì)調(diào)查局(USGS)官網(wǎng)(https://www.usgs.gov/)，其中包括1998—2011年的Landsat 4-5 TM數(shù)據(jù)和2013—2019年的Landsat 8 OLI數(shù)據(jù)。數(shù)據(jù)空間分辨率為30 m，地圖投影坐標(biāo)為UTM-WGS 84投影坐標(biāo)系。

為獲取研究區(qū)植被覆蓋度的結(jié)果，選擇植物生長期的7～9月Landsat圖像作為源數(shù)據(jù)，并且通過對(duì)比2009年7～9月共3個(gè)月的影像處理結(jié)果，發(fā)現(xiàn)該3個(gè)月植被覆蓋度結(jié)果較為一致，因此，對(duì)青島地區(qū)植被覆蓋度的計(jì)算以7～9月為主，選擇云量較少的影像進(jìn)行計(jì)算。然而，由于Landsat影像時(shí)間分辨率較低，為減少季節(jié)變化和云量對(duì)結(jié)果的影響，本次研究選取的影像主要集中在8月和9月(除2005年7月)且時(shí)間跨度在一個(gè)月內(nèi)(表1)。所選取的衛(wèi)星影像云量大多在5%以下，多于5%的影像中云量主要集中在海域或研究區(qū)范圍外，不會(huì)影響研究區(qū)內(nèi)植被覆蓋度的計(jì)算。利用ENVI 5.3和ArcGIS 10.2軟件作為處理平臺(tái)，對(duì)獲取的Landsat影像進(jìn)行預(yù)處理，在ENVI中利用快速大氣校正工具(QUAC)進(jìn)行大氣校正操作，去除部分云的影響。

表1 青島市1998—2019年各級(jí)植被覆蓋度組成

同時(shí)，為了探討研究區(qū)植被覆蓋度變化對(duì)氣候因子的響應(yīng)過程，從國家氣象科學(xué)數(shù)據(jù)中心(http://data.cma.cn/)獲取了青島市氣象站(區(qū)站號(hào)54857)記錄的1998—2019年的逐月氣象數(shù)據(jù)，包括月降水量、月平均氣溫及月最低溫等。另外，為了分析人為因素對(duì)城市環(huán)境植被覆蓋度的潛在影響，從青島市統(tǒng)計(jì)局(http://qdtj.qingdao.gov.cn/)獲取了1998—2018年間青島市人口、國內(nèi)生產(chǎn)總值(GDP)、城市建成區(qū)面積和鋪裝道路面積數(shù)據(jù)。

2.2 研究方法

遙感圖像信息普遍受到大氣、土壤背景與植被類型等因素的影響，植被覆蓋度的提取還受到地形、植物群落結(jié)構(gòu)、太陽高度角以及大氣噪音等方面的干擾，利用像元二分模型和歸一化植被指數(shù)(NDVI)可以有效地降低以上因素帶來的誤差，提高植被覆蓋度計(jì)算精度[21]。NDVI具有植被檢測靈敏度高、植被覆蓋度檢測范圍寬的特點(diǎn)，可以消除地形和群落結(jié)構(gòu)的干擾，削弱太陽高度角和大氣噪音的影響[21]。像元二分模型主要是利用NDVIsoil與NDVIveg兩個(gè)調(diào)節(jié)因子，僅保留植被覆蓋度信息，降低大氣、土壤背景及植被類型等因素的影響。已有研究利用該方法對(duì)青島市嶗山區(qū)植被覆蓋度進(jìn)行分析，獲取的模型結(jié)果與實(shí)測結(jié)果呈顯著相關(guān)，表明利用該方法可以有效地提取區(qū)域植被覆蓋度結(jié)果[29]。

2.2.1 歸一化植被指數(shù)提取 歸一化植被指數(shù)(NDVI，又稱為標(biāo)準(zhǔn)化植被指數(shù))，是根據(jù)遙感傳感器接收的地物光譜信息計(jì)算得出的地表植被狀況定量值，用以區(qū)分植被與其他地物，是目前廣泛應(yīng)用的一種植被指數(shù)[30]。其計(jì)算公式為

N=(NIR-R)/(NIR+R)，

(1)

其中R為紅波段，NIR為近紅外波段。具體應(yīng)用到Landsat 4-5 TM和Landsat 8 OLI遙感影像的波段中時(shí)，分別對(duì)應(yīng)第3與第4通道，第4與第5通道，波長范圍分別為0.63～0.69 μm和0.76～0.90 μm，0.630～0.680 μm和0.845～0.885 μm[31]。

2.2.2 植被覆蓋度的計(jì)算——像元二分模型法 像元二分模型法假設(shè)像元只由植物覆蓋的區(qū)域和無植物覆蓋的區(qū)域兩部分構(gòu)成，這兩部分的面積在像元中所占的比例就是這兩部分的比重，地面有植物覆蓋的部分所占像元的比例即為影像的植被覆蓋度[30]。本文在改進(jìn)后的像元二分模型法[31]基礎(chǔ)上利用歸一化植被指數(shù)來計(jì)算青島市的植被覆蓋度：

f=(N-Nsoil)/(Nveg-Nsoil),

(2)

其中，Nsoil是裸土或沒有植物覆蓋區(qū)域的N數(shù)值，Nveg代表全部被植物所覆蓋區(qū)域的N數(shù)值[31]。

根據(jù)前人的研究結(jié)果，本文在統(tǒng)計(jì)影像 NDVI 值、區(qū)域 NDVI 累積概率分布表的基礎(chǔ)上，查閱相關(guān)文獻(xiàn)[32,33]，依次選取NDVI累積概率分布為5%和95%的數(shù)值作為Nsoil和Nveg并代入式(2)，計(jì)算青島市不同時(shí)期的植被覆蓋度。

2.2.3 植被覆蓋度分級(jí)制圖 根據(jù)青島市的區(qū)域特點(diǎn)，查閱已有研究文獻(xiàn)[21,24]，將青島市植被覆蓋度按等間距方法進(jìn)行分級(jí)：0

3 結(jié)果分析

3.1 青島市近21年植被覆蓋度變化

基于1998—2019年期間8、9月Landsat數(shù)據(jù)，利用像元二分模型法獲取了青島市近21年的植被覆蓋度變化和各級(jí)植被覆蓋度的空間分布特征(表1，圖2)。從中可以發(fā)現(xiàn)，1998—2019年間青島市植被覆蓋度總體處于一個(gè)較高的水平，以高覆蓋為主，除2009,2011和2016年外，高覆蓋所占比例均高于50%；其次是較高覆蓋，然后依次是中覆蓋、較低覆蓋和低覆蓋。對(duì)不同等級(jí)的植被覆蓋度對(duì)比分析也發(fā)現(xiàn)，近21年來低覆蓋、較低覆蓋、中覆蓋三個(gè)等級(jí)的變化幅度較小，較高覆蓋、高覆蓋兩個(gè)等級(jí)的變化幅度較大，其中高覆蓋的變化幅度最為顯著(表1)，可以有效地用來指示研究區(qū)植被的變化。從時(shí)間上來看，近21年來青島市植被覆蓋度持續(xù)波動(dòng)，總體呈增長趨勢，高覆蓋度區(qū)域所占比例從1998年的55.41%增長到2019年的67.48%(表1)。從空間分布上來看，青島市植被覆蓋度呈現(xiàn)中間低，兩邊高的特點(diǎn)，北部的平度市、萊西市以及南部的黃島區(qū)高植被覆蓋度所占比例較大(圖2)。

圖2 青島市1998—2019年植被覆蓋度空間分布特征

3.2 青島市植被覆蓋度變化對(duì)氣候因子的響應(yīng)

城市植被覆蓋度的變化既受到自然因素如氣候條件等因素影響，又受到人為因素的潛在影響[34]。改革開放以來，青島市發(fā)展迅速，城市人口的增加、經(jīng)濟(jì)的發(fā)展、城市建設(shè)等可以直接影響到植被覆蓋度的變化；另一方面，城市產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整，實(shí)行生態(tài)建設(shè)改善城市生態(tài)環(huán)境，一定程度上也可以提高城市植被覆蓋度[34]。本研究主要從人口、GDP、城市建成區(qū)面積、鋪裝道路面積等方面對(duì)植被覆蓋度的影響進(jìn)行分析。統(tǒng)計(jì)分析結(jié)果表明，青島市高植被覆蓋度變化與人口、GDP、城市建成區(qū)面積、鋪裝道路面積的相關(guān)系數(shù)分別為-0.069 7，-0.081 5，-0.004 6和-0.138 3，均呈弱負(fù)相關(guān)關(guān)系，表明隨著城市擴(kuò)張的進(jìn)行，植被覆蓋度呈小幅的下降趨勢。由于人為因素對(duì)植被覆蓋度變化的影響相對(duì)有限，本文主要針對(duì)氣候因子和植被覆蓋度的關(guān)系進(jìn)行討論。

3.2.1 青島市近21年植被覆蓋度年際變化影響因素 影響植被覆蓋度的氣候因子范圍比較廣泛，當(dāng)月氣溫、降水，本季度的平均氣溫、降水，年平均氣溫、降水，都會(huì)對(duì)植被的生長產(chǎn)生影響?？紤]到氣候因子對(duì)植被覆蓋度的影響可能存在滯后效應(yīng)，因此本文對(duì)前一月的氣溫、降水，前兩月的氣溫、降水和上一年冬季的平均氣溫、降水以及年最低溫度也進(jìn)行了相關(guān)性分析，選擇中等程度及以上相關(guān)性的各氣候因子進(jìn)行統(tǒng)計(jì)分析(圖3)。

從結(jié)果中可以發(fā)現(xiàn)，青島地區(qū)與植被覆蓋度變化呈中等程度及以上相關(guān)的氣候因子包括當(dāng)月氣溫、夏季(7～9月)平均氣溫、前一月降水及年最低氣溫，相關(guān)性大小依次為：年最低氣溫、當(dāng)月氣溫、前一月降水和夏季平均氣溫(圖3)?？傮w上，植被覆蓋度變化與氣溫的相關(guān)性高于降水，其中與當(dāng)月平均氣溫的相關(guān)系數(shù)為0.55，略高于前一月降水(相關(guān)系數(shù)0.42)。在年尺度上，植被覆蓋度變化與年最低氣溫的相關(guān)系數(shù)為0.62，呈現(xiàn)較強(qiáng)的正相關(guān)關(guān)系，說明在年尺度上低溫對(duì)研究區(qū)植物生長具有較強(qiáng)的影響。

圖3 青島市植被覆蓋度年際變化與各氣候因子關(guān)系圖(虛線代表高植被覆蓋度，實(shí)線代表不同氣候因子)

3.2.2 月尺度植被覆蓋度與氣候因子關(guān)系 由以上分析可知，青島市植被覆蓋度年際變化與當(dāng)月氣溫和前一月降水的相關(guān)性較高，為了進(jìn)一步探討區(qū)域植被覆蓋度對(duì)氣候因素的響應(yīng)，尤其是植被覆蓋度變化對(duì)氣候因子影響的滯后性，我們?cè)谀觌H變化數(shù)據(jù)結(jié)果中增加4～6月份Landsat圖像處理結(jié)果(除2012年因數(shù)據(jù)缺失，每年均有植被覆蓋度結(jié)果)，在此基礎(chǔ)上，對(duì)植被生長與前一月以及當(dāng)月氣候因子的相關(guān)系數(shù)進(jìn)行統(tǒng)計(jì)分析(圖4)。

青島市高植被覆蓋度與當(dāng)月平均氣溫的相關(guān)系數(shù)為0.87，與前一月氣溫的相關(guān)系數(shù)為0.76，與當(dāng)月降水量的相關(guān)系數(shù)為0.38，與前一月降水量的相關(guān)系數(shù)為0.74(圖4)。其中，高植被覆蓋度與當(dāng)月平均氣溫為極強(qiáng)相關(guān)，與前一月降水和當(dāng)月氣溫為強(qiáng)相關(guān)關(guān)系，表明青島市植被生長與同期氣溫的相關(guān)性更加顯著，與降水的關(guān)系存在一定的滯后性。

圖4 植被覆蓋度與前一月及當(dāng)月氣候因子關(guān)系圖(虛線代表植被覆蓋度，實(shí)線代表各氣候因子)

4 結(jié)論

基于Landsat-NDVI遙感數(shù)據(jù)，本文利用像元二分模型法提取了青島市1998—2019年的植被覆蓋度，同時(shí)結(jié)合氣象觀測數(shù)據(jù)，分析了植被覆蓋度變化和氣候因子的相關(guān)關(guān)系，得出以下主要結(jié)論：

(1)青島市植被覆蓋度在1998—2019年期間總體處于較高水平，高覆蓋度占50%以上，近21年來植被覆蓋度不斷波動(dòng)變化，總體呈增長趨勢。

(2)植被覆蓋度年際變化與氣溫變化的相關(guān)性總體高于降水，與氣候因子在月尺度上的相關(guān)性高于年尺度，且植被變化同樣對(duì)溫度變化更為敏感。在年際、月際尺度上植被覆蓋度與氣溫的相關(guān)性均高于降水，表明氣溫是影響研究區(qū)植被生長的主要因子。

(3)植被覆蓋度對(duì)氣溫變化的最大響應(yīng)月為當(dāng)月, 無明顯的滯后性；植被覆蓋度對(duì)降水變化則存在1個(gè)月左右的時(shí)間滯后。

5 討論

基于Landsat數(shù)據(jù)，利用像元二分模型法獲取的青島市各級(jí)植被覆蓋度組成與之前用MODIS數(shù)據(jù)建立的植被覆蓋組成情況基本一致[24]，但高覆蓋度的占比較MODIS數(shù)據(jù)結(jié)果偏高，可能與Landsat影像的空間分辨率更高有關(guān)。近21年來青島市植被覆蓋度總體呈增長的趨勢，對(duì)山東省2005—2015年及2000—2014年期間研究中也發(fā)現(xiàn)植被覆蓋度呈現(xiàn)增長的趨勢，與本文的研究結(jié)論相一致[25,35]。然而，趙鑫等人[24]基于MODIS數(shù)據(jù)發(fā)現(xiàn)近十年來青島市植被覆蓋度呈下降的趨勢，可能是由于該研究中所選取的研究時(shí)段以及數(shù)據(jù)時(shí)間間隔不同所致。在空間上，青島市植被覆蓋度呈現(xiàn)中間低、兩邊高的特征，這主要是由于青島市南北都屬于農(nóng)業(yè)種植區(qū)，植被覆蓋度高于城區(qū)和未開發(fā)的裸地。相應(yīng)地，城市化率高的地區(qū)植被覆蓋度要明顯低于城市化率較低的地區(qū)，二者在空間上呈顯著的負(fù)相關(guān)關(guān)系，植被覆蓋度的整體空間分布特征與基于MODIS數(shù)據(jù)所得到的結(jié)果基本一致[24]。此外，研究發(fā)現(xiàn)2009年植被覆蓋度出現(xiàn)低值，與華北平原1981—2013年植被覆蓋度時(shí)空變化特征研究中所發(fā)現(xiàn)的極低值一致[34]。資料記錄表明，2007年我國發(fā)生了大范圍、持續(xù)性的嚴(yán)重干旱，2008 年“中國大旱”拉開序幕，隨后在2009年發(fā)生特大旱情，對(duì)華北地區(qū)影響較大，造成植被覆蓋度偏低[34]。另外，2016年中國自然災(zāi)害頻發(fā)，風(fēng)雹、洪澇和地質(zhì)災(zāi)害災(zāi)情偏重且主要集中在6、7月，尤其是華北地區(qū)災(zāi)情嚴(yán)重，導(dǎo)致2016年植被覆蓋度迅速下降，在一定程度上反映了區(qū)域植被覆蓋度對(duì)氣候變化響應(yīng)的敏感性。

青島市植被覆蓋度變化與氣候因子的相關(guān)性分析表明，其與氣溫的相關(guān)性大于降水。陳歡等發(fā)現(xiàn)中國大陸植被NDVI與溫度和降水具有明顯的空間響應(yīng)，并且與溫度的相關(guān)性更加明顯[36]。青島市植被生長與月尺度氣候因子關(guān)系比較密切，且與同期氣溫的相關(guān)性更加顯著，與降水關(guān)系存在滯后性。山東省2005—2015 年植被狀況研究同樣發(fā)現(xiàn)研究區(qū)內(nèi)NDVI的重心遷移軌跡與降水、溫度在方向上具有明顯的相關(guān)性[25]。溫度和降水分別通過土壤養(yǎng)分有效性和水分有效性影響植被生長， 其中植被對(duì)氣溫的滯后響應(yīng)與土壤溫度變化、土壤有機(jī)質(zhì)分解和營養(yǎng)元素釋放有關(guān), 而植被對(duì)降水的滯后響應(yīng)則反映了降水事件與降水到達(dá)植物根部，影響植物生長的時(shí)間間隔[37]。有研究表明，植被生長對(duì)溫度和降水變化的時(shí)間滯后隨植被類型、土壤特征、地形條件等因素而變化[38]。青島地區(qū)植被類型主要為溫帶落葉闊葉林，水土光熱條件較好，屬于豐水地區(qū)，相對(duì)于干旱地區(qū)而言，植被對(duì)降水的響應(yīng)速度更慢，滯后時(shí)間更長。對(duì)1982—2012年中國植被覆蓋時(shí)空變化特征研究中也得出我國北方地區(qū)植被對(duì)降水有約1個(gè)月的響應(yīng)時(shí)間[8]，阿多等人則發(fā)現(xiàn)華北平原植被覆蓋度與區(qū)域前一個(gè)月的月均溫的相關(guān)性低于當(dāng)月均溫，與前一個(gè)月降水量的相關(guān)性高于當(dāng)月降水量[34]。青島市植被生長受溫度的影響高于降水，可能是因?yàn)榍鄭u市位于沿海地區(qū)，降水豐富，因此植被生長對(duì)降水變化敏感性弱于溫度。

從上述分析結(jié)果中發(fā)現(xiàn)，人為因素對(duì)青島市植被覆蓋度變化的影響較為有限，但均呈現(xiàn)負(fù)相關(guān)關(guān)系，表明隨著人類活動(dòng)的干預(yù)及城市的擴(kuò)張，城市植被覆蓋度存在小幅下降的趨勢。因此，政府應(yīng)當(dāng)將環(huán)境保護(hù)納入國民經(jīng)濟(jì)與社會(huì)發(fā)展計(jì)劃，在經(jīng)濟(jì)發(fā)展中防治環(huán)境污染和生態(tài)破壞，提高人們保護(hù)環(huán)境的意識(shí)，保障區(qū)域生態(tài)環(huán)境平衡。