中蒙俄經(jīng)濟(jì)走廊荒漠化時空格局變化及其驅(qū)動因子

范澤孟,李賽博

1 中國科學(xué)院地理學(xué)與資源研究所,資源與環(huán)境信息系統(tǒng)國家重點實驗室, 北京 100101 2 中國科學(xué)院大學(xué)資源與環(huán)境學(xué)院, 北京 100049 3 江蘇省地理信息資源開發(fā)與利用協(xié)同創(chuàng)新中心, 南京 210023

土地荒漠化作為一個重要的全球性生態(tài)環(huán)境及社會經(jīng)濟(jì)問題,不僅造成生態(tài)環(huán)境的惡化和全球糧食產(chǎn)量的降低,甚至?xí)：Φ饺祟惿婧蜕鐣€(wěn)定[1-3]。2015年《聯(lián)合國防止荒漠化公約》統(tǒng)計顯示,全球約有30億人口居住在土地荒漠化區(qū)；而且到2045年,大約1.35億人將會因荒漠化問題而流離失所。由于荒漠化和干旱,全球農(nóng)用地正以每年12萬km2的速度轉(zhuǎn)為荒地,而通過修復(fù)生態(tài)系統(tǒng)退化區(qū)的土壤,有望每年儲存30億t碳[4]。中蒙俄經(jīng)濟(jì)走廊作為“一帶一路”倡議中的六大經(jīng)濟(jì)走廊之一,而中蒙俄跨境鐵路、公路和原油管道作為支撐中蒙俄經(jīng)濟(jì)走廊區(qū)域合作發(fā)展的關(guān)鍵樞紐和通道,其沿線區(qū)域的荒漠化問題直接關(guān)系到跨境鐵路、公路和原油管道建設(shè)及安全運營重大生態(tài)環(huán)境[5]。因此,對中蒙俄經(jīng)濟(jì)走廊開展土地荒漠化程度識別與驅(qū)動機(jī)理研究,不僅可以作為有效防止和治理荒漠化的基本前提,而且有助于開展生物多樣性、自然保護(hù)區(qū)及其他生態(tài)環(huán)境保護(hù)領(lǐng)域的技術(shù)交流合作。

自20世紀(jì)70年代以來,國內(nèi)外學(xué)者[6-8]利用遙感數(shù)據(jù)覆蓋范圍廣,獲取速度快等優(yōu)勢,將不同平臺或傳感器的影像數(shù)據(jù)反演的植被指數(shù)[9-11]與土壤[12]、氣候[13]、地形[14]等多種數(shù)據(jù)源結(jié)合,對土地荒漠化程度識別。譬如,Lamchin等基于NDVI、表土粒徑指數(shù)(TGSI)和地表反照率(Albedo),并結(jié)合決策樹方法識別了蒙古國的霍格諾汗自然保護(hù)區(qū)荒漠化[15]。在此基礎(chǔ)上,大量學(xué)者嘗試將遙感影像反演參數(shù)作為輸入數(shù)據(jù),運用分類回歸樹(CART)[16]、支持向量機(jī)(SVM)[17]和隨機(jī)森林(RF)[18]等機(jī)器學(xué)習(xí)方法對荒漠化狀況進(jìn)行識別。其中,CART方法是Breiman在1984年提出的基于決策樹構(gòu)建的既可用于分類,也可用于連續(xù)變量預(yù)測的算法[19]；SVM模型是一種基于統(tǒng)計學(xué)習(xí)理論的重要機(jī)器學(xué)習(xí)算法,即將非線性可分?jǐn)?shù)據(jù)映射到高維空間,完成線性可分,已被廣泛應(yīng)用于遙感圖像的分類[20]；RF方法是Breiman在2001年提出的一種基于分類回歸樹的分類器,能夠處理較為復(fù)雜的遙感圖像,適用于多類別與多特征數(shù)據(jù)[21]。另外,還有學(xué)者在對地表反照率(Albedo)和歸一化植被指數(shù)(NDVI)進(jìn)行集成并構(gòu)建Albedo-NDVI特征空間方法,有效實現(xiàn)了荒漠化程度的等級劃分[22]。

