三峽庫區(qū)香溪河流域降雨侵蝕力的時(shí)空分布特征

劉惠英，任洪玉，張長(zhǎng)偉，張平倉?

(1.南昌工程學(xué)院，330099，南昌； 2.長(zhǎng)江科學(xué)院水土保持研究所，430012，武漢)

三峽庫區(qū)香溪河流域降雨侵蝕力的時(shí)空分布特征

劉惠英1，任洪玉2，張長(zhǎng)偉2，張平倉2?

(1.南昌工程學(xué)院，330099，南昌； 2.長(zhǎng)江科學(xué)院水土保持研究所，430012，武漢)

為研究流域降雨侵蝕力變化規(guī)律，利用三峽庫區(qū)香溪河流域內(nèi)10個(gè)雨量站1971—2010年的日降雨資料，采用降雨侵蝕力日降雨簡(jiǎn)易模型，分析該流域降雨侵蝕力的年內(nèi)分配和年際變化規(guī)律。在ArcGIS軟件支持下，采用克里格插值研究流域降雨和降雨侵蝕力時(shí)空變化特征。結(jié)果表明：香溪河流域降雨侵蝕力多年變化范圍為2 465.26～7 419.29 MJ·mm/(hm2·h)，多年平均均值為4 535.63 MJ·mm/(hm2·h)，降雨侵蝕力R值的年際分配差異明顯，最大年R值為最小年R值的3倍；流域侵蝕力空間變化趨勢(shì)為從西向東逐漸遞減；流域近40多年的降雨和降雨侵蝕力系列比較平穩(wěn)，經(jīng)Mann-kendall檢驗(yàn)無顯著的變化趨勢(shì)。流域降雨量、侵蝕性降雨量和降雨侵蝕力年內(nèi)分布較集中，汛期降雨量、汛期侵蝕性降雨量、汛期降雨侵蝕力占全年的比例分別為85.4%、92.4%和94.0%。

降雨侵蝕力； 日降雨量； 時(shí)空分布； 香溪河流域

降雨是引起土壤侵蝕的動(dòng)力因素之一，雨滴擊濺造成土壤顆粒分離、搬運(yùn)最終導(dǎo)致土壤流失。降雨侵蝕力因子R是降雨能量的定量指標(biāo)，土壤侵蝕模型USLE、RUSLE及其他眾多流域農(nóng)業(yè)非點(diǎn)源污染模型均把R值作為主要參數(shù)之一，對(duì)于R的計(jì)算方法和空間分布規(guī)律的研究越來越被人們重視。正確評(píng)估降雨侵蝕力，對(duì)預(yù)測(cè)土壤侵蝕及優(yōu)化水土保持措施具有重要意義。降雨侵蝕力的計(jì)算方法由W. H. Wischmeier等[1]在1978年提出，定義為一場(chǎng)降雨的降雨總動(dòng)能E與最大30min降雨強(qiáng)度I30的乘積，該指標(biāo)已在世界范圍內(nèi)得到廣泛應(yīng)用。我國(guó)對(duì)降雨侵蝕力指標(biāo)的研究始于20世紀(jì)八十年代，王萬忠等[2]通過研究代表性小區(qū)降雨侵蝕資料，證實(shí)EI30在我國(guó)同樣適用；但實(shí)際工作中因較難獲得E和I30資料，許多學(xué)者[3-6]嘗試采用氣象站點(diǎn)常規(guī)降雨數(shù)據(jù)建立侵蝕力簡(jiǎn)易算法，構(gòu)建了多種經(jīng)驗(yàn)估算模型，常見的是利用氣象站整編各時(shí)間尺度的降雨量等常規(guī)資料來計(jì)算降雨侵蝕力。C. W. Richardson等[5]的冪函數(shù)結(jié)構(gòu)的日雨量侵蝕力簡(jiǎn)易模型，得到許多學(xué)者的進(jìn)一步分析驗(yàn)證。章文波等[6-7]對(duì)冪函數(shù)日降雨量侵蝕力模型做了修正并提出了適用我國(guó)的參數(shù)估算方法，并用其分析我國(guó)降雨侵蝕力及其季節(jié)分布。謝云等[8]利用陜北團(tuán)山溝徑流場(chǎng)數(shù)據(jù)定量判斷出我國(guó)的侵蝕性降雨標(biāo)準(zhǔn)。

