龍訓(xùn)建， 翁薛柔， 葉 琰， 張逸軒， 徐廷兵， 葉 勇

(1.西南大學(xué) 資源環(huán)境學(xué)院， 重慶 400715； 2.重慶市人工影響天氣辦公室，重慶 400074； 3.西南大學(xué) 后勤保障部， 重慶 400715)

1 研究背景

全球氣候變化及生態(tài)環(huán)境問題是聯(lián)合國政府間氣候變化專門委員會(huì)(Iternational Panel on Climate Change, IPCC)的重要關(guān)注點(diǎn)，水土流失是其多年評(píng)估報(bào)告均涉及的生態(tài)環(huán)境問題之一[1]。降雨是產(chǎn)生陸域水土流失的能量來源[2-4]。水土流失不僅會(huì)引起土壤及其營養(yǎng)因素的大量流失，還會(huì)造成陸域土地退化、河海塘庫泥沙淤積、蓄水及防洪干旱能力受限[5-7]。導(dǎo)致水土流失的侵蝕性降雨通常只包括會(huì)形成地表徑流并引發(fā)土壤侵蝕的降雨[8]，這已被不同區(qū)域分布、時(shí)間尺度和土壤類型的研究成果驗(yàn)證[4,9-10]。最具代表性和普適性的侵蝕性降雨臨界值是Wischmeier等[11]確定的次降雨量≥12.7 mm或一次降雨的15 min雨量≥6.4 mm，但我國學(xué)者針對(duì)紅壤、黃土、紫壤等多種土壤類型，應(yīng)用頻率分析法[12]、侵蝕力偏差系數(shù)[13]、徑流小區(qū)實(shí)測(cè)[14]方法建立了不同的適合各自研究區(qū)域的侵蝕性降雨標(biāo)準(zhǔn)，其中，應(yīng)用最為廣泛的是次降雨量≥12 mm[8,15-16]。

次降雨的界定標(biāo)準(zhǔn)對(duì)深入了解土壤水蝕的過程和機(jī)理有重要意義，但對(duì)于以行政區(qū)劃為管理單元的水土保持政策措施以及管理工作，還需有侵蝕性降雨長時(shí)間序列的變化特點(diǎn)、時(shí)空分布特征的理論支撐，尤其應(yīng)掌握侵蝕性降雨特征參數(shù)中雨量、雨強(qiáng)、歷時(shí)、發(fā)生次數(shù)等特征值的時(shí)空變化。重慶屬于我國水土保持區(qū)劃中的西南紫色土區(qū)，是國家級(jí)水土流失重點(diǎn)防治區(qū)[17]?！吨貞c市水土保持公報(bào)》結(jié)果表明[18]，截至2018年，重慶幅員面積的31.32%仍然因水力侵蝕和重力侵蝕產(chǎn)生水土流失。另外，受數(shù)據(jù)采集難度、時(shí)間序列尺度限制影響，重慶市侵蝕性降雨和降雨侵蝕力的已有研究成果多集中于小流域的實(shí)驗(yàn)性分析或單個(gè)行政區(qū)域的初步研究[9,19-20]，基于多站點(diǎn)的侵蝕性降雨特征因子與土壤侵蝕之間的統(tǒng)計(jì)學(xué)相關(guān)分析研究還未見報(bào)道。因此，基于重慶市32個(gè)國家級(jí)氣象站點(diǎn)2009-2018年的逐分鐘降雨量數(shù)據(jù)，將數(shù)理統(tǒng)計(jì)與空間插值相結(jié)合，采用聚類分析方法，分析探討區(qū)域侵蝕性降雨的雨量、歷時(shí)、雨強(qiáng)、發(fā)生次數(shù)、降雨侵蝕力以及侵蝕性雨情的時(shí)空分布特征，可為區(qū)域水土流失預(yù)測(cè)、生態(tài)環(huán)境保護(hù)治理和災(zāi)害防治工作提供理論參考。

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

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

