基于MODIS的寧夏植被覆蓋時空變化特征分析

吳霞 張永宏 王長軍 樊麗琴 李磊

摘要：以2000—2015年寧夏回族自治區NDVI數據及高程數據為基礎，采用趨勢分析、變異系數、GIS空間統計等方法，分析寧夏近年來植被覆蓋時空分布及變化特征。結果表明，從年際變化情況看，近16年來寧夏植被覆蓋整體呈緩慢增長趨勢，NDVI值增速為0.008/10年，其中固原市增速最高;從月際變化情況看，寧夏植被覆蓋表現為冬春季低，夏秋季高，呈先增加后減少的趨勢。近16年來寧夏全區年均NDVI值為0.37，植被覆蓋整體處于中低水平，中低覆蓋地區面積占比達到70.39%;從地域分布來看，表現為南北植被覆蓋高，中部地區植被覆蓋低;高程方面，植被覆蓋隨著高程增加呈現“降低—升高—降低”的變化趨勢，其中2 400～<2 700 m海拔范圍內植被覆蓋最高。近年來寧夏植被覆蓋整體維持穩定且有所改善，植被覆蓋呈改善趨勢的區域占比為50.17%，其中南部黃土高原地區植被覆蓋改善最明顯，城鎮地區植被退化較嚴重。寧夏植被覆蓋變異系數介于0.002～0.651，整體呈現中高波動態勢，中高波動區域占比達到57.88%。

關鍵詞：MODIS;NDVI;寧夏;植被覆蓋;時空變化特征

中圖分類號： S127文獻標志碼： A文章編號：1002-1302（2021）02-0204-05

收稿日期：2020-02-06

基金項目：寧夏自治區自然科學基金（編號：2019AAC03290）;寧夏青年科技人才托舉項目（編號：TJGC2018055）;寧夏農林科學院科技創新全產業鏈項目（編號：NKYZ-16-0907）。

作者簡介：吳霞（1984—），女，寧夏吳忠人，碩士，助理研究員，主要從事土地資源研究。E-mail：wuxia-xia@163.com。

通信作者：張永宏，研究員，主要從事作物栽培與耕作研究。E-mail：zyh8401@163.com。

植被指地球表面某一地區所覆蓋的植物群落，它與氣候、土壤、地形及水資源等自然環境要素密切相關，是生態系統的重要組成部分[1-2]。植被作為研究全球氣候變化的“指示器”，在全球能量交換、生物地球化學和水文循環方面起著重要的作用，能直接反映生態環境的總體狀況，因此通過監測植被覆蓋狀況來掌握生態環境變化特征具有重要的科學意義和現實價值[3-5]。歸一化植被指數（normal difference vegetation index，NDVI）別稱標準化植被指數，在使用遙感圖像進行植被研究以及植物物候研究中得到廣泛應用，它是植物生長狀態以及植被空間分布密度的最佳指示因子，與植被分布密度呈線性相關[6-7]。隨著遙感技術的發展，越來越多國內外學者采用NDVI數據對植被覆蓋變化進行科學研究，從而探明研究區植被覆蓋的空間分布特征及其隨時間延長的演變規律，以及地形、氣候、人類活動等因素對植被覆蓋的影響[8-11]。寧夏生態本底相對脆弱，屬于典型的生態脆弱區，多年來受干旱和水資源短缺的困擾，植被覆蓋度低，水土流失和土地荒漠化現象嚴重。因此，研究該區地表植被覆蓋時空變化，對當地的生態建設與環境保護具有很好的指示意義。目前，已有學者對寧夏植被覆蓋變化情況進行分析[12-13]，范錦龍等使用PATHFINDER 8 km NDVI數據進行研究，發現寧夏全區植被具有明顯的1年1季的季節特征[14];裴志方等基于MODIS（中分辨率成像光譜儀）NDVI數據進行研究，發現寧夏植被覆蓋整體向良好方向發展[15]。但是目前有關寧夏植被覆蓋穩定性方面的研究尚未見相關報道，因此本研究采用2000—2015年寧夏NDVI數據，在分析寧夏植被覆蓋時空演變規律的基礎上，對植被覆蓋穩定性進行定量化分析，為進一步分析寧夏植被動態變化機制、預測其變化趨勢提供科學依據。

