基于ＲＷＥＱ模型的石渠縣土壤風(fēng)蝕評價及影響因素分析

楊雪琴　丁銳　謝駿鍇　劉恩勤

摘要：應(yīng)用修正風(fēng)蝕模型（RWEQ）對四川省甘孜藏族自治州石渠縣進行土壤風(fēng)蝕模數(shù)估算，分析土壤風(fēng)蝕在石渠縣的空間分布特征及各個因素對石渠縣土壤風(fēng)蝕作用的影響。結(jié)果表明，石渠縣北部抗土壤風(fēng)蝕能力優(yōu)于南部，西南部受土壤風(fēng)蝕影響最大;氣象因子為土壤風(fēng)蝕最主要的影響因子，在6月汛期時土壤風(fēng)蝕模數(shù)達到最大;植被指數(shù)的空間分布特征與土壤風(fēng)蝕模數(shù)空間分布特征契合，對其影響較深。

關(guān)鍵詞：RWEQ模型;土壤風(fēng)蝕;防風(fēng)固沙;石渠縣

中圖分類號：S157.1? ? ? ? ?文獻標(biāo)識碼：A

文章編號：0439-8114（2019）15-0026-04

DOI：10.14088/j.cnki.issn0439-8114.2019.15.006? ? ? ? ? ?開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

Evaluation of wind erosion and analysis of its influencing factors

in Shiqu county based on REWQ model

YANG Xue-qin，DING Rui，XIE Jun-kai，LIU En-qin

（College of Earth Science，Chengdu University of Technology/Key Laboratory of Information Technology & Application of Land and Resource，Chengdu 610059，China）

Abstract： The modified wind erosion model （RWEQ） was used to estimate the soil wind erosion modulus of Shiqu county， Ganzi prefecture， Sichuan province， and the spatial distribution characteristics of soil wind erosion and the effects of various factors on the soil wind erosion in Shiqu county were analyzed. The results showed that， the resistance to soil wind erosion in northern Shiqu county was better than that in southern， and the southwest was the most affected by soil wind erosion. Meteorological factors were the most important factors affecting soil wind erosion. The soil wind erosion modulus reached the maximum in the flood season in June. The spatial distribution of vegetation index was consistent with the spatial distribution characteristics of soil wind erosion modulus， and its influence was deep.

Key words： REWQ model; wind erosion; sand-fixing; Shiqu county

土壤侵蝕是土壤退化的一種表現(xiàn)形式。土壤退化的兩個主要因素分別為水蝕和風(fēng)蝕，二者導(dǎo)致的全球退化土地面積約占84%[1]，使土壤過度侵蝕成為目前世界上最重要的環(huán)境問題之一。土壤侵蝕過程可能非常緩慢，也可能以較快的速度造成表土的嚴重流失。隨著中國基礎(chǔ)建設(shè)腳步的加快，人類活動導(dǎo)致土壤侵蝕速率也開始加快。土壤侵蝕使?fàn)I養(yǎng)豐富的上層土壤流失，將最終導(dǎo)致土地荒漠化。農(nóng)業(yè)泥沙流入河流還會導(dǎo)致水體富營養(yǎng)化，進一步污染水體和土壤。在土壤風(fēng)蝕過程中產(chǎn)生的懸浮土壤顆粒是空氣污染的主要來源之一，當(dāng)土壤顆粒被農(nóng)藥、化工燃料等有毒化學(xué)品污染，再著陸或被人體吸入時會造成嚴重的生態(tài)和公共健康問題。中國西部山區(qū)生態(tài)脆弱，水土流失現(xiàn)象較為嚴重，威脅著國家的生態(tài)安全，當(dāng)?shù)仫L(fēng)蝕模數(shù)的估算工作對維護其生態(tài)環(huán)境可持續(xù)發(fā)展具有一定意義。

