王麗霞 張珈瑋 張雙成 隋立春 楊耘 劉招

摘要??以陜西省作為研究區(qū)域，首先基于MODIS NDVI中國月合成數(shù)據(jù)計算研究區(qū)2000、2005、2010和2015年度的植被覆蓋度并進(jìn)行等級劃分;而后利用CA-Markov模型，以植被覆蓋度等級作為元胞類型，計算不同時期各植被覆蓋度等級的轉(zhuǎn)移矩陣，由此模擬2010和2015年的NDVI值;對比模擬NDVI結(jié)果和原始影像數(shù)據(jù)，評價模擬精度，并預(yù)測2020年NDVI的空間分布狀況。結(jié)果表明，利用CA-Markov模型進(jìn)行植被覆蓋空間分布的模擬，得到2010和2015年模擬結(jié)果的Kappa系數(shù)分別為0.797 5、0.853 2，符合精度要求，可以用于植被覆蓋空間分布的預(yù)測。陜西省植被覆蓋存在明顯的空間差異性，呈現(xiàn)出陜北—關(guān)中—陜南地區(qū)植被覆蓋度逐漸遞增的緯度地帶性規(guī)律。2000—2020年各級植被覆蓋區(qū)均有向更高一級植被覆蓋區(qū)變化的趨勢?？傮w上植被覆蓋度上升趨勢明顯，尤以陜北地區(qū)北部變化最為顯著，關(guān)中地區(qū)和陜南秦巴山區(qū)植被覆蓋度增幅較小。

關(guān)鍵詞??歸一化植被指數(shù);植被覆蓋度;CA-Markov模型;模擬預(yù)測

中圖分類號??P??237文獻(xiàn)標(biāo)識碼??A

文章編號??0517-6611（2020）04-0053-04

doi：10.3969/j.issn.0517-6611.2020.04.016

開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

Simulation and Prediction of Vegetation Coverage in Shaanxi Province Based on CA-Markov Model

WANG Li-xia1， ZHANG Jia-wei2， ZHANG Shuang-cheng1 et al??（1.School of Geology Engineering and Geomatics， Changan University， Xian， Shaanxi 710054;2. School of Earth Science and Resources， Changan University， Xian， Shaanxi 710054）

Abstract??Taking Shaanxi Province as the research area. Firstly， based on MODIS NDVI China monthly synthetic data， the vegetation coverage in the region of 2000， 2005， 2010 and 2015 was calculated and graded. Then using the CA-Markov model and taking the vegetation coverage level as the cell type， the transfer matrix of various vegetation coverage levels in different periods was calculated， and NDVI in 2010 and 2015 was simulated. Comparing simulated NDVI with the original image data， the accuracy of simulation was evaluated and the spatial distribution of NDVI in 2020 was predicted. The results indicated that using CA-Markov model to simulate the spatial distribution of vegetation coverage， the Kappa coefficient of the simulation results in 2010 was 0.797 5， and in 2015 was 0.853 2， which met the accuracy requirements and could be used to predict the spatial distribution of vegetation coverage.There were obvious spatial differences in vegetation cover of Shaanxi Province， showing the latitudinal zonal law of increasing vegetation coverage from the northern of Shaanxi to Guanzhong to the southern of Shaanxi area. Between 2000 and 2020， vegetation coverage areas at all levels had a tendency to change to higher-level vegetation coverage areas. In general， the vegetation coverage increased obviously， especially in the northern part of northern Shaanxi. The vegetation coverage in Guanzhong and Qinba Mountains in southern Shaanxi increased little.

Key words??Normalized difference vegetation index;Vegetation coverage;CA-Markov model;Simulation and prediction

