李曉東，閆守剛

(白城師范學院，吉林 白城 137000)

地表植被的生長是以年為周期的一種自然現(xiàn)象。植物的生長既是對地球表層生態(tài)環(huán)境變化的反饋，同時也是對一個時間階段以來地區(qū)生態(tài)環(huán)境累積過程的反應(yīng)，是連接生物圈各生態(tài)因子的關(guān)鍵鏈，還是地球表面自然地理因子發(fā)生全球變化的“指示器”[1-2]。在農(nóng)業(yè)生產(chǎn)上，農(nóng)作物春播的預測和農(nóng)作物防災(zāi)減災(zāi)都與農(nóng)作物物候研究關(guān)系緊密，所以農(nóng)作物物候研究已經(jīng)成為地理學研究的熱點之一[3-11]。國志興等[12]以東北地區(qū)地表植被為研究對象，選擇分段式Logistic函數(shù)擬合法模擬了研究區(qū)20年的植被物候期，分析了在全球環(huán)境變化背景下，東北地區(qū)地表植被對氣候變化的響應(yīng)。楊永民等[13]以隴西黃土高原地表植被為研究對象，利用Logistic和Gaussian兩種非線性的方法模型提取了研究區(qū)的物候信息，為大范圍植被監(jiān)測提供了新思路。武永鋒等[14]以中國大陸地表植被為研究對象，建立了陸地地表植被綠度期的遙感監(jiān)測模型，即Logistic擬合模型，并驗證了該模型結(jié)果的可靠性和優(yōu)越性。

筆者參考常見的幾種遙感監(jiān)測方法，針對吉林西部的植被生長情況，由實驗室內(nèi)的人工判定取代傳統(tǒng)的回歸方程分析，對最終數(shù)學運算結(jié)果進行插值運算，得到研究區(qū)域地表植被生長期的起始與終止日期(即植被綠度始期和綠度末期)以及生長期的持續(xù)天數(shù)，并結(jié)合趨勢分析研究本地區(qū)植被物候期在最近10年的變化特點，以期為該區(qū)域農(nóng)業(yè)生產(chǎn)決策和農(nóng)牧業(yè)種植方式的調(diào)整提供科學依據(jù)，并對該地區(qū)在全球變化背景下，分析植被對氣候變化的響應(yīng)提供一定借鑒。

1 材料與方法

1.1研究區(qū)概況 以松嫩平原為研究對象。主要研究區(qū)包括白城、松原兩個地級市，東邊是東北平原腹地，西鄰呼倫貝爾草原，南接科爾沁沙區(qū)，北依大興安嶺林地，位于農(nóng)牧過渡帶;屬中溫帶大陸性季風氣候，春季干燥多風、升溫比較快，夏季炎熱、降水集中，秋季涼爽、變溫快、溫差大、冬季漫長、干燥寒冷。

1.2數(shù)據(jù)獲取與處理分析

1.2.1地表植被指數(shù)數(shù)據(jù) 采用的NDVI數(shù)據(jù)為國家自然科學基金委員會“中國西部環(huán)境與生態(tài)科學數(shù)據(jù)中心”(Data source:Environment & Ecology Scientific Data Center of Western China, National Natural Science Foundation of China)制作的1998-2008年每15 d最大值合成(MVC)的SPOT-NDVI數(shù)據(jù)集，該數(shù)據(jù)集主要用于地表植被和生態(tài)環(huán)境的研究[15]。

1.2.2植被生長物候期的計算 將1998-2008年36旬的NDVI時間序列數(shù)據(jù)轉(zhuǎn)化為頻率累積方程，其頻率累積方程曲線表現(xiàn)出與植被生物量曲線相接近的“S”型生長變化特點。根據(jù)這一特點，構(gòu)建基于NDVI頻率累積曲線方程的Logistic擬合模型來求算植被生長的物候期。構(gòu)建過程主要包括頻率累積轉(zhuǎn)化、曲線擬合、擬合優(yōu)度檢驗以及植被生長物候期(包括始期和終期)的確定[16-17]。

