基于相對降水利用率的內(nèi)蒙古風(fēng)蝕區(qū)沙漠化研究

泉濤 王靜璞 黨文碩 劉新秋 沈瑞華

摘要 以MODIS NDVI數(shù)據(jù)和降水量站點數(shù)據(jù)為基礎(chǔ)，計算內(nèi)蒙古風(fēng)蝕區(qū)2001—2018年的相對降水利用率（RUEr），通過分析RUEr的時間變化、空間分布及變異系數(shù)，綜合評價內(nèi)蒙古風(fēng)蝕區(qū)18年間的沙漠化演變情況。結(jié)果表明，18年間內(nèi)蒙古風(fēng)蝕區(qū)內(nèi)的RUEr呈現(xiàn)出微弱的下降趨勢，沙漠化程度加劇;空間分布上，62.98%的地區(qū)為正斜率，植被呈恢復(fù)狀態(tài)，降水為影響沙漠化的主導(dǎo)因素;37.02%的地區(qū)為負斜率，植被出現(xiàn)退化現(xiàn)象，人類活動為影響沙漠化的主導(dǎo)因素。研究區(qū)內(nèi)89.09%的地區(qū)為輕度沙漠化，重度沙漠化地區(qū)集中在內(nèi)蒙古西北部，中度沙漠化和非沙漠化地區(qū)在研究區(qū)內(nèi)零星分布。18年間內(nèi)蒙古風(fēng)蝕區(qū)的沙漠化變異程度較大，87.91%的地區(qū)表現(xiàn)為脆弱到非常脆弱，僅有12.09%的地區(qū)表現(xiàn)為穩(wěn)定到非常穩(wěn)定。

關(guān)鍵詞 沙漠化;相對降水利用率;時空變化;內(nèi)蒙古風(fēng)蝕區(qū)

中圖分類號 X144文獻標識碼 A

文章編號 0517-6611（2021）08-0072-05

doi：10.3969/j.issn.0517-6611.2021.08.019

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

Study on the Desertification in Wind Erosion Area of Inner Mongolia Based on RUEr

QUAN Tao，WANG Jing-pu，DANG Wen-shuo et al （School of Resources and Environmental Engineering，Ludong University，Yantai，Shandong 264025）

Abstract Based on MODIS NDVI data and precipitation station data，this paper calculated the relative rainfall utilization efficiency （RUEr） in the wind erosion area of Inner Mongolia from 2001 to 2018.By analyzing the temporal variation，spatial distribution and coefficient of variation of RUEr，the evolution of desertification in the wind erosion area of Inner Mongolia in the past 18 years was evaluated comprehensively.The results showed that during the past 18 years，the RUEr in the wind erosion area of Inner Mongolia had shown a weak downward trend，which means the degree of desertification had increased.In terms of spatial distribution，62.98% of the study areas had a positive slope，vegetation was in a recovery state，precipitation was the leading factor affecting desertification，37.02% of the study areas had a negative slope，vegetation was in a degraded state，human activity was the leading factor affecting desertification.89.09% of the study areas were light desertification.The areas with heavy desertification were concentrated in the northwest of Inner Mongolia，while the areas with moderate and no desertification are scattered in the study areas.In the past 18 years，the degree of desertification variation in the wind erosion area of Inner Mongolia was relatively large，87.91% of the study areas were fragile to very fragile，and only 12.09% of the study areas were stable to very stable.

Key words Desertification;Relative rainfall utilization efficiency;Spatial-temporal variation;Wind erosion area of Inner Mongolia

近年來，由于全球變暖和人類對環(huán)境的過度干預(yù)，沙漠化成為一個影響全球的環(huán)境問題[1]。內(nèi)蒙古自治區(qū)面積廣大，氣候干旱多風(fēng)沙，易造成水土流失，加之近年來經(jīng)濟的快速發(fā)展帶來過度放牧、資源開采不合理等問題，使內(nèi)蒙古成為我國沙漠化形勢最嚴峻、治理任務(wù)最艱巨的地區(qū)之一[2-3]。