植被是陸地生態(tài)系統(tǒng)的主體，起著連接土壤、大氣和水分等要素，參與生態(tài)系統(tǒng)物質(zhì)循環(huán)的重要作用。植被生長和覆蓋狀況存在明顯的時空差異性，在一定程度上反映了氣候變化的趨勢，是生態(tài)環(huán)境研究的核心問題之一[1]。遙感技術(shù)以其信息量大、觀測范圍廣、精度高、速度快以及實時性和動態(tài)性等特點為開展植被監(jiān)測和評估提供了有效的手段[2]。目前研究主要體現(xiàn)為通過遙感影像提取各類植被指數(shù)，在不同空間尺度下對植被覆蓋進(jìn)行長時間序列的動態(tài)監(jiān)測[3-8]。其中歸一化植被指數(shù)（normalized difference vegetation index，NDVI）在生態(tài)環(huán)境評價、旱情預(yù)警監(jiān)測、植被生長力評估和土地利用評價等方面應(yīng)用最為廣泛。CA-Markov模型綜合了馬爾科夫模型（Markov）長期預(yù)測的優(yōu)點以及元胞自動機(jī)（CA）模擬復(fù)雜系統(tǒng)空間變化的能力，在分析和模擬土地利用變化方面取得了良好的效果[9-15]。

陜西省位于105°29′～111°15′E、31°42′～39°35′N，地處我國西北地區(qū)邊緣，與甘肅、寧夏、山西等8個省、區(qū)、市毗鄰，總面積為20.58 萬km2。陜西省的地貌形態(tài)復(fù)雜，中部低、南北高，南北方向分化較為明顯，由西向東地勢也有明顯傾斜。由于地貌和大氣環(huán)流的影響，陜西省季風(fēng)氣候的性質(zhì)明顯，全省年平均氣溫13.7 ℃，年平均降水量為676.4 mm，降雨量集中在7—10月，從南向北逐漸遞減，氣候類型從陜南的濕潤性氣候分別過渡到關(guān)中陜北的半濕潤性和半干旱性氣候。陜西省地貌類型多樣，植被覆蓋存在明顯的地域差異性[16]，由北向南依次為陜北防護(hù)林生態(tài)區(qū)、關(guān)中平原農(nóng)業(yè)生態(tài)區(qū)以及陜南落葉闊葉和常綠闊葉混交林生態(tài)區(qū)。因此探究該區(qū)植被覆蓋的時空分布規(guī)律，對于陜北黃土高原生態(tài)修復(fù)、關(guān)中平原農(nóng)業(yè)種植標(biāo)準(zhǔn)化建設(shè)以及陜南秦巴山區(qū)天然林木保護(hù)具有重要意義。

1??資料與方法

1.1??影像數(shù)據(jù)來源及處理

該研究采用2000—2015年MODIS NDVI中國月合成數(shù)據(jù)，由地理空間數(shù)據(jù)云平臺下載（http：//www.gscloud.cn）。中國500 m NDVI月合成產(chǎn)品是由MODND1D計算得到，計算方法為取月內(nèi)每天最大值?？臻g分辨率為500 m。除2015年10月數(shù)據(jù)缺失外，總共191期NDVI數(shù)據(jù)。為研究陜西省植被覆蓋變化特征，統(tǒng)計2000—2015年NDVI月平均值，并對每年1—12月的數(shù)據(jù)進(jìn)行合成，得到年際NDVI數(shù)據(jù)。將影像的每一個柵格作為一個元胞進(jìn)行分析，其大小為500 m×500 m，使用的濾波器為5×5濾波器。

1.2??其他數(shù)據(jù)來源及處理??中國1∶400萬基礎(chǔ)數(shù)據(jù)、植被覆蓋度劃分標(biāo)準(zhǔn)、模型計算公式等資料來自其他地理信息網(wǎng)站。中國1∶400萬基礎(chǔ)數(shù)據(jù)為Krasovsky_1940坐標(biāo)，將其更換為MODIS NDVI數(shù)據(jù)所用的WGS84坐標(biāo)系，以進(jìn)行地理配準(zhǔn)。并利用陜西省界SHP文件對原數(shù)據(jù)進(jìn)行裁剪操作。

1.3??植被覆蓋度計算及分級

NDVI是植物生長、覆蓋、生物量和植被種類情況等的綜合反映，與植被分布密度呈線性相關(guān)[17]。為詳細(xì)統(tǒng)計植被覆蓋的空間情況，可由NDVI提取植被覆蓋度。植被覆蓋度的計算方法為：