雖然利用遙感數(shù)據(jù)和荒漠化程度評估模型可以識別荒漠化狀況的動態(tài),但區(qū)分并測度荒漠化進(jìn)程中氣候變化和人類活動相對作用強(qiáng)度仍然是一個具有挑戰(zhàn)性的工作[23]?，F(xiàn)有的基于實際與潛在的植被狀況[24-25]、殘差趨勢法[26]、降水利用率指數(shù)[27]雖然能揭示土地荒漠化過程中人類活動的潛在驅(qū)動作用,但是無法定量評估氣候變化和人類活動在荒漠化過程中的相對作用強(qiáng)度[28]?；趦舫跫壣a(chǎn)力的人類占用 (Human Appropriation of NPP, HANPP) 部分[29],即氣候和遙感數(shù)據(jù)模擬的潛在凈初級生產(chǎn)力(Potential NPP,PNPP)和實際存留在生態(tài)系統(tǒng)的實際NPP(Actual NPP,ANPP)之間的差異,可分別用PNPP和HANPP表征氣候變化和人類活動的影響[30],并利用線性趨勢分析方法定量評估荒漠化過程中氣候變化和人類活動所占比例。

論文旨在針對中蒙俄經(jīng)濟(jì)走廊區(qū)域的荒漠化格局及程度的定量分析問題,利用Landsat、Sentinel-2等多源遙感影像數(shù)據(jù)構(gòu)建荒漠化程度判別標(biāo)準(zhǔn)和指標(biāo)體系,基于Google Earth Engine(GEE)云平臺,在CART、SVM、RF和Albedo-NDVI四個模型對荒漠化動態(tài)識別精度對比的基礎(chǔ)上,篩選最優(yōu)模型方法,實現(xiàn)中蒙俄經(jīng)濟(jì)走廊荒漠化格局及程度的定量識別。并在此基礎(chǔ)上,引入線性趨勢法計算 PNPP 和 HANPP 在2000—2015年的斜率,實現(xiàn)該時段內(nèi)PNPP和HANPP二者斜率的正負(fù)組合的情景分析,進(jìn)而實現(xiàn)氣候變化和人文因素對中蒙俄經(jīng)濟(jì)走廊荒漠化驅(qū)動系數(shù)的定量測度。

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

1.1 數(shù)據(jù)收集與處理

1.2 荒漠化程度的判別指標(biāo)

鑒于荒漠化程度判別指標(biāo)應(yīng)具備可反演性、可靠性和實用性的原則[35],選取植被、土壤、地形和氣候等參數(shù)作為荒漠化程度的判別指標(biāo)。根據(jù)提取的2000年和2015年的地表反射率影像數(shù)據(jù),采用各個指標(biāo)求算模型(表2),分別獲得：表征地表植被狀態(tài)生物物理特征的歸一化植被指數(shù)(NDVI)、植被覆蓋度(FVC)和修正的土壤調(diào)節(jié)植被指數(shù)(MSAVI)；描述土壤信息的表土粒徑指數(shù)(TGSI)、地表反照率(Albedo)、歸一化水體指數(shù)(NDWI)和裸土指數(shù)(BSI)。另外,地形指標(biāo)和氣候指標(biāo)包括：高程、坡度、坡向和陸表溫度(LST),坡度和坡向是基于SRTM 30m數(shù)據(jù),經(jīng)GEE平臺的SRTM數(shù)字高程模型函數(shù)計算得出。

表2 荒漠化程度判別的植被及土壤參數(shù)

表2中,Red為傳感器紅色波段,NIR為近紅外波段,MIR為中紅外波段,Blue為藍(lán)色波段,Green為綠色波段,SWIR1為短波紅外1波段,SWIR2為短波紅外2波段；FVC中NDVIS和NDVIv分別表示裸土和植被的NDVI值,為排除研究區(qū)NDVI的噪音現(xiàn)象,選取區(qū)域內(nèi)NDVI頻率分布直方圖的5%和95%的NDVI值分別作為NDVIS和NDVIv值。