植被破壞是沙漠化最為顯著的特征之一，因此植被信息成為評價沙漠化的一個重要指標。遙感技術(shù)可以快速準確地獲取地面植被信息，加之全天候、多角度、大范圍的優(yōu)勢，其在土地沙漠化監(jiān)測中的應(yīng)用已經(jīng)非常普遍[4]。近年來，許多學(xué)者從不同角度對內(nèi)蒙古沙漠化進行了研究，其中大多分析了內(nèi)蒙古地區(qū)的荒漠化指數(shù)、濕潤指數(shù)、植被指數(shù)、植被覆蓋度和降水利用率等指標的時空變化特征，并與氣候因子、地形因子等進行了相關(guān)性分析[5-8]。

降水利用率（rainfall utilization efficiency，RUE）是凈初級生產(chǎn)力（net primary productivity，NPP）與降水量的比值，用來反映植被生長與耗水特性之間的關(guān)系，同時干旱、半干旱地區(qū)的植被生產(chǎn)力與降水量具有密切的關(guān)系，因此，RUE的變化對干旱、半干旱地區(qū)荒漠草原生態(tài)系統(tǒng)的健康具有顯著的指示作用，能科學(xué)地反映沙漠化的演變特征[8-12]。該研究利用2001—2018年生長季4—10月的MODIS NDVI數(shù)據(jù)和氣象站點逐日降水量數(shù)據(jù)，通過計算各年相對降水利用率來反演內(nèi)蒙古風(fēng)蝕區(qū)的沙漠化演變情況，科學(xué)評價內(nèi)蒙古風(fēng)蝕區(qū)的生態(tài)環(huán)境，以期為內(nèi)蒙古風(fēng)蝕區(qū)沙漠化監(jiān)測預(yù)警及治理工作提供數(shù)據(jù)及理論支持。

1 資料與方法

1.1 研究區(qū)概況

內(nèi)蒙古自治區(qū)（37°24′～53°23′N，97°12′～126°04′E）是我國北部重要的生態(tài)屏障，緯度偏高，經(jīng)度范圍跨度較大（圖1），自東向西表現(xiàn)出極強的經(jīng)度地帶性，植被類型相應(yīng)表現(xiàn)為溫帶森林、溫帶草原和溫帶荒漠。為更加科學(xué)準確地研究內(nèi)蒙古沙漠化的演變特征，該研究將內(nèi)蒙古境內(nèi)的無土壤風(fēng)蝕發(fā)生的溫帶森林區(qū)去除，將剩余易發(fā)生風(fēng)蝕的溫帶草原區(qū)、溫帶荒漠區(qū)作為研究區(qū)，綜合研究該區(qū)域的土地沙漠化演變情況。

內(nèi)蒙古風(fēng)蝕區(qū)約占全國面積的1/10，是氣候變化的敏感脆弱區(qū)，也是易發(fā)生沙漠化的典型地區(qū)[13]。該區(qū)地勢較高，平均海拔在1 000 m左右，有大小河流千余條，大小湖泊近千個。主要氣候類型為溫帶大陸性氣候，冬春漫長，干寒多風(fēng)沙;夏季短促，雨熱同期，日照充足，風(fēng)能資源豐富。在脆弱的自然環(huán)境和以畜牧業(yè)為主的國民經(jīng)濟影響下，該區(qū)土地沙漠化的形勢較為嚴峻。

1.2 數(shù)據(jù)源及預(yù)處理 該研究采用的遙感影像選自NASA/MODIS系列數(shù)據(jù)的MOD13Q1植被指數(shù)數(shù)據(jù)集（來源：https：//modis.gsfc.nasa.gov），下載了2001—2018年生長季4—10月覆蓋內(nèi)蒙古風(fēng)蝕區(qū)的6個瓦片（h25v03、h25v04、h25v05、h26v04、h26v05、h27v04）的影像，共 1 428 景。利用MRT（MODIS reprojection tool）對同期數(shù)據(jù)進行拼接、投影，并提取出各年NDVI數(shù)據(jù)。在ENVI軟件中利用最大值合成法計算各年各像元的NDVI最大值，作為計算RUE的基礎(chǔ)數(shù)據(jù)。采用的氣象數(shù)據(jù)為中國地面氣候資料日值數(shù)據(jù)集（V3.0）中的降水量數(shù)據(jù)（來源：http：//data.cma.cn）。利用Excel軟件篩選出內(nèi)蒙古風(fēng)蝕區(qū)及周邊共80個氣象站點的數(shù)據(jù)，計算出各年生長季的月平均降水量，并在ArcGIS軟件中進行克里金插值，得到2001—2018年生長季內(nèi)蒙古風(fēng)蝕區(qū)空間分辨率250 m（與遙感影像保持一致）的月平均降水量的柵格圖像。