FVC=NDVI-NDVIminNDVImax-NDVImin×100%（1）

ArcGIS的柵格計算器可對數(shù)據(jù)進(jìn)行運算，由年NDVI數(shù)據(jù)得出年植被覆蓋數(shù)據(jù)。

根據(jù)水利部1996年頒布的《土壤侵蝕分類分級標(biāo)準(zhǔn)》（SL190—96），按植被覆蓋度大小分為5個等級區(qū)：≤30%為低度植被覆蓋區(qū);>30%～45%為較低植被覆蓋區(qū);>45%～60%為中度植被覆蓋區(qū);>60%～75%為較高植被覆蓋區(qū);>75%為高度植被覆蓋區(qū)。將陜西省2000、2005、2010、2015年NDVI數(shù)據(jù)分為上述5個等級，結(jié)果如圖1所示。

1.4??研究方法

自然界許多事物都沒有固定而精確的變化模式，受多種因素的干擾，常表現(xiàn)出很強的不確定性。為了研究方便，常采用一些隨機(jī)過程的灰箱模型來研究變化方式，其中Markov模型以概率論作為理論基礎(chǔ)，對隨機(jī)事件的變化規(guī)律和趨勢進(jìn)行分析，可以較好地分析這個過程，并預(yù)測短期的變化結(jié)果[18]。但因為Markov模型無后效性，在空間預(yù)測上就遇到了較大問題。CA-Markov模型綜合了Markov模型對于長時間序列的預(yù)測和元胞自動機(jī)（CA）模擬復(fù)雜系統(tǒng)空間變化的能力，能夠更加準(zhǔn)確地從時間和空間上模擬和預(yù)測植被覆蓋的空間分布情況。

此次試驗將進(jìn)行2010、2015年植被覆蓋的模擬和2020年植被覆蓋的預(yù)測。轉(zhuǎn)移概率矩陣的計算在IDRISI軟件中進(jìn)行。IDRISI支持ASCII格式文件的導(dǎo)入，于ArcGIS中將2000、2005、2010、2015年未進(jìn)行重分類的植被覆蓋數(shù)據(jù)轉(zhuǎn)換為ASCII格式，并導(dǎo)入IDRISI，即可進(jìn)行后續(xù)操作。重分類步驟于IDRISI中進(jìn)行，按植被覆蓋度分級標(biāo)準(zhǔn)進(jìn)行分級。

首先需要計算Markov轉(zhuǎn)移矩陣，以獲得5年期的轉(zhuǎn)移方向和轉(zhuǎn)移程度。使用IDRISI中的Markov模塊，輸入2000、2005年數(shù)據(jù)，將時間間隔與推演間隔均設(shè)置為5年，比例誤差設(shè)置為0.1，獲得2000—2005年的Markov轉(zhuǎn)移矩陣（表1）。同理可獲得2005—2010和2010—2015年的轉(zhuǎn)移矩陣。

2??結(jié)果與分析

2.1??植被覆蓋模擬

基于2000和2005年區(qū)域NDVI數(shù)據(jù)，在IDRISI平臺下，獲得2000—2005年各植被覆蓋等級概率轉(zhuǎn)移矩陣，運用CA-Markov模型，生成適宜性圖譜，采用5×5濾波器進(jìn)行運算，分別可以得到2010和2015年的模擬數(shù)據(jù)。將2010和2015年植被覆蓋的解譯結(jié)果和模擬結(jié)果進(jìn)行對比（圖2），

模擬結(jié)果表明，陜西省植被覆蓋存在明顯的緯度地帶性，呈現(xiàn)陜北—關(guān)中—陜南地區(qū)植被覆蓋度逐漸遞增的規(guī)律。陜北地區(qū)處于風(fēng)沙區(qū)，降水較少，年平均溫度低，植被覆蓋度較低;關(guān)中平原地區(qū)以農(nóng)田作物為主，雨熱條件次之，其植被覆蓋度屬中等水平;陜南地區(qū)分布有秦嶺山地和大巴山區(qū)，降水豐沛、氣溫較高，其植被覆蓋度對應(yīng)較高。