求出NDVI頻率累積時序曲線方程的曲率最大值。曲率最大值處彎曲程度最大，說明在該點處對應(yīng)日期之后地表植被開始迅速增長，可以判定此極值點對應(yīng)日期為地表植被生長的初始日期,同樣可以確定植被生長終期與生長期的持續(xù)天數(shù)。處理分析都是在遙感軟件平臺ERDAS與地信軟件平臺Surfer和微軟Excel軟件中完成。

2 結(jié)果與分析

2.1數(shù)據(jù)結(jié)果的驗證 在研究區(qū)內(nèi)，選取180個樣點，利用1998、2003和2008年(分別代表前期、中期和后期)的研究區(qū)內(nèi)地表植被累計分布函數(shù)對應(yīng)研究區(qū)的同期擬合數(shù)值，使用Logistic函數(shù)對研究區(qū)內(nèi)的NDVI頻率累積方程進行擬合驗證(圖1-3)。

1998、2003和2008年研究區(qū)內(nèi)地表植被累計分布函數(shù)與擬合曲線的相關(guān)性數(shù)值分別為，0.995 666 652，0.995 529 646和0.997 939 054。圖中擬合優(yōu)度(R2)越大，表明相關(guān)越好。NDVI頻率累積曲線與擬合數(shù)值之間的擬合優(yōu)度均在0.994以上，說明整體擬合效果良好。

圖1 1998年松嫩平原地表植被累計分布函數(shù)與擬合曲線圖Fig.1 NDVI cumulative frequency and fitting curve of Songnen Plain, in 1998

圖2 2003年松嫩平原地表植被累計分布函數(shù)與擬合曲線圖Fig.2 NDVI cumulative frequency and fitting curve of Songnen Plain, in 2003

圖3 2008年松嫩平原地表植被累計分布函數(shù)與擬合曲線圖Fig.3 NDVI cumulative frequency and fitting curve of Songnen Plain, in 2008

2.2松嫩平原地表植被生長始期10年來的變化情況 對生長季初始日期進行趨勢分析得到變化斜率數(shù)值。對趨勢分析結(jié)果數(shù)據(jù)作負值處理，最后得到的數(shù)據(jù)正值代表生長初期提前的天數(shù)，負值代表延后的天數(shù)。基于180個采樣點上的地表植被生長初始期數(shù)據(jù)計算得到，研究區(qū)內(nèi)地表植被生長平均初始期的趨勢變化斜率為0.86，即整個研究區(qū)地表植被生長初始期平均提前0.86 d。在空間分布上，研究區(qū)以北的地表植被生長平均初始期變化斜率為0，與正常年份持平；長春以西，松嫩平原以南地區(qū)的地表植被生長初始期的變化斜率數(shù)值以正值為主，其中部分地區(qū)生長初始期的變化斜率≥10，即提前10 d以上(圖4)。

2.3松嫩平原地表植被生長終止期10年來的變化情況 由生長季終止日期得到的變化斜率的實際意義就是延后的天數(shù)，所以變化斜率正值代表了生長終止期延后的天數(shù)，負值代表了提前結(jié)束的天數(shù)。基于180個采樣點上的地表植被數(shù)據(jù)計算得到，研究區(qū)內(nèi)地表植被生長平均終止期的趨勢變化斜率為3.17，即整個研究區(qū)地表植被生長終止期平均延后3.17 d。在空間分布上，研究區(qū)內(nèi)的中部地區(qū)地表植被生長平均終止期變化斜率為0，與正常年份持平，嫩江流域和松嫩平原南部的地表植被生長平均終止期的變化斜率≥5，即有5 d以上的延后趨勢(圖5)。

圖4 松嫩平原地表植被生長初始期趨勢變化圖Fig.4 Changes trend of the beginning date of greenness period, Songnen Plain

2.4松嫩平原地表植被生長期總天數(shù)10年來的變化情況 在研究區(qū)內(nèi)地表植被生長期總天數(shù)的趨勢分析中，生長季持續(xù)天數(shù)變化斜率的正值代表了生長期增加的天數(shù)，負值代表了減少的天數(shù)?；?80個采樣點上的地表植被的生長季持續(xù)天數(shù)，利用統(tǒng)計軟件對其進行趨勢分析，得出研究區(qū)內(nèi)的地表植被生長期總天數(shù)趨勢變化斜率為2.61，即生長期的天數(shù)有2～3 d的增加幅度。在空間分布上，研究區(qū)內(nèi)的中、西部地區(qū)地表植被生長期總天數(shù)的變化斜率為0，變化穩(wěn)定，長春以西，松嫩平原以南地區(qū)和嫩江流域以北地區(qū)的趨勢變化斜率≥20，即生長季持續(xù)天數(shù)延長20 d以上(圖6)。