1.3 研究方法

1.3.1 相對降水利用率。姬盼盼等[14]研究表明，在干旱、半干旱地區(qū)NDVI對NPP的作用相對較強;林志東等[15]研究指出在自然狀態(tài)下NDVI與NPP的相關(guān)系數(shù)較高。因此該研究用NDVI替代NPP與降水量（P）求商來計算RUE，計算公式如下：

RUE=NDVIP（1）

由于RUE在數(shù)學(xué)意義上僅反映2個指數(shù)的比值關(guān)系，因此當研究大范圍地區(qū)沙漠化演變情況時，在降水量較少的干旱地區(qū)，會出現(xiàn)2個指數(shù)均較低但RUE值反而高的現(xiàn)象，故對研究區(qū)內(nèi)RUE數(shù)據(jù)進行頻率統(tǒng)計，以[5%，95%]為置信區(qū)間，分別求出各年潛在降水利用率最值RUEmin和RUEmax，進而計算研究區(qū)各年的相對降水利用率（RUEr），計算公式如下：

RUEr=RUE-RUEminRUEmax-RUEmin（2）

RUEr在一定程度上消除了降水帶來的影響，其值反映了該地區(qū)內(nèi)自然植被的相對狀態(tài)。同時，年平均相對降水利用率RUEmer可反映該土地長期的生產(chǎn)能力，年極大相對降水利用率RUEexr可反映該土地受到破壞后的恢復(fù)彈力，且潛在RUE的最值是根據(jù)RUE的概率分布的5%和95%分位數(shù)確定的，RUEr的計算結(jié)果會超出0～1的范圍，因此可根據(jù)RUEr的數(shù)值來評價土地生產(chǎn)力[16]，計算公式如下：

RUEmer=1nni=1xi（3）

RUEexr=xt（4）

式中，n為年份總數(shù)，即n為18;xi為第i年的RUEr;xt為2001—2018年NDVI值最大的年份對應(yīng)的RUEr。

1.3.2 趨勢分析方法。

用一元線性回歸方程的斜率slope來分析RUEr的年際變化趨勢，以此來分析內(nèi)蒙古風(fēng)蝕區(qū)沙漠化的年際變化特征。若slope>0，說明RUEr呈增加趨勢，植被呈恢復(fù)狀態(tài)，沙漠化形勢轉(zhuǎn)好;slope=0，說明RUEr保持不變;slope<0，說明RUEr呈下降趨勢，植被呈退化狀態(tài)，沙漠化形勢加劇。計算公式如下：

slope=n×ni=1（i×xi）-ni=1ini=1xin×ni=1i2-（ni=1i）2（5）

式中，n為年份總數(shù)，即n為18;i為1～18的年序號;xi為第i年的RUEr。

1.3.3 變異系數(shù)。

變異系數(shù)（CV）是反映一組單個觀測數(shù)中各觀測量變化程度的統(tǒng)計量，也是反映時間和空間變異程度的重要指標，數(shù)值越大說明其變異程度越大。為進一步分析研究區(qū)沙漠化發(fā)展的穩(wěn)定程度，該研究計算了研究區(qū)內(nèi)RUEr的變異系數(shù)，并將其劃分為非常穩(wěn)定（CV≤0.15）、穩(wěn)定（0.150.45）共4個等級。計算公式如下：

CV=ni=1（xi-）2n-1（6）

式中，n為年份總數(shù)，即n為18;xi為第i年的RUEr;為RUEmer。

2 結(jié)果與分析