土壤侵蝕模型的構(gòu)建能夠更好地探尋侵蝕的原因，以便對土壤侵蝕進行一系列預(yù)測工作，減輕其對生態(tài)環(huán)境的影響。由于侵蝕過程較為復(fù)雜，因此建模需要充分考慮各項因素，例如氣象、水文、地質(zhì)環(huán)境、土壤條件等。對于風(fēng)蝕模型研究，Woodruff等[2]提出了風(fēng)蝕方程（Wind erosion equation，WEQ），其被美國農(nóng)業(yè)部用于對農(nóng)田管理措施土壤風(fēng)蝕的潛在影像評估。在充分考慮土壤結(jié)皮及土壤可蝕性因素后，F(xiàn)ryrear等[3]在WEQ基礎(chǔ)上得到修正風(fēng)蝕方程（Revised wind erosion equation，RWEQ）。國內(nèi)學(xué)者也利用RWEQ模型進行了諸多研究，并取得了較好的成果，如：遲文峰等[4]基于RWEQ模型對內(nèi)蒙古高原進行風(fēng)蝕研究，反演了20世紀90年代以來內(nèi)蒙古高原土壤風(fēng)蝕模數(shù);江凌等[5]研究發(fā)現(xiàn)RWEQ模型估算結(jié)果與137Cs推算結(jié)果具有較高擬合性，在國內(nèi)區(qū)域尺度應(yīng)用具有一定可行性;鞏國麗等[6]針對模型中的土壤結(jié)皮因子和可蝕性因子進行了改進和應(yīng)用。這些研究結(jié)果表明，在理想數(shù)據(jù)情況下，RWEQ模型能很好完成對土壤風(fēng)蝕情況的預(yù)測，為防止土壤荒漠化提供決策依據(jù)。

本研究根據(jù)RWEQ方程構(gòu)建四川省石渠縣地理空間數(shù)據(jù)庫，對模型各個因子進行計算，進而評估石渠縣土壤風(fēng)蝕模數(shù)。根據(jù)侵蝕強度對結(jié)果進行分級，并分析其時間和空間分布特征，探尋土壤風(fēng)蝕影響因素，為當(dāng)?shù)厣鷳B(tài)保護和防風(fēng)固沙工作提供決策依據(jù)和經(jīng)驗指導(dǎo)。

1? 數(shù)據(jù)來源與研究方法

1.1? 研究區(qū)概況

石渠縣隸屬于四川省甘孜藏族自治州，位于四川省西北部，青藏高原東南緣，介于北緯32°19′—34°20′，東經(jīng)97°20′—99°16′，是四川省最偏遠、交通最不方便的縣之一。石渠縣幅員面積25 191 km2，全縣平均海拔4 000 m，為大陸性季風(fēng)高原氣候，日照時間長，氣溫偏低，晝夜溫差大，無絕對無霜期，暴風(fēng)雪災(zāi)害易發(fā)。該地區(qū)地質(zhì)結(jié)構(gòu)復(fù)雜，又因大規(guī)模的資源開采導(dǎo)致當(dāng)?shù)厣鷳B(tài)環(huán)境更加脆弱，崩塌、泥石流、滑坡等地質(zhì)災(zāi)害頻發(fā)。

1.2? 數(shù)據(jù)來源

構(gòu)建RWEQ模型所使用的數(shù)據(jù)主要分為氣候數(shù)據(jù)、太陽輻射數(shù)據(jù)、雪蓋數(shù)據(jù)、地形數(shù)據(jù)、遙感數(shù)據(jù)和土壤數(shù)據(jù)集。其中，氣候數(shù)據(jù)和太陽輻射數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)（http：//data.cma.cn/）;雪蓋數(shù)據(jù)和土壤數(shù)據(jù)集來源于寒區(qū)旱區(qū)科學(xué)數(shù)據(jù)中心（http：//westdc.westgis.ac.cn/data/）;30 m分辨率DEM數(shù)據(jù)來源于地理空間數(shù)據(jù)云（http：//www.gscloud.cn/），被用于計算地形起伏度;遙感數(shù)據(jù)來源于美國航空航天局（https：//ladsweb.modaps.eosdis.nasa.gov/）。

1.3? 研究方法與數(shù)據(jù)處理

1.3.1? RWEQ模型? 生態(tài)系統(tǒng)中防風(fēng)固沙服務(wù)功能是指在地表無植被狀況下的土壤風(fēng)蝕量與植被覆蓋條件下的土壤風(fēng)蝕量的差值[6，7]。在眾多模型中，RWEQ模型充分考慮了氣候條件因素、土壤可蝕性及結(jié)皮因素、地表粗糙度、植被覆蓋狀況等，其因子覆蓋全面，數(shù)據(jù)易于獲取，常用于在小于1年的時間尺度上進行土壤風(fēng)蝕量預(yù)報，適用區(qū)域范圍更廣，時間尺度更靈活[8]。

