周炎廣，武子豐，胡日娜，李紅悅，王卓然，哈斯額爾敦

毛烏素沙地新墾地土壤風(fēng)蝕特征

周炎廣，武子豐，胡日娜，李紅悅，王卓然，哈斯額爾敦※

（北京師范大學(xué)地理科學(xué)學(xué)部自然資源學(xué)院，北京 100875）

為了理解近期沙地農(nóng)業(yè)開發(fā)的土壤風(fēng)蝕風(fēng)險及其對區(qū)域沙漠化的影響。該文對毛烏素沙地灌叢沙堆地、沙質(zhì)草地以及耕地不同位置進行了多組風(fēng)沙同步觀測，并對所收集的沙物質(zhì)進行粒度分析。結(jié)果表明：草地與灌叢地被開墾為耕地后，土壤風(fēng)蝕強度成倍數(shù)增長。在可蝕性顆粒充足的耕地，輸沙率順風(fēng)向趨于增加的同時，隨高度的增加，輸沙率在0～10 cm高度內(nèi)以指數(shù)函數(shù)規(guī)律遞減，在10～20 cm高度內(nèi)以冪函數(shù)規(guī)律遞減；總輸沙率與風(fēng)速呈較好的冪函數(shù)關(guān)系，侵蝕模數(shù)約為9 657 t/(km2·a)。該文認(rèn)為，新墾地的出現(xiàn)對區(qū)域沙漠化具有重要影響，其面積應(yīng)作為評價區(qū)域沙漠化的重要指標(biāo)；減小裸露耕地面積，建立防風(fēng)固沙林，保留一定高度的留茬，可有效控制新墾地的土壤風(fēng)蝕問題。

風(fēng)；侵蝕；粒度；輸沙率；風(fēng)沙流結(jié)構(gòu)；毛烏素沙地

0 引 言

土壤風(fēng)蝕是干旱、半干旱區(qū)沙漠化的首要環(huán)節(jié)[1-4]，并受下墊面及氣候環(huán)境要素等[5-7]的共同影響。強烈的土壤風(fēng)蝕加速區(qū)域沙漠化的發(fā)展，影響人類生活質(zhì)量，進而引起嚴(yán)重的生態(tài)與環(huán)境問題[8-9]。毛烏素沙地位于中國北方農(nóng)牧交錯帶中部，對氣候變化和人類活動具有較高的敏感性和脆弱性[10-12]。在歷史時期曾經(jīng)歷了多次沙丘活化、沙漠?dāng)U張與沙丘固定、沙漠收縮。在過去幾十年前該區(qū)風(fēng)沙活動強烈、沙漠化比較嚴(yán)重。后來尤其是2000年以來，隨著政府實施多項生態(tài)治理工程（諸如京津風(fēng)沙源治理、三北防護林、退耕還林等），植被覆蓋度得到大幅度增加，大部分沙丘得以固定，沙漠化出現(xiàn)逆轉(zhuǎn)，整體的生態(tài)環(huán)境得到有效改善[13-15]。與此同時，隨著快速城市化和工礦交通建設(shè)占用農(nóng)田，以及農(nóng)業(yè)技術(shù)的發(fā)展，地下淡水資源較為豐富的毛烏素沙地成為農(nóng)業(yè)開發(fā)的主要場所。近年來，在毛烏素沙地東南部與西南部，大面積草地與灌叢地開墾為耕地的現(xiàn)象極為普遍，致使新墾地的面積快速增長[16]。這些開墾后的耕地在干旱多風(fēng)的冬、春季節(jié)多處于裸露狀態(tài)，易遭受強烈風(fēng)蝕。因此，研究新墾地的土壤風(fēng)蝕問題，并找出控制方法具有重要意義。本文以毛烏素沙地南部為研究區(qū)域，在該區(qū)的灌叢沙堆地、沙質(zhì)草地與新墾地的不同位置上進行氣流與集沙的同步觀測，試圖闡明新墾地的土壤風(fēng)蝕特征，揭示新墾地對區(qū)域沙漠化的影響。同時，通過分析耕地輸沙率與高度、風(fēng)速的關(guān)系，結(jié)合前人的研究成果，探討出減輕新墾地土壤風(fēng)蝕的有效方法與措施。

1 研究區(qū)概況與研究方法

1.1 研究區(qū)概況

研究區(qū)位于毛烏素沙地南部的鄂爾多斯市鄂托克前旗，地理坐標(biāo)為：38°05′27″～38°13′58″N，107°24′28″～107°37′30″E（圖1），地貌景觀以NW-SE向的帶狀沙地與平坦灘地相間分布為特色。該區(qū)屬溫帶半干旱氣候，冬春季干冷多風(fēng)，夏季溫暖濕潤。根據(jù)鄂托克前旗氣象站多年（1959－2015年）觀測數(shù)據(jù)的統(tǒng)計分析，年均溫7.9oC，平均年降水量254.3 mm，年蒸發(fā)量2 456.10 mm，年均風(fēng)速2.4 m/s，尤以4－5月風(fēng)速最高。起沙風(fēng)風(fēng)向主要集中于西北，且來自WNW方向所占比例最大。土壤類型以風(fēng)沙土、栗鈣土和草甸土為主。植被主要為沙地植被和一些草甸植被，以油蒿（）、北沙柳（）、藏錦雞兒（）、紫花苜蓿（）與芨芨草（）等為主。

注：圖中小框為觀測區(qū)域。

1.2 研究方法

1.2.1 觀測點位置

在研究區(qū)中，選取一片地勢平坦且具有多種下墊面類型的小塊區(qū)域。依據(jù)下墊面類型，依次分為灌叢地、草地與耕地3類。灌叢地的主要植被類型為油蒿（），平均高度約60 cm，蓋度達40%，且地表無明顯結(jié)皮現(xiàn)象；草地分布較為均勻，植被高度在10 cm左右，蓋度約50%，地表有近2 mm厚的物理-生物結(jié)皮；耕地為2011年后開墾種植玉米，地表平坦且均為松軟沙物質(zhì)。觀測期間，耕地裸露且無灌溉裝置。