基礎(chǔ)數(shù)據(jù)來源于重慶市32個(gè)國家級(jí)氣象站點(diǎn)2009-2018年的逐分鐘降雨量數(shù)據(jù)。降雨事件中降雨間歇時(shí)間6 h以上或連續(xù)6 h降雨量不足1.2 mm，則視為2次降雨事件，否則為1次降雨事件[21]。次降雨事件并不都會(huì)引起水土流失，結(jié)合研究區(qū)內(nèi)相關(guān)文獻(xiàn)的研究結(jié)果[19]，將次降雨量高于12 mm的降雨事件劃分為侵蝕性降雨。依據(jù)重慶市氣候中心統(tǒng)計(jì)本地30 a以上的氣溫平均值進(jìn)行四季劃分，其中3-5月為春季，6-9月為夏季，10-12月為秋季，1-2月為冬季。同時(shí)，在地理區(qū)域上，將研究區(qū)劃分為渝西、渝東南和渝東北3個(gè)片區(qū)，重慶市地理區(qū)域劃分及氣象站點(diǎn)位置見圖1。

圖1 重慶市地理區(qū)域劃分及氣象站點(diǎn)位置

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

2.2.1 降雨侵蝕力 單場(chǎng)降雨侵蝕力是雨強(qiáng)和動(dòng)能的綜合作用結(jié)果[3]，定義為：

R=E·I30

(1)

式中：I30為單場(chǎng)次降雨最大30 min雨強(qiáng)，mm/h，根據(jù)分鐘降雨量統(tǒng)計(jì)得到；E為單場(chǎng)次降雨動(dòng)能，MJ/hm2，是降雨-歷時(shí)函數(shù)，其計(jì)算方法見公式(2)[22]。

(2)

式中：D為次降雨歷時(shí)，h；v為次降雨量，mm；i為降雨強(qiáng)度，mm/h。

2.2.2 空間插值 采用精確性插值法的反距離權(quán)重法(inverse distance weighted，IDW)實(shí)現(xiàn)空間插值。基于插值區(qū)域內(nèi)部樣本點(diǎn)的相似性，計(jì)算與到鄰近區(qū)域樣點(diǎn)的加權(quán)平均值來估算出單元格的值，進(jìn)而插值得到面值[23]。計(jì)算公式如下：

(3)

式中：v0為未知點(diǎn)的估計(jì)值；vi為第i(i=1,2,…,n)個(gè)樣本點(diǎn)的值；di為采樣點(diǎn)與未知點(diǎn)之間的距離；k為距離的冪值，默認(rèn)值為2，本次取0.5～3。

2.2.3 聚類分析 聚類分析按照一定規(guī)則進(jìn)行數(shù)據(jù)分類，判別出相似性較高的同一類對(duì)象和差異較大的不同類對(duì)象[24]。系統(tǒng)聚類法、K均值聚類法(K-means)和主成分聚類法是目前最常用的3種方法[25-26]?？紤]到雨情分類過程中變量有限、數(shù)據(jù)較多的特點(diǎn)，結(jié)合各種聚類分析方法的適用范圍，本文應(yīng)用K均值聚類法將研究區(qū)內(nèi)2009-2018年發(fā)生的所有侵蝕性降雨進(jìn)行雨情分類，其中分類特征變量包括侵蝕性降雨雨量、降雨歷時(shí)和最大30 min降雨強(qiáng)度，降雨事件發(fā)生次數(shù)為標(biāo)記分析區(qū)域內(nèi)侵蝕性降雨的雨情特征[19]。其次，應(yīng)用K均值聚類法對(duì)32個(gè)氣象站點(diǎn)的侵蝕性降雨特征值進(jìn)行分類，獲得研究區(qū)空間上的雨情分區(qū)。

2.2.4 趨勢(shì)變化參數(shù) 時(shí)間序列的年際變化特征利用離差系數(shù)CV和趨勢(shì)系數(shù)r來描述，其中CV越大，表明年降雨侵蝕力波動(dòng)幅度越大；r>0表明降雨侵蝕力在研究時(shí)段內(nèi)有線性增加的趨勢(shì)，反之亦然，∣r∣大小表示隨時(shí)間變化快慢幅度[27]。CV和r的計(jì)算公式如下：

(4)

(5)

3 結(jié)果與分析

3.1 侵蝕性降雨時(shí)間變化特征