本研究采用RWEQ修正風(fēng)蝕方程對石渠縣的土壤風(fēng)蝕狀況進行定量評價，該方程基本表達式如下：

SR=SL潛-SL? ?（1）

SL潛=■×Qmax潛×e■

Qmax潛=109.8×（WF×EF×SCF×K′）S潛=150.71×（WF×EF×SCF×K′）-0.371 1? （2）

SL=■×Qmax×e■

Qmax=109.8×（WF×EF×SCF×K′×C）S=150.71×（WF×EF×SCF×K′×C）-0.371 1 （3）

式中，SR是土壤風(fēng)蝕模數(shù)[t/（hm2·年）];SL潛是潛在土壤風(fēng)蝕模數(shù)[t/（hm2·年）];SL是實際風(fēng)蝕模數(shù)[t/（hm2·年）];Z為下風(fēng)向最大風(fēng)蝕出現(xiàn)距離，Z=50 m;S潛、S是關(guān)鍵地塊的長度（m）;Qmax潛、Qmax為最大轉(zhuǎn)移量（kg/m）;WF為氣候因子（kg/m）;EF為土壤可蝕性因子（無量綱）;SCF為土壤結(jié)皮因子（無量綱）;K′為地表粗糙度因子;C為植被覆蓋度因子。

1.3.2? 評價指標(biāo)的獲取

1）氣象因子。風(fēng)速、氣溫、降雨、日照輻射以及雪蓋天數(shù)等氣候狀況都會影響土壤風(fēng)蝕模數(shù)，其共同組成了氣候因子，該因子計算表達式：

WF=WE×■×SW×SD? （4）

WE=μ2×（μ2-μ1）2×Nd? （5）

？籽=348.01（■）（6）

SW=■? ?（7）

ETP=0.016 2×■×（DT+17.8）? （8）

SD=1-P? ?（Hsnow>25.4 mm）? （9）

式中，WF為氣象因子（kg/m）;WE為風(fēng)場強度因子（m3/s3），由監(jiān)測風(fēng)速μ2（m/s）、起沙風(fēng)速μ1（假定為5 m/s）和觀察周期天數(shù)Nd計算得到;？籽為空氣密度（kg/m3），由海拔高度EL（km）和絕對溫度T（K）計算得到;g為重力加速度（m/s2）;SW為土壤濕度因子（無量綱）;R為降雨量（mm）;I為灌溉量（mm）;Rd為降雨次數(shù)和（或）灌溉天數(shù);ETP為地表潛在相對蒸發(fā)量（mm），由太陽輻射SR（cal/cm2）和平均溫度DT（℃）計算得到;SD為雪覆蓋因子（無量綱）;P為計算時段內(nèi)積雪覆蓋深度（Hsnow）大于25.4 mm的概率。

2）土壤可蝕性因子。土壤可蝕性是指土壤對侵蝕的敏感性。對于不同的機械組成和理化性質(zhì)的土壤種類，粒度越小，有機質(zhì)含量越低，其土壤可蝕性越大，越容易被侵蝕;反之粒度越粗，有機質(zhì)含量越高，其可蝕性越小，越不容易被侵蝕。對于土壤可蝕性因子的計算公式如下[9]：

EF=■

（10）

土壤結(jié)皮是指某些低等生物與土表相互作用或降水滴濺在土表上形成的微層，一般按產(chǎn)生機理可分為生物性結(jié)皮和物理性結(jié)皮。其中，生物性結(jié)皮有利于抵抗土壤風(fēng)蝕;物理性結(jié)皮易碎，反而加速了土壤被風(fēng)蝕的過程。其計算公式如下[10]：

SCF=■? （11）

在兩式中，所有參量均為百分制，sa為土壤沙粒含量;si為土壤粉沙含量;cl為黏土含量;OM為有機質(zhì)含量;Caco為碳酸鈣含量，可不予考慮。

