王英敏

(本溪市水利工程質(zhì)量與安全監(jiān)督站，遼寧 本溪 117022)

隨著人類活動(dòng)的增加，全球氣候在近百年來(lái)發(fā)生了顯著的變化，并造成極端氣候事件的增加，進(jìn)而對(duì)全球的水資源、土壤侵蝕以及生態(tài)系統(tǒng)產(chǎn)生了潛在影響[1]。遼寧省西部地區(qū)是遼寧省乃至我國(guó)水土流失嚴(yán)重的地區(qū)之一，而氣候變化也對(duì)當(dāng)?shù)氐耐寥狼治g特征造成顯著影響[2]?；诖耍疚囊赃|西地區(qū)朝陽(yáng)市為例，研究氣候變化背景下土壤侵蝕力在時(shí)間序列上的變化規(guī)律，為當(dāng)?shù)氐乃帘３忠?guī)劃和綜合防治工作提供參考。

1 研究區(qū)概況

朝陽(yáng)市位于我國(guó)遼寧省西部，南鄰河北省，北接內(nèi)蒙古自治區(qū)，東連遼寧省中部的工業(yè)地帶，具有十分優(yōu)越的地理位置。朝陽(yáng)市下轄2個(gè)市轄區(qū)、2個(gè)縣和1個(gè)自治縣，總面積19736km2，人口334.9萬(wàn)人。朝陽(yáng)市為典型的溫帶大陸性氣候，四季分明、雨熱同期，降水偏少，年均降雨量450～580mm。朝陽(yáng)市境內(nèi)的主要河流有大凌河、小凌河、青龍河和老哈河等，多年平均徑流量13.22億m3。由于朝陽(yáng)市干旱少雨，大多數(shù)地區(qū)為丘陵地形，植被覆蓋率較低，大部分土地的地表土壤不能有效截蓄雨水、分散徑流、降低流速，加上局部的干旱和暴雨頻發(fā)，因此水土流失主要以水蝕為主且比較嚴(yán)重[3]。改革開(kāi)放以來(lái)，隨著經(jīng)濟(jì)發(fā)展水平的不斷提升，朝陽(yáng)市加大了在環(huán)境治理領(lǐng)域的投入力度，水土流失治理取得了顯著成績(jī)[4]。經(jīng)過(guò)多年治理，朝陽(yáng)市的水土流失面積減少了1.22萬(wàn)km2，但是仍有水土流失面積6553km2，治理任務(wù)仍舊十分艱巨。因此，研究當(dāng)?shù)亟涤昵治g力的分布特征，對(duì)采取正確的水土保持策略，提高水土流失防治效果具有重要的借鑒意義。

2 降雨侵蝕力的計(jì)算模型

2.1 算法的選擇

研究中所使用的氣象數(shù)據(jù)來(lái)自于遼寧省氣象信息中心，從其數(shù)據(jù)庫(kù)中提取朝陽(yáng)市13個(gè)主要?dú)庀笳军c(diǎn)1958—2019年的降雨和地理信息數(shù)據(jù)，作為本次研究的基礎(chǔ)氣象數(shù)據(jù)。通過(guò)查閱和整理該領(lǐng)域的研究文獻(xiàn)，結(jié)合遼寧省西部地區(qū)的實(shí)際特點(diǎn)，選擇其中比較有代表性的朱強(qiáng)模型[5]、陳世發(fā)模型[6]、宮冰模型[7]以及劉慧英模型[8]，分別記作模型1、模型2、模型3和模型4，以朝陽(yáng)市太平房和中三家兩地的降雨侵蝕力精確值為基準(zhǔn)，利用有效系數(shù)和相對(duì)偏差分析方法，展開(kāi)上述4個(gè)模型之間的差異分析，并對(duì)各個(gè)模型的優(yōu)劣進(jìn)行必要的評(píng)價(jià)。根據(jù)四個(gè)模型的具體計(jì)算要求，分別計(jì)算出太平房在1997年至2003年以及中三家在2001年至2005年的降雨侵蝕力變化特征，結(jié)果分別如圖1、2所示。由圖可知，4個(gè)模型在降雨侵蝕量計(jì)算結(jié)果方面的變幅較大。從三家地的結(jié)果來(lái)看，4個(gè)模型的R值的變化趨勢(shì)與基準(zhǔn)值的變化趨勢(shì)基本一致，除了模型1在1998年的數(shù)值略大于基準(zhǔn)值之外，其他模型和年份的數(shù)值均小于基準(zhǔn)值，且基本呈現(xiàn)出基準(zhǔn)值>模型1>模型2>模型3>模型4的基本特征。從圖2來(lái)看，4個(gè)模型的R值和基準(zhǔn)值之間存在比較明顯的交叉現(xiàn)象。

圖1 太平房4種模型計(jì)算成果對(duì)比

圖2 三家地4種模型計(jì)算成果對(duì)比

模型的優(yōu)劣主要通過(guò)有效系數(shù)和相對(duì)偏差兩個(gè)指標(biāo)進(jìn)行衡量，文中選擇的4個(gè)模型的計(jì)算成果見(jiàn)表1。由表中的計(jì)算結(jié)果來(lái)看，太平房和三家地在模型1的計(jì)算結(jié)果上均有較高的有效系數(shù)以及較低的偏差系數(shù)，有效系數(shù)分別為0.86和0.76，偏差系數(shù)分別為0.05和0.10。從太平房的計(jì)算結(jié)果來(lái)看，模型4的有效系數(shù)最小，為0.63，R值也相對(duì)偏小，但是在三家地的有效系數(shù)最高。總體來(lái)看，模型4的穩(wěn)定性明顯不足。此外，模型2和模型3的模擬能力與模型1相比明顯偏低。綜合上述分析，模型1的整體模擬效果相對(duì)較好，因此在后續(xù)研究中選擇模型1進(jìn)行研究區(qū)的降雨侵蝕力計(jì)算。