在確定以上荒漠化程度提取指標(biāo)及參數(shù)計算方法的基礎(chǔ)上,基于Google Earth Engine云計算平臺上,分別實現(xiàn)分類回歸樹(CART)、支持向量機(jī)(SVM)、隨機(jī)森林(RF)和Albedo-NDVI(表3)的模型算法,提取中蒙俄經(jīng)濟(jì)走廊荒漠化程度的空間分布。

表3 荒漠化程度的分類模型

Albedo-NDVI理論公式：

Albedo=k×NDVI+b

α×k=-1

DDI=α×NDVI-Albedo

式中,DDI為荒漠化差值指數(shù)。

1.3 荒漠化信息提取精度驗證與評價

在分別運用CART、SVM、RF和Albedo-NDVI四個模型獲取中蒙俄經(jīng)濟(jì)走廊荒漠化程度空間分布的基礎(chǔ)上,選取120個樣點對荒漠化程度的分類和精度進(jìn)行驗證。其中包括研究區(qū)內(nèi)24個全球通量觀測塔站點和GEE平臺隨機(jī)生成的96個樣點。在四種分類模型中,基于GEE平臺的隨機(jī)算子抽取30%的樣本點用于建立混淆矩陣,得出總體精度(Overall Accuracy, OA)、生產(chǎn)者精度(Producer′s Accuracy, PA)、用戶精度(User′s Accuracy, UA)及Kappa系數(shù)評價指標(biāo),從不同的側(cè)面對比荒漠化識別的精度,其中生產(chǎn)者精度是一個樣本點被正確分類的概率,它是對丟失誤差的測量。而用戶精度是每一行中正確分類的樣本點的個數(shù)與該行中所有樣本點數(shù)的比值,它表示了包含誤差[47-42]。其理論公式可表達(dá)如下:

PAi=Xii/Xi+

UAi=Xii/X+i

式中,N代表樣本點數(shù)的總和,Xii代表混淆矩陣中第i種荒漠化類型反演結(jié)果與實際結(jié)果相同的樣本點數(shù),Xi+代表混淆矩陣中第i種荒漠化類型的列總和,X+i代表混淆矩陣中第i種荒漠化類型的行總和。

通過建立混淆矩陣,得出各個模型的識別精度估量指標(biāo)(表4)。結(jié)果表明,在荒漠化分類結(jié)果中,CART模型是研究區(qū)荒漠化識別精度最高的分類方法,Kappa系數(shù)為0.754,其次是Albedo-NDVI模型,而SVM模型的精度最低。通過對生產(chǎn)者精度對比分析,Albedo-NDVI模型對中度、重度和極重度荒漠化具有較高的生產(chǎn)者精度,說明對這三種荒漠化級別的漏測現(xiàn)象較小,但是對輕度荒漠化的生產(chǎn)精度極低。CART模型對各個荒漠化級別的生產(chǎn)者精度都較高,表明總體漏測現(xiàn)象較小。另外,CART模型對各個荒漠化級別的用戶精度同樣保持較高值,且相比其他分類模型對非荒漠化區(qū)域的用戶精度最高,表明總體的多測誤差較小,而Albedo-NDVI模型對輕度荒漠化的用戶精度極低,說明在識別輕度荒漠化區(qū)域時,具有很高的多測誤差。SVM和RF的用戶精度和生產(chǎn)者精度都較低,表明兩者在研究區(qū)的荒漠化識別效果并不好。

表4 各種分類算法的精度評價

1.4 荒漠化驅(qū)動機(jī)理綜合分析模型

基于MODIS數(shù)據(jù)和利用BIOME-BGC模型與光能利用率模型(CASA)模擬獲得的NPP(ANPP)數(shù)據(jù)產(chǎn)品,提取中蒙俄經(jīng)濟(jì)走廊區(qū)域的500m分辨率的NPP(ANPP)空間數(shù)據(jù)。在GEE云平臺上,實現(xiàn)Thornthwaite memorial 模型[43]的模型算法設(shè)計,進(jìn)而求解PNPP和HANPP。并結(jié)合降水和氣溫數(shù)據(jù)的空間數(shù)據(jù),在GEE平臺上構(gòu)建并實現(xiàn)了中蒙俄經(jīng)濟(jì)走廊荒漠化驅(qū)動機(jī)理綜合分析模型,其理論公式可表達(dá)如下：