1研究區概況與研究方法

1.1研究區概況

寧夏回族自治區位于我國西北部，黃河中上游地區，疆域輪廓南北長、東西短，呈“十”字形，南北相距約456 km，東西相距約250 km，總面積為 6.64萬km2 [16]。其地勢南高北低，西部高差較大，東部起伏較緩，從西面、北面至東面，由騰格里沙漠、烏蘭布和沙漠和毛烏素沙地相圍，南面與黃土高原相連[17]。寧夏屬典型的大陸性氣候，為溫帶半干旱區和半濕潤地區，年均氣溫5～9 ℃，年降水量180～680 mm，均由南向北遞減，降水多集中于6—9月，且年際變化率大。全區平均年水面蒸發量 1 250 mm，變幅在800～1 600 mm，年均徑流量為 9.493億m3，年均徑流深18.3 mm，人均水量均遠低于全國和黃河流域平均值，干旱威脅嚴重[18-19]。受水熱條件尤其是水分因素的制約，寧夏植被自南向北呈現森林草原—干草原—荒漠草原—草原化荒漠的水平分布規律，其中以草原和荒漠草原為主。2000年以來，寧夏依托三北防護林、天然林保護、退耕還林等國家重大林業工程，組織實施了封山禁牧、防沙治沙、濕地保護、生態修復等重點工程，植被覆蓋率不斷提升，山川面貌不斷美化優化，生態環境得到了恢復和進一步提升[20-21]。

1.2數據來源

本研究所需NDVI數據及數字高程模型（DEM）數據均來源于中國科學院資源環境數據云平臺，NDVI數據集是基于連續時間序列的SPOT/VEGETATION NDVI衛星遙感數據，采用最大值合成法生成年度及月度最大NDVI數據，該數據集的空間分辨率為1 km，時間分辨率為10 d;DEM數據是基于最新的SRTM V4.1數據經整理拼接生成，空間分辨率為90 m。本研究選取經裁剪得到的2000—2015年寧夏NDVI年度和月度數據以及DEM數據，對寧夏植被覆蓋時空變化特征進行分析。

1.3研究方法

1.3.1趨勢分析法一元線性回歸分析可以模擬出影像中每個柵格的變化趨勢，綜合反映研究區域NDVI時空格局演變特征[22]。本研究對2000—2015年寧夏NDVI數據在像元尺度上進行一元線性回歸分析，計算NDVI的變化趨勢，具體計算過程如下。

slope=n×∑ni=1i×NDVIi-（∑ni=1i）×（∑ni=1NDVIi）n×∑ni=1i2-（∑ni=1i）2。（1）

式中：slope表示像元NDVI回歸方程斜率;n表示研究期總年數，本研究中n=16 年;i表示年序號;NDVIi表示第i年NDVI的最大值。若slope>0，表示NDVI隨時間變化呈增長趨勢，slope值越大，表明增長越快;若slope<0，則表示NDVI隨時間延長呈下降趨勢，slope絕對值越大表明下降越明顯[23]。

1.3.2變異系數法變異系數（CV）別稱標準差率，是衡量資料中各觀測值變異程度的一個統計量，本研究通過變異系數（CV）反映NDVI隨時間序列變化的波動情況，其計算過程如下。

CV=12n∑ni=1（NDVIi-NDVI）2NDVI×100%。（2）

式中：n表示研究期總年數，本研究中n = 16年;i表示年序號;NDVIi表示第i年NDVI的最大值;NDVI表示研究區內16年的NDVI平均值。CV值越大，表明NDVI年際變化越大，時序波動較大;反之，則表明NDVI年際變化不大，時序較穩定[24-25]。

