基于土壤指數(shù)的中國(guó)荒漠化時(shí)空演變分析

摘要：研究以?xún)蓚€(gè)土壤指數(shù)：亮荒漠指數(shù)和暗荒漠指數(shù)，確定了2001～2017年的中國(guó)荒漠化范圍，并結(jié)合重點(diǎn)防治區(qū)域分析時(shí)空演變特征，探討了生態(tài)工程的投入與荒漠化的聯(lián)系。揭示了該提取方法可以用于大尺度的荒漠化年度監(jiān)測(cè)，總體精度達(dá)到94%。結(jié)果表明：近17年中國(guó)荒漠化減少超過(guò)12×10⁴km2，相當(dāng)于一個(gè)福建省的面積，集中分布在半干旱荒漠治理區(qū);生態(tài)工程的投入與荒漠化減少量呈現(xiàn)顯著的正相關(guān)，表明政策投入對(duì)土地退化起到明顯的逆轉(zhuǎn)作用。

關(guān)鍵詞：荒漠化;土壤指數(shù);時(shí)空演變;生態(tài)工程;年度監(jiān)測(cè)

中圖分類(lèi)號(hào)：P208 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1674-9944（ 2020） 2-0025-03

1 引言

荒漠化是指包含氣候變異和人類(lèi)活動(dòng)在內(nèi)的種種因素造成的干旱、半干旱和亞濕潤(rùn)干旱地區(qū)的土地退化。為了掌握荒漠化和沙化狀況，我國(guó)每五年組織開(kāi)展一次荒漠化和沙化土地監(jiān)測(cè)工作，以原國(guó)家林業(yè)局第五次全國(guó)荒漠化和沙化監(jiān)測(cè)結(jié)果分析，截止2014年底，全國(guó)荒漠化土地面積為261×10⁴km2[1]。最近有學(xué)者開(kāi)發(fā)了兩個(gè)土壤指數(shù)：亮荒漠指數(shù)RIBS（ the Product In-dex for Bright Soil）和暗荒漠指數(shù)PIDS（ the ProductIndex for Dark Soil）[2]，在對(duì)于土壤的光譜分離性和分類(lèi)精度方面具有可靠的性能，但是所提出的這兩個(gè)土壤指數(shù)目前未應(yīng)用到荒漠化的監(jiān)測(cè)中，其有效性需要在大尺度的荒漠化提取中進(jìn)一步評(píng)估。因此本文以中國(guó)區(qū)域作為研究對(duì)象來(lái)評(píng)價(jià)土壤指標(biāo)，基于亮荒漠指數(shù)和暗荒漠指數(shù)判斷規(guī)則的分類(lèi)方法，得到2001～2017年的中國(guó)范圍的荒漠化增加和減少區(qū)域，并在荒漠化分區(qū)的基礎(chǔ)上，進(jìn)一步探索荒漠化的時(shí)空演變和重心點(diǎn)的遷移特征，為生態(tài)環(huán)境建設(shè)和荒漠化防治策略提供成果參考。

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

2.1 數(shù)據(jù)獲取及預(yù)處理