圖5 松嫩平原地表植被生長終止期趨勢變化圖Fig.5 Changes trend of the end date of greenness period, Songnen Plain

圖6 松嫩平原地表植被生長期總天數(shù)趨勢變化圖Fig.6 Changes trend of the total days of greenness period, Songnen Plain

3 討論與結(jié)論

研究區(qū)內(nèi)地表植被生長初始期在第126天前后(即5月初)開始。本研究表明,研究區(qū)內(nèi)的大部分地區(qū)地表植被生長的平均初始期與正常年份相當，但部分地區(qū)有提前趨勢，提前10 d以上。Zhang等[18]認為,氣候變暖使得中國北部地區(qū)植被物候期在春季提前。趙秀蘭等[19]對東北地區(qū)近50年的氣溫變化研究表明，東北地區(qū)平均氣溫呈上升趨勢，春季升溫明顯，降水量呈逐年減少趨勢，秋季減少明顯。王宗明等[20]對東北地區(qū)植被對氣候因子的響應(yīng)研究表明，春季氣溫對植被的影響高于其余季節(jié)。春季氣溫變暖和穩(wěn)定降水是研究區(qū)內(nèi)地表植被生長初始期提前的原因之一。

研究區(qū)內(nèi)地表植被生長終期在第266天前后(9月下旬至10月上旬)結(jié)束。本研究表明，研究區(qū)內(nèi)大部分地區(qū)地表植被生長的平均終止期與正常年份相比有推后趨勢，平均推后3～4 d，但長春以西，松嫩平原以南地區(qū)的地表植被生長的平均終止期與正常年份持平。李波等[21]對吉林西部近10年氣溫變化研究表明，研究區(qū)內(nèi)秋季氣溫變化穩(wěn)定，波動較小，9月與10月的平均氣溫有增溫趨勢。暖干為主的秋季氣候，再加上該研究區(qū)內(nèi)近年來農(nóng)田灌溉設(shè)施的推廣應(yīng)用，使得研究區(qū)內(nèi)地表植被生長終期后延。

從10年地表植被生長期總天數(shù)平均變化情況看，研究區(qū)內(nèi)的地表植被生長期總天數(shù)存在增加趨勢，變化幅度為2～3 d，在長春以西，松嫩平原以南地區(qū)和嫩江流域以北地區(qū)有20 d以上的延長情況存在。沈柏竹等[22]和涂鋼等[23]對東北地區(qū)幾十年來的水熱變化特點研究表明，研究區(qū)的生長季氣溫變化穩(wěn)定，在10年尺度上基本沒有振蕩特征出現(xiàn)，同樣研究區(qū)內(nèi)地表降水變化趨勢穩(wěn)定，這一特點有利于研究區(qū)內(nèi)農(nóng)作物的規(guī)律性耕種。整體而言，本研究區(qū)域內(nèi)地表植被生長期總天數(shù)基本沒有變化，但在長春以西，松嫩平原以南地區(qū)和嫩江流域以北地區(qū)有延長趨勢。

本研究中的植被生長期是包括研究區(qū)農(nóng)作物在內(nèi)的所有地表植被生命活動的始期與終期，并非專指農(nóng)作物的生長特點，所以與農(nóng)業(yè)上的農(nóng)作物開始與終止期意義是不同的。實驗室得出的數(shù)據(jù)運算結(jié)果本身存在的實際應(yīng)用價值與意義需要在農(nóng)業(yè)生產(chǎn)領(lǐng)域進行驗證，判斷是否符合實際的地表生長狀況。另外，研究區(qū)地表植被除受氣候因子影響之外，人為因素也是不可忽略的，這一方面進行深入研究的意義較高。

致謝：在數(shù)據(jù)的整理過程中得到了侯旭、婁思文、姚希奎、杜保佳、趙英杰、孫繼梅、盧瑤、孫晶和王蔚悅等同學的大力幫助，特此感謝。