2結果與分析

2.1寧夏植被覆蓋時間變化特征

2.1.1NDVI年際變化特征由圖1可知，2000—2015年寧夏植被覆蓋呈緩慢增長趨勢，NDVI值增速僅為0.008/10年，最大值出現在2013年，為047;最小值出現在2000年，為0.29。各地區植被覆蓋均呈增長趨勢，固原市植被覆蓋增長最明顯，NDVI值從2000年的0.41增加到2015年的0.58，增長率為41.46%;中衛市植被覆蓋增長最緩慢，16年間NDVI值僅增長21.74%;銀川市和石嘴山市植被覆蓋均呈增長趨勢，增速基本與全區水平持平。

2.1.2NDVI月際變化特征由圖2可知，1年內寧夏植被覆蓋呈先增加再減少的趨勢，這與干旱半干旱地區氣候特征和降水量月際分布密切相關。寧夏1—12月NDVI值介于0.11～0.37，1—4月屬于冬春季節，干旱少雨，且耕地處于冬閑狀態，NDVI值較低且增長緩慢，最大值僅為4月的0.138;5—9月是寧夏的多雨季節，植被生長迅速，NDVI值快速增長，至7月達到最大，為0.374;10—12月NDVI快速下降，下降率達到48.28%。

2.2寧夏植被覆蓋空間分布特征

2.2.1NDVI空間分布特征利用ArcGIS的像元統計工具計算得到寧夏2000—2015年平均年度NDVI空間分布圖（圖3），并根據寧夏實際情況將植被覆蓋劃分為極低覆蓋度（NDVI<0.15）、低覆蓋度（0.15≤NDVI<0.30）、中低覆蓋度（0.30≤NDVI<0.45）、中覆蓋度（0.45≤NDVI<0.60）、高覆蓋度（NDVI≥0.60）5類。由圖3可知，寧夏平均年度NDVI值介于0.10～0.86，全區平均年度NDVI值為

0.37，總體處于中低覆蓋水平。由表1可知，寧夏極低覆蓋度的區域面積為693 km2，占1.33%，主要集中分布于中衛市西北部騰格里沙漠邊緣;低覆蓋度的區域面積為23 176 km2，占44.60%，主要分布于寧夏中部干旱少雨地區，包括中衛市海原縣、吳忠市同心縣、鹽池縣等;中低覆蓋度和中覆蓋度的區域面積分別為12 711、9 695 km2，合計占比為4311%，主要分布在位于黃土高原北部的固原地區;高覆蓋度的區域面積為5 694 km2，占10.96%，分布于固原市六盤山森林地區以及沿黃河灌溉地帶。

2.2.2不同高程NDVI變化特征寧夏海拔高度介于1 086～345 7 m，按照300 m等間距劃分為8個區間。由圖4可知，隨著海拔高度的增加，寧夏NDVI值呈現“降低—升高—降低”的變化態勢，NDVI低值區出現在海拔1 200～<1 800 m之間，主要分布于寧夏中部地區，該地區降水量稀少、蒸發量大且缺乏灌溉水源，不利于植被生長;NDVI高值區出現在海拔2 400～<3 000 m之間，分布于寧夏南部六盤山地區，該區域降水量充沛且溫度適宜，植被生長旺盛。

2.3寧夏植被覆蓋時空格局演變

2.3.1空間趨勢特征采用一元線性回歸趨勢分析法計算得到2000—2015年寧夏植被覆蓋變化趨勢，參照NDVI變化類型將趨勢線斜率（slope）劃分為顯著退化（slope<-0.01）、輕微退化（-0.01≤slope<-0.005）、基本不變（-0.005 ≤slope<0005）、輕微改善（0.005 ≤slope<0.01）和顯著改善（slope ≥ 0.01）5類， 并繪制NDVI變化趨勢空間