2.1 內(nèi)蒙古風(fēng)蝕區(qū)沙漠化程度時間變化特征

由內(nèi)蒙古風(fēng)蝕區(qū)2001—2018年生長季月平均降水量年際變化趨勢（圖2）可知，18年間內(nèi)蒙古風(fēng)蝕區(qū)內(nèi)降水量呈增加趨勢。研究區(qū)內(nèi)生長季月平均降水量的平均值為379.64 mm，最大值出現(xiàn)在2018年，為509.18 mm，最低值出現(xiàn)在2005年，為294.64 mm。從生長季月平均降水量的距平值變化來看，2001、2004、2005、2006、2007、2009、2010、2011、2015年呈現(xiàn)出負距平，其余年份均呈現(xiàn)出正距平。

由于降水量存在很大的年際變化，僅計算年際比例不足以達到沙漠化監(jiān)測的預(yù)期效果，因此以時間為自變量，RUEr為因變量，對時間序列RUEr進行線性回歸分析，其回歸斜率主要取決于NDVI與降水量的相關(guān)性，斜率為正，表明植被正在恢復(fù)，斜率為負，表明植被正在退化，以此來監(jiān)測長期以來植被對降水的響應(yīng)。

從內(nèi)蒙古風(fēng)蝕區(qū)2001—2018年相對降水利用率的回歸斜率及生長季月平均降水量的空間分布（圖3）可看出，研究區(qū)內(nèi)RUEr的回歸斜率具有明顯的空間差異性，大致表現(xiàn)為由東部向西部逐漸遞增，其值在-0.15～0.26，平均為-0.01，表現(xiàn)為微弱的下降趨勢，說明植被有退化現(xiàn)象，沙漠化形勢在日趨加重，需予以重視。從RUEr時空趨勢來看，全區(qū)內(nèi)大部分地區(qū)RUEr呈增加趨勢，面積為62.02萬km2，占全區(qū)面積的62.98%，集中分布在內(nèi)蒙古中西部地區(qū);該地區(qū)植被類型主要為溫帶草原及溫帶荒漠，大陸性氣候顯著，冬季嚴寒，夏季炎熱，氣溫日較差及年較差較大，日照充足，阿拉善高原的西部地區(qū)日照時間更可達3 400 h以上，多風(fēng)沙天氣，降水稀少，地表徑流不發(fā)育，因此降水是影響該地區(qū)沙漠化的主導(dǎo)因素。此外，還有37.02%的地區(qū)RUEr呈下降趨勢，面積為36.46萬km2，主要分布在東部大興安嶺周圍地區(qū)以及中部的河套平原地區(qū)，前者植被類型主要為溫帶草原，季風(fēng)性氣候

顯著，降水較中西部地區(qū)明顯增多，后者為經(jīng)濟發(fā)展迅速的城市地區(qū)，因此降水對沙漠化的影響作用減弱，人類活動成為導(dǎo)致植被退化、加劇沙漠化的主導(dǎo)因素。

2.2 內(nèi)蒙古風(fēng)蝕區(qū)沙漠化程度空間分布特征

由內(nèi)蒙古風(fēng)蝕區(qū)2001—2018年RUEmer和RUEexr空間分布（圖4）可知，2001—2018年大興安嶺周圍地區(qū)以及河套平原地區(qū)為RUEmer值較高的地區(qū)，具有較好的生產(chǎn)能力，東部的呼倫貝爾草原、科爾沁草原次之，西部阿拉善盟地區(qū)為RUEmer值較低的地區(qū)，部分地區(qū)不具備生產(chǎn)能力;RUEexr值與RUEmer值的分布大致相同，其差異在于呼倫貝爾草原、科爾沁草原屬于較高地區(qū)，具備較好的恢復(fù)彈力，而中部地區(qū)屬于較低區(qū)域，恢復(fù)彈力較差。