[1] 張峰,吳炳方,劉成林,等.利用時序植被指數(shù)監(jiān)測作物物候的方法研究[J].農(nóng)業(yè)工程學報,2004,20(1):155-159.

[2] 吳文斌,楊鵬,唐華俊,等.基于NDVI數(shù)據(jù)的華北地區(qū)耕地物候空間格局[J].中國農(nóng)業(yè)科學,2009,42(2):552-560.

[3] 陳效逑,王林海.遙感物候?qū)W研究[J].地理科學進展,2009,28(1):33-40.

[4] 張學霞,葛全勝,鄭景云.遙感技術(shù)在植物物候研究中的應(yīng)用綜述[J].地球科學進展,2003,18(4):534-544.

[5] 林忠輝,莫興國.NDVI時間序列諧波分析與地表物候信息獲取[J].農(nóng)業(yè)工程學報,2006,22(12):138-144.

[6] Karlsen S R,Solheim I,Beck P S A,etal.Variability of the start of the growing season in Fennoscandia,1982-2002[J].International Journal of Biometeorology,2007,51(6):513-524.

[7] 于信芳,莊大方.基于MODIS NDVI數(shù)據(jù)的東北森林物候期監(jiān)測[J].資源科學,2006,28(4):111-117.

[8] 余振,孫鵬森,孫世榮.中國東部南北樣帶主要植被類型物候期的變化[J].植物生態(tài)學報,2010,34(3):316-329.

[9] 宮攀,陳仲新.基于MODIS數(shù)據(jù)的東北地區(qū)植被物候參數(shù)提取[J].土壤通報,2009,40(2):213-217.

[10] 宋開山,劉殿偉,王宗明,等.1954年以來三江平原土地利用變化及驅(qū)動力[J].地理學報,2008,63(1):93-104.

[11] 楊淑杰,李曉東,王蘭蕓.吉林西部地表植被指數(shù)時空分布特點研究[J].東北林業(yè)大學學報,2010(5):223-334.

[12] 國志興,張曉寧,王宗明,等.東北地區(qū)植被物候?qū)夂蜃兓捻憫?yīng)[J].生態(tài)學雜志,2010,29(3):578-58.

[13] 楊永民,馮兆東,龍愛華,等.基于MODIS時序植被指數(shù)的植被物候空間格局[J].蘭州大學學報(自然科學版),2010,46(6):14-19.

[14] 武永峰,李茂松,李京.中國植被綠度期遙感監(jiān)測方法研究[J].遙感學報,2008,12(1):92-103.

[15] 陳懷亮.黃淮海地區(qū)植被覆蓋變化及其對氣候與水資源影響研究[D].北京:中國氣象科學研究院,2007:39-59.

[16] Stockli R,Vidale P L.European plant phenological and climate as seen in a 20-year AVHRR land-surface parameter dataset[J].International Journal of Remote Sensing,2004,25(17):3303-3330.

[17] Chen J,Jonsson P,Tamura M,etal.A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay fliter[J].Remote Sensing of Environment,2004,91:332-344.

[18] Zheng J Y,Ge Q S,Hao Z X,etal.Impacts of climate warming on plants phenophases in China for the last 40 years[J].Chinese Science Bulletin,2002,21:1826-1831.

[19] 趙秀蘭.近 50 年中國東北地區(qū)氣候變化對農(nóng)業(yè)的影響[J].東北農(nóng)業(yè)大學學報,2010,41(9):144-149.

[20] 王宗明,國志興,宋開山,等.中國東北地區(qū)植被NDVI對氣候變化的響應(yīng)[J].生態(tài)學雜志,2009,28(6):1041-1048.

[21] 李波,劉晶如,李曉東.白城市近35年來的氣溫變化趨勢及對農(nóng)業(yè)生產(chǎn)的影響[J].白城師范學院學報,2009,6(3):71-74.

[22] 沈柏竹,廉毅,楊涵洧,等.中國東北夏季溫度年代際變化特征[J].地理科學,2012,32(6):94-100.

[23] 涂鋼,劉波,王淑瑜.基于陸面模式NCAR/CLM3.5的東北區(qū)域生長季地表干濕狀況時空分布特征模擬研究[J].地理科學,2012,32(6):101-106.