從各級植被覆蓋區(qū)空間分布特征來看，低度植被覆蓋區(qū)主要分布在陜北地區(qū)，關(guān)中地區(qū)城市坐落地及其輻射地區(qū)亦有極少分布;較低植被覆蓋區(qū)分布在陜北地區(qū)北部和中部;中度植被覆蓋區(qū)分布在陜北地區(qū)南部和關(guān)中地區(qū)中部;較高植被覆蓋區(qū)主要分布在陜北地區(qū)南部，關(guān)中地區(qū)東部和陜南地區(qū)東部;高度植被覆蓋區(qū)分布在陜北地區(qū)南部和陜南地區(qū)。

從時間變化特征來看，2000—2015年陜西省植被覆蓋度整體呈現(xiàn)顯著增加的趨勢，其中陜北地區(qū)變化最為顯著：低度植被覆蓋區(qū)轉(zhuǎn)化為較低植被覆蓋區(qū)，較低植被覆蓋區(qū)面積占比增加，由陜北地區(qū)中部逐漸擴(kuò)展到東北部;而低度植被覆蓋區(qū)域面積逐漸縮減，到2015年主要集中在陜北地區(qū)西北部長城沿線風(fēng)沙區(qū)—毛烏素沙地;中度植被覆蓋區(qū)轉(zhuǎn)變?yōu)檩^高植被覆蓋區(qū)，并由陜北地區(qū)南部轉(zhuǎn)移到中部;較高植被覆蓋區(qū)主要分布在陜北地區(qū)南部，逐漸轉(zhuǎn)化為高度植被覆蓋區(qū)。對于關(guān)中地區(qū)和陜南地區(qū)來說，中度植被覆蓋區(qū)和較高植被覆蓋區(qū)都有向高度植被覆蓋區(qū)域明顯變化的跡象?？傮w上較低植被覆蓋區(qū)和高度植被覆蓋區(qū)的面積占比逐年增加，而低度植被覆蓋區(qū)、中度植被覆蓋區(qū)和較高植被覆蓋區(qū)面積逐漸減少。陜北地區(qū)北部植被覆蓋度增加非常顯著，關(guān)中地區(qū)和陜南秦巴山區(qū)植被覆蓋增幅較小。

2.2??植被覆蓋預(yù)測??為了評價模型的實用性，對模擬的精度進(jìn)行了分析，得到2010和2015年模擬結(jié)果的Kappa系數(shù)分別為0.797 5、0.853 2，精度較高符合要求。因此基于此參數(shù)設(shè)置，運用CA-Markov模型，利用2010—2015年的植被覆蓋等級轉(zhuǎn)移矩陣，得到2020年預(yù)測的結(jié)果，如圖3所示，其植被覆蓋度面積分布如表2所示。

預(yù)測結(jié)果表明，陜西省在2015—2020年各級植被覆蓋區(qū)面積發(fā)生了較大變化，總體上植被覆蓋度上升明顯。低度植被覆蓋區(qū)面積減少了2.35%，向較低植被覆蓋區(qū)轉(zhuǎn)變;而大部分較低植被覆蓋區(qū)轉(zhuǎn)化為中度植被覆蓋區(qū)，面積減少2.9%，這3級植被覆蓋區(qū)主要分布在陜北地區(qū)北部;陜北地區(qū)南部大部分斑塊狀的中度植被覆蓋區(qū)向較高植被覆蓋區(qū)轉(zhuǎn)變，關(guān)中地區(qū)少數(shù)較高植被覆蓋區(qū)向高度植被覆蓋區(qū)擴(kuò)張;而陜南地區(qū)基本沒有明顯變化。因此較高植被覆蓋區(qū)增幅最大，增加了3.15%，其次是高度植被覆蓋區(qū)，增加了1.2%。

3??結(jié)論與討論

該研究利用2000—2015年陜西省NDVI月時序數(shù)據(jù)，合成年度NDVI數(shù)據(jù)，在IDRISI軟件中采用CA-Markov模型對2000、2005、2010、2015年4期陜西省植被覆蓋數(shù)據(jù)進(jìn)行了運算分析，模擬結(jié)果與原影像分級結(jié)果具有較高的一致性，以此預(yù)測2020年陜西省植被覆蓋的空間分布，得出以下結(jié)論：