交際的雙方要達(dá)到預(yù)期的交際效果，就必須擁有共同的語言背景知識或語用前提，交流時省去一些不言而喻，不言自明的內(nèi)容，這部分被省略的共享的背景知識稱為“情景缺省”。被缺省的部分與語篇之外的文化背景相關(guān)，就是“文化缺省”?？缥幕浑H中，不屬于該文化的接受者，接觸到這些文化缺省時常常遇到“意義真空”，無法將語篇內(nèi)信息與語篇之外的背景知識和經(jīng)驗聯(lián)系起來，難以建立起理解話語所必需的語義連貫和情景連貫。（王東風(fēng)，1997：55）因此對于文化缺省的部分，譯者應(yīng)該采用恰當(dāng)?shù)姆绞絹砜朔幕町?，幫助譯文讀者付出最少的努力，來獲取最大或者足夠的語境效果。

PNPP=3000(1-e-0.0009695(v-20))

L=3000+25t+0.05t3

HANPP=PNPP-ANPP

式中,r表示年降水量(mm),t為年平均氣溫(℃),x表示年份,y為PNPP時得到SlopePNPP(Sp),為HANPP時得到SlopeHANPP(SH),n為時間間隔。

對于研究區(qū)域,首先基于給定時間段內(nèi)荒漠化發(fā)展?fàn)顩r,即荒漠化程度加重的區(qū)域定義為“退化區(qū)”,荒漠化程度減輕的區(qū)域定義為“恢復(fù)區(qū)”,將Sp和SH的正負(fù)組合分為8種情景(表5),以測度氣候變化和人類活動所占比例。在驅(qū)動因子定性方面,當(dāng)Sp為正值、或SH為負(fù)值時,表明在2000—2015年,氣候變化或者人類活動有利于荒漠化區(qū)域的恢復(fù)；而當(dāng)Sp為負(fù)值、或SH為正值時,表明2000—2015年,氣候變化或者人類活動會加重土地荒漠化；在驅(qū)動因子定量方面,如在荒漠化恢復(fù)區(qū)域,情景1表示恢復(fù)完全歸因于氣候變化,情景2表示恢復(fù)完全歸因于人類活動,情景3表示荒漠化恢復(fù)由氣候變化和人類活動共同引起,情景4表示這種狀況超出模型的合理范圍,氣候變化和人類活動都無法解釋荒漠化恢復(fù)現(xiàn)象,此時的“Error”區(qū)域是指被標(biāo)記為錯誤的區(qū)域,定義為氣候變化和人類活動對植被的影響,與研究區(qū)植被的實際生長狀況并不符。即研究區(qū)氣候變化和人類活動的驅(qū)動作用將加重土地荒漠化,但實際上研究區(qū)的荒漠化現(xiàn)象被逐漸改善。

表5 氣候變化和人類活動定量評估模型

2 結(jié)果與分析

2.1 荒漠化時空分布格局及其轉(zhuǎn)換趨勢

根據(jù)荒漠化程度指標(biāo)體系和最佳荒漠化程度評估模型獲得2000年和2015年研究區(qū)荒漠化現(xiàn)狀(圖1)。通過對比兩期荒漠化現(xiàn)狀,發(fā)現(xiàn)研究區(qū)荒漠化空間分布具有很強(qiáng)的地帶性特征,極重度荒漠化區(qū)域主要分布在內(nèi)蒙古自治區(qū)西部,蒙古國的南戈壁盟、東戈壁盟和中戈壁盟,重度和中度荒漠化區(qū)域主要分布在蒙古國的東戈壁盟、及其與內(nèi)蒙古自治區(qū)交界處,輕度荒漠化主要分布在蒙古國的東北部,內(nèi)蒙古自治區(qū)的南部,而非荒漠化區(qū)域主要分布在俄羅斯區(qū)域和中國東北區(qū)域。另外,在2000年到2015年間,中蒙俄經(jīng)濟(jì)走廊區(qū)域的輕度荒漠化面積從71.97萬km2增加到94.8萬km2,以1.52萬km2/a的速度擴(kuò)張,極重度、重度和中度荒漠化面積均呈現(xiàn)減小趨勢,其中極重度荒漠化面積減少最快,以4873km2/a速度遞減(表6)。其中,中國區(qū)域的荒漠化呈恢復(fù)趨勢,尤其是極重度、重度和中度荒漠化區(qū)域呈縮小趨勢；蒙古國區(qū)域的荒漠化狀況則呈加劇趨勢,尤其是極重度、中度和輕度荒漠化區(qū)域呈逐漸擴(kuò)張趨勢；而俄羅斯區(qū)域荒漠化的加重和減輕兩個方向的狀況均最為緩慢。