研究選取了從2001年1月到2017年12月份、8d最大化合成、500 m空間分辨率的MOD09A1數(shù)據(jù)，來(lái)自美國(guó)宇航局（ http：//ladsweb. nascom. nasa. gov/da-ta/search. html）。在MATLAB的操作環(huán)境下，首先對(duì)獲取得到遙感數(shù)據(jù)幾何糾正、鑲嵌和投影變換，在此基礎(chǔ)上，采用基于最新測(cè)繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下的標(biāo)準(zhǔn)地圖對(duì)投影后的影像進(jìn)行裁剪，得到用于荒漠化監(jiān)測(cè)的預(yù)處理產(chǎn)品。為了更好地分析大范圍的荒漠化空間格局特征，將研究區(qū)分為4個(gè)治理區(qū)[3]：（I）;干旱荒漠地帶沙丘活化及流沙入侵沙漠化治理發(fā)展區(qū)（Ⅱ）;半濕潤(rùn)地帶沙漠化治理零星分布區(qū)（Ⅳ）;西部高寒地帶治理區(qū)（Ⅲ）。在這四大區(qū)下，包含35個(gè)更小的重點(diǎn)防治亞區(qū)[4]。通過(guò)無(wú)人機(jī)實(shí)地低空航拍影像、GPS定位野外景觀照片、不同時(shí)期Google Earth影像及相關(guān)專(zhuān)題地圖等手段共獲得4733個(gè)點(diǎn)位。本文將荒漠概括為亮荒漠和暗荒漠，亮荒漠主要是指地面完全被沙所覆蓋、植被非常稀少、雨水稀少且空氣干燥的荒蕪地區(qū)。暗荒漠指地面幾乎被粗砂、礫石所覆蓋，植物稀少的荒漠地帶，多數(shù)地區(qū)不是沙漠而是裸巖，這兩種均為土地退化點(diǎn)位。

2.2 荒漠化監(jiān)測(cè)方法

基于亮度一暗度一綠度模型框架建立起亮荒漠指數(shù)（ RIBS）和暗荒漠指數(shù)（PIDS），其中歸一化積雪指數(shù)（ NDSI）、經(jīng)過(guò)纓帽變化的亮度分量（TCI）和歸一化植被指數(shù)（ NDVI）三個(gè)指標(biāo)分別對(duì)應(yīng)亮度、暗度和綠度。首先計(jì)算上述三個(gè)分量的空間歸一化結(jié)果，然后逐像元計(jì)算RIBS和PIDS，基本公式如下：

式（2）中：RIBS為像元亮荒漠指數(shù)值，N—NDSI為采用公式（1）空間歸一化后的積雪指數(shù)，同樣，N一TCI為空間歸一化后的纓帽變化亮度分量。

RIDS= TCI×NDVI（3）

式（3）中：PIDS為像元暗荒漠指數(shù)值，TC1為纓帽變化中的亮度分量;NDVI為歸一化植被指數(shù)。

在所有地物類(lèi)型中，亮荒漠的RIBS值域范圍最小，且內(nèi)部異質(zhì)性較小;暗荒漠的PIDS同樣是最小值。由于在PIDS中，暗荒漠、亮荒漠與水體的分離度較低，因而首先采用RIBS提取亮荒漠，并同時(shí)提取水體作為掩膜，然后再使用PIDS獲得暗荒漠，最后綜合亮荒漠和暗荒漠空間分布得到荒漠化分布產(chǎn)品。該研究監(jiān)測(cè)方法流程如圖1所示。

其中，θ1、θ2和θ3是閾值常數(shù)。

依據(jù)上述流程圖，首先獲得年際上的荒漠指數(shù)空間分布圖，并使用地物類(lèi)型穩(wěn)定點(diǎn)位訓(xùn)練得到3個(gè)閾值常數(shù)項(xiàng)，其中2001年的θ1、θ2和θ3依次為1.30、0.036、3.1，2017年的θ1、θ2和03依次為1. 37、0.041、2.9。將上述閾值代人荒漠化的判斷規(guī)則中（圖1），以2017年為例，即亮荒漠指數(shù)小于閾值θ1-1. 37得出了亮荒漠，然后暗荒漠指數(shù)小于閾值θ2-0. 041且亮荒漠指數(shù)小于閾值θ3-2.9得到暗荒漠，最終獲得了2017年荒漠化分布圖（圖2a），同理獲得2001年荒漠分布圖。基于上述兩幅年度荒漠化范圍，獲得2001～2017年荒漠化增加和減少區(qū)域（圖2b）。

3 結(jié)果與分析

3.1 精度驗(yàn)證