（1）陜西省植被覆蓋存在明顯的空間差異性，呈現(xiàn)出陜北—關(guān)中—陜南地區(qū)植被覆蓋度逐漸遞增的緯度地帶性規(guī)律，這與當(dāng)?shù)氐臍夂驐l件有著密不可分的關(guān)系。陜北地區(qū)各級植被覆蓋區(qū)均有分布，關(guān)中地區(qū)分布有低度植被覆蓋區(qū)、中度植被覆蓋區(qū)和較高植被覆蓋區(qū)，而陜南地區(qū)主要分布有較高植被覆蓋區(qū)和高度植被覆蓋區(qū)。

（2）利用CA-Markov模型進(jìn)行植被覆蓋空間分布的模擬，得到2010和2015年模擬結(jié)果的Kappa系數(shù)分別為0.797 5、0.853 2，符合精度要求，可以用于植被覆蓋空間分布的預(yù)測。

（3）在2000—2015年較長時間尺度上，陜西省植被覆蓋整體呈現(xiàn)顯著增加的趨勢。其中陜北地區(qū)北部植被覆蓋度變化最為明顯，各級植被覆蓋區(qū)均有向更高一級植被覆蓋區(qū)變化的趨勢;對于關(guān)中地區(qū)和陜南地區(qū)來說，中度植被覆蓋區(qū)和較高植被覆蓋區(qū)都有向高度植被覆蓋區(qū)域明顯轉(zhuǎn)變的跡象，植被覆蓋度增幅較小。

（4）預(yù)測到2020年，各級植被覆蓋區(qū)面積將總體呈現(xiàn)上升趨勢。陜北地區(qū)北部低度植被覆蓋區(qū)與較低植被覆蓋區(qū)面積均減少，分別轉(zhuǎn)化為較低植被覆蓋區(qū)和中度植被覆蓋區(qū);陜北地區(qū)南部和關(guān)中地區(qū)同樣向更高一級的植被覆蓋區(qū)轉(zhuǎn)化;陜南地區(qū)保持穩(wěn)定。

上述植被覆蓋顯著增加的地區(qū)主要分布在陜北黃土高原和陜南秦巴山區(qū)，與退耕還林工程實施區(qū)域一致;陜北地區(qū)南部、關(guān)中地區(qū)和秦巴山中低山區(qū)植被覆蓋穩(wěn)中有增，與天然林保護(hù)工程分布區(qū)一致。充分說明退耕還林工程和天然林保護(hù)工程在陜西省取得了較好的成效。

傳統(tǒng)的空間模擬及預(yù)測依據(jù)大量數(shù)據(jù)進(jìn)行回歸分析，以回歸方程估算屬性特征，其分析過程及模擬結(jié)果存在較多不確定性。該研究基于IDRISI平臺的CA-Markov模型，其工作原理是以預(yù)測基期的植被覆蓋空間分布為初始狀態(tài)，以基期和之后各級植被覆蓋區(qū)轉(zhuǎn)移面積及適宜性圖集表述的像元適宜植被覆蓋度為依據(jù)，對植被覆蓋空間分布進(jìn)行重新匹配，從而得到預(yù)測的各級植被覆蓋區(qū)面積，有效地模擬和預(yù)測了植被的時空分布特征。

目前CA-Markov模型的運用主要在于土地利用格局的模擬預(yù)測，對于植被覆蓋模擬預(yù)測的研究案例較少。植被作為覆蓋的一種類型，其動態(tài)變化實質(zhì)上也是土地覆蓋的動態(tài)變化，并且相對于不同土地利用類型之間的相互轉(zhuǎn)化，單一的植被變化模擬和預(yù)測更為簡單。由于此次研究首次嘗試?yán)肅A-Markov模型進(jìn)行植被覆蓋分布的模擬預(yù)測，著重模擬了自然演變條件下的時空規(guī)律，而植被覆蓋變化是一個復(fù)雜的過程，受到自然、人類活動以及土地利用發(fā)展政策等多種不確定因素的影響。因此，如何在全面分析和綜合考慮多種因素作用的基礎(chǔ)上調(diào)整模型參數(shù)也有待進(jìn)一步深入探究。