圖1 研究區(qū)荒漠化時空分布格局Fig.1 The spatiotemporal distribution pattern of desertification in the study area

表6 研究區(qū)荒漠化類型面積統(tǒng)計

為更好地理解研究區(qū)2000—2015年荒漠化空間動態(tài)變化,將荒漠化空間格局變化分為五類(圖2,表7),分別為強(qiáng)烈退化(荒漠化加劇兩級及其以上)、退化(荒漠化加劇到相鄰一級)、穩(wěn)定(荒漠化程度不變)、恢復(fù)(荒漠化降低到相鄰一級)和顯著恢復(fù)(荒漠化降低兩級及其以上)?？臻g統(tǒng)計結(jié)果顯示,在2000—2015年間,中蒙俄經(jīng)濟(jì)走廊的荒漠化程度保持穩(wěn)定的區(qū)域主要分布在中國東部區(qū)域和俄羅斯區(qū)域,退化和嚴(yán)重退化區(qū)主要分布在蒙古國的東部、中國區(qū)域與蒙古國區(qū)域的西部交界處,恢復(fù)和顯著恢復(fù)區(qū)域主要分布在內(nèi)蒙古自治區(qū)和蒙古國的中部地區(qū)。另外,在2000—2015年間,中蒙俄經(jīng)濟(jì)走廊大部分區(qū)域的荒漠化程度未發(fā)生改變(占研究區(qū)面積的80.91%),而荒漠化程度加重區(qū)域和恢復(fù)區(qū)域分別占總面積的9.68%和9.41%。其中,在中國區(qū)域,荒漠化恢復(fù)區(qū)域的面積比荒漠化程度加重區(qū)域的面積多15.18萬km2,而在蒙古國區(qū)域和俄羅斯區(qū)域,恢復(fù)區(qū)域的面積則分別比加重區(qū)域的總面積少13.20萬km2和3.76萬km2。

圖2 研究區(qū)荒漠化時空格局變化 Fig.2 Changes of spatiotemporal pattern of desertification in the study area

表7 研究區(qū)荒漠化時空格局變化面積統(tǒng)計

同時,基于荒漠化轉(zhuǎn)移矩陣(表8)可以發(fā)現(xiàn),2000—2015年研究區(qū)荒漠化恢復(fù)區(qū)域主要來源于極重度荒漠化轉(zhuǎn)為重度荒漠化和輕度荒漠化轉(zhuǎn)為非荒漠化,而荒漠化加重區(qū)域主要來源于重度荒漠化轉(zhuǎn)為極重度荒漠化。

表8 研究區(qū)荒漠化轉(zhuǎn)移矩陣/(×104km2)

2.2 荒漠化驅(qū)動力定量評估

結(jié)合荒漠恢復(fù)區(qū)和退化區(qū)荒漠化程度數(shù)據(jù),運行荒漠化驅(qū)動機(jī)理綜合分析模型,獲得氣候變化和人類活動在中蒙俄經(jīng)濟(jì)走廊荒漠化“恢復(fù)區(qū)”和“退化區(qū)”的影響作用系數(shù)(圖3,圖4,表9)。結(jié)果顯示,在2000—2015年間,氣候變化和人類活動在荒漠化恢復(fù)區(qū)的驅(qū)動作用在空間上具有明顯的空間分異性。其中,氣候變化對蒙古國東部和內(nèi)蒙古自治區(qū)的大部分區(qū)域荒漠化恢復(fù)起主導(dǎo)作用,而且氣候變化驅(qū)動作用>50%的荒漠化恢復(fù)區(qū)域占總恢復(fù)面積的68.8%,而由人類活動驅(qū)動作用>50%的恢復(fù)區(qū)域僅占整個恢復(fù)區(qū)面積的27.36%。氣候變化對荒漠恢復(fù)的驅(qū)動作用,在中國區(qū)域、蒙古國區(qū)域和俄羅斯區(qū)域分別占荒漠恢復(fù)區(qū)面積的77%、58%和38%。另外,在2000—2015年間,在中蒙俄經(jīng)濟(jì)走廊的荒漠化程度加重區(qū)域,人類活動的驅(qū)動效應(yīng)則占主導(dǎo)作用,其中,人類活動驅(qū)動作用>50%的荒漠化加重區(qū)占所有荒漠化加重區(qū)域面積的比例高達(dá)69.68%,主要分布在蒙古國的東部和內(nèi)蒙古的中部地區(qū)。