3.1.1 年際變化特征 對(duì)2009-2018年重慶市侵蝕性降雨的雨量、歷時(shí)、降雨強(qiáng)度、發(fā)生次數(shù)及占總降雨的百分比進(jìn)行統(tǒng)計(jì)，結(jié)果見表1。

表1 2009-2018年重慶市侵蝕性降雨特征參數(shù)統(tǒng)計(jì)表

由表1可知，重慶市年均侵蝕性降雨量為807.0 mm，占年均降雨量的75%。其中，渝東北片區(qū)年平均侵蝕性降雨量最大，達(dá)865.5 mm，占年均降雨量的78%；渝西片區(qū)年平均侵蝕性降雨量最小，為768.3 mm，占年均降雨量的73%。降雨量和侵蝕性降雨在研究區(qū)的分布上，均呈現(xiàn)出由東向西遞減的特征。這說明渝東北片區(qū)不僅年降雨總量相對(duì)豐富，侵蝕性降雨發(fā)生概率也較大；渝西地區(qū)年降雨總量相對(duì)偏少，產(chǎn)生侵蝕性降雨的總雨量也相對(duì)較小。

研究區(qū)年均侵蝕性降雨歷時(shí)相差較小，在413～425 h之間，其中歷時(shí)最長的425 h為渝西片區(qū)，歷時(shí)最短的413 h在渝東南片區(qū)；在歷時(shí)占比上，則是渝東北最長，達(dá)48%，渝西最少，為42%。研究區(qū)年均侵蝕性降雨強(qiáng)度為71.2 mm/h，占年總降雨強(qiáng)度的55%；高值區(qū)分布在渝東南片區(qū)，為73.5 mm/h；低值區(qū)分布于渝西片區(qū)，僅69.8 mm/h。對(duì)2009-2018年期間發(fā)生侵蝕性降雨次數(shù)的統(tǒng)計(jì)結(jié)果表明，研究區(qū)內(nèi)發(fā)生侵蝕性降雨的平均次數(shù)差距不大，為24～26次，僅占年均總降雨發(fā)生次數(shù)的17%～20%，但此部分侵蝕性降雨的雨量占年均降雨量的73%～78%。由此說明，研究區(qū)內(nèi)侵蝕性降雨的發(fā)生次數(shù)雖遠(yuǎn)低于全年降雨次數(shù)，歷時(shí)不及總降雨歷時(shí)的一半，但侵蝕性雨量的占比達(dá)到了3/4，對(duì)全年水力侵蝕起主要貢獻(xiàn)作用。

3.1.2 年內(nèi)變化特征 2009-2018年重慶市侵蝕性降雨特征平均值的年內(nèi)變化曲線見圖2。由圖2可知，侵蝕性雨量年內(nèi)變化呈雙峰特征(2(a))，峰值集中發(fā)生在5-6月和9月，在此期間的累計(jì)侵蝕性雨量占比約為全年侵蝕性雨量的47%，占全年降雨量的35%；侵蝕性降雨歷時(shí)年內(nèi)變化也呈雙峰特征(2(b))，峰值月份為5和9月，其累計(jì)歷時(shí)約占全年侵蝕性降雨歷時(shí)的30%；侵蝕性雨強(qiáng)呈單峰特征(2(c))，峰值出現(xiàn)在8月，該時(shí)段內(nèi)侵蝕性降雨歷時(shí)僅為全年的16%，但侵蝕性雨量占全年的27%，即該時(shí)段侵蝕性降雨特征表現(xiàn)為雨量大、歷時(shí)短，這與區(qū)域在該時(shí)段內(nèi)易發(fā)生小范圍強(qiáng)對(duì)流天氣有關(guān)。每年5-6月發(fā)生侵蝕性降雨次數(shù)最多(2(d))，約占全年侵蝕性降雨發(fā)生次數(shù)的31%。雨強(qiáng)同時(shí)受雨量和歷時(shí)影響，結(jié)合張弛等[28]對(duì)西南地區(qū)夏季極端降水的水汽來源分析可知，每年的5-9月西太平洋副熱帶高壓西伸北移，西南季風(fēng)增強(qiáng)，下沉氣流控制重慶大部分地區(qū)，理論上降雨可能性較小，但該地區(qū)大氣層結(jié)不穩(wěn)定，絕對(duì)水汽含量高，易產(chǎn)生短歷時(shí)的強(qiáng)對(duì)流性降雨，從而表現(xiàn)出較高的強(qiáng)降雨占比[29]。另外，12月至次年2月研究區(qū)氣溫低，大氣層結(jié)穩(wěn)定，絕對(duì)水汽含量低，降雨量少，因而侵蝕性降雨的雨量、歷時(shí)、雨強(qiáng)及發(fā)生次數(shù)在12月至次年2月期間分別僅占全年的2%、9%、3%和5%。

