林建琳

(本溪市水利電力勘測設計有限責任公司，遼寧 本溪 117022)

土壤侵蝕對我國自然環(huán)境和社會發(fā)展造成惡劣影響。一方面，水土流失使我國喪失了大面積的水生態(tài)鄰近區(qū)域良好質(zhì)量土地，優(yōu)質(zhì)耕地面積與林地減少，農(nóng)業(yè)經(jīng)濟大受影響[1- 2]；另一方面，不論是土壤水力侵蝕或風力侵蝕均給區(qū)域生態(tài)環(huán)境帶來惡劣的影響[3- 4]。

現(xiàn)有研究大多是關于水力作用下土壤的受侵蝕狀況、因素、影響以及防治措施等[5- 7]，而對于風力影響下土壤的侵蝕研究較為不足。梁會等[8]通過野外土壤采樣及風洞試驗修正模型參數(shù)，運用遙感、地理信息系統(tǒng)技術，對2000—2015年土壤風蝕強度時空變化特征進行分析，運用變異系數(shù)法確定各影響因子權重，明確了區(qū)域影響風蝕的主要因子為植被覆蓋和土壤濕度；韓柳等[9]利用1971—2011年北方風蝕區(qū)146個氣象站的逐時風速資料，將5m/s以上風速劃分為和風、強風、大風、狂風4個等級，采用線性趨勢分析、小波分析、M-K突變檢驗分析和空間相似度分析等方法分析了該地區(qū)不同等級土壤風蝕力時空分布特征，并討論了與平均風速的時空相關性；張安東[10]等以寧夏中部干旱帶不同自然地貌作為試驗樣地，應用誘捕法收集風蝕沙粒，通過激光粒度分析儀分析沙粒粒徑后，判斷各試驗點的土壤組成，分析評價沙粒粒徑與地表風蝕相關性，發(fā)現(xiàn)風蝕沙源多外來沉積，風蝕強度降低;耕作農(nóng)田中黃綿土粉粒比例一般高于風沙土，抗風蝕能力較強。

我國北方地區(qū)存在大范圍的風蝕區(qū)域，尤其以覆沙黃土水力侵蝕流失最為嚴重，因此，本文將以某區(qū)域為例，基于多年監(jiān)測數(shù)據(jù)建立地形因子數(shù)據(jù)庫，并深入分析區(qū)域風蝕主要影響因素，進一步建立區(qū)域風蝕量預測模型，以期能夠為區(qū)域風蝕防治措施的開展提供一定的依據(jù)。

1 流域現(xiàn)狀

北方某流域土地利用類型以草地、耕地、林地為主，屬于多沙流沙區(qū)，是典型的風蝕-水蝕復合型侵蝕區(qū)域。據(jù)相關數(shù)據(jù)統(tǒng)計，流域多年平均降雨量為330mm，且主要集中在河流汛期7—10月；流域每年3—5月有大風天氣，最大風速在25m/s左右。受風蝕-水蝕綜合影響下，流域總水土流失面積達到389.26km2，土壤侵蝕情況嚴重。圖1為基于ArcGIS的區(qū)域衛(wèi)星遙感圖，由圖1可見，區(qū)域北部以農(nóng)田種植區(qū)域為主，綠化程度高，水土保持情況良好，片區(qū)土壤輕度侵蝕；相對之下，流域中部為大面積的風蝕區(qū)域，南部支流較多，形成了廣闊的水蝕區(qū)域。

圖1 區(qū)域土壤侵蝕分布圖

圖2 區(qū)域土壤累計風蝕量變化趨勢

根據(jù)上述分析可知，該流域主要土壤侵蝕類型為風蝕，因此本文將著重于區(qū)域強力風蝕條件下土壤流失情況，根據(jù)相關數(shù)據(jù)顯示，區(qū)域在一年內(nèi)土壤累計風蝕量情況如圖2所示。區(qū)域風蝕主要發(fā)生在3—5月強風期，因此累計曲線自此階段增速較快。

2 基于地形因子的流域風蝕模型

2.1 地形因子

選定某一片區(qū)單元進行研究，選取不同樁點的監(jiān)測數(shù)據(jù)，運用tecplot軟件對點數(shù)據(jù)進行光滑拼接，形成圖3所示該區(qū)域沙丘形態(tài)變化趨勢3D云圖，由圖3可見，在風力侵蝕的作用下，該沙丘逐漸由北向南推移，坡度、表面粗糙度、起伏度等均發(fā)生了較大的變化。