為了直觀地看出荒漠化空間格局的準(zhǔn)確性，本研究選取了1292個(gè)和1733個(gè)植被覆蓋類(lèi)型穩(wěn)定的點(diǎn)位對(duì)兩年的結(jié)果進(jìn)行精度驗(yàn)證。2001年和2017年總體精度均達(dá)到了94%，Kappa系數(shù)均在85以上，表明該算法適用于大范圍的荒漠化信息提取。發(fā)現(xiàn)在荒漠化混淆的點(diǎn)位中，由于光譜特征相似造成識(shí)別成草地或者草地誤分成荒漠化。

從另一方面，依據(jù)第二次和第五次全國(guó)荒漠化和沙化狀況，即從1999年底到2014年底，全國(guó)荒漠化面積由267. 54×10⁴km2減少至261.16×10⁴km2，本研究基于2001～2017年的荒漠化分類(lèi)結(jié)果，全國(guó)荒漠化面積由272. 04×10⁴km2減少至259. 81×10⁴km2，凈減少12.23×10⁴km2.相較于荒漠化公報(bào)中的凈減少面積6. 38×10⁴km2，并兼顧到時(shí)間范圍的不一致，該算法還是取得了預(yù)期良好的精度效果。

3.2 荒漠時(shí)空演變特征

總體來(lái)看，荒漠化的特征是“整體減少，局部增加”的態(tài)勢(shì)，即4個(gè)分區(qū)荒漠化均表現(xiàn)出凈減少的趨勢(shì)，面積達(dá)到11.8×10⁴km2，占全國(guó)凈減少面積的96. 8%，說(shuō)明四大區(qū)域基本上覆蓋了絕大部分荒漠化區(qū)域，其中半干旱荒漠治理區(qū)（I）的面積凈減少比重最大，凈減少比例近56%。在更小的35個(gè)分區(qū)中，有29個(gè)分區(qū)是面積減少大于面積增加，即呈現(xiàn)出凈減少的趨勢(shì)，總量達(dá)‘到12.7×10⁴km2，而在其余的6個(gè)分區(qū)中，荒漠化面積凈增加量超過(guò)0.9×10d km2。

空間格局上，在呈現(xiàn)凈減少趨勢(shì)的29個(gè)亞區(qū)中，鄂爾多斯地區(qū)（I -2）最為顯著，荒漠化面積凈減少面積（單位：×10⁴km2）達(dá)到1.6×l0⁴km2，其次是錫林郭勒地區(qū)（I -7）是1.5×104 km2、呼倫貝爾草原（I-9）是1.4 × 10⁴km2、河北壩上地區(qū)（I-5）是1.3×10⁴km2、塔里木盆地（Ⅱ-4）是0.8×10⁴km2和科爾沁草原（I-8）是0.7×10⁴km2。干旱治理區(qū)（Ⅱ）和西部高寒治理區(qū)（Ⅲ）的部分地區(qū)存在荒漠化凈增加的趨勢(shì)，尤其以吐一哈盆地（Ⅱ-3）荒漠化面積占絕對(duì)優(yōu)勢(shì)超過(guò)0.8×10⁴km2，其次是柴達(dá)木盆地（Ⅲ-2）是0.07×10⁴km2、昆侖山脈（Ⅲ-1）是0.06×10⁴km2等，除此之外的半濕潤(rùn)零星分布區(qū)中的部分區(qū)域的荒漠化面積也呈現(xiàn)微弱增加的態(tài)勢(shì)（圖3b）。