圖2 2009-2018年重慶市侵蝕性降雨特征平均值的年內(nèi)變化曲線

2009-2018年重慶市降雨侵蝕力平均值的年內(nèi)變化曲線見圖3。

圖3 2009-2018年重慶市降雨侵蝕力平均值的年內(nèi)變化曲線

由圖3可看出，降雨侵蝕力年內(nèi)變化曲線為單峰型，各片區(qū)峰值均出現(xiàn)在7月，且重慶市、渝西片區(qū)、渝東南片區(qū)和渝東北片區(qū)在該月的降雨侵蝕力分別占全年降雨侵蝕力的24%、24%、26%和22%；5-9月各區(qū)域累計(jì)降雨侵蝕力分別占全年降雨侵蝕力的89%、90%、87%和89%，此時(shí)段屬于侵蝕性降雨的集中發(fā)生時(shí)段。此外，這一時(shí)期集中的高強(qiáng)度傳統(tǒng)農(nóng)業(yè)種植方式所必須的土壤翻耕也是加劇區(qū)域土壤侵蝕的原因之一。

3.2 侵蝕性降雨空間變化特征

3.2.1 空間變化特征 基于反距離權(quán)重法得到的研究區(qū)多年平均侵蝕性降雨各參數(shù)空間分布如圖4所示。由圖4可見，侵蝕性降雨量(圖4(a))高值區(qū)主要分布在渝東南的酉陽、渝東北的城口、開州和云陽；次高值區(qū)主要分布在渝東南的彭水，渝東北的巫溪和梁平。整體上，長江沿線河谷區(qū)域的降雨量大多低于兩側(cè)山區(qū)，即侵蝕性降雨量低值區(qū)位于渝西片區(qū)的潼南-大足-榮昌-永川一帶，多年平均值介于680～710 mm之間；侵蝕性降雨量高值區(qū)主要位于渝東北的城口-開州-云陽一帶，多年平均值超過了1 000 mm，是全市年平均侵蝕性降雨量的1.25倍，是低值區(qū)的1.47倍。方德賢等[30]利用地面站點(diǎn)數(shù)據(jù)和CMORPH融合逐小時(shí)降水產(chǎn)品對(duì)重慶市2008-2016年的年降雨量空間分布特征進(jìn)行了分析，其結(jié)果與本次得到的多年平均侵蝕性降雨量的空間分布具有一致性。侵蝕性降雨歷時(shí)(圖4(b))高值區(qū)主要分布在渝東北片區(qū)的城口和渝西片區(qū)的萬盛，多年平均侵蝕性降雨歷時(shí)達(dá)530 h，低值區(qū)主要分布在渝西片區(qū)的潼南-大足-榮昌-永川一帶和渝東北片區(qū)的奉節(jié)、巫山和巫溪，年均降雨歷時(shí)為370 h，僅為高值區(qū)侵蝕性降雨歷時(shí)的70%。研究區(qū)侵蝕性降雨強(qiáng)度(圖4(c))介于23.5～30.6 mm/h之間，其中，最大值在渝東南片區(qū)的酉陽，其次是渝東北片區(qū)的巫溪和開州。多年平均侵蝕性降雨發(fā)生次數(shù)(圖4(d))最多的地區(qū)為渝西片區(qū)的萬盛、渝東北片區(qū)的城口、渝東南片區(qū)的酉陽，低值中心則集中在渝西片區(qū)的潼南-大足-榮昌-永川一帶，且渝西片區(qū)從東南到西北發(fā)生侵蝕性降雨的次數(shù)呈下降趨勢(shì)。