參考文獻(xiàn)

[1]

馬明國，董立新，王雪梅，過去21a中國西北植被覆蓋動態(tài)監(jiān)測與模擬[J].冰川凍土，2003，25（2）：232-236.

[2] 顧祝軍，曾志遠(yuǎn).遙感植被蓋度研究[J].水土保持研究，2005（2）：18-21.

[3] SCHUCKNECHT A，MATSCHULLAT J，ERASMI S.Spatial and temporal variability of vegetation status in Paraíba， Northeastern Brazil[C]//IEEE international geoscience and remote sensing symposium.Munich，Germany：IEEE，2012.

[4] TESTA S，SOUDANI K，BOSCHETTI L，et al.MODIS-derived EVI， NDVI and WDRVI time series to estimate phenological metrics in French deciduous forests[J]. International journal of applied earth observations and geoinformation，2018，64：132-144.

[5] 郭廣猛，謝高地，甄霖.涇河上游固原地區(qū)的NDVI變化與降水的相關(guān)性研究[J].資源科學(xué)，2007，29（2）：178-182.

[6] 馬娜，胡云鋒，莊大方，等.基于遙感和像元二分模型的內(nèi)蒙古正藍(lán)旗植被覆蓋度格局和動態(tài)變化[J]. 地理科學(xué)，2012，32（2）：251-256.

[7] 陳京華，賈文雄，趙珍，等.1982—2006 年祁連山植被覆蓋的時空變化特征研究[J].地球科學(xué)進(jìn)展，2015，30（7）：834-845.

[8] 王光鎮(zhèn)，王靜璞，鄒學(xué)勇，等.基于像元三分模型的錫林郭勒草原光合植被和非光合植被覆蓋度估算[J]. 生態(tài)學(xué)報，2017，37（17）：5722-5731.

[9] MEMARIAN H，BALASUNDRAM S K，TALIB J B，et al.Validation of CA-markov for simulation of land use and cover change in the Langat basin， Malaysia[J].Journal of geographic information system，2012，4：542-554.

[10] DEEP S，SAKLANI A.Urban sprawl modeling using cellular automata[J]. The Egyptian journal of remote sensing and space sciences，2014，17（2）：179-187.

[11] MONDAL M S，SHARMA N，GARG P K，et al.Statistical independence test and validation of CA Markov land use land cover （LULC） prediction results[J]. The Egyptian journal of remote sensing and space sciences，2016，19：259-272.

[12] 楊國清，劉耀林，吳志峰.基于CA-Markov模型的土地利用格局變化研究[J].武漢大學(xué)學(xué)報（信息科學(xué)版），2007，32（5）：414-418.

[13] 鄭青華，羅格平，朱磊，等.基于CA-Markov模型的伊犁河三角洲景觀格局預(yù)測[J]. 應(yīng)用生態(tài)學(xué)報，2010，21（4）：873-882.

[14] 趙冠偉，陳穎彪，陳健飛，等.CA-Markov模型的空間尺度敏感性研究[J].地理科學(xué)，2011，31（8）：897-902.

[15] 趙建軍，張洪巖，喬志和，等.基于CA-Markov模型的向海濕地土地覆被變化動態(tài)模擬研究[J]. 自然資源學(xué)報，2009，24（12）：2178-2186.

[16] 秦超. 陜西省植被覆蓋時空演變及其驅(qū)動因子研究[D].西安：陜西師范大學(xué)，2015.

[17] 陳云浩，李曉兵，陳晉，等. 1983—1992年中國陸地植被NDVI演變特征的變化矢量分析[J].遙感學(xué)報，2002，6（1）：12-18.

[18] 韓文權(quán)，常禹.景觀動態(tài)的Markov模型研究：以長白山自然保護(hù)區(qū)為例[J].生態(tài)學(xué)報，2004，24（9）：1958-1965.