3）植被覆蓋度因子。不同的植被有著不同的根系，也就有不同的固水固沙能力。植被覆蓋度因子表示一定植被覆蓋條件下，對土壤風(fēng)蝕的抑制作用大小。依據(jù)研究區(qū)LUCC分類圖，將植被分為林地、灌木、草地、農(nóng)田、裸地5個植被類型，根據(jù)不同的系數(shù)計算各植被覆蓋度因子。其表達式[11]：

C=e■? （12）

式中，ai為不同植被類型的系數(shù)，其中，林地為0.153 5，灌木為0.092 1，草地為0.151 1，農(nóng)田為0.043 8，裸地為0.076 8;SC為植被覆蓋度（無量綱），由NDVI數(shù)據(jù)集計算得到。

4）地表粗糙度因子。地表粗糙度是指由地形引起的土地表面粗糙程度對土壤風(fēng)蝕的影響。其表達式[11]：

K′=e■? ?（13）

Kr=0.2×■? ?（14）

式中，Kr為因地形起伏產(chǎn)生的地形粗糙度長度（cm）;Crr為隨機糙度因子，取0;ΔH為距離L范圍內(nèi)的海拔高程差（m），根據(jù)不同的起伏地形情況，L具有不同的值。

2? 結(jié)果與分析

為了分析該地區(qū)的土壤風(fēng)蝕變化，參照水利部《土壤侵蝕分類分級標(biāo)準(zhǔn)》（SL 190-2007）并結(jié)合研究區(qū)的實際情況采用自然斷點法將土壤風(fēng)蝕模數(shù)分為潛在風(fēng)蝕（0，5.45]、微度風(fēng)蝕（5.45，10.65]、輕度風(fēng)蝕（10.65，15.62]、中度風(fēng)蝕（15.62，24.50]、重度風(fēng)蝕（24.50，55.24]5類。

2.1? 土壤風(fēng)蝕等級空間分布特征

根據(jù)上述方法得到2015年石渠縣土壤風(fēng)蝕等級空間分布圖（圖1）。整體上看，石渠縣土壤風(fēng)蝕情況不容樂觀，南部較北部侵蝕狀況更為嚴重。中部呷衣鄉(xiāng)以及西南部真達鄉(xiāng)、奔達鄉(xiāng)、正科鄉(xiāng)、鄧柯鄉(xiāng)和麻呷鄉(xiāng)風(fēng)蝕程度最為嚴重;東南部大部分地區(qū)如蝦扎鄉(xiāng)、起塢鄉(xiāng)、長沙干馬鄉(xiāng)、長須貢馬鄉(xiāng)和瓦須鄉(xiāng)存在輕度風(fēng)蝕;以北部集中的其余地區(qū)風(fēng)蝕狀況得到有效控制，均為微度或潛在侵蝕。

結(jié)合石渠縣當(dāng)?shù)財?shù)字高程模型（DEM）等地形資料，發(fā)現(xiàn)石渠縣整體海拔較高，在西南部以及中部呷衣鄉(xiāng)存在較大的高程差，地勢起伏大，屬于高原山地，風(fēng)蝕狀況更為嚴重;而北部雖然海拔較高，但高程差較小，屬于高原地形，風(fēng)蝕狀況相較而言較輕。故分析認為地表粗糙度因子對土壤風(fēng)蝕存在較大的影響，地表粗糙度因子越大，土壤風(fēng)蝕模數(shù)越高。

2.2? 氣象變化對土壤風(fēng)蝕模數(shù)的影響

利用2015年全年各月的月均氣象數(shù)據(jù)，結(jié)合石渠縣2015年植被覆蓋、土地利用、土壤機械組成及理化性質(zhì)等數(shù)據(jù)，計算該時間尺度內(nèi)的月均土壤風(fēng)蝕模數(shù)，以探究氣象變化對土壤風(fēng)蝕模數(shù)的影響。