沿順風(fēng)向，選取5個典型位置作為階梯式集沙儀的布設(shè)點位，分別用字母A、B、C、D及E表示（圖2）。A點選在距草地邊界一定距離的灌叢地中部；B點在耕地與灌叢地邊界等間距的草地中部；C、D、E代表了耕地上的不同位置，C點與耕地上邊界相距5m，D點與C點相距45m，E點與C點相距85m。風(fēng)沙觀測期間，在A、B、C、D、E位置及其附近進行多組重復(fù)。

注：A、B、C、D和E為階梯式集沙儀的布設(shè)點位。

1.2.2 氣流觀測

觀測儀器為已校準(zhǔn)的 EC 9-1型三杯風(fēng)速儀，與Campell公司生產(chǎn)的CR1000數(shù)據(jù)采集儀相連，該系統(tǒng)可每5 s自動采集并儲存觀測數(shù)據(jù)。選定風(fēng)沙活動強烈的大風(fēng)日，在D點位置0.2、0.3、0.6、1.0、1.2、1.5、2.0、3.0及4.0 m共9個高度處架設(shè)杯狀風(fēng)速儀，并在4.0 m高度處架設(shè)校準(zhǔn)后的風(fēng)向標(biāo)。本次野外觀測時間為2018年5月，共進行12組風(fēng)沙觀測。其中5組為順風(fēng)向A、B、C、D、E 5個位置同步進行的風(fēng)沙觀測，意在對比灌叢地、草地與耕地在同流情況下的輸沙情況，觀測時段分別為12:25－12:55、12:56－13:25、14:43－15:13、15:24－15:54及15:55－16:24。剩余7組為耕地內(nèi)部3個位置（C、D、E）的風(fēng)沙同步觀測，觀測時段分別為09:37－09:57、09:58－10:17、10:18－10:37、10:38－10:47、10:48－11:07、11:08－11:27、11:28－11:47。

1.2.3 輸沙率觀測與粒度測量

本文使用的階梯式集沙儀總高度為20 cm，共10個進沙口，進沙單口面積為2 cm×2 cm，集沙效率為90%[17-18]。集沙儀的進沙口中心距地面高度分別為1、3、5、7、9、11、13、15、17和19 cm。風(fēng)沙觀測時，所有集沙儀的底部與地面齊平，集沙與測風(fēng)同步進行。為保證集沙量的充足性與準(zhǔn)確性，集沙口的方向需依據(jù)當(dāng)時風(fēng)向做出相應(yīng)的調(diào)整，且每次集沙觀測時間至少20 min，同步觀測記錄多組數(shù)據(jù)。集沙儀所采集的沙粒樣品經(jīng)風(fēng)干處理后用千分之一電子天平（精度0.001 g）稱質(zhì)量記錄，并以觀測高度內(nèi)單位時間單位寬度的總輸沙量表示輸沙率。觀測開始前，在所設(shè)的5個觀測位置（A、B、C、D、E）采集表層沉積物采樣，且在每個位置上采取3組樣品，用以平均。觀測結(jié)束后，通過Malvern公司生產(chǎn)的Mastersizer2000型激光粒度儀，分別對研究區(qū)地表所采樣品與集沙儀所收集物質(zhì)進行粒度分析。

2 研究結(jié)果

圖3是觀測期間距地表2m高度處的風(fēng)向與風(fēng)速。風(fēng)向集中于NNW（326.26°～348.75°）與N（348.76°～11.25°）之間，且絕大多數(shù)在NNW方向，平均風(fēng)速集中于6～8 m/s。

注：各觀測時段以虛線劃分，并用羅馬數(shù)字表示。

耕地表層與對應(yīng)風(fēng)蝕物的粒徑頻率曲線（圖4）均呈逼近正態(tài)分布的單峰。耕地不同位置（C、D、E）表層沉積物粒徑級配均以中、細砂（64～500m）組分為主（占總組分的85%～95%），且眾數(shù)對應(yīng)的粒徑均介于150～250m之間，頻率在12%～16%之間。與其相對應(yīng)的風(fēng)蝕物粒徑級配也以中、細砂（64～500m）組分為主（占總組分的90%～96%），眾數(shù)對應(yīng)的粒徑介于150～210m之間，頻率在15%～18%之間。

圖4 耕地（C、D、E）表層及對應(yīng)風(fēng)蝕物粒徑頻率曲線

2.1 輸沙率的順風(fēng)向變化

圖5是同流條件下順風(fēng)向不同地表的輸沙率變化。耕地表面的輸沙率較高，且總體隨風(fēng)速增大而增大；而草地和灌叢地的輸沙率極低，與風(fēng)速的大小無明顯關(guān)系。這說明草地和灌叢地幾乎不遭受風(fēng)蝕，而裸露的耕地較易遭受侵蝕。在耕地，順風(fēng)向由C點至E點，輸沙率整體呈增大趨勢。在平均風(fēng)速為5.78、7.38與7.56 m/s時，耕地輸沙率由順風(fēng)向依次增加，且增幅較大；在平均風(fēng)速為6.46與7.20 m/s時，E點的輸沙率低于D點，但高于C點，這可能是與復(fù)雜多變的野外環(huán)境有關(guān)[19]。

表1為各地表輸沙率之間的對比關(guān)系。耕地與灌叢地輸沙率的比值為2.409～41.792，平均18.63；耕地與草地的比值2.605～41.952，平均17.95；草地與灌叢地的比值趨近于1。耕地與灌叢地和耕地與草地的輸沙率比值，隨著風(fēng)速的增大，呈指數(shù)函數(shù)規(guī)律增長（2：0.921，0.857）。

圖5 順風(fēng)向不同地表的輸沙率變化