表9 氣候變化和人為因素在退化和恢復(fù)區(qū)占主導(dǎo)驅(qū)動作用的面積/(×10km2)

2000—2015年間,研究區(qū)荒漠化驅(qū)動機(jī)理結(jié)果中含有部分“Error”區(qū)域(圖3,圖4),其中恢復(fù)區(qū)的“Error”區(qū)域定義為氣候變化和人類活動不利于荒漠化土地恢復(fù),但實際上有利于植被生長的區(qū)域,這些“Error”區(qū)域主要分布在俄羅斯境內(nèi),占所有恢復(fù)區(qū)域的3.83%。另外,退化區(qū)的“Error”區(qū)域定義為氣候變化和人類活動有利于荒漠化土地恢復(fù),但實際上抑制了植被生長的區(qū)域,這些“Error”區(qū)域占所有退化區(qū)域的28.87%。事實上,產(chǎn)生這種現(xiàn)象主要是因為遙感反演荒漠化結(jié)果存在一定的不確定性,并且基于NPP計算驅(qū)動機(jī)理,往往也受氣候因子的噪音影響,因此這些超出模型解釋范疇的“Error”區(qū)域難以避免。

圖3 氣候變化和人為因素在恢復(fù)區(qū)的驅(qū)動作用 Fig.3 The driving strength of climate change and human factors in desertification restoration area

圖4 氣候變化和人為因素在退化區(qū)的驅(qū)動作用 Fig.4 The driving role of climate change and human factors in desertification degradation area

3 討論

自21世紀(jì)以來,中蒙俄經(jīng)濟(jì)走廊區(qū)域的極重度、重度和中度荒漠化面積減少,輕度荒漠化面積增加；中國區(qū)域荒漠化土地恢復(fù)效果顯著,蒙古國區(qū)域荒漠化狀況在加劇等結(jié)果與現(xiàn)有研究成果[15, 18, 44-46]一致。譬如,中蒙俄經(jīng)濟(jì)走廊區(qū)域的49.12%為草地,裸地僅占14.59%。在氣候變化和人類活動干擾的共同驅(qū)動下,草地退化現(xiàn)象嚴(yán)重,其中新增裸地面積的89.94%均自于草地退化,但是裸地面積凈減少近800km2,表明極重度荒漠化狀況呈恢復(fù)趨勢。

為了進(jìn)一步解釋線性趨勢分析方法對氣候變化和人類活動的荒漠化驅(qū)動作用測度的不足,分別選取降水量、平均氣溫和年蒸散發(fā)量3個植被生長的關(guān)鍵氣候參數(shù)[38,40,47]和能夠綜合反映人類活動狀態(tài)的人均GDP[48],對模型中未能解釋的“Error”區(qū)域進(jìn)行分析?？臻g統(tǒng)計分析顯示,在2000—2015年間,中蒙俄經(jīng)濟(jì)走廊區(qū)域的年降水量減少29.7mm,年平均氣溫增加1.57℃,年蒸散發(fā)量減少17.3mm,人均GDP增加6369＄。年降水量減少區(qū)域主要分布 “非荒漠化”區(qū)(圖1),尤其是俄羅斯區(qū)域的中部地區(qū)和中國區(qū)域的長白山脈部分地區(qū)。在中國區(qū)域和蒙古國區(qū)域的年降水量增加區(qū)域主要分布在荒漠化區(qū),其中中國區(qū)域的荒漠化“顯著恢復(fù)”和“恢復(fù)”區(qū)與年降水量增加量大于50mm區(qū)的空間分布一致,尤其是中國區(qū)域的東北平原、大興安嶺南側(cè)地區(qū)和鄂爾多斯高原(圖1,圖5)。