近17年來(lái)荒漠化重心大致位于在甘肅和青海的接壤處附近，重心向西南方向偏移，向西偏移經(jīng)度1.5。，向南偏移緯度0.4。，遷移量為201 km。重心偏移是由于內(nèi)蒙古荒漠化減少顯著，因而重心向新疆和西藏地區(qū)偏移，荒漠化整體邊界向西南方向移動(dòng)。但是各個(gè)區(qū)內(nèi)荒漠化重心點(diǎn)坐標(biāo)遷移方向不一致，半干旱治理區(qū)（I）重心點(diǎn)位于錫林郭勒草原附近，逐漸向東北方向遷移，揭示了烏蘭察布草原和鄂爾多斯地區(qū)荒漠化面積減少幅度相較于呼倫貝爾草原等地區(qū)較大，與全區(qū)荒漠化土地重心遷移相反。半干旱治理區(qū)（Ⅱ）重心向西北方向移動(dòng);磋部高寒治理區(qū)和半濕潤(rùn)零星分布區(qū)荒漠化重心點(diǎn)均向西南方向遷移，與全區(qū)荒漠化土地重心遷移一致。因而四個(gè)地區(qū)的荒漠化重心都由半干旱中西部地區(qū)向四周遷移，由于該地區(qū)是鄂爾多斯附近，荒漠化減少量最為顯著，并且干旱地區(qū)的微弱增加和該地區(qū)以東的減少量相對(duì)較小導(dǎo)致重心由中部地區(qū)向邊緣擴(kuò)散。從荒漠化重心空間位置遷移量來(lái)看，半干旱治理區(qū)遷移量最大，為40.2 km;半濕潤(rùn)零星分布區(qū)和西部高寒地帶荒漠化遷移量較大，分別為13.3 km和15.2 km;干旱治理區(qū)遷移量最小，為9.9 km，表明該區(qū)荒漠化面積變化幅度最小（圖3a）。

3.3 生態(tài)工程與荒漠化的關(guān)系

本研究重點(diǎn)分析了中國(guó)政府在荒漠化斗爭(zhēng)中生態(tài)工程上的努力，涉及到的生態(tài)工程包括三北防護(hù)林體系建設(shè)工程（SDP TN）、中國(guó)一全球環(huán)境基金干旱生態(tài)系統(tǒng)土地退化防治伙伴關(guān)系（PRE-GEF）、津京風(fēng)沙源治理工程（ SCP- BT）和草原生態(tài)保護(hù)補(bǔ)助獎(jiǎng)勵(lì)項(xiàng)目（GE-SA），整理分析得到這四個(gè)生態(tài)工程在基于省域單元上的投資金額和涉及面積，需要說(shuō)明的是在統(tǒng)計(jì)投資金額時(shí)，由于時(shí)間跨度較長(zhǎng)，因而需要考慮居民消費(fèi)價(jià)格指數(shù)來(lái)獲得總投資額。生態(tài)工程主要涉及17個(gè)省市，其中內(nèi)蒙古自治區(qū)的投資數(shù)額達(dá)到567億元，并且投資涉及面積比例最高，達(dá)到77×10⁴kmz，因而內(nèi)蒙古荒漠化改善效果最為顯著。可以發(fā)現(xiàn)（圖3c）內(nèi)蒙古、新疆等地區(qū)，三種曲線(xiàn)趨勢(shì)以及變化都相對(duì)較為一致，最不一致的出現(xiàn)在京津冀地區(qū)，投資金額與投資涉及面積比值較高，而青海、西藏等地區(qū)，投資金額與投資涉及面積比值較低，以上與政策的具體投入力度相關(guān)。但是總體表現(xiàn)出荒漠化減少面積與生態(tài)工程投資面積呈現(xiàn)顯著的正相關(guān)，相關(guān)性系數(shù)達(dá)到0. 84，并且與生態(tài)工程投資金額同樣呈現(xiàn)顯著的正相關(guān)，相關(guān)系數(shù)同樣高達(dá)0. 76，均通過(guò)了顯著性檢驗(yàn)。其中草原生態(tài)保護(hù)補(bǔ)助獎(jiǎng)勵(lì)項(xiàng)目對(duì)我國(guó)的荒漠化減少起到了顯著的積極作用，揭示了近些年生態(tài)工程效益顯著，治理工程的投入對(duì)于荒漠化的減少具有直接的驅(qū)動(dòng)作用。

4 結(jié)論