三峽庫區(qū)是長(zhǎng)江上游四大重點(diǎn)水土流失片區(qū)之一，我國(guó)政府已將三峽庫區(qū)列為全國(guó)水土保持重點(diǎn)防治區(qū)。庫區(qū)土壤侵蝕產(chǎn)生的泥沙直接入庫，庫區(qū)水土流失治理對(duì)減緩庫區(qū)泥沙淤積，延長(zhǎng)水庫使用壽命具有重大的實(shí)際意義。庫區(qū)降雨侵蝕力研究不僅對(duì)該區(qū)水土流失治理、水土保持效益評(píng)價(jià)具有重要意義，而且對(duì)三峽水庫區(qū)間經(jīng)過“長(zhǎng)治”工程、“天保”工程和退耕還林(草)工程等治理后流域水沙關(guān)系變化原因分析提供一定的理論支持。從庫區(qū)目前研究現(xiàn)狀來看，至少存在采用數(shù)據(jù)周期過短和采用的測(cè)站數(shù)量過少2方面問題。W. H. Wischmeier等[1]指出，考慮到氣候變化的周期波動(dòng)，多年平均降雨侵蝕力的計(jì)算一般要求至少有20年以上的降雨過程資料；同時(shí)閆業(yè)超等[9]研究表明，降雨資料的時(shí)間序列長(zhǎng)度對(duì)R平均值的估計(jì)置信度有顯著影響。目前三峽水庫庫區(qū)侵蝕力研究在這2方面都不盡如人意。吳昌廣等[10]、汪言在等[11]和張照錄等[12]在研究整個(gè)三峽庫區(qū)侵蝕力時(shí)分別只用了有7、4和5個(gè)站的日降雨資料；花利忠等[13]研究庫區(qū)流域大寧河時(shí)也只用了8年的資料；史東梅等[14]、繆馳遠(yuǎn)等[15]和張革等[16]各分別在研究涪陵、重慶主城區(qū)和香溪河流域時(shí)也只用了1個(gè)站的資料。單個(gè)站點(diǎn)對(duì)于分析降雨侵蝕力時(shí)空分布，指導(dǎo)流域水土保持措施規(guī)劃遠(yuǎn)遠(yuǎn)不夠，對(duì)流域治理的指導(dǎo)意義有限。

筆者以三峽庫區(qū)支流香溪河流域?yàn)槔捎昧饔騼?nèi)10個(gè)雨量站點(diǎn)1971—2010年的日降雨量資料對(duì)流域降雨侵蝕力進(jìn)行系統(tǒng)研究，增加了R值分布密度，提高了其空間插值精度，并分析其時(shí)空分布規(guī)律，為定量評(píng)價(jià)流域水土流失提供了重要背景參數(shù)，同時(shí)對(duì)進(jìn)一步開展水土流失調(diào)查、土壤侵蝕評(píng)價(jià)、水土保持效益分析和開展生態(tài)恢復(fù)工作提供科學(xué)依據(jù)。

1 研究區(qū)概況

香溪河流域地處長(zhǎng)江流域上游(圖1)，發(fā)源于神農(nóng)架南麓，系長(zhǎng)江三峽大壩壩首的第一大支流，位于湖北省西北部，流域總面積3 183 km2，其中神農(nóng)架林區(qū)865 km2，興山縣2 106 km2，秭歸縣212 km2，為峽谷型河流。香溪河由北向南縱貫興山縣全境，至游家河入秭歸縣后又匯屈平河之水，于香溪鎮(zhèn)東注入長(zhǎng)江，匯長(zhǎng)江之水流入西陵峽，下泄至三峽大壩[17]。流域處于亞熱帶大陸性季風(fēng)氣候區(qū)。三峽工程蓄水后，形成了典型的香溪河庫灣。流域土地資源短缺，因多年盲目墾荒、毀林造田、亂砍濫伐，20世紀(jì)80年代末森林植被破壞嚴(yán)重，土壤、巖石裸露加劇，造成水土流失，大量泥沙入河。遙感影像顯示，1988年流域水土流失面積2 244.72 km2，占流域總面積的70.5%，其中強(qiáng)烈及以上水土流失程度面積占流域面積的25.97%，占水土流失總面積的36.82%[18]。