根據(jù)上述方法得到2015年石渠縣各月土壤風(fēng)蝕模數(shù)統(tǒng)計圖（圖2）。由圖2可知，1—3月和10—11月處于雪期，土壤風(fēng)蝕模數(shù)較小。從3月之后迅速上升，在6月達到最高峰，期間處于春夏季節(jié)交替，風(fēng)速較大，風(fēng)力強勁，6月進入汛期，在風(fēng)力和降雨的共同作用下，土壤風(fēng)蝕模數(shù)也達到了全年最大。9月同樣對應(yīng)一個小高峰，該時期當(dāng)?shù)丶磳⑦M入雪期，植被枯萎，風(fēng)力強勁，致使土壤侵蝕模數(shù)在此攀升。由此可知，氣象對土壤侵蝕模數(shù)有較大影響，且因子復(fù)雜，不同的時期有不同的主力因子，視當(dāng)?shù)貧庀髼l件決定。

2.3? 植被覆蓋變化對土壤風(fēng)蝕模數(shù)的影響

植被覆蓋能有效地抑制風(fēng)蝕速度。作為風(fēng)蝕的阻抗因子，不同的植被類型具有不同強弱的固沙能力，同一種植被不同程度的覆蓋度也會造成不同強弱的固沙能力。將有植被覆蓋的地表固沙量（Qsa）與無植被覆蓋（如裸地）的地表風(fēng)蝕量（Qsoil）相比，就能消除氣候等因素的影響，得到防風(fēng)固沙服務(wù)功能的保有率（F）。

其計算公式如下：

F=■? （15）

保有率反映了植被變化對防風(fēng)固沙服務(wù)能力的影響，其值越高，土壤風(fēng)蝕模數(shù)越小，固沙能力越強;反之值越低，土壤風(fēng)蝕模數(shù)越大，固沙能力越弱。

根據(jù)上述方法得到2015年石渠縣防風(fēng)固沙保有率空間分布圖（圖3）。石渠縣保有率范圍較大，最小為0.56，最高為4.92。整體依然是北部高于南部，由北到南逐漸減少，與土壤風(fēng)蝕模數(shù)空間分布特征相一致，說明植被覆蓋程度對固沙量有較大影響，且植被覆蓋度因子較低的區(qū)域固沙量更大，土壤風(fēng)蝕程度較微弱，植被覆蓋度因子高的區(qū)域固沙量偏小，土壤風(fēng)蝕程度較為嚴重。

3? 小結(jié)

本研究利用修正土壤風(fēng)蝕模型來定量分析土壤風(fēng)蝕模數(shù)，評價石渠縣風(fēng)蝕等級，并對氣象因子和植被因子分別與土壤風(fēng)蝕模數(shù)之間的關(guān)系進行了深入的分析與探究。具體結(jié)論如下：

1）石渠縣的各區(qū)域抗土壤風(fēng)蝕能力有差異，北部整體優(yōu)于南部，北部較為集中的以潛在和微度風(fēng)蝕為主，東南部的蝦扎鄉(xiāng)、起塢鄉(xiāng)、長沙干馬鄉(xiāng)、長須貢馬鄉(xiāng)和瓦須鄉(xiāng)存在輕度風(fēng)蝕，西南部的真達鄉(xiāng)、奔達鄉(xiāng)、正科鄉(xiāng)、鄧柯鄉(xiāng)和麻呷鄉(xiāng)等則以中度、重度風(fēng)蝕為主。

2）土壤風(fēng)蝕模數(shù)隨著一年中時間的推進，氣象的變化而在曲線中表現(xiàn)出兩個小高峰，最高峰由降雨和風(fēng)力共同導(dǎo)致，次高峰由植被狀況和風(fēng)力共同導(dǎo)致。氣象因素在不同的時期有不同的主力因子，視當(dāng)?shù)貧庀髼l件具體分析。

3）防風(fēng)固沙保有率能較好地定量描述由植被引起的土壤風(fēng)蝕量及其空間分布特征。植被對防風(fēng)固沙服務(wù)影響程度較深，整體上保有率良好，其空間分布特征契合土壤風(fēng)蝕模數(shù)空間分布特征。

參考文獻：

[1] BLANCO-CANQUI H， LAL R. Principles of soil conservation and management[M/OL].Springer，2010：http：//www.doc88.com/p-31975794580.html，2008-08/2019-02-18.

[2] WOODRUFF N P，SIDDOWAY F H. A wind erosion equation[J].Proceedings of the soil science society of america，1965，29（5）：602-608.