表1 各地表輸沙率之間的對比關(guān)系

2.2 輸沙率的垂線變化

圖6是不同風(fēng)速條件下的耕地（C、D、E）輸沙率隨高度的變化圖。隨高度增加，輸沙率趨于減少。但是，以距地表10 cm高度為界，出現(xiàn)兩種遞減形式：以上顯示輸沙率與高度均呈現(xiàn)較好的直線關(guān)系，以下則呈明顯的曲線形態(tài)。這說明在10 cm高度上下，耕地輸沙率與高度之間的關(guān)系可能有兩種函數(shù)關(guān)系。

2.2.1 距地表10 cm高度內(nèi)的輸沙率變化

在耕地0～10 cm高度范圍內(nèi)的輸沙通量占總輸沙通量（0～20 cm高度內(nèi)）的90%以上。圖7是耕地（C、D、E）輸沙率與高度（范圍0～10 cm）的關(guān)系圖。在10 cm高度內(nèi)，輸沙率隨高度的增加而迅速減少，減少率范圍在80%～96%。近地表0～2 cm高度內(nèi)的風(fēng)蝕物含量變化較大，且無明顯規(guī)律，這與近地表沙粒復(fù)雜移動方式（蠕移和躍移）相關(guān)。在0～10 cm區(qū)間，耕地（C、D、E）的輸沙率與高度呈良好的指數(shù)函數(shù)關(guān)系。指數(shù)函數(shù)公式為

=ab（1）

式中為輸沙率，g/(cm·min)；為高度，cm；、為回歸系數(shù)。隨風(fēng)速的增大，氣流中的總輸沙量整體呈增大趨勢；在0～10 cm高度區(qū)間內(nèi)，輸沙率與高度的指數(shù)函數(shù)相關(guān)系數(shù)均達0.99以上，說明相關(guān)性極高。由此，指數(shù)函數(shù)關(guān)系式能夠充分且準(zhǔn)確反映10cm高度內(nèi)輸沙率與高度的關(guān)系。

圖6 耕地（C、D、E）輸沙率與高度的關(guān)系

圖7 耕地（C、D、E）輸沙率在0～10 cm高度內(nèi)的垂線變化

2.2.2 距地表10-20cm高度內(nèi)的輸沙率變化

圖8是高度在10～20cm范圍內(nèi)耕地（C、D、E）輸沙率與高度之間的關(guān)系。在10～20cm高度范圍內(nèi)，風(fēng)蝕物含量隨高度的增加而減少，但減少程度較??；隨風(fēng)速的增加，各層輸沙率呈增加的趨勢。在10～20 cm高度內(nèi)，風(fēng)蝕物含量隨高度的分布更趨近于冪函數(shù)形式。冪函數(shù)公式為

=ax（2）

式中為輸沙率，g/(cm·min)；為高度，cm；、為回歸系數(shù)。在10～20 cm高度內(nèi)，輸沙率與高度的冪函數(shù)相關(guān)系數(shù)（2）均達0.93以上，相關(guān)性極高。由此，冪函數(shù)關(guān)系式能夠充分反映10～20cm高度范圍內(nèi)的風(fēng)蝕物含量隨高度分布的實際。

圖8 耕地（C、D、E）輸沙率在10～20 cm高度內(nèi)的垂線變化

2.2.3 總輸沙率的計算

利用上述經(jīng)驗公式（1）和（2）模擬出的輸沙率值與實測值之間的對比結(jié)果顯示（圖9），模擬公式可以較為準(zhǔn)確地模擬耕地風(fēng)沙流輸沙率與高度的關(guān)系。至于近地表1 cm處的輸沙率值與函數(shù)模擬值具有較大的差異，說明貼近地表處沙物質(zhì)移動方式的復(fù)雜性。模擬后的輸沙率與高度的關(guān)系在10 cm高度處上下具有明顯的改變。在0～10 cm以內(nèi)，輸沙率與高度呈指數(shù)函數(shù)規(guī)律；而在10～20 cm高度內(nèi)，輸沙率與高度呈冪函數(shù)規(guī)律。函數(shù)公式所算出的模擬值與集沙儀所測的實際數(shù)據(jù)高度吻合，說明依據(jù)模擬函數(shù)可以準(zhǔn)確計算出不同高度上的輸沙率值。

注：C-1指在觀測點C上進行的第一次測量時，所得輸沙率實測值與模擬值的對比關(guān)系，余同。v為平均風(fēng)速。

本文所用集沙儀高度為20 cm，低于風(fēng)沙層高度。因此，為得出耕地近地面層的總輸沙率，需通過上述模擬結(jié)果進行計算，計算公式與步驟如下

式中()是高度在時的輸沙率，g/(cm·min)，為高度，cm，數(shù)值0.9為集沙儀的集沙效率，為回歸系數(shù)。通過公式（3），可計算出耕地不同位置（C、D、E）上每組積沙的總輸沙率。

2.3 總輸沙率與風(fēng)速的關(guān)系

圖10是耕地（C、D、E）上總輸沙率與風(fēng)速之間的關(guān)系圖。其中輸沙率是通過模擬函數(shù)公式（3），計算所得的總輸沙率，風(fēng)速為每組觀測期間2 m處的平均風(fēng)速。耕地不同位置的總輸沙率與平均風(fēng)速呈較好的冪函數(shù)關(guān)系，相關(guān)系數(shù)（2）均在0.86以上，說明隨風(fēng)速的增加，輸沙率顯著增高。同時，順風(fēng)向距耕地上邊界越遠，輸沙率隨風(fēng)速增長的增幅更大?；诳傒斏陈逝c平均風(fēng)速間的關(guān)系公式，可通過平均風(fēng)速而反推出耕地各位置在此時期的總輸沙率。

圖10 耕地（C、D、E）上總輸沙率與風(fēng)速的關(guān)系

2.4 耕地風(fēng)蝕模數(shù)的估算