分布圖 （圖5）。由表2可知，2000—2015年寧夏植被覆蓋整體維持穩定且有所好轉，植被狀況呈現顯著改善和輕微改善的面積分別為11 546 、14 524? km2，合計占比為50.17%，主要分布于引黃灌溉地區和南部黃土高原地區;植被覆蓋基本不變區域面積為24 739 km2，占比為47.60%，分布于中衛市、吳忠市等中部干旱少雨地區;植被覆蓋呈退化趨勢的區域面積合計為1 160 km2，占2.23%，呈點狀集中分布于城鎮地區，其中銀川市周圍植被退化最明顯。

2.3.2空間波動特征經計算得到2000—2015年寧夏NDVI變異系數（CV）， 參照自然間斷法劃分為低波動、較低波動、中等波動、較高波動和高波動5類，并繪制變異系數空間分布圖（圖6）。據統計，寧夏16年間植被覆蓋變異系數介于0.002～0.651，整體呈現中等波動和高波動變化態勢，區域波動差異性較大。低波動和較低波動區域面積合計占比為4212%，主要分布于黃河沿岸灌水便利地區以及南部山區降水相對充沛區域;中高波動區域主要分布于中部干旱帶、荒漠地區以及北部賀蘭山地區等，其中植被覆蓋較高波動和高波動區域面積合計為10 721 km2，占比達到20.63%，與寧夏降水量年度波動情況密切相關。

3結論

2000—2015年寧夏植被覆蓋隨時間變化特征、空間分布情況、發展趨勢以及穩定性方面表現如下：第一，從年際變化來看，寧夏植被覆蓋呈緩慢增長趨勢，NDVI值增速為0.008/10年，其中固原地區植被覆蓋增長最快，增長率達到41.46%;從月際變化來看，植被覆蓋冬春季低，夏秋季高，NDVI最低值出現在2月，為0.11，最高值出現在7月，為037。第二，近16年寧夏平均年度NDVI值為037，整體處于中低覆蓋水平，中低覆蓋地區面積占比達到70.39%，地域分布上呈現南北兩頭植被覆蓋較高，中部地區植被覆蓋偏低的態勢;植被覆蓋隨著高程增加呈現“降低—升高—降低”的變化趨勢，其中2 400～2 700 m海拔范圍內是寧夏水熱條件較好的地區，植被覆蓋最高。第三，從植被覆蓋發展趨勢來看，2000—2015年寧夏植被覆蓋整體維持穩定且有所改善，50.17%的地區植被覆蓋呈現改善趨勢，其中南部黃土高原地區植被覆蓋改善最明顯，城鎮地區特別是銀川市周圍植被退化較嚴重。第四，植被覆蓋穩定性方面，寧夏植被覆蓋變異系數介于0.002～0.651，整體呈現中高波動變化態勢，中高波動區域占比達到57.88%，區域差異表現為“南北波動小，中部波動大”。

本研究利用NDVI數據對寧夏2000—2015年植被覆蓋時空變化情況進行具體分析，發現寧夏植被覆蓋在空間分布和時間變化趨勢方面均與寧夏地形地貌、氣候特征、降水量分布以及人類活動等密切相關。本研究并未就溫度、降水、經濟以及人口等因素對植被覆蓋的影響作深入探討，在今后的研究中有必要綜合考慮多因素對寧夏植被覆蓋的影響，明確各因素對植被覆蓋的影響因子，進一步預測植被覆蓋的發展趨勢。

參考文獻：

[1]賈艷紅，趙傳燕，南忠仁. 西北干旱區黑河下游植被覆蓋變化研究綜述[J]. 地理科學進展，2007，26（4）：64-74.

[2]張春森，胡艷，史曉亮. 基于AVHRR和MODIS NDVI數據的黃土高原植被覆蓋時空演變分析[J]. 應用科學學報，2016，34（6）：702-712.

[3]杜艷秀. 四川省植被覆蓋時空變化及驅動力定量研究[D]. 成都：成都理工大學，2016.

[4]史海萍. 山西氣候變化及植被生態系統的響應研究[D]. 蘭州：蘭州大學，2006.