圖4 2009-2018年重慶市多年平均侵蝕性降雨特征值空間分布

圖5為2009-2018年研究區(qū)多年平均降雨侵蝕力的空間分布。圖5顯示，多年平均降雨侵蝕力高值區(qū)分布在渝東北片區(qū)的開州、渝西片區(qū)的銅梁和巴南，達(dá)2 500 MJ·mm/(hm2·h)；渝東北片區(qū)的梁平、渝東南片區(qū)的酉陽為次高值區(qū)，為2 330 MJ·mm/(hm2·h)；渝西片區(qū)的南川區(qū)最小，僅為1 220 MJ·mm/(hm2·h)。結(jié)合研究區(qū)域地形地貌、土壤類型和人類活動(dòng)等因素綜合分析多年平均降雨侵蝕力空間分布差異顯著的原因，高值分布區(qū)的開州、梁平和酉陽主要為山地丘陵區(qū)，地形復(fù)雜，坡度較大，又多黃壤、石灰?guī)r土、黃棕壤等易于發(fā)生降雨侵蝕的土壤類型，因此更易產(chǎn)生土壤流失；渝西的高值區(qū)主要分布在土地復(fù)種指數(shù)較高的區(qū)縣。這些區(qū)域的降雨侵蝕潛在危險(xiǎn)性較高，降雨導(dǎo)致的水土流失、土壤侵蝕等問題需持續(xù)關(guān)注。

3.2.2 侵蝕性降雨雨情分區(qū) 侵蝕性降雨的年際和年內(nèi)分布均表現(xiàn)出不同的閾值區(qū)域，不同參數(shù)在不同時(shí)間不同區(qū)域出現(xiàn)極值。將侵蝕性降雨量、降雨歷時(shí)、最大30 min降雨強(qiáng)度(I30)、降雨發(fā)生次數(shù)以及產(chǎn)生的降雨侵蝕力作為分類因子，利用主成分聚類法對(duì)32個(gè)氣象站點(diǎn)進(jìn)行侵蝕性雨情分區(qū)，以探討研究區(qū)內(nèi)侵蝕性降雨在空間分布上的共性特點(diǎn)，結(jié)果見圖6。

圖5 2009-2018年重慶市多年平均降雨侵蝕力的空間分布 圖6 基于主成分聚類的重慶市侵蝕性降雨雨情特征分區(qū)

主成分聚類法將研究區(qū)劃分為3個(gè)分區(qū)。第1類分區(qū)為雨量大、歷時(shí)長、雨強(qiáng)較高的強(qiáng)烈侵蝕區(qū)，包括銅梁、巴南、沙坪壩、合川、酉陽、彭水、梁平、巫溪、云陽、開州、城口11個(gè)行政區(qū)域。這些區(qū)域潛在土壤侵蝕危險(xiǎn)性相對(duì)最高，可劃為重度侵蝕區(qū)；第2類涵蓋5個(gè)行政區(qū)域，包括長壽、南川、石柱、忠縣、豐都，屬于雨強(qiáng)最小區(qū)域，因而這些區(qū)降雨侵蝕力低于其他兩類分區(qū)，劃為輕度侵蝕區(qū)；其余區(qū)域?yàn)榈?類，歸為中度侵蝕區(qū)。將分區(qū)結(jié)果與水力侵蝕危險(xiǎn)性相對(duì)應(yīng)，可將第1、2、3類設(shè)為重點(diǎn)水土流失防治區(qū)、一般水土流失防治區(qū)和中等水土流失防治區(qū)。

3.3 侵蝕性降雨雨情特征