通過風(fēng)速與總輸沙率的關(guān)系（圖10），結(jié)合鄂托克旗氣象站的風(fēng)速數(shù)據(jù)，可初步推算出該片新墾地的風(fēng)蝕模數(shù)。因區(qū)域氣象站所測得的風(fēng)速數(shù)據(jù)為10 m處曠野風(fēng)速，而文中各類分布關(guān)系均是以2 m高度處平均風(fēng)速計算，故在使用氣象站風(fēng)速數(shù)據(jù)前，需將10 m處風(fēng)速轉(zhuǎn)換為2 m處風(fēng)速。本次野外連續(xù)觀測了9個不同高度處的平均風(fēng)速，可通過下式推算10 m高處風(fēng)速

式中*是摩阻風(fēng)速，0.2、2、10分別表示高度為0.2、2及10 m高度處的平均風(fēng)速。由公式（4），結(jié)合0.2 m與2 m處實際測量的多組數(shù)據(jù)，推出對應(yīng)10 m高處的平均風(fēng)速。將所得10 m處平均風(fēng)速與2 m處平均風(fēng)速進行擬合，可得公式（5），擬合系數(shù)2等于0.988。

因此，研究區(qū)的2 m處風(fēng)速與10 m處風(fēng)速可由公式（5）進行轉(zhuǎn)換。

耕地在5月底至10月底常種有玉米作物，在此時期耕地幾乎不發(fā)生風(fēng)蝕。因而估算耕地的風(fēng)蝕模數(shù)時，應(yīng)重點考慮耕地處于裸露狀態(tài)時的風(fēng)蝕量，即冬春季期間。耕地下風(fēng)向（C點）與耕地上風(fēng)向（E點）之間相距85m，用E點位置的風(fēng)蝕量減去C點位置的風(fēng)蝕量，便可得出兩位置間風(fēng)蝕物的流失量，計算公式如下。

式中()是平均風(fēng)速在時的輸沙率，g/(cm·min)；為平均風(fēng)速。若該片耕地長期保持此地表狀態(tài)，那么便可通過公式（6），結(jié)合該區(qū)氣象站的風(fēng)速數(shù)據(jù)，初步推算出該片耕地一年的風(fēng)蝕量。結(jié)果顯示，該片新墾地風(fēng)蝕模數(shù)為9 657 t/(km2·a)。

3 討 論

耕地風(fēng)蝕強度的大小受區(qū)域風(fēng)速、土壤含水量、不可蝕性顆粒含量、耕作方式等[9,20-21]多種因素的影響。本文單一考慮了區(qū)域風(fēng)速對耕地風(fēng)蝕強度的影響，而未考慮其他因素，故文中所用風(fēng)蝕模型具有較大的局限性，僅能反映所研究耕地的風(fēng)蝕狀況。盡管如此，通過該模型計算所得的風(fēng)蝕模數(shù)已能充分反映出新墾地風(fēng)蝕強烈。近年來，對不同區(qū)域的農(nóng)田風(fēng)蝕量問題開展了各類模型研究[22-26]。王仁德等[24]通過粒度對比法，測得河北壩上無覆蓋物農(nóng)田2013年的風(fēng)蝕模數(shù)為960～5 700 t/(km2·a)，平均為2 852.14 t/(km2·a)。李勝龍等[25]通過修正風(fēng)蝕方程（RWEQ）估算得出,東北地區(qū)的壟作無覆蓋農(nóng)田風(fēng)蝕模數(shù)高達181.7～86 582.9 t/(km2·a)。遲文峰等[26]同樣使用（RWEQ）方法估算得出，內(nèi)蒙古西部2015年土壤風(fēng)蝕模數(shù)為5 946 t/(km2·a)。綜上可知，無覆蓋耕地在不同區(qū)域和狀態(tài)條件下均表現(xiàn)為中強度侵蝕，其土壤風(fēng)蝕問題極為嚴(yán)重。

近幾十年來，毛烏素沙地的耕地面積大幅度增長[11,27]，且仍有加劇增長的趨勢。據(jù)研究[16]，從2010年到2015年，毛烏素沙地的耕地面積共增加了411.29 km2。究其原因，主要包括以下3個方面：一是區(qū)域政府的占補平衡政策[28-29]，即城市建設(shè)所占用的耕地，需通過增加新墾地得到補充；二是農(nóng)業(yè)生產(chǎn)技術(shù)與水平的提高，隨著各類噴灌、滴灌技術(shù)[30-31]的普及，半干旱區(qū)耕種已成為一種高效而經(jīng)濟的選擇；三是禁牧政策的間接影響，草場被封育后，牧民需開墾新耕地用于種植畜牧飼料。觀測結(jié)果表明，灌叢地與草地被開墾為耕地后，風(fēng)蝕強度增強，輸沙率呈倍數(shù)增長。大面積開墾耕地不僅直接導(dǎo)致土壤風(fēng)蝕問題的加重，同時還會間接引發(fā)眾多環(huán)境問題，如地下水位下降及水污染問題等[32]。因此，耕地面積的大小對區(qū)域沙漠化程度的強弱具有重要影響。目前，國內(nèi)學(xué)者對區(qū)域沙漠化監(jiān)測評價的指標(biāo)做過較多研究[33-34]，構(gòu)建了多種沙漠化評價體系，但所選指標(biāo)往往集中在裸沙面積與植被蓋度等，并未考慮耕地面積的影響，故而可能低估了區(qū)域沙漠化程度。同時，在進行遙感解譯和監(jiān)測時，常忽略了影像時間的問題。若影像在夏、秋季期間，此時耕地上具有作物覆蓋，對區(qū)域沙漠化的影響較低；而在冬、春季期間，耕地上幾乎無作物覆蓋，對區(qū)域沙漠化的影響巨大。由此，該文認(rèn)為在評價區(qū)域沙漠化程度時，應(yīng)將耕地面積作為重要評價指標(biāo)，并將各季節(jié)分開進行評價。