[5]Song W，Mu X，Ruan G，et al.Estimating fractional vegetation cover and the vegetation index of bare soil and highly dense vegetation with a physically based method[J]. International Journal of Applied Earth Observations and Geoinformation，2017，58：168-176.

[6]李卓，劉淑亮，孫然好，等. 黃淮海地區耕地復種指數的時空格局演變[J]. 生態學報，2018，38（12）：4454-4460.

[7]Carlson T N，Ripley D A. On the relation between NDVI，fractional vegetation cover，and leaf area index[J]. Remote Sensing of Environment，1997，62（3）：241-252.

[8]Ghazal M，Khalil Y A，Hajjdiab H. UAV-based remote sensing for vegetation cover estimation using NDVI imagery and level sets method[C]//International Symposium on Signal Processing & Information Technology，2016.

[9]王傳勝，李秋秋. 2000年以來西秦嶺地區植被覆蓋變化特征——以隴南市為例[J]. 水土保持研究，2016，23（2）：308-312.

[10]穆少杰，李建龍，陳奕兆，等. 2001—2010年內蒙古植被覆蓋度時空變化特征[J]. 地理學報，2012，67（9）：1255-1268.

[11]劉哲，邱炳文，王壯壯，等. 2001—2014年間黃土高原植被覆蓋狀態時空演變分析[J]. 國土資源遙感，2017，29（1）：192-198.

[12]簡太敏，王小麗，沙金燕. 關于寧夏植被覆蓋變化與氣候因子相關性分析[J]. 北京農業，2015（34）：136-136，137.

[13]強英云，張慧娟. 1987—2010年銀川平原植被蓋度變化規律研究[J]. 寧夏工程技術，2016，15（4）：398-400.

[14]范錦龍，張曉煜. 寧夏近18年植被動態變化遙感時序分析[J]. 干旱區研究，2009，26（1）：53-58.

[15]裴志方，楊武年，吳彬，等. 2000—2016年寧夏植被覆蓋景觀格局遙感動態分析[J]. 水土保持研究，2018，25（1）：215-219.

[16]牛志男. 從南到北走“塞上”——寧夏回族自治區采訪手記[J]. 中國民族，2008（9）：28-31.

[17]馬啟智. 《走進寧夏》序言[J]. 寧夏史志，2004（2）：5-6.

[18]丁永紅，王文，陳曉光，等. 寧夏近44年暴雨氣候特征和變化規律分析[J]. 高原氣象，2007，26（3）：630-636.

[19]徐利崗，湯英，杜歷，等. 寧夏降水空間分布格局及其多時間尺度變化特征分析[J]. 安徽農業科學，2011，39（34）：21303-21306，21310.

[20]吳素琴，溫淑萍，楊瑞全，等. 寧夏干草原、荒漠草原及其治理[J]. 寧夏農林科技，2005（3）：35-37.

[21]陳華轅. 美麗寧夏綠色隨行——自治區人大常委會審議全區生態林業建設情況報告側記[J]. 寧夏人大，2015（4）：20-21.

[22]李雙雙，延軍平，萬佳. 近10年陜甘寧黃土高原區植被覆蓋時空變化特征[J]. 地理學報，2012，67（7）：960-970.

[23]曹磊. 江蘇省植被NDVI動態變化及其與氣候因子的關系[D]. 南京：南京農業大學，2014.

[24]李杰，張軍，劉陳立，等. 基于MODIS-NDVI的云南怒江流域植被覆蓋時空變化特征研究[J]. 林業資源管理，2018（4）：90-99.

[25]Barbosa H A，Kuma T V L，Ertürk A G. Using the satellite-derived NDVI-OLR feedbacks over West Sahel Africa to assess Land-Atmosphere responses to environmental change[C]//American Institute of Physics，2009：357-360.李忠良，張余慶，陳成. 西藏農業產值影響因子敏感性分析[J]. 江蘇農業科學，2021，49（2）：209-217.