3.3.1 侵蝕性降雨事件聚類結(jié)果 按照次降雨和侵蝕性降雨統(tǒng)計(jì)標(biāo)準(zhǔn)，2009-2018年期間，研究區(qū)共發(fā)生44 645場(chǎng)降雨，其中8 016場(chǎng)為侵蝕性降雨，占總降雨場(chǎng)數(shù)的18%，但侵蝕性雨量累計(jì)值占降雨總量的75%。單場(chǎng)侵蝕性降雨的雨量平均值為34.6 mm，平均歷時(shí)為16.5 h，平均產(chǎn)生82.16 MJ·mm/(hm2·h)的降雨侵蝕力。單場(chǎng)降雨侵蝕力變化劇烈，存在極端暴雨事件，以2009年銅梁站記錄的一場(chǎng)特大暴雨為例，該次雨量累計(jì)達(dá)396.9 mm，歷時(shí)為46 h，產(chǎn)生的降雨侵蝕力達(dá)3 953.13 MJ·mm/(hm2·h)，為研究時(shí)段內(nèi)銅梁平均降雨侵蝕力的38倍，區(qū)域平均值的48倍，是研究時(shí)段內(nèi)發(fā)生的最嚴(yán)重的一場(chǎng)侵蝕性降雨。

通常水土流失的嚴(yán)重程度受侵蝕性降雨雨情特征影響明顯[31-32]。因此，以侵蝕性降雨的平均降雨量P、平均降雨歷時(shí)T、最大30 min降雨強(qiáng)度I30和降雨次數(shù)為分類特征變量，采用K均值聚類法對(duì)研究區(qū)內(nèi)的雨情特征進(jìn)行分類，得到不同雨情分類相應(yīng)降雨侵蝕力的特點(diǎn)，結(jié)果見表2。研究時(shí)段內(nèi)的8 016場(chǎng)侵蝕性降雨可劃分為3類，第Ⅰ類屬于小雨量、短歷時(shí)、低強(qiáng)度、高頻率的雨情類型，其產(chǎn)生的平均降雨侵蝕力僅為24.96 MJ·mm/(hm2·h)，共計(jì)發(fā)生5 897次，占侵蝕性降雨總次數(shù)的74%，單次侵蝕量少但累計(jì)占比較大，是區(qū)域內(nèi)主要的侵蝕降雨類型；第Ⅱ類雨情的分類特征變量值均介于第Ⅰ類和第Ⅲ類之間，產(chǎn)生的降雨侵蝕力約為147.21 MJ·mm/(hm2·h)。第Ⅲ類雨情發(fā)生次數(shù)較少，但雨量大、歷時(shí)長、雨強(qiáng)高，2009-2018年總計(jì)發(fā)生325次，僅占侵蝕性雨量場(chǎng)次的4%，對(duì)應(yīng)產(chǎn)生的降雨侵蝕力達(dá)600.88 MJ·mm/(hm2·h)，是第Ⅰ類雨情的24倍，第Ⅱ類雨情的4倍，這類雨情具有單次影響巨大，發(fā)生頻率相對(duì)較低，總量占比不高，但易造成嚴(yán)重災(zāi)害損失的特點(diǎn)。

表2 2009-2018年重慶市侵蝕性降雨事件K均值聚類結(jié)果

3.3.2 各類侵蝕性降雨的時(shí)空變化特征 將K均值聚類的各類侵蝕性降雨對(duì)應(yīng)到年內(nèi)各月份，統(tǒng)計(jì)出各類侵蝕性降雨次數(shù)占相應(yīng)月份降雨事件的比例，統(tǒng)計(jì)結(jié)果見圖7。由圖7可看出，第Ⅰ類雨情主要發(fā)生在5-9月，其中6月發(fā)生比例最大，達(dá)20%；第Ⅱ類雨情各月分布相對(duì)均勻，4-8月發(fā)生比例均超過10%；第Ⅲ類雨情在年內(nèi)分布極為不均，5-9月占全年發(fā)生第Ⅲ類雨情次數(shù)的92.62%，而在12-次年3月期間幾乎不發(fā)生第Ⅲ類侵蝕性降雨。由此可見，5-9月是各類侵蝕性降雨的高發(fā)時(shí)期。結(jié)合水文學(xué)上對(duì)汛期的劃分，研究區(qū)侵蝕性降雨高發(fā)期與汛期同步，也表明汛期的侵蝕性降雨不僅次數(shù)多，且雨情復(fù)雜。

圖7 重慶市各類侵蝕性降雨次數(shù)的年內(nèi)分布比例