圖3 區(qū)域沙丘形態(tài)變化趨勢圖

基于流域2015—2016年監(jiān)測數(shù)據(jù)建立起流域風蝕情況變遷數(shù)據(jù)庫，經(jīng)過篩選分析，結合流域?qū)嶋H情況，得出流域風蝕因素主要地形因子包括坡度、坡度變率、坡向變率、粗糙度以及起伏度。

2.2 模型建立

根據(jù)這5種地形因子建立起數(shù)據(jù)庫，得出主要參數(shù)見表1。

通過以上數(shù)據(jù)建立起各因素與風蝕量的一元回歸模型(如圖4所示)，發(fā)現(xiàn)區(qū)域風蝕量與坡度、坡變化率、坡向變率、粗糙度以及起伏度都具有一定的相關性，但相對于坡度、坡變化率以及坡向變率，粗糙度與起伏度相關性更高，R2分別達到0.605與0.749。由此可見，區(qū)域風蝕量與粗糙度及起伏度關系更大，因此，文章將主要圍繞粗糙度與起伏度兩個地形因子指標建立起區(qū)域土壤風蝕模型。

根據(jù)區(qū)域監(jiān)測體系所提供的數(shù)篩選后建立起區(qū)域土壤風蝕數(shù)據(jù)庫，基于多元逐步回歸模型基本理論建立區(qū)域土壤風蝕多因素模型，為突出多元模型的優(yōu)勢，并與所建立的一元回歸模型進行對比，以說明多元模型更高的科學性與準確性。模型的基本公式參數(shù)如下：

(1)模型A：影響因子為起伏度l，R=0.865，R2=0.749，修正R2=0.736，標準估計誤差m=1.623；

(2)模型B：影響因子為粗糙度k與起伏度l，R=0.922，R2=0.850，修正后R2=0.834，標準估計誤差m=1.288。

利用SPSS軟件進一步求取模型基本參數(shù)的平方差、自由度、均方以及F的回歸、殘差和總計，經(jīng)過參數(shù)的輸入、篩選及比對，最終輸出結果見表2。

根據(jù)上述分析，得出最終風蝕量與地形因子的A、B模型基本函數(shù)關系式如下：

模型A：Y=43.345k-0.468

(1)

模型B：Y=40.807k+6.453l+1.25

(2)

為精確表述地形因子與風蝕量模型的正確性，對排除的變量所帶來的紊亂性與不穩(wěn)定也進行分析，計算結果見表3。

由以上不穩(wěn)定性中模型A計算結果可知，對于模型A，坡度m、坡度變率ma以及坡向變率mb都有較大的共線不穩(wěn)定性，共線紊亂系數(shù)分別達到了0.720、0.766和0.975，由此可見，坡度m、坡度變率ma以及坡向變率mb三個地形因子無法加入模型A中建立更為完善的多元分析模型；相較于上述3個地形因子，粗糙度K在模型A中帶來的紊亂系數(shù)則較小，僅達到0.313，顯然可以根據(jù)模型A基本算式帶入粗糙度k建立更為復雜的多元模型，即所述模型B。

表1 不同地形因子與風蝕量變化表

圖4 風蝕量與各地形因子相關性

表2 方差分析

表3 排除變量共線不穩(wěn)定性分析

對于模型B，坡度m、坡度變率ma以及坡向變率mb三個地形因子的共線紊亂系數(shù)均在0.5以上，分別達到0.698、0.684和0.881，顯然不能進一步帶入模型B中建立起更多因素的復雜模型，因此我們可以認為，基于粗糙度k與起伏度l建立起的區(qū)域風蝕為最優(yōu)解。

3 結論

通過篩選區(qū)域多年地表微地形因子監(jiān)測數(shù)據(jù)建立起區(qū)域土壤風蝕影響因子數(shù)據(jù)庫，并利用tecplot等數(shù)值軟件建立地表覆沙黃土走勢傾向圖，發(fā)現(xiàn)區(qū)域土壤在風力影響下整體由西北向東南移動。進一步對地形監(jiān)測數(shù)據(jù)進行發(fā)掘，建立起區(qū)域土壤風蝕量與多微地形因子之間的經(jīng)驗模型，能夠較好的模擬預測區(qū)域風蝕量，能夠為區(qū)域土壤侵蝕防治提供一定的依據(jù)。