圖1 香溪河流域地理位置及雨量站分布圖Fig.1 Location and distribution of rain gauges in the Xiangxi River Watershed

2 資料和方法

2.1 數(shù)據(jù)來源

本研究采用的1971—2010年各雨量站日降雨數(shù)據(jù)均來自長(zhǎng)江水利委員會(huì)編撰的水文年鑒。數(shù)據(jù)的可靠性和精度滿足計(jì)算要求。各站點(diǎn)基本信息見表1。

2.2 侵蝕力計(jì)算方法

降雨侵蝕力最初提出是針對(duì)一場(chǎng)降雨中的最大30 min降雨強(qiáng)度和降雨動(dòng)能的乘積EI30而言的；但是很多國(guó)家和地區(qū)都缺少場(chǎng)次降雨資料，因此在利用EI30計(jì)算降雨侵蝕力時(shí)受到很大限制。同時(shí)日降雨和次降雨并非一一對(duì)應(yīng)，直接采用日雨量估算日降雨侵蝕力一方面存在困難，另一方面也不準(zhǔn)確。

表1 香溪河流域雨量站位置信息和年降雨特征

在土壤侵蝕模型USLE 和RUSLE 中，降雨侵蝕力季節(jié)變化均以半月為步長(zhǎng)。本文采用章文波等[8]提出的半月時(shí)段模型計(jì)算。

(1)

式中：Rk為第k個(gè)半月的降雨侵蝕力，MJ·mm/(hm2·h)，k=1，2，3，…，24；Pdk為第k個(gè)半月的日侵蝕性降雨量,mm；j為第k個(gè)半月的時(shí)間，d，j=13，14，15，16；a和b為反映當(dāng)?shù)亟涤晏卣鞯哪Ｐ蛥?shù)。計(jì)算公式如下：

(2)

a=21.586b-7.189 1。

(3)

式中：Pd12為日降雨量≥12 mm的日平均降雨量，mm；Py12為日降雨量≥12 mm的年平均降雨量，mm。

2.3 分析方法

2.3.1 氣候傾向率 氣候傾向率在1988年被英國(guó)學(xué)者P. D. Jones[17]用來研究地球表面溫度變化特征之后，在我國(guó)也得到了廣泛的關(guān)注和推廣，被用來研究降水、溫度等氣象要素的變化趨勢(shì)[19-21]。氣候傾向率指每10年平均變化的絕對(duì)值：把氣候要素表示為時(shí)間的線性函數(shù)，即在計(jì)算侵蝕力R變化趨勢(shì)時(shí)，采用最小二乘法，計(jì)算降雨侵蝕力隨時(shí)間t的線性回歸系數(shù)α，α為降雨侵蝕力變化速率，正值表示增加趨勢(shì)，負(fù)值表示減小趨勢(shì)，要素的變化則可用一次線性方程表示:

R=αt+β。

(4)

式中β為截距。文中以線性回歸系數(shù)α的10倍作為降雨侵蝕力的氣候傾向率。

2.3.2 Mann-Kendall趨勢(shì)檢驗(yàn) Mann-Kendall趨勢(shì)檢驗(yàn)法是世界氣象組織推薦并已廣泛使用的非參數(shù)檢驗(yàn)方法。非參數(shù)檢驗(yàn)亦稱為無分布檢驗(yàn)，其優(yōu)點(diǎn)是樣本不需要遵從一定的分布，也不受少數(shù)異常值的干擾，更適合用于類型變量和順序變量，廣泛應(yīng)用于水文和氣象時(shí)間序列的變化趨勢(shì)分析中，詳細(xì)方法參見文獻(xiàn)[22-23]。

降雨侵蝕力的平均值、方差、變異系數(shù)及相關(guān)性等描述性統(tǒng)計(jì)特征采用統(tǒng)計(jì)軟件SPSS計(jì)算，在ArcGIS10.2軟件中采用克里格插值方法對(duì)年降雨量和年降雨侵蝕力進(jìn)行插值。

3 結(jié)果與分析

3.1 年降雨量和年降雨侵蝕力的空間分布

圖2 香溪河流域1971—2010年多年平均降雨量和多年平均降雨侵蝕力空間分布Fig.2 Spatial distribution of annual precipitation and annual rainfall erosivity in the Xiangxi River Watershed during 1971—2010