沙漠化作為土地退化的一種形式，其發(fā)展程度可以通過土地生產(chǎn)力的評價結(jié)果來反演。根據(jù)李春娥等[16-17]的研究，對研究區(qū)的土地生產(chǎn)力進行等級劃分，并通過土地生產(chǎn)力的評價結(jié)果來劃分出不同等級的沙漠化程度（圖5）：①當RUEmer和RUEexr均小于0時，表明該區(qū)域既沒有生產(chǎn)能力，又不具備恢復(fù)彈力，屬于重度沙漠化地區(qū)。該類土地面積為1.72萬km2，占全區(qū)面積的1.75%，主要集中在內(nèi)蒙古西北部與蒙古國接壤的戈壁地帶，該地區(qū)長期處于低植被覆蓋狀態(tài)，生態(tài)環(huán)境惡劣。

②當RUEmer<0且 0

④當RUEmer>1或RUEexr>1時，表示該區(qū)域具備較好的生產(chǎn)能力或恢復(fù)彈力，屬于非沙漠化地區(qū)。該類土地面積為5.99萬km2，占全區(qū)面積的6.08%，零星分布于研究區(qū)。

2.3 內(nèi)蒙古風(fēng)蝕區(qū)沙漠化穩(wěn)定性分析

為進一步研究內(nèi)蒙古風(fēng)蝕區(qū)內(nèi)沙漠化的演變情況，分析其穩(wěn)定程度，該研究逐像元計算了研究區(qū)內(nèi)RUEr的變異系數(shù)CV，并進行等級劃分，結(jié)果表明，研究區(qū)內(nèi)RUEr的變異程度較大（圖6）。

非常穩(wěn)定地區(qū)面積為1.94萬km2，占全區(qū)面積的1.97%，其分布與重度沙漠化地區(qū)大致相同，主要為內(nèi)蒙古西部植被覆蓋度極低、不具備生產(chǎn)能力的戈壁地區(qū)。

穩(wěn)定地區(qū)和脆弱地區(qū)面積分別為9.97萬、25.99萬km2，分別占全區(qū)面積的10.13%、26.39%，兩者交錯分布于內(nèi)蒙古風(fēng)蝕區(qū)的東部地區(qū)，該地區(qū)主要植被類型為溫帶草原，錫林郭勒草原、烏蘭察布草原等分布于此，植被覆蓋度相對較高，加之近年來退耕還草的落實與普及，該地區(qū)沙漠化形勢有了一定的控制。

非常脆弱地區(qū)占據(jù)研究區(qū)大部分地區(qū)，再次表明研究區(qū)內(nèi)土地沙漠化形勢嚴峻，這類地區(qū)面積為60.58萬km2，占全區(qū)面積的61.51%，主要分布在中西部溫帶荒漠區(qū)，主要有巴丹吉林沙漠、庫布齊沙漠、烏蘭布和沙漠、毛烏素沙地等，該地區(qū)植被覆蓋較低，沙漠化形勢嚴峻，但RUEr有增加趨勢。此外，還有少數(shù)分布在東部的呼倫貝爾沙地和科爾沁沙地，該地區(qū)雖然植被覆蓋較高，但RUEr有下降趨勢，需要引起重視，及時采取措施。

整體而言，內(nèi)蒙古風(fēng)蝕區(qū)18年間有87.91%的地區(qū)RUEr變化較大，表現(xiàn)為脆弱到非常脆弱，僅有12.09%的地區(qū)表現(xiàn)為穩(wěn)定到非常穩(wěn)定。

3 結(jié)論

該研究通過計算內(nèi)蒙古風(fēng)蝕區(qū)2001—2018年生長季4—10月的相對降水利用率，分析其時間變化、空間分布及變異系數(shù)，綜合評價了

18年間內(nèi)蒙古風(fēng)蝕區(qū)的沙漠化演變情況，得出以下結(jié)論：

（1）整體上內(nèi)蒙古風(fēng)蝕區(qū)沙漠化呈現(xiàn)出微弱的加重趨勢，其空間分布上具有地域差異性，中西部地區(qū)RUEr的斜率為正，植被呈恢復(fù)狀態(tài)，降水為影響沙漠化的主導(dǎo)因素;東部地區(qū)斜率為負，植被出現(xiàn)退化現(xiàn)象，人類活動為影響沙漠化的主導(dǎo)因素。