風(fēng)蝕物含量隨高度的變化規(guī)律，在不同的高度區(qū)間具有一定差異，這與風(fēng)蝕物質(zhì)的運動方式轉(zhuǎn)變具有密切聯(lián)系。Fryrear等[35-37]在美國得克薩斯州大斯普林土壤風(fēng)蝕研究站長期觀測研究后認(rèn)為，躍移層輸沙量隨高度的分布可用指數(shù)函數(shù)表示，而懸移層則可用冪函數(shù)表示。哈斯[38]在河北壩上高原研究發(fā)現(xiàn)，風(fēng)蝕物含量隨高度的增加，在0～20 cm高度內(nèi)以指數(shù)規(guī)律遞減，反映了以躍移質(zhì)為主的風(fēng)沙流結(jié)構(gòu)；而在20～100 cm高度內(nèi)以冪函數(shù)規(guī)律遞減，反映了以懸移質(zhì)為主的風(fēng)沙流結(jié)構(gòu)。本文研究中，輸沙率與高度的關(guān)系同樣具有分段性，在0～10 cm高度區(qū)間內(nèi)服從指數(shù)函數(shù)分布，在10～20 cm高度區(qū)間服從冪函數(shù)規(guī)律。文中的分段處在10 cm高度處，與哈斯[38]研究的20 cm高度處分段具有一定的差異，這可能是由于所研究區(qū)域下墊面性質(zhì)的差異而導(dǎo)致。僅使用以躍移運動為主的指數(shù)函數(shù)對總輸沙率與高度的關(guān)系進行模擬，將無法準(zhǔn)確統(tǒng)計到懸移物質(zhì)的含量。文中通過分段函數(shù)計算耕地上的總輸沙率值是更接近實際情況的一種方式。

新墾地土壤風(fēng)蝕問題的根本原因在于干旱且大風(fēng)日多的冬春季，耕地常處在裸露狀態(tài)，沒有基本的防風(fēng)固沙措施。因此，對于新墾地本身的防護和管理也應(yīng)引起足夠的重視。本文通過所得結(jié)果分析，結(jié)合前人的研究[39-43]，提出以下幾點基本的防護措施。首先，減少對草地與灌叢地的破壞，控制新墾地的開墾面積。區(qū)域政府應(yīng)高度重視并把控住耕地的開墾面積，不應(yīng)以犧牲環(huán)境來促進發(fā)展。其次，在耕地周圍建立科學(xué)配置后的防風(fēng)固沙林[39-41]。因在可蝕性物質(zhì)充足的耕地上，風(fēng)蝕物含量隨風(fēng)速的增大呈倍數(shù)增長，故降低耕地周圍風(fēng)速是解決風(fēng)蝕問題的必要手段。最后，保留適當(dāng)高度的留茬和殘存物[42-43]。耕地留茬和殘存物可增加耕地表面的粗糙度，改變近地表的氣流狀況，同時阻攔沙物質(zhì)的流失，進而減輕土壤侵蝕危害。

4 結(jié) 論

根據(jù)分析結(jié)果及討論，得出以下結(jié)論：

1）沙質(zhì)草地和灌叢沙丘地被開墾后，土壤風(fēng)蝕強度成倍增長。

2）耕地上，輸沙率隨高度的增加而減小，在0～10 cm高度內(nèi)服從指數(shù)分布，在10～20 cm高度內(nèi)服從冪函數(shù)規(guī)律；沙源充足的耕地上，隨風(fēng)速的增大，總輸沙率呈冪函數(shù)形式增長。

3）毛烏素沙地新墾地的出現(xiàn)對區(qū)域沙漠化具有重要影響，其面積應(yīng)被考慮到區(qū)域沙漠化的評價體系中。

4）控制耕地面積，建立防風(fēng)防沙林，保留適宜高度的留茬，可有效減輕土壤風(fēng)蝕危害。

[1]董光榮，李長治，金炯，等. 關(guān)于土壤風(fēng)蝕風(fēng)洞模擬實驗的某些結(jié)果[J]. 科學(xué)通報，1987，32(4)：297－301.

Dong Guangrong, Li Changzhi, Jin Jiong, et al. Some results of wind tunnel simulation of soil erosion[J]. Science Bulletin, 1987, 32(4): 297－301. (in Chinese with English abstract)

[2]陳智，麻碩士，趙永來，等. 保護性耕作農(nóng)田地表風(fēng)沙流特性[J]. 農(nóng)業(yè)工程學(xué)報，2010，26(1)：118－122.