根據(jù)香溪河流域10個(gè)雨量站的逐日降雨數(shù)據(jù)，采用簡(jiǎn)易算法計(jì)算各個(gè)站點(diǎn)降雨侵蝕力。計(jì)算結(jié)果表明：流域年降雨量在688.12～1 318.78 mm之間，多年均值為985.86 mm，降雨量呈現(xiàn)由西北向東南遞減趨勢(shì)(圖2(a))，與流域高程變化趨勢(shì)基本一致。流域西北部的紅花和九沖站多年平均降雨量均超過1 250 mm，比流域均值多近30%。多年降雨最少區(qū)域?yàn)榍嗌秸?、水月寺、鄭家坪、興山和峽口站附近。

流域降雨侵蝕力的高值區(qū)位于西北部的紅花、九沖河站一帶，年均降雨侵蝕力均超過6 000 MJ·mm/(hm2·h)；侵蝕力低值區(qū)在峽口鎮(zhèn)、水月寺和青山站附近，多年均值均小于4 500 MJ·mm/(hm2·h)，峽口站最小，降雨侵蝕力多年均值僅為3 999.4 MJ·mm/(hm2·h )。10個(gè)站點(diǎn)中，降雨侵蝕力年際變化最大的是九沖站，其極值比達(dá)到了4.43，中陽埡極值比最小為2.47，其余8個(gè)站點(diǎn)的極值比在3～4之間。用克里格插值法對(duì)站點(diǎn)多年平均侵蝕力進(jìn)行空間插值，繪制出流域降雨侵蝕力的等值線圖(圖2(b))。香溪河流域降雨侵蝕力的變化趨勢(shì)為從西向東逐漸遞減。流域以興山和鄭家坪所在位置為分界線，其以西區(qū)域年降雨侵蝕力值大且變化劇烈，年侵蝕力值隨高程的降低而顯著減小；以東區(qū)域年降雨侵蝕值小，但侵蝕力變化不大，基本分布在4 000～4 500 MJ·mm/(hm2·h)。本文采用“泰森多邊形法”確定每個(gè)雨量站對(duì)應(yīng)的權(quán)重，流域內(nèi)10個(gè)雨量站較好地反映了流域降雨侵蝕力空間異質(zhì)性，根據(jù)每個(gè)雨量站權(quán)重和該雨量站的降雨侵蝕力數(shù)據(jù)，計(jì)算得到香溪河流域多年侵蝕力變化范圍為2 465.26～7 419.29 MJ·mm/(hm2·h)，多年均值為4 535.63 MJ·mm/(hm2·h·a)。

流域年降雨量和年降雨侵蝕力分布基本一致，表明該區(qū)域降雨量和降雨侵蝕力之間有較好的協(xié)同性。具體表現(xiàn)為：流域西北部降雨量和降雨侵蝕力最大；流域東南部降雨量和降雨侵蝕力最小。鄭家坪和興山站一帶年降雨量較小，但是降雨侵蝕力卻相對(duì)較大，表明流域中部鄭家坪、興山一帶的年降雨中侵蝕性降雨量比例較大。青山站和水月寺站降雨量屬中等降雨程度，但是年降雨侵蝕力卻最小。表明流域東南部青山、水月寺一帶降雨量雖大，但侵蝕性降雨量相對(duì)較少；與此相同情況也發(fā)生在流域中北部的水果園區(qū)域：水果園多年平均降雨量比中陽埡多30 mm，但多年降雨侵蝕力均值差異很小。

3.2 降雨侵蝕力年際變化

通過計(jì)算香溪河流域各站點(diǎn)的降雨侵蝕力氣候傾向率，發(fā)現(xiàn)1971—2010年間，流域南陽河、興山和峽口3個(gè)站點(diǎn)表現(xiàn)正氣候傾向率，每10年增加量均在200 MJ·mm/(hm2·h)以上，均值為227.3 MJ·mm/(hm2·h)，其中南陽河的正氣候傾向率最大，每10年增加量為354.8 MJ·mm/(hm2·h)；而其余6個(gè)站點(diǎn)降雨侵蝕力年際變化表現(xiàn)為不同程度的負(fù)趨勢(shì)，即降雨引起土壤侵蝕的能力在減小，平均每10年減小313.8 MJ·mm/(hm2·h)；流域侵蝕力總體上為減小趨勢(shì)，但不顯著。降雨侵蝕力的這種年際變化特征(圖3)與降雨量和降雨強(qiáng)度密切相關(guān)，但在某些年份二者關(guān)系并不一致：1979年由于侵蝕性降雨量大而降雨總量較小使得侵蝕力較大；但1980年和1994年降雨總量較大，而降雨侵蝕力較小，原因在于這2年降雨量雖大但侵蝕性降雨量相對(duì)較小。