（2）從沙漠化程度空間分布來看，研究區(qū)內(nèi)89.09%的地區(qū)為輕度沙漠化，重度沙漠化地區(qū)集中在內(nèi)蒙古西北部，中度沙漠化和非沙漠化地區(qū)在研究區(qū)內(nèi)零星分布。

（3）18年間內(nèi)蒙古風(fēng)蝕區(qū)的沙漠化變異程度較大，87.91% 的地區(qū)表現(xiàn)為脆弱到非常脆弱，僅有12.09%的地區(qū)表現(xiàn)為穩(wěn)定到非常穩(wěn)定。

參考文獻

[1]

郭瑞霞，管曉丹，張艷婷.我國荒漠化主要研究進展[J].干旱氣象，2015，33（3）：505-513.

[2]孫浩然，趙志根.基于SPOT NDVI數(shù)據(jù)的內(nèi)蒙古地區(qū)土地退化監(jiān)測研究[J].黑龍江工程學(xué)院學(xué)報，2019，33（3）：30-34.

[3]要永在.內(nèi)蒙古荒漠化驅(qū)動力及水土保持研究[J].中國農(nóng)業(yè)資源與區(qū)劃，2018，39（2）：13-17，87.

[4]李亞云，楊秀春，朱曉華，等.遙感技術(shù)在中國土地荒漠化監(jiān)測中的應(yīng)用進展[J].地理科學(xué)進展，2009，28（1）：55-62.

[5]陳寶瑞，辛?xí)云?，朱玉霞，?內(nèi)蒙古荒漠化年際動態(tài)變化及與氣候因子分析[J].遙感信息，2007，22（6）：39-44，104.

[6]朱玉霞，覃志豪，徐斌.基于MODIS數(shù)據(jù)的草原荒漠化年際動態(tài)變化研究：以內(nèi)蒙古自治區(qū)為例[J].中國草地學(xué)報，2007，29（4）：2-8.

[7]李云鵬，那日蘇，烏日娜，等.內(nèi)蒙古草原荒漠化遙感監(jiān)測分析[J].中國草地學(xué)報，2011，33（3）：79-86.

[8]殷賀，李正國，王仰麟，等.基于時間序列植被特征的內(nèi)蒙古荒漠化評價[J].地理學(xué)報，2011，66（5）：653-661.

[9]殷賀，李正國，王仰麟.荒漠化評價研究進展[J].植物生態(tài)學(xué)報，2011，35（3）：345-352.

[10]LAUENROTH W K，BURKE I C，PARUELO J M.Patterns of production and precipitation-use efficiency of winter wheat and native grasslands in the central great plains of the United States[J].Ecosystems，2000，3（4）：344-351.

[11]PARUELO J M，LAUENROTH W K，BURKE I C，et al.Grassland precipitation-use efficiency varies across a resource aradient[J].Ecosystems，1999，2（1）：64-68.

[12]VARNAMKHASTI A S，MILCHUNAS D G，LAUENROTH W K，et al.Production and rain use efficiency in short-grass steppe：Grazing history，defoliation and water resource[J].Journal of vegetation science，1995，6（6）：787-796.

[13]孫艷玲，郭鵬，延曉冬，等.內(nèi)蒙古植被覆蓋變化及其與氣候、人類活動的關(guān)系[J].自然資源學(xué)報，2010，25（3）：407-414.

[14]姬盼盼，高敏華，付曉紅，等.中國西北半干旱區(qū)凈初級生產(chǎn)力驅(qū)動因子空間計算分析[J].生態(tài)學(xué)報，2019，39（24）：9023-9032.

[15]林志東，武國勝.基于MODIS的大田縣均溪谷地NPP與NDVI相關(guān)性的時空變化特征[J].亞熱帶資源與環(huán)境學(xué)報，2015，10（1）：27-33.

[16]李春娥.新疆土地荒漠化時空變化特征分析[J].測繪科學(xué)，2018，43（9）：33-39.

[17]王紅巖.基于NPP和植被降水利用效率土地退化遙感評價與監(jiān)測技術(shù)研究[D].北京：中國林業(yè)科學(xué)研究院，2013：82-85.