研究基于MODIS遙感影像數(shù)據(jù)，采用荒漠化指數(shù)的閾值判斷規(guī)則獲取2001～2017年中國(guó)的荒漠化變化范圍，并對(duì)荒漠化時(shí)空演變進(jìn)行研究，簡(jiǎn)要分析了當(dāng)前生態(tài)工程的效益，得出如下結(jié)論。

（1）整體來(lái)看，中國(guó)荒漠化特征是“整體呈現(xiàn)改善狀態(tài)，部分地區(qū)呈現(xiàn)退化狀態(tài)”，荒漠化重心整體向西南方向偏移，重心由內(nèi)蒙古中西部地區(qū)向四周遷移。

（2）從荒漠化動(dòng)態(tài)變化熱點(diǎn)分區(qū)來(lái)看，在4大治理區(qū)域中，荒漠化均表現(xiàn)出減少的趨勢(shì)，尤其在半干旱地區(qū)最為突出，在其劃分的亞區(qū)中，鄂爾多斯地區(qū)的荒漠化凈減少程度最高，達(dá)到1.5×10⁴km2。但是依然存在少部分亞區(qū)荒漠化出現(xiàn)增加的趨勢(shì)，如吐一哈盆地凈增加面積達(dá)到0.8×10⁴km2以上。

17年間，中國(guó)荒漠化得到了有效的改善。近些年大力投入的生態(tài)工程與荒漠化面積的減少是緊密相關(guān)的。綜合各地區(qū)荒漠化演變特征，對(duì)于土地改善的地區(qū)，如鄂爾多斯草原和錫林郭勒草原等，繼續(xù)加強(qiáng)植被的保護(hù)和生態(tài)工程投入[5]。對(duì)于惡化的地區(qū)，如吐一哈盆地，需要加強(qiáng)生態(tài)恢復(fù)與重建，對(duì)于西部高寒地帶生態(tài)敏感脆弱的地區(qū)，避免過(guò)度農(nóng)墾、過(guò)度放牧，從而確保荒漠化防治的順利進(jìn)行[6]。本文對(duì)于荒漠化信息提取方法的研究，主要基于前面學(xué)者的研究成果繼而推廣到大范圍的研究區(qū)域上，從而得到各熱點(diǎn)研究區(qū)的定量荒漠變化面積，因此該荒漠化變化結(jié)果分析與驗(yàn)證有待于進(jìn)一步研究。

參考文獻(xiàn)：

[1]屠志方，李夢(mèng)先，孫濤.第五次全國(guó)荒漠化和沙化監(jiān)測(cè)結(jié)果及分析[J].林業(yè)資源管理，2016，（1）：1-5.

[2] Qiu B，Zhang K， Tang Z，et al. Developing soil indices based onbrightness，darkness，and greenness to improve land surface map-ping accuracy [Jl. Giscience&Remote Sensing， 2017：1—19.

[3]周日平.中國(guó)荒漠化分區(qū)與時(shí)空演變[J].地球信息科學(xué)學(xué)報(bào)，2019， 21（5）：675-687.

[4]王濤.我國(guó)沙漠化研究的若干問(wèn)題——3.沙漠化研究和防治的重點(diǎn)區(qū)域[J].中國(guó)沙漠，2004，24（1）：1—9.

[5]李春蘭，朝魯門(mén)，包玉海。21世紀(jì)初期氣候波動(dòng)下渾善達(dá)克沙地荒漠化動(dòng)態(tài)變化分析[J].干旱區(qū)地理（漢文版），2015， 38（3）：556-564.

[6]孫建國(guó)，王濤，顏長(zhǎng)珍，氣候變化和人類(lèi)活動(dòng)在榆林市荒漠化過(guò)程中的相對(duì)作用[J].中國(guó)沙漠，2012，32（3）：625-630.

作者簡(jiǎn)介：燕雄飛（1995-），男，碩士研究生，研究方向?yàn)榛哪O(jiān)測(cè)相關(guān)研究。