同時(shí)還可采用Mann-kendall趨勢(shì)檢驗(yàn)對(duì)流域內(nèi)10個(gè)站點(diǎn)降雨和降雨侵蝕力做趨勢(shì)分析，結(jié)果見表2。Z為Mann-kendall檢驗(yàn)的統(tǒng)計(jì)量，在0.05顯著性水平時(shí)，|Z1-0.05/2|=1.96和|Z1-0.10/2|=1.64，可判斷序列的變化趨勢(shì)及其顯著性。Mann-kendall對(duì)站點(diǎn)侵蝕力趨勢(shì)檢驗(yàn)結(jié)果與氣候傾向率結(jié)果一致，但所有站點(diǎn)的趨勢(shì)變化均沒有通過0.05顯著性水平的檢驗(yàn)，僅紅花站通過了0.10顯著性水平的檢驗(yàn)。對(duì)流域各站點(diǎn)降雨Mann-kendall檢驗(yàn)結(jié)果發(fā)現(xiàn)：只有南陽河站是微弱上升趨勢(shì)，峽口和興山站點(diǎn)都是下降趨勢(shì)，這個(gè)侵蝕力變化趨勢(shì)不一致。九沖站降雨減小趨勢(shì)超過了0.05顯著性水平的檢驗(yàn)，但紅花站沒有發(fā)生顯著的變化趨勢(shì)。說明流域年侵蝕力變化存在一定的隨機(jī)波動(dòng)，但無顯著地氣候趨勢(shì)。

3.3 降雨侵蝕力年內(nèi)分布

4—10月為長(zhǎng)江流域的汛期，10個(gè)站點(diǎn)汛期降雨量占了全年降雨量的84.0%～87.0%，汛期侵蝕性降雨量占了全年侵蝕性降雨量的91.0%～94.2%，汛期降雨侵蝕力占了全年降雨侵蝕力的92.6%～95.5%，三者的均值分別是85.4%、92.4%和94.0%，表明香溪河流域降雨量、侵蝕性降雨量和降雨侵蝕力在年內(nèi)分布集中。各雨量站點(diǎn)降雨侵蝕力的年內(nèi)變化和降雨量年內(nèi)分配特征一致。因?yàn)榻涤昵治g力年內(nèi)分配的集中程度對(duì)于研究土壤侵蝕的季節(jié)變化十分重要，故以連續(xù)6個(gè)半月R值和的最大值與年平均R值的比值來表示侵蝕力年內(nèi)分配的集中程度。香溪河流域10個(gè)站集中程度均值為57.8%，表明降雨侵蝕力R值年內(nèi)分布集中度較大。相對(duì)而言鄭家坪站的集中程度最高為61.8%，而峽口站降雨侵蝕力分布最均勻，集中程度為53.9%，其余站點(diǎn)都分布在55%～60%之間。年降雨侵蝕力的分布最集中的時(shí)段是從6月下半月到9月上半月(圖4)，該時(shí)段的農(nóng)耕活動(dòng)和開發(fā)建設(shè)要遵循侵蝕力分布規(guī)律，避免侵蝕加劇。

圖3 香溪河流域1971—2010年降雨侵蝕力和年降雨量變化Fig.3 Annual precipitation and annual rainfall erosivity in the Xiangxi River Watershed from 1971 to 2010

表2 香溪河流域各站點(diǎn)降雨量和降雨侵蝕力MK趨勢(shì)檢驗(yàn)結(jié)果

圖4 流域各站點(diǎn)半月降雨侵蝕力的年內(nèi)分布Fig.4 Inter-annual distribution of semi-monthly rainfall erosivity for rain gauges in the Xiangxi River Watershed

4 結(jié)論