Chen Zhi, Ma Shuoshi, Zhao Yonglai, et al. Characteristics of drifting sand flux over conservation tillage field[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(1): 118－122. (in Chinese with English abstract)

[3]董治寶，李振山，嚴(yán)平. 國外土壤風(fēng)蝕的研究歷史與特點[J]. 中國沙漠，1995，15(1)：100－104.

Dong Zhibao, Li Zhengshan, Yan Ping. An outline of the wind erosion research history in the world[J]. Journal of Desert Research, 1995, 15(1): 100－104. (in Chinese with English abstract)

[4]楊彩紅，馮福學(xué)，柴強，等. 小麥玉米田耕作模式的防風(fēng)蝕效果[J]. 中國沙漠，2019，39(4)：9－15.

Yang Caihong, Feng Fuxue, Chai Qiang, et al. Soil wind erosion in wheat/maize intercropping with zero tillage and straw retention[J]. Journal of Desert Research, 2019, 39(4): 9－15. (in Chinese with English abstract)

[5]劉玉璋，董光榮，李長治. 影響土壤風(fēng)蝕主要因素的風(fēng)洞實驗研究[J]. 中國沙漠，1992，12(4)：41－49.

Liu Yuzhang, Dong Guangrong, Li Changzhi. Study on some factors influencing soil erosion by wind tunnel experiment[J]. Journal of Desert Research, 1992, 12(4): 41－49. (in Chinese with English abstract)

[6]孫悅超，麻碩士，陳智，等. 植被蓋度和殘茬高度對保護性耕作農(nóng)田防風(fēng)蝕效果的影響[J]. 農(nóng)業(yè)工程學(xué)報，2010，26(8)：156－159.

Sun Yuechao, Ma Shuoshi, Chen Zhi, et al. Influences of vegetation coverage and residual stubble height on preventing resistance to wind erosion effect of conservation tillage farmland[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(8): 156－159. (in Chinese with English abstract)

[7]王萍，王燕. 民勤荒漠-綠洲過渡帶不同下墊面條件的土壤風(fēng)蝕特征[J]. 農(nóng)業(yè)工程學(xué)報，2012，28(11)：138－145.

Wang Ping, Wang Yan. Characteristics of aeolian sediment under different underlying surfaces in oasis-desert transitional region of Minqin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(11): 138－145. (in Chinese with English abstract)

[8]南嶺，董治寶，肖鋒軍. 農(nóng)牧交錯帶農(nóng)田土壤風(fēng)蝕PM10釋放特征[J]. 中國沙漠，2017，37(6)：1079－1084.

Nan Ling, Dong Zhibao, Xiao Fengjun. Characteristic of PM10 emission of cultivated soil in ecotone between agriculture and animal husbandry[J]. Journal of Desert Research, 2017, 37(6): 1079－1084. (in Chinese with English abstract)

[9]哈斯，陳渭南. 耕作方式對土壤風(fēng)蝕的影響：以河北壩上地區(qū)為例[J]. 水土保持學(xué)報，1996，2(1)：10－16.

Ha Si，Chen Weinan. Effect of tillage praetices on soil wind erosion: Taking Bashang area of Hebei province as an example[J]. Journal of Soil and Water Conservation, 1996, 2(1): 10－16. (in Chinese with English abstract)

[10]張新時. 毛烏素沙地的生態(tài)背景及其草地建設(shè)的原則與優(yōu)化模式[J]. 植物生態(tài)學(xué)報，1994，18(1)：1－16.

Zhang Xinshi. Principles and optimal models for development of Mu Us sandy grassland[J]. Acta Phytoecologica Sinica, 1994, 18(1): 1-16. (in Chinese with English abstract)

[11]閆峰，吳波. 近40 a毛烏素沙地沙漠化過程研究[J]. 干旱區(qū)地理，2013，36(6)：987－996.

Yan Feng, Wu Bo. Desertification progress in Mu Us Sandy Land over the past 40 years[J]. Arid Land Geography, 2013, 36(6): 987－996. (in Chinese with English abstract)

[12]徐小玲，延軍平. 毛烏素沙地的脆弱性與可持續(xù)發(fā)展研究[J]. 干旱區(qū)研究，2004，21(3)：286－289.

Xu Xiaoling, Yan Junping. Study on the vulnerability and the Sustainable development in the Mu Us sandy land[J]. Arid Zone Research, 2004, 21(3): 286－289. (in Chinese with English abstract)

[13]李曉嵐. 毛烏素沙地沙漠化逆轉(zhuǎn)過程及成因分析[D]. 西安：陜西師范大學(xué)，2017.

Li Xiaolan. The process and cause of desertification reversal in Maowusu sandy land[D]. Xi'an: Shaanxi Normal University, 2017. (in Chinese with English abstract)

[14]吳徽. 近50年來毛烏素沙地的沙漠化過程研究[J]. 中國沙漠，2001，21(2)：164－169.

Wu Hui. Study on proeess of desertiflcation in Mu Us sandy land for last 50 years , China[J]. Journal of Desert Research, 2001, 21(2): 164－169. (in Chinese with English abstract)

[15]郭堅，王濤，韓邦帥，等. 近30a來毛烏素沙地及其周邊地區(qū)沙漠化動態(tài)變化過程研究[J]. 中國沙漠，2008，28(6)：1017－1021.

Guo Jian, Wang Tao, Han Bangshuai, et al. Dynamic process of aeolian desertification land variation in Mu Us sandy land and its surrounding area in recent 30 years[J]. Journal of Desert Research, 2008, 28(6): 1017－1021. (in Chinese with English abstract)

[16]Li Sen, Yan Changzhen, Wang Tao, et al. Monitoring grassland reclamation in the Mu Us Desert using remote sensing from 2010 to 2015[J]. Environmental Earth Sciences, 2019, 78(10): 1－9.

[17]Wang Xunming, Dong Zhibao, Qu Jianjun, et al. Dynamic processes of a simple linear dune: A study in the Taklimakan Sand Sea, China[J]. Geomorphology, 2003, 52(3/4): 233－241.

[18]Chen W, Yang Z, Dong Z, et al. Vertical distribution of windblown sand flux in the surface layer during sand storms in the Taklimakan desert, Central Asia[J]. Physical Geography, 1996, 17(3): 193－218.

[19]黃曉琪. 戈壁地表風(fēng)蝕過程研究[D]. 北京：北京師范大學(xué)，2019.

Huang Xiaoqi. Research on Gobi Surface Wind Erosion Process[D]. Beijing: Beijing Normal University, 2019. (in Chinese with English abstract)

[20]董治寶，陳渭南，李振山，等. 風(fēng)沙土水分抗風(fēng)蝕性研究[J]. 水土保持通報，1996，16(2)：17－23.

Dong Zhibao, Chen Weinan, Li Zhenshan, et al. On resistance of eolian sandy moisture to wind erosion[J]. Bulletin of Soil and Water Conservation,1996, 16(2): 17－23. (in chinese with English abstract)

[21]蒙仲舉，高永，王淮亮，等. 基于土壤粒度和大風(fēng)日數(shù)的風(fēng)蝕風(fēng)險預(yù)報[J]. 農(nóng)業(yè)工程學(xué)報，2015，31(6)：186－192.

Meng Zhongju, Gao Yong, Wang Huailiang, et al. Wind erosion prediction based on soil particle size characteristics and strong wind days[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(6): 186－192. (in Chinese with English abstract)

[22]Zhang Chunlai, Zou Xueyong, Yang Ping, et al. Wind tunnel test and137Cs tracing study on wind erosion of several soils in Tibet[J]. Soil and Tillage Research, 2007, 94(2): 269－282.

[23]Chi Wenfeng, Zhao Yuanyuan, Kuang Wenhui, et al. Impacts of anthropogenic land use/cover changes on soil wind erosion in China[J]. Science of The Total Environment, 2019, 688(2): 204－215.

[24]王仁德，常春平，彭帥，等. 基于粒度對比法的壩上農(nóng)田風(fēng)蝕與粉塵釋放量估算[J]. 農(nóng)業(yè)工程學(xué)報，2013，29(21)：108－114.

Wang Rende, Chang Chunping, Peng Shuai, et al. Estimation on farmland wind-erosion and dust emission amount in Bashang of Hebei province by grain composition contrast[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(21): 108－114. (in Chinese with English abstract)

[25]李勝龍，李和平，林藝，等. 東北地區(qū)不同耕作方式農(nóng)田土壤風(fēng)蝕特征[J]. 水土保持學(xué)報，2019，33(4)：110－118.

Li Shenglong, Li Heping, Lin Yi, et al. Effects of tillage methods on wind erosion in farmland of northeastern China[J]. Journal of Soil and Water Conservation, 2019, 33(4): 110－118. (in Chinese with English abstract)

[26]遲文峰，白文科，劉正佳，等. 基于RWEQ模型的內(nèi)蒙古高原土壤風(fēng)蝕研究[J]. 生態(tài)環(huán)境學(xué)報，2018，27(6)：38－47.

Chi Wenfeng, Bai Wenke, Liu Zhengjia, et al. Wind erosion in Inner Mongolia Plateau using the Revised Wind Erosion Equation[J]. Ecology and Environmental Sciences, 2018, 27(6): 38－47. (in Chinese with English abstract)

[27]劉娟，劉華民，卓義，等. 毛烏素沙地1990－2014年景觀格局變化及驅(qū)動力[J]. 草業(yè)科學(xué)，2017，34(2)：255－263.

Liu Juan, Liu Huamin, Zhuo Yi, et al. Dynamics and driving forces of landscape patterns in Mu Us sandy land, from1990 to 2014[J]. Pratacultural Science, 2017, 34(2): 255－263. (in Chinese with English abstract)

[28]Liu Luo, Liu Zhenjie, Gong Jianzhou, et al. Quantifying the amount, heterogeneity, and pattern of farmland: Implications for China’s requisition-compensation balance of farmland[J]. Land Use Policy, 2019, 81: 256－266.

[29]Liu Y, Fang F, Li Y. Key issues of land use in China and implications for policy making[J]. Land Use Policy, 2014, 40: 6－12.

[30]Sui J, Wang J, Gong S, et al. Assessment of maize yield-increasing potential and optimum N level under mulched drip irrigation in the Northeast of China[J]. Field Crops Research, 2018, 215: 132－139.

[31]Peng Yaoqi, Xiao Yingxin, Fu Zetian, et al. Precision irrigation perspectives on the sustainable water-saving of field crop production in China: Water demand prediction and irrigation scheme optimization[J]. Journal of Cleaner Production, 2019, 230: 365－377.

[32]石輝，劉秀花，陳占飛，等. 陜北榆林毛烏素沙地大規(guī)模土地整治開發(fā)的生態(tài)環(huán)境問題及其對策[J]. 生態(tài)學(xué)雜志，2019，38(7)：2228－2235.

Shi Hui, Liu Xiuhua, Chen Zhanfei, et al. Ecoenvironmental problems and their solution strategy for largescale land consolidation and development in Mu Us Sandy Land of Yulin in North Shaanxi[J]. Chinese Journal of Ecology, 2019, 38(7): 2228－2235. (in Chinese with English abstract)

[33]高尚武，王葆芳，朱靈益，等. 中國沙質(zhì)荒漠化土地監(jiān)測評價指標(biāo)體系[J]. 林業(yè)科學(xué)，1998，34(2)：1－10.

Gao Shangwu, Wang Baofang, Zhu Lingyi, et al. Monitoring and evaluation indicator system on sandy desertification of china[J]. Scientia Silvae Sinicae, 1998, 34(2): 1－10. (in Chinese with English abstract)

[34]王濤，吳薇. 沙質(zhì)荒漠化的遙感監(jiān)測與評估：以中國北方沙質(zhì)荒漠化區(qū)內(nèi)的實踐為例[J]. 第四紀(jì)研究，1998(2)：108－118.

Wang Tao, Wu Wei. Remote sensing monitoring and assessing sandy desertification: An example from the sandy desertification region of north China [J]. Quaternary Sciences, 1998(2): 108－118. (in Chinese with English abstract)

[35]Fryrear D W, Saleh A. Field wind erosion-vertical- distribution[J]. Soil Science, 1993, 155: 294－300.

[36]Vories E D, Fryrear D W. Vertical-distribution of wind-eroded soil over a smooth, bare field[J]. Transactions of the ASAE, 1991, 34: 1763－1768.

[37]吳正等. 風(fēng)沙地貌與治沙工程學(xué)[M]. 北京：科學(xué)出版社，2003：63－65.

[38]哈斯. 河北壩上高原土壤風(fēng)蝕物垂直分布的初步研究[J]. 中國沙漠，1997，17(1)：9－14.

Ha Si. Preliminary study on the vertical distributions of wind dust over Bashang Plateau, Hebei Province[J]. Journal of Desert Research, 1997, 17(1): 9－14. (in Chinese with English abstract)

[39]孫悅超，麻碩士，陳智. 保護性耕作農(nóng)田和檸條帶狀配置草地防風(fēng)蝕效果的風(fēng)洞測試[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(11)：140－146.

Sun Yuechao, Ma Shuoshi, Chen Zhi. Wind tunnels test on anti-wind erosion effect of conservation tillage farmland and caragana collocated with grassland in strip[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(11): 140－146. (in Chinese with English abstract)

[40]左忠，潘占兵，張安東，等. 干旱風(fēng)沙區(qū)農(nóng)田防護林網(wǎng)空間風(fēng)速與地表風(fēng)蝕特征[J]. 農(nóng)業(yè)工程學(xué)報，2018，34(2)：135－141.

Zuo Zhong, Pan Zhanbing, Zhang Andong, et al. Spatial wind speed and surface wind erosion characteristics of farm-shelter forest network in arid sandy area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(2): 135－141. (in Chinese with English abstract)

[41]厲靜文，劉明虎，郭浩，等. 防風(fēng)固沙林研究進展[J]. 世界林業(yè)研究，2019，32(5)：28－33.

Li Jingwen, Liu Minghu, Guo Hao, et al. Research progress in windbreak and sand fixation forest[J]. World Forestry Research, 2019, 32(5): 28－33. (in Chinese with English abstract)

[42]Agen L J. Crop residue effects on aerodynamic processes and wind erosion[J]. Theoretical & Applied Climatology, 1996, 54(1/2): 39－46.

[43]Wang E, Harman W L, Williams J R, et al. Simulated effects of crop rotations and residue management on wind erosion in Wuchuan, West-Central Inner Mongolia, China[J]. Journal of Environmental Quality, 2002, 31(4): 1240－1247.

Characteristics of soil wind erosion in new reclaimation land of Mu Us sandy land, China

Zhou Yanguang, Wu Zifeng, Hu Rina, Li Hongyue, Wang Zhuoran, Hasi Eerdun※

(,,,100875)

Soil wind erosion is the main cause of land degradation in arid and semi-arid areas in northern China, which is mainly affected by surface type, climate and so on. Strong soil wind-erosion will not only accelerate regional desertification and influence the quality of human life, but also cause serious ecological and environmental problems. Since the 2000 s, our country has gradually strengthened the control of land desertification and carried out a series of sand control works, which effectively improved the overall ecological environment of the sandy land. But in recent years, there are some common phenomenon that extensive grassland and brushland was turned into cultivated land in the southeast and southwest of Mu Us sandy land. However, in the dry and windy winter or spring, these new reclaimation land is usually so exposed that it is vulnerable to suffer from strong wind erosion. Therefore, it is great significance to study the characteristics of soil wind erosion in new reclaimation land and find out the control methods. Based on the above question, this paper selected the southern part of Mu Us sandy land (38°05′27″-38°13′58″N, 107°24′28″-107°37′30″E) as the research area, as well as set up 5 observation points on different surfaces, which were represented respectively by A, B, C, D and E. Point A was on the brushland, point B was on the middle of the grassland, and other three points were respectively on different parts of the cultivated land. On the day with strong wind activity, 12 groups simultaneous observation of wind and sand were made at these 5 points (A, B, C, D, E), with each observation time at least 20 minutes. After these observations, the material of the sand samplers were weighed with an electronic balance of 1/1 000, and the material of surface and sand samplers on cultivated land were analyzed by Mastersizer2000 Laser Particle Size Analyzer. The results showed that the soil wind erosion intensity increased exponentially after grassland and shrubland were cultivated; on cultivated land, with the increase of height, the sediment transport rate decreases with exponential function within 0-10 cm height and power function within 10-20 cm height. Due to the height of the sand collector used in this paper is 20cm, lower than the height of the sand flow layer. In order to obtain the total sediment transport rate of the surface of the cultivated land, this paper must calculate it through the function simulation obtained from the above result. After calculating the total sediment transport rate, fitting it with the wind speed. We found that the total sediment transport rate has a good power function relationship with the wind speed in the cultivated land with sufficient erodibility particles. Based on the formula of the relationship between wind speed and total sediment transport rate, as well as the wind speed data of the meteorological stations in the region, the wind erosion modulus of the new reclaimation land can be preliminarily calculated. The result shows that the wind erosion modulus of the new reclaimed land is 9 657t/(km2·a). Through the results obtained in this paper, combined with the former research results, this paper argues that the emergence of new reclaimed land has an important impact on regional desertification, and its area should be used as an important index when evaluating regional desertification. Meanwhile, Controlling cultivated land area, establishing windbreak and sand fixation forest, as well as retaining stubble at a certain height can effectively control soil wind erosion in newly reclaimed land.

wind; erosion; particle size; sediment transport rate; structure of drifting sand flux; Mu Us sandy land

周炎廣，武子豐，胡日娜，李紅悅，王卓然，哈斯額爾敦. 毛烏素沙地新墾地土壤風(fēng)蝕特征[J]. 農(nóng)業(yè)工程學(xué)報，2020，36(1)：138－147.doi：10.11975/j.issn.1002-6819.2020.01.016 http://www.tcsae.org

Zhou Yanguang, Wu Zifeng, Hu Rina, Li Hongyue, Wang Zhuoran, Hasi Eerdun. Characteristics of soil wind erosion in new reclaimation land of Mu Us sandy land, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(1): 138－147. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.01.016 http://www.tcsae.org

2019-08-27

2019-11-21

國家重點研發(fā)計劃課題（2016YFC0500805）；國家自然科學(xué)基金項目（41671002）

周炎廣，博士生，主要從事干旱區(qū)地貌研究。Email：201521190011@mail.bnu.edu.cn

哈斯額爾敦，教授，博士生導(dǎo)師，主要從事干旱區(qū)地貌與環(huán)境研究。Email：hasi@bnu.edu.cn

10.11975/j.issn.1002-6819.2020.01.016

S157.1; S157.3

A

1002-6819(2020)-01-0138-10