圖5 研究區(qū)2000—2015氣候和人類活動因子變化量Fig.5 The change of climate and human activity factors in the study area from 2000 to 2015

蒙古國區(qū)域的“極重度荒漠化”和“重度荒漠化”區(qū)的年降水量雖呈增加(0—50mm)趨勢,但其升溫幅度(>2.8℃)和年蒸散發(fā)量增加幅度(>25mm),尤其是荒漠化“嚴(yán)重退化”和“退化”區(qū)的升溫幅度在3.0℃以上,再加上該地區(qū)的人均GDP增加幅度(>8500＄)也大于其他區(qū)域。因此,由于該區(qū)域升溫和蒸散發(fā)增加在空間上的疊加,使得荒漠化程度(尤其在“Error”區(qū)內(nèi))呈加重趨勢。同時,位于“Error”區(qū)的蒙古國南戈壁盟,由于對該區(qū)域豐富礦產(chǎn)資源(如世界最大的銅礦奧尤陶勒蓋金銅礦)大量開采[49],在致使人均GDP快速增長的同時,也加重了土地沙化的程度。另外,雖然中國區(qū)域的人均GDP也呈快速增長趨勢,但由于京津風(fēng)沙源治理工程、退牧還草工程等系列國家重大生態(tài)工程的實施,使得中國區(qū)域的荒漠化呈恢復(fù)趨勢。

4 結(jié)論

基于GEE的云計算平臺和多源數(shù)據(jù),結(jié)合植被、土壤、氣候和地形等荒漠化程度的判別指標(biāo)模型,實現(xiàn)CART、SVM、RF和Albedo-NDVI 4個模型對荒漠化動態(tài)識別精度對比分析的基礎(chǔ)上,構(gòu)建荒漠化驅(qū)動機(jī)理綜合分析模型,實現(xiàn)中蒙俄經(jīng)濟(jì)走廊荒漠化程度的時空變化及驅(qū)動效應(yīng)分析。

精度驗證表明,CART模型對中蒙俄經(jīng)濟(jì)走廊荒漠化程度信息提取和分類的精度高于SVM、RF和Albedo-NDVI模型,其Kappa系數(shù)為0.754,且對各個荒漠化程度分類的用戶精度和生產(chǎn)者精度均比較穩(wěn)定。雖然Albedo-NDVI模型的總體精度也比較高,但是對于輕度荒漠化識別效果較差。

研究結(jié)果顯示,在2000—2015年間,中蒙俄經(jīng)濟(jì)走廊區(qū)域的荒漠化程度整體呈減弱趨勢,尤其是極重度荒漠化、重度和中度荒漠化的區(qū)域呈縮減趨勢,但輕度荒漠化區(qū)域卻呈增加趨勢。中國區(qū)域的荒漠化呈恢復(fù)趨勢,凈恢復(fù)面積達(dá)15.18萬km2,而蒙古國區(qū)域和俄羅斯區(qū)域的荒漠化呈加重趨勢,其新增的荒漠化面積分別為13.2萬km2和3.76萬km2。

荒漠化程度時空變化及驅(qū)動作用分析表明,中蒙俄經(jīng)濟(jì)走廊荒漠化程度減弱主要表現(xiàn)為極重度荒漠化轉(zhuǎn)為重度荒漠化、輕度荒漠化轉(zhuǎn)為非荒漠化兩種狀態(tài),而荒漠化程度加重區(qū)域則主要為重度荒漠化轉(zhuǎn)為極重度荒漠化。氣候變化在中蒙俄經(jīng)濟(jì)走廊荒漠化程度減弱和恢復(fù)種起到主導(dǎo)作用,而局部的荒漠化加重區(qū)域,人類活動的驅(qū)動作用超過了氣候變化驅(qū)動作用。