1)香溪河流域多年降雨侵蝕力變化范圍為2 465.26～7 419.29 MJ·mm/(hm2·h)，多年均值為4 535.63 MJ·mm/(hm2·h)。流域降雨侵蝕力最大值位于其西北部，降雨侵蝕力總的變化趨勢(shì)從西北向東南逐漸遞減，流域以興山和鄭家坪所在位置為分界線，以西區(qū)域年降雨侵蝕力大，年降雨侵蝕力值隨高程的降低而顯著減小；以東區(qū)域年降雨侵蝕力小，且變化緩慢。

2)各站點(diǎn)年降雨量和年降雨侵蝕力年際變化分別以氣候傾向率和Mann-kendall 2種方法表示時(shí)，除了九沖站降雨量有顯著性減小趨勢(shì)(α=0.05)、紅花站降雨侵蝕力有明顯減小趨勢(shì)(α=0.10)以外，其余站點(diǎn)均沒有顯著的變化趨勢(shì)。表明香溪河流域近40多年的降雨量和降雨侵蝕力比較平穩(wěn)，無大的波動(dòng)。

3)汛期降雨量占全年降雨總量的84%～87%，汛期侵蝕性降雨量占全年侵蝕性降雨總量的 91%～94.2%，汛期降雨侵蝕力占全年降雨侵蝕力的92.6%～95.5%，三者的均值分別是85.4%、92.4%和94%，表明香溪河流域侵蝕性降雨和降雨侵蝕力年內(nèi)分布較集中。年降雨侵蝕力分布最集中的時(shí)段是從6月下半月到9月上半月，此時(shí)段的農(nóng)耕活動(dòng)和開發(fā)建設(shè)要遵循侵蝕力分布規(guī)律，避免侵蝕加劇。

[1] Wischmeier W H，Smith D D．Predicting rainfall erosion losses： a guide to conservation planning[M]．Washington：Agriculture Handbank，1978

[2] 王萬忠，焦菊英，郝小品，等．中國(guó)降雨侵蝕力R值的計(jì)算與分布：I[J]．水土保持學(xué)報(bào)，1995，9(4)：5-18

[3] Renard K G, Freimund J R. Using monthly precipitation data to estimate the R-factor in the revised USLE[J]. Journal of Hydrology, 1994，174：287-306

[4] Richardson C W，F(xiàn)oster G R，Wright D A．Estimation of erosion index from daily rainfall amount [J]．Transactions of the ASAE，1983，26(1)：153-156

[5] Richardson C W，F(xiàn)oster G R，Wright D A. Estimation of erosion index from daily rainfall amount [J]．Transactions of the ASAE，1983，26(1)：153-156

[6] 章文波，謝云，劉寶元．利用日雨量計(jì)算降雨侵蝕力的方法研究[J]．地理科學(xué)，2002，22(6)：705-711

[7] 章文波，付金生．不同類型雨量資料估算降雨侵蝕力[J].資源科學(xué)，2003，25(1):35-41

[8] 謝云，劉寶元，章文波．侵蝕性降雨標(biāo)準(zhǔn)研究[J]．水土保持學(xué)報(bào)，2000，14(4)：6-11

[9] 閆業(yè)超，岳書平，張樹文．降雨資料時(shí)間序列長(zhǎng)度對(duì)降雨侵蝕力平均值置信度的影響[J]．自然資源學(xué)報(bào)，2012，28(2)：21-27

[10] 吳昌廣，林德生，肖文發(fā)，等．三峽庫區(qū)降雨侵蝕力時(shí)空分布特征[J]．應(yīng)用生態(tài)學(xué)報(bào)，2011，22(1)：151-158

[11] 汪言在，閏冬，李月臣，等．日降雨對(duì)降雨侵蝕力年雨量簡(jiǎn)易算法的影響分析[J]．水土保持通報(bào)，2012，32(3)：150-155

[12] 張照錄，薛重生．湖北三峽庫區(qū)降雨侵蝕力的計(jì)算及其特征分析[J]．中國(guó)水土保持，2005(8)：10-12

[13] 花利忠，賀秀斌，顏昌宙，等．三峽庫區(qū)大寧河流域降雨侵蝕力時(shí)空分布特征[J]．水土保持通報(bào)，2008，28(4)：22-25

[14] 史東梅，盧喜平，蔣光毅，等．紫色丘陵區(qū)降雨侵蝕力簡(jiǎn)易算法的模擬[J]．農(nóng)業(yè)工程學(xué)報(bào)，2010，26(2)：116-123