各類侵蝕性降雨雨情在研究區(qū)各片區(qū)的發(fā)生頻率統(tǒng)計(jì)結(jié)果見表3。由表3可見，各類侵蝕性降雨在各片區(qū)的分布比例差異較小，分別介于21%～24%、71%～75%和3%～5%之間。

表3 各類侵蝕性降雨雨情在重慶市各片區(qū)的分布比例 %

4 討 論

已有研究表明[19-20, 33]，侵蝕性降雨的雨量、歷時(shí)、雨強(qiáng)等特征因子可對(duì)土壤侵蝕產(chǎn)生不同程度的影響。Tao等[34]認(rèn)為大雨量、中歷時(shí)、大雨強(qiáng)的降雨雨型是喀斯特地區(qū)土壤侵蝕的主要雨型。常松果等[35]基于紫色丘陵區(qū)小流域，觀測(cè)研究得出該區(qū)域主要降雨類型為小雨量、短歷時(shí)、低雨強(qiáng)的降雨，但土壤流失的主要?jiǎng)恿s來自于大雨量、長歷時(shí)、高雨強(qiáng)的雨型。前人研究結(jié)果表明，無論是基于區(qū)域大尺度還是小流域的研究成果，驅(qū)動(dòng)紫色土區(qū)域土壤侵蝕的主要降雨為雨量大、雨強(qiáng)高的侵蝕性降雨，與本次研究結(jié)果具有高度一致性。

史東梅等[36]根據(jù)重慶涪陵區(qū)1952-2001年日降雨數(shù)據(jù)，對(duì)區(qū)域降雨侵蝕力的年際和年內(nèi)變化特征進(jìn)行了分析，本次研究時(shí)段涉及涪陵區(qū)降雨侵蝕力變化特征與史東梅等[36]的結(jié)論總體具有較高的一致性，但也存在局部差異，這說明整體趨勢(shì)尚且一致，但局部差異也體現(xiàn)了降雨隨機(jī)性特征。

5 結(jié) 論

本文以重慶市行政區(qū)劃為研究單元，以逐分鐘降雨觀測(cè)數(shù)據(jù)為基礎(chǔ)，開展侵蝕性降雨的雨量、歷時(shí)、降雨強(qiáng)度及降雨侵蝕力特征值分析，探討研究區(qū)侵蝕性降雨時(shí)空分布及雨情特征，得到以下結(jié)論：

(1)時(shí)間尺度上，年平均雨量、歷時(shí)、雨強(qiáng)及發(fā)生次數(shù)分別為807.0 mm、422 h、71.2 mm/h和25次，分別占多年平均總降雨的75%、44%、55%和18%，年均降雨侵蝕力為1 895.46 MJ·mm/(hm2·h)；各參數(shù)年際變幅劇烈，均呈增加趨勢(shì)。侵蝕性降雨年內(nèi)主要集中于5-9月，其中雨量、歷時(shí)的年內(nèi)變化為雙峰型，雨強(qiáng)、發(fā)生次數(shù)及降雨侵蝕力則為單峰型。

(2)基于反距離權(quán)重插值法得到侵蝕性降雨空間特征：同一區(qū)域侵蝕性雨量大，歷時(shí)長，次數(shù)多，雨強(qiáng)??；同一參數(shù)在年際的空間分布受地形地貌、氣團(tuán)運(yùn)動(dòng)影響表現(xiàn)出空間差異性。

(3)研究區(qū)侵蝕性降雨由主成分聚類分析法分為雨量大、歷時(shí)長、雨強(qiáng)較高的第1類強(qiáng)烈侵蝕區(qū)、第2類輕度侵蝕區(qū)以及第3類中度侵蝕區(qū)。

(4)研究區(qū)3種侵蝕性雨情由K均值聚類法得到，小雨量、短歷時(shí)、低強(qiáng)度、高頻率的第Ⅰ類雨情發(fā)生次數(shù)最多，雨量、歷時(shí)、最大I30及出現(xiàn)次數(shù)中等的第Ⅱ類雨情介于第Ⅰ類和第Ⅲ類之間，雨量大、歷時(shí)長、雨強(qiáng)高、易造成嚴(yán)重災(zāi)害損失的第Ⅲ類雨情發(fā)生次數(shù)較少。