表2 各模型R值、有效系數(shù)和相對(duì)偏差計(jì)算成果

2.2 計(jì)算模型的構(gòu)建

根據(jù)上文篩選的月降雨侵蝕力計(jì)算方法，計(jì)算研究時(shí)段內(nèi)的月降雨侵蝕力，并匯總年降雨侵蝕力。相關(guān)性分析結(jié)果顯示，年降雨侵蝕力和降雨量之間存在較強(qiáng)的相關(guān)性，因此將逐年降雨侵蝕力的R值作為y值，逐年降雨量作為x值，對(duì)所有的站點(diǎn)利用線性函數(shù)、對(duì)數(shù)函數(shù)、冪函數(shù)和指數(shù)函數(shù)等四種模型進(jìn)行回歸分析[9]，從計(jì)算結(jié)果來(lái)看，冪函數(shù)模型的決定系數(shù)最大，為0.930，所占比例為97.3%，也最大，因此回歸擬合效果最佳。因此，利用冪函數(shù)進(jìn)行回歸分析，獲得整個(gè)朝陽(yáng)市的回歸擬合方程，其表達(dá)式為：

y=0.046x1.5887

(1)

根據(jù)上式，可以獲得朝陽(yáng)市降雨侵蝕力的R值計(jì)算公式：

R=0.046P1.5887

(2)

式中，R—年降雨侵蝕力，MJ·hm-2h-1a-1；P—逐年降雨量，mm。

2.3 模型的檢驗(yàn)

利用公式(2)對(duì)朝陽(yáng)市另外8個(gè)站點(diǎn)的逐年降雨侵蝕力進(jìn)行計(jì)算，并對(duì)其進(jìn)行線性回歸驗(yàn)證，結(jié)果顯示，8個(gè)氣象站的回歸決定系數(shù)為0.908～0.963，均值為0.937；與基準(zhǔn)值的線性回歸斜率為0.934～1.099，均值為1.000，其中斜率在0.9以上的占100%，斜率在0.95至1.05之間的占70.56%。由此可見(jiàn)，上述模型在計(jì)算降雨侵蝕力方面精度較高，計(jì)算方法較為簡(jiǎn)便，且受降雨資料的限制較少，可以用于后續(xù)研究。

3 計(jì)算結(jié)果與分析

3.1 年內(nèi)分布特征分析

根據(jù)上文構(gòu)建的模型，對(duì)朝陽(yáng)市各月的降雨侵蝕力R值的多年均值以及降雨侵蝕力的比例進(jìn)行計(jì)算，根據(jù)計(jì)算結(jié)果繪制出降雨侵蝕力年內(nèi)分布特征如圖3所示。由圖3可知，朝陽(yáng)市的降雨侵蝕力R值與降雨的季節(jié)之間存在較為一致的變化特征，呈現(xiàn)出明顯的單峰型特點(diǎn)，在1—7月逐漸增加，在8—12月逐漸減小。從降雨侵蝕比例上來(lái)看，主要集中于降雨量較大的4—10月，占到年R值的97.98%，特別是汛期的7月份的降雨侵蝕力占到了全年的36.98%。

3.2 年際分布特征分析

根據(jù)上文構(gòu)建的模型，對(duì)朝陽(yáng)市1984—2019年間的降雨侵蝕力R值進(jìn)行計(jì)算，結(jié)果見(jiàn)表2。由表2中的結(jié)果可知，朝陽(yáng)市多年平均降雨侵蝕力為1898.28 MJ·hm-2h-1a-1，其最大值出現(xiàn)在1998年，為5448.39 MJ·hm-2h-1a-1，最小值出現(xiàn)在1986年，為990.47MJ·hm-2h-1a-1，振幅為4457.92 MJ·hm-2h-1a-1。由此可見(jiàn)，朝陽(yáng)市年降雨侵蝕力的差別較大。

圖3 朝陽(yáng)市降雨侵蝕力年內(nèi)分布特征

表2 年降雨侵蝕力計(jì)算成果 單位：MJ·hm-2h-1a-1

4 結(jié)論

此次研究以遼西地區(qū)朝陽(yáng)市為例，構(gòu)建降雨侵蝕力計(jì)算模型，對(duì)降雨侵蝕力時(shí)間序列分布特征進(jìn)行研究，并獲得如下主要結(jié)論：

(1)借鑒朱強(qiáng)的研究成果，利用冪函數(shù)擬合，構(gòu)建起朝陽(yáng)市降雨侵蝕力計(jì)算模型，驗(yàn)證結(jié)果顯示模型具有較高的計(jì)算精度。

(2)朝陽(yáng)市的降雨侵蝕力R值與季節(jié)之間存在較為一致的變化特征，呈現(xiàn)出明顯的單峰型特點(diǎn)；從降雨侵蝕比例上來(lái)看，主要集中于降雨量較大的4—10月，占到年R值的97.98%。

(3)朝陽(yáng)市年降雨侵蝕力的差別較大。