[15] 繆馳遠(yuǎn)，徐霞，魏欣，等．重慶市主城區(qū)降雨侵蝕力特征分析[J]．資源科學(xué)，2007，29(4)：54-61

[16] 張革，劉德富，宋林旭，等．香溪河流域降雨侵蝕力特征及簡(jiǎn)易算法初步研究[J]．水土保持研究，2013，20(5)：7-12

[17] Jones P D. Hemispheric surface air temperature variations: recent trend and an updata to 1978[J]. Journal of Climate,1988,1:654-660

[18] 孫長(zhǎng)安.香溪河流域土地利用與水土流失的關(guān)系研究[D].北京：北京林業(yè)大學(xué)，2008：22-56

[19] 徐國(guó)昌，姚輝，李珊.我國(guó)干旱半干旱地區(qū)現(xiàn)代降水量和歷史干旱頻率的變化[J].氣象學(xué)報(bào)，1992,50(3)：378-382

[20] 任國(guó)玉，吳虹，陳正宏.我國(guó)降水變化趨勢(shì)的空間特征[J].應(yīng)用氣象學(xué)報(bào)，2000,11(3):322-330

[21] 章文波，劉寶元.北方農(nóng)牧交錯(cuò)帶降水極值變化空間特征[J].自然資源學(xué)報(bào)，2003,18(3):274-280

[22] 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京: 氣象出版社，2007

[23] Wan ZawishWan Zin, Suhaila Jamaludin, Sayang Mohd Deni, et al. Recent changes in extreme rainfall events in Peninsular Malaysia:1971—2005[J]. Theor Appl Climatol, 2010, 99:303-314

(責(zé)任編輯：程 云)

Temporal and spatial variations of rainfall erosivity in the Xiangxi River watershed of the Three Gorges Reservoir region

Liu Huiying1，Ren Hongyu2，Zhang Changwei2，Zhang Pingcang2

(1. Nanchang Institute of Technology，330099， Nanchang，China; 2.Yangtze River Scientific Research Institute，430012，Wuhan，China)

Based on the daily rainfall data of ten rain gauging stations from 1971 to 2010 in the Xiangxi River watershed of the Three Gorges Reservoir region，we studied the spatio-temporal distribution characteristics of rainfall erosivity estimated with our focus on the annual and inter-annual trends of the rainfall erosivity (R-factor) for the watershed through Mann Kendall nonparametric tests and Kriging space interpolation of Arcgis software. The results are shown as follows. 1) From 1971 to 2010, the annual rainfall erosivity of the watershed changed from 2 465.26 to 7 419.29 MJ·mm/(hm2·h), with an average annual rainfall erosivity of 4 535.63 MJ·mm/(hm2·h). 2) There was a great variation of inter-annual R-factor with the maximum three times greater than the minimum. The intra-annual distribution of rainfall，erosive rainfall and rainfall erosivity were highly concentrated with a single peak, mainly during the period from April to October，accounting for 85.4%, 92.4% and 94% of whole year values, respectively. 3) For the whole Xiangxi River watershed, annual precipitation and rainfall erosivity had no significant tendencies of change over 40 years. 4) The results also showed that the spatial distribution of R in the Xiangxi River Watershed decreased rapidly from the west to the east.

rainfall erosivity； daily rainfall； spatio-temporal distribution； Xiangxi river watershed

2014-09-22

2015-04-07

劉惠英(1973—)，女，講師，博士研究生。主要研究方向：坡面土壤侵蝕和流域水體監(jiān)測(cè)。E-mail：jlfx7401@163.com

?通信作者簡(jiǎn)介： 張平倉(1961—)，男，教授，博士生導(dǎo)師。主要研究方向：土壤侵蝕。E-mail：zhangpc@mail.crsri.cn

S157.1； TP79

A

1672-3007(2015)03-0001-07

項(xiàng)目名稱： 水利部公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目“長(zhǎng)江流域山洪災(zāi)害區(qū)域特征及防御體系研究”(201301059)；江西省土壤侵蝕與防治重點(diǎn)實(shí)驗(yàn)室開放基金“贛江上游水沙時(shí)空演變及其植被恢復(fù)的響應(yīng)”(JXSB201303)；江西省自然科學(xué)基金“贛江流域產(chǎn)匯流參數(shù)及水沙變化對(duì)LUCC的響應(yīng)研究”(20132BAB203032)