邱美娟, 穆 佳, 郭春明, 王冬妮, 任景全, 曲思邈, 李建平

(吉林省氣象科學(xué)研究所, 吉林 長(zhǎng)春 130062)

1961-2015年吉林省降雨侵蝕力的時(shí)空變化特征

邱美娟, 穆 佳, 郭春明, 王冬妮, 任景全, 曲思邈, 李建平

(吉林省氣象科學(xué)研究所,吉林長(zhǎng)春130062)

[目的] 分析吉林省1961—2015年降雨侵蝕力的時(shí)空變化趨勢(shì)，為該省的農(nóng)業(yè)和生態(tài)保護(hù)、水土保持等工作提供科學(xué)依據(jù)。 [方法] 利用吉林省46個(gè)自動(dòng)氣象站1961—2015年逐日降雨量資料估算吉林省逐氣象站的降雨侵蝕力，并采用相關(guān)系數(shù)、氣候傾向率和反距離空間插值方法分析吉林省降雨侵蝕力的時(shí)空變化趨勢(shì)。 [結(jié)果] 吉林省年平均降雨侵蝕力在空間分布上從集安開始呈向西北和東北逐漸遞減的變化趨勢(shì)，其空間分布特征與年平均降水量的空間分布特征基本一致。時(shí)間分布上與多年平均降水量的時(shí)間分布特征具有高度一致性，在7月達(dá)到峰值。有34.8%的氣象站點(diǎn)降雨侵蝕力呈上升趨勢(shì)，中、西部大部分地區(qū)呈下降趨勢(shì)，東部有1/2以上呈上升趨勢(shì)，但只有長(zhǎng)白站的下降趨勢(shì)通過顯著性檢驗(yàn)。不同地區(qū)各年代平均降雨侵蝕力變化也不一致，具有波動(dòng)性。不同年代各降雨侵蝕力等值線在空間分布上總體變化不大。 [結(jié)論] 吉林省降雨侵蝕力在時(shí)空變化上與降水量一致，不同地區(qū)降雨侵蝕力變化趨勢(shì)不一樣，幾乎沒有通過顯著性檢驗(yàn)。

降雨侵蝕力; 時(shí)空變化; 氣象; 相關(guān)系數(shù); 氣候傾向率

文獻(xiàn)參數(shù)： 邱美娟, 穆佳, 郭春明, 等.1961—2015年吉林省降雨侵蝕力的時(shí)空變化特征[J].水土保持通報(bào)，2017,37(4)：295-301.DOI:10.13961/j.cnki.stbctb.2017.04.050； Qiu Meijuan, Mu Jia, Guo Chunming, et al. Temporal-spatial variation characteristics of rainfall erosivity in Jilin Province during 1961 to 2015[J]. Bulletin of Soil and Water Conservation, 2017,37(4)：295-301.DOI:10.13961/j.cnki.stbctb.2017.04.050

降雨是導(dǎo)致土壤侵蝕的主要?jiǎng)恿σ蛩刂?，在美?guó)通用土壤流失方程USLE[1]以及國(guó)內(nèi)改進(jìn)的一些土壤流失方程[2]，比如中國(guó)坡面土壤流失方程CSLE(chinese soil loss equation)[3]等眾多土壤流失方程中土壤侵蝕的降雨因素用降雨侵蝕力(rainfall erosivity)表示[4]。降雨侵蝕力是指降雨引起土壤侵蝕的潛在能力[5]，降雨侵蝕力已經(jīng)成為土壤侵蝕、產(chǎn)沙和水環(huán)境建模的主要參數(shù)之一[6]。準(zhǔn)確評(píng)估降雨侵蝕力，認(rèn)識(shí)其時(shí)空分布規(guī)律，對(duì)于區(qū)域土壤水蝕的評(píng)估與預(yù)報(bào)以及制定有效的水土保持措施具有重大意義，因此降雨侵蝕力研究已經(jīng)成為一項(xiàng)重要的研究課題。國(guó)內(nèi)外很多學(xué)者對(duì)降雨侵蝕力R值進(jìn)行了大量系統(tǒng)的研究[7-9]。其中美國(guó)學(xué)者Wischmeier等[10]于1965年建立了著名的土壤流失方程(USEL， universal soil loss equation)，并首次提出以次降雨總動(dòng)能E與30 min最大雨強(qiáng)I30的乘積EI30作為降雨侵蝕能力指標(biāo)[11]。國(guó)內(nèi)有更多的研究者[12-16]則用自己的試驗(yàn)資料統(tǒng)計(jì)分析得出了適應(yīng)一定地區(qū)的降雨侵蝕力最佳計(jì)算組合，如賈志軍等[12-13]和王萬忠等[14-15]都認(rèn)為，黃土高原降雨侵蝕力最佳計(jì)算組合是一次暴雨總動(dòng)能E與一次暴雨的最大10 min最大雨強(qiáng)I10的乘積的累加。而周伏建[17]等認(rèn)為適合福建的降雨侵蝕力最佳計(jì)算組合是一次暴雨總能量E與一次暴雨的60 min最大雨強(qiáng)I60的乘積的累加。但是由于計(jì)算EI30，I10或者I60等指標(biāo)需要連續(xù)記錄的降雨過程資料，而我國(guó)常規(guī)氣象觀測(cè)尚不具備這種條件，因此一般利用氣象站常規(guī)降雨統(tǒng)計(jì)資料估算降雨侵蝕力的簡(jiǎn)易方法[18]。目前，國(guó)內(nèi)常用的降雨侵蝕力簡(jiǎn)易算法有基于年降雨量的簡(jiǎn)易算法[18-21]、基于月降雨量的簡(jiǎn)易算法[22]、基于日降雨量的簡(jiǎn)易算法[23-25]以及基于小時(shí)雨量的簡(jiǎn)易算法等[26]4種。但是小時(shí)降雨量數(shù)據(jù)獲取較困難，而年、月雨量信息對(duì)降雨侵蝕力的估算精度有局限性[27]，相比之下，用日雨量數(shù)據(jù)估算降雨侵蝕力精度可靠。國(guó)內(nèi)利用日雨量數(shù)據(jù)估算降雨侵蝕力的研究已有不少研究報(bào)道。如馬良等[28]利用江西省16個(gè)氣象站的降雨資料分析了該省降雨侵蝕力變化的時(shí)空特征。賴成光等[29]根據(jù)珠江流域43個(gè)氣象站1960—2012年逐日降雨資料計(jì)算各站點(diǎn)降雨侵蝕力，對(duì)流域降雨侵蝕力進(jìn)行了時(shí)空變化分析。Xin等[30]利用1956—2008年日降雨資料分析了黃土高原的降雨侵蝕力時(shí)空變化特征，指出該地區(qū)的降雨侵蝕力在過去50 a多呈下降趨勢(shì)，且在黃河中游的河口鎮(zhèn)至龍門段區(qū)域下降尤為明顯；劉斌濤等[31]利用中國(guó)590個(gè)氣象站1960—2009年逐日降雨量資料估算中國(guó)降雨侵蝕力時(shí)空變化特征，計(jì)算結(jié)果發(fā)現(xiàn)青藏高原中東部、東緣、天山山脈以及東南沿海區(qū)域年降雨侵蝕力顯著增加，青海省的諾木洪—都蘭—曲麻萊—伍道梁一帶增加趨勢(shì)最為顯著。以上關(guān)于降雨侵蝕力等的研究均取得顯著的研究成果，但是，關(guān)于吉林省近幾十年來降雨侵蝕力的空間分布特征和時(shí)間變化規(guī)律的詳細(xì)研究尚未見報(bào)道，吉林省降雨侵蝕力變化的時(shí)空差異還不清晰，需要開展相應(yīng)的研究。因此，本文利用吉林省1961—2015年的逐日降雨量資料通過計(jì)算分析吉林省近幾十年來降雨侵蝕力時(shí)空變化趨勢(shì)，以期為政府相關(guān)部門制定水土保持工作方針提供參考。

1 資料與方法

1.1 資料及處理

考慮資料的完整性和連續(xù)性的基礎(chǔ)上，篩選出吉林省46個(gè)自動(dòng)氣象站1961—2015年的逐日降雨量數(shù)據(jù)。利用Fortran計(jì)算機(jī)編程語言計(jì)算各站逐月、逐年、以及多年平均降雨侵蝕力。利用數(shù)理統(tǒng)計(jì)分析方法分析其時(shí)空變化特征，采用地理信息系統(tǒng)軟件ArcGIS對(duì)其進(jìn)行空間表達(dá)。

1.2 降雨侵蝕力計(jì)算

采用第一次全國(guó)水利普查水土保持專項(xiàng)普查使用的降雨侵蝕力算法，該算法由章文波等[24,28]提出降雨侵蝕力估算模型改進(jìn)而來，其計(jì)算公式為：

(1)

α=21.239β-7.3967

(2)

(3)

(4)

(5)

1.3 相關(guān)系數(shù)

相關(guān)系數(shù)計(jì)算公式為：

(6)

1.4 氣候傾向率

將氣象要素的趨勢(shì)變化用一次線性方程表示，即

(7)

2 結(jié)果與分析

2.1 降雨侵蝕力的時(shí)空分布特征

圖1為吉林省1961—2015年多年平均降雨侵蝕力空間分布圖。從圖1可見，吉林省年平均降雨侵蝕力從集安開始呈向西北、向東北逐漸遞減的變化趨勢(shì)。年平均降雨侵蝕力的空間分布特征與年平均降水量的空間分布特征(圖1)基本一致。其中集安地區(qū)降雨侵蝕力最高達(dá)到7 023 MJ·mm/(hm2·h·a)，其次是通化地區(qū)遞減到6 500 MJ·mm/(hm2·h·a)以下，繼續(xù)向西北、東北一直到柳河—輝南東部—樺甸南部—東崗西部一帶，降雨侵蝕力降低到4 500～5 500，4 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線大致與700 mm年均降水量等值線相近；3 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線大致與600 mm年均降水量等值線相近，東部在敦化、二道一帶，西北部在四平—公主嶺—長(zhǎng)春—九臺(tái)—榆樹一帶，兩條等值線幾乎重合；2 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線大致與500 mm年均降水量等值線幾乎重合，在雙遼—長(zhǎng)嶺—前郭—扶余西部一帶。

從吉林省平均及代表站點(diǎn)(洮南、長(zhǎng)春、敦化、通化)1961—2015年多年平均降雨侵蝕力時(shí)間分布情況看(圖2)，降雨侵蝕力在1—3月非常低，有的站點(diǎn)甚至為0，之后隨著月份逐漸升高，到7月達(dá)到最高值，隨后又逐漸降低，11，12月非常低。年平均降雨侵蝕力的時(shí)間分布特征與年平均降水量的時(shí)間分布特征具有高度一致性。

圖1 吉林省1961-2015年年均降雨侵蝕力和年均降雨量空間分布

2.2 降雨侵蝕力的時(shí)空變化趨勢(shì)

2.2.1 年際變化 為分析吉林省1961—2015年55 a來降雨侵蝕力的年際變化特征，計(jì)算了各氣象站的降雨侵蝕力1961—2015年的氣候傾向率并對(duì)其進(jìn)行顯著性檢驗(yàn)，并利用地理信息系統(tǒng)軟件ArcGIS的反距離插值方法對(duì)降雨侵蝕力氣候傾向率和顯著性檢驗(yàn)結(jié)果進(jìn)行空間表達(dá)(圖3—4)，并繪制代表站點(diǎn)年際變化曲線圖(圖5)。

氣候傾向率反映了降雨侵蝕力的變化趨勢(shì)，由圖3可見，近55 a來吉林省中、西部大部分地區(qū)降雨侵蝕力呈下降趨勢(shì)，即氣候傾向率小于0 〔MJ·mm/(hm2·h·a)〕/10 a，其中中部的四平、梨樹、榆樹、遼源和西部的長(zhǎng)嶺、農(nóng)安、通榆等地區(qū)降雨侵蝕力的下降趨勢(shì)相對(duì)最大，氣候傾向率均在-70 〔MJ·mm/(hm2·h·a)〕/10 a以下；其中四平站降雨侵蝕力的氣候傾向率最大，為-135 〔MJ·mm /(hm2·h·a)〕/10 a，其次是梨樹，為-95 〔MJ·mm /(hm2·h·a)/10 a〕；其他地區(qū)降雨侵蝕力的氣候傾向率-64～-4 〔MJ·mm /(hm2·h·a)〕/10 a，且表現(xiàn)為西北部小于東南部地區(qū)(圖3)。

圖2 吉林省及代表站年平均降雨侵蝕力和降水量時(shí)間分布

中部的雙遼、長(zhǎng)春和西部的鎮(zhèn)賚、大安等地區(qū)降雨侵蝕力呈上升趨勢(shì)，氣候傾向率分別為7,8,7,32 〔MJ·mm/(hm2·h·a)〕/10 a。東部地區(qū)有1/2以上站點(diǎn)降雨侵蝕力呈上升趨勢(shì)，其中通化站的上升趨勢(shì)最大，為105 〔MJ·mm/(hm2·h·a)〕/10 a；龍井、靖宇、東崗、吉林城郊和安圖的氣候傾向率在45～70 〔MJ·mm/(hm2·h·a)〕/10 a；臨江、煙筒山、磐石、和龍、延吉和樺甸等站的氣候傾向率在0.08～25 〔MJ·mm/(hm2·h·a)〕/10 a。東部其它地區(qū)降雨侵蝕力均呈下降趨勢(shì)，除了長(zhǎng)白站下降趨勢(shì)最大，為-165 〔MJ·mm /(hm2·h·a)〕/10 a，其他各站的氣候傾向率在-71～-18 〔MJ·mm/(hm2·h·a)〕/10 a。

圖3 吉林省年降雨侵蝕力氣候傾向率空間分布

在研究的46個(gè)氣象站中，有16個(gè)氣象站年降雨侵蝕力呈上升趨勢(shì)，占34.8%，呈上升趨勢(shì)的氣象站主要集中的東部地區(qū)。對(duì)各站降雨侵蝕力的變化趨勢(shì)進(jìn)行顯著性檢驗(yàn)發(fā)現(xiàn)(圖4)，吉林省大部降雨侵蝕力的上升或者下降趨勢(shì)幾乎均未通過任何水平的顯著性檢驗(yàn)，只有長(zhǎng)白站的下降趨勢(shì)通過0.05水平(R=0.265)的顯著性檢驗(yàn)。

圖4 吉林省年降雨侵蝕力顯著性檢驗(yàn)相關(guān)系數(shù)值空間分布

圖5 吉林省不同代表站點(diǎn)年降雨侵蝕力年際變化

2.2.2 年代際變化 各個(gè)代表站以及吉林省平均20世紀(jì)60，70，80，90，2000年代以及最近5 a(2011—2015年)降雨侵蝕力變化特征詳見表1，圖5—6。從圖5和表1可以看出，吉林省平均年降雨侵蝕力呈波動(dòng)變化趨勢(shì)，其中80年代最高，90年代最低，最高最低相差7.9%。洮南站近5 a最高，2000年代最低，最高最低相差65.8%，差異較大。長(zhǎng)春站80年代最高，90年代最低，最高最低相差18.3%。敦化站80年代最高，70年代最低，最高最低相差10.6%。通化站2000年代最高，70年代最低，最高最低相差21.9%。因此，不同氣象站點(diǎn)各年代平均降雨侵蝕力變化不一致，有一定的差異性。

表1 吉林省平均及代表站點(diǎn)降雨侵蝕力年代際變化 MJ·mm/(hm2·h·a)

70年代與60年代相比，多數(shù)呈下降趨勢(shì)，變化最大的是琿春，相對(duì)下降了22.5%，其次是雙陽，相對(duì)上升了16.1%；而鎮(zhèn)賚、長(zhǎng)春和靖宇變化在±1%以內(nèi)，幾乎無變化；其他均變化-17.0%～11.9%。80年代與70年代相比，大多數(shù)呈上升趨勢(shì)，洮南、鎮(zhèn)賚、大安、乾安、舒蘭、蛟河、汪清、和龍、琿春和延吉上升趨勢(shì)相對(duì)較大，在20.1%～28.5%之間，扶余幾乎無變化在±1%以內(nèi)，其他均變化-5.3%～18.5%。90年代與80年代相比，大多數(shù)呈下降趨勢(shì)，農(nóng)安、德惠、公主嶺、蛟河和遼源下降趨勢(shì)相對(duì)較大，相對(duì)下降21.1%～29.3%；龍井上升趨勢(shì)最大，相對(duì)上升23.4%；白城幾乎無變化；其它均變化-18.7%～9.5%。2000年代與90年代相比，差不多一半的站點(diǎn)呈下降趨勢(shì)，白城、洮南、鎮(zhèn)賚、乾安和長(zhǎng)嶺下降趨勢(shì)較大，相對(duì)下降23.0%～40.4%；煙筒山和遼源上升趨勢(shì)較大，分別相對(duì)上升23.3%和22.8%；前郭、伊通、安圖和磐石幾乎無變化；其他變化在-16.2%～17.7%之間。近5 a與2000年代相比，中、西部多數(shù)呈上升趨勢(shì)，而東部多數(shù)呈下降趨勢(shì)；安圖和梅河口下降趨勢(shì)最大，分別相對(duì)下降23.8%和25.9%；四平、伊通和東豐幾乎無變化；而白城、洮南、鎮(zhèn)賚、乾安、前郭、通榆和汪清上升趨勢(shì)較大，相對(duì)上升21.4%～76.1%，其他變化-16.1%～19.6%(表略)。

圖6 吉林省年均降雨侵蝕力〔MJ·mm/(hm2·h·a)〕年代際變化空間分布

從吉林省年均降雨侵蝕力年代際變化空間分布上看，60年代2 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線在雙遼—長(zhǎng)嶺—前郭—扶余一帶，另外還有和龍、龍井延吉零星站點(diǎn)；3 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線在梨樹—伊通—雙陽—吉林城郊—榆樹一帶，以及東部的二道和敦化一帶；4 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線在柳河—輝南南部—樺甸南部—靖宇—東崗西部—臨江一帶；5 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線在通化站以北；6 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線在通化集安兩站。70年代與60年代相比，2 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線中西部幾乎無變化，而東部擴(kuò)大到汪清地區(qū)；3 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線略微向西北部移動(dòng)；4 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線向南移動(dòng)到靖宇南部；5 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線略微向南移動(dòng)；6 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線已經(jīng)移出吉林省地界。80年代與70年代相比2 500和3 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線變化非常??；4 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線有向西、向北移動(dòng)到梅河口—輝南樺甸—東崗一帶；5 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線幾乎無變化；6 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線又移動(dòng)至集安地區(qū)。90年代與80年代相比，2 500，3 500和4 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線均向東南移動(dòng)1～2個(gè)站點(diǎn)左右；而6 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線移動(dòng)到通化站地區(qū)。2000年代與90年代相比，各等值線又向西北移動(dòng)，但移動(dòng)幅度較小，未達(dá)到80年代的位置；6 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線擴(kuò)大到通化和集安兩地。近5 a與2000年代相比，4 500和6 500 MJ·mm/(hm2·h·a)年均降雨侵蝕力等值線向南移動(dòng)了1個(gè)站點(diǎn)左右，其他等值線移動(dòng)幅度很小。

3 討論與結(jié)論

(1) 吉林省年平均降雨侵蝕力在空間分布上從集安開始呈向西北、向東北逐漸遞減的變化趨勢(shì)，其空間分布特征與年平均降水量的空間分布特征基本一致。集安地區(qū)降雨侵蝕力最高達(dá)到7 023 MJ·mm/(hm2·h·a)。

(2) 從吉林省平均及代表站點(diǎn)1961—2015年多年平均降雨侵蝕力時(shí)間分布情況看，降雨侵蝕力隨著月份逐漸升高，到7月達(dá)到峰值，隨后又逐漸降低。年平均降雨侵蝕力的時(shí)間分布特征與年平均降水量的時(shí)間分布特征具有高度一致性。

(3) 在研究的46個(gè)氣象站中，有16個(gè)氣象站年降雨侵蝕力呈上升趨勢(shì)，占34.8%。中、西部大部分地區(qū)降雨侵蝕力呈下降趨勢(shì)，個(gè)別站降雨侵蝕力呈上升趨勢(shì)；東部地區(qū)有1/2以上站點(diǎn)降雨侵蝕力呈上升趨勢(shì)，其中通化站的上升趨勢(shì)最大，為105 MJ·mm/(hm2·h·a)/10 a，長(zhǎng)白站下降趨勢(shì)最大，為-165 MJ·mm/(hm2·h·a)/10 a。只有長(zhǎng)白站的下降趨勢(shì)通過了0.05水平的顯著性檢驗(yàn)。

(4) 吉林省平均年降雨侵蝕力呈波動(dòng)變化趨勢(shì)，其中80年代最高，90年代最低，最高最低相差7.9%。不同氣象站點(diǎn)各年代平均降雨侵蝕力變化也不一致，有一定的差異性。70年代與60年代相比，多數(shù)呈下降趨勢(shì)；80年代與70年代相比，大多數(shù)呈上升趨勢(shì)；90年代與80年代相比，大多數(shù)呈下降趨勢(shì)；2000年代與90年代相比，差不多1/2的站點(diǎn)呈下降趨勢(shì)；近5 a與2000年代相比，中、西部多數(shù)呈上升趨勢(shì)，而東部多數(shù)呈下降趨勢(shì)。不同年代各降雨侵蝕力等值線在空間分布上總體變化不大。

[1] 于興修,楊桂山.通用水土流失方程因子定量研究進(jìn)展與展望[J].自然災(zāi)害學(xué)報(bào),2003,12(3):14-18.

[2] 胡續(xù)禮,潘劍君,楊樹江,等.幾種降雨侵蝕力模型的比較研究[J].水土保持通報(bào),2006,26(1):68-70.

[3] 鄭粉莉,王占禮,楊勤科.我國(guó)水蝕預(yù)報(bào)模型研究的現(xiàn)狀、挑戰(zhàn)與任務(wù)[J].中國(guó)水土保持科學(xué),2005,3(1):7-14

[4] 沈照偉,田剛,李鋼.浙江省降雨侵蝕力變化特征分析[J].水土保持通報(bào),2013,33(4):119-124.

[5] 章文波,謝云,劉寶元.降雨侵蝕力研究進(jìn)展[J].水土保持學(xué)報(bào),2002,16(2):43-46.

[6] 王愛娟,李智廣,劉峰.長(zhǎng)江上游水蝕區(qū)降雨侵蝕力的時(shí)空分布特征[J].水土保持通報(bào),2013,33(1):8-11.

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

[8] 段文明,穆興民,王飛,等.嘉陵江流域降雨侵蝕力時(shí)空變化分析[J].水土保持通報(bào),2012,32(5):182-185.

[9] 許月卿,周巧富,李雙成.貴州省降雨侵蝕力時(shí)空分布規(guī)律分析[J].水土保持通報(bào),2005,25(4):11-14.

[10] Wischmeier W H, Smith D D. A Universal Soil-loss Equation to Guide Conservation Farm Planning[C]. US：Transactions 7th Interenational Congress Soil Science, 1960.

[11] 章文波,謝元,劉寶元.降雨侵蝕力研究進(jìn)展[J].水土保持學(xué)報(bào),2002,16(2):43-46.

[12] 賈志軍,王小平,李俊義.晉西黃土丘陵溝壑區(qū)降雨侵蝕力指標(biāo)R值的確定[J].中國(guó)水土保持,1987(6):18-20.

[13] 賈志軍,王小平,李俊義,等.晉西黃土高原降雨侵蝕力研究[J].中國(guó)水土保持,1991(1):43-46.

[14] 王萬忠.黃土地區(qū)降雨特性與土壤流失關(guān)系的研究[J].水土保持通報(bào),1983,3(5):62-64.

[15] 王萬忠.黃土地區(qū)降雨侵蝕力R指標(biāo)的研究[J].中國(guó)水土保持,1987(12):34-38.

[16] 范瑞瑜.黃河中游地區(qū)小流域土壤流失量計(jì)算方程的研究[J].中國(guó)水土保持,1985(2):12-18.

[17] 周伏建,明華,林福興,等.福建省降雨侵蝕力指標(biāo)的初步探討[J].福建水土保持,1989(2):58-60.

[18] 伍育鵬,謝云,章文波.國(guó)內(nèi)外降雨侵蝕力建議計(jì)算方法的比較[J].水土保持學(xué)報(bào),2001,15(3):31-34.

[19] 繆馳遠(yuǎn),徐霞,魏欣,等.重慶市主城區(qū)降雨侵蝕力特征分析[J].資源科學(xué),2007,29(4):55-60.

[20] Lee J H, Heo J H. Evaluation of estimation methods for rainfall erosivity based on annual precipitation in Korea[J]. Journal of Hydrology, 2011, 409(1/2):30-48.

[21] Renard K G, Freimund J R. Using monthly precipitation data to estimate the R-factor in the revised USLE[J]. Journal of Hydrology, 1994, 157(1/4):287-306.

[22] Richardson C W, Forster G R, Wright D A. Estimation of erosion index from daily rainfall amount[J]. Transaction of the ASAE, 1983,26(1):153-156.

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

[24] Angulo M M, Begueria S. Estimating rainfall erosivity from daily precipitation records:A comparison among methods using data from the Ebro Basin(NE Spain)[J]. Journal of Hydrology, 2009,379(1/2):111-121.

[25] 楊軒,梁音,方繼青,等.基于日降雨信息的月降雨侵蝕力模型[J].土壤學(xué)報(bào),2010,47(2):216-222.

[26] 殷水清,謝云,王春剛.用小時(shí)降雨資料估算降雨侵蝕力的方法[J].地理研究,2007,26(3):541-547.

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

[28] 馬良,姜廣輝,左長(zhǎng)清,等.江西省50余年來降雨侵蝕力變化的時(shí)空分布特征[J].農(nóng)業(yè)工程學(xué)報(bào),2009,25(10):61-68.

[29] 賴成光,陳曉宏,王兆禮,等.珠江流域1960—2012年降雨侵蝕力時(shí)空變化特征[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(8):159-167.

[30] Xin Zhongbao, Yu Xinxiao, Li Qingyun, et al. Spatiotemporal variation in rainfall erosivity on the Chinese Loess Plateau during the period 1956—2008[J]. Regional Environmental Change, 2011, 11(1):149-159.

[31] 劉斌濤,陶和平,宋春風(fēng),等.1960—2009年中國(guó)降雨侵蝕力的時(shí)空變化趨勢(shì)[J].地理研究,2013,32(2):245-256.

Temporal-spatial Variation Characteristics of Rainfall Erosivity in Jilin Province During1961to2015

QIU Meijuan, MU Jia, GUO Chunming, WANG Dongni, REN Jingquan, QU Simiao, LI Jianping

(Institute of Meteorological Sciences of Jilin Province, Changchun, Jilin 130062, China)

[Objective] The temporal and spatial variation trend of rainfall erosivity in Jilin Province from 1961 to 2015 were analyzed to provide scientific bases for agricultural and ecological protection, as well as soil and water conservation. [Methods] Based on the daily rainfall data from 46 meteorological stations in Jilin Province from 1961 to 2015, rainfall erosivity was estimated and the temporal and spatial distribution was analyzed with correlation coefficient, climate tendency rate and the inverse distance interpolation methods, respectively. [Results] The spatial distribution of rainfall erosivity has a decreasing trend from the station of Jian to the northwest and northeast region. The spatial distribution of average annual rainfall erosivity is similar to that of average annual precipitation. The temporal distribution of yearly averaged rainfall erosivity is similar to that of multi-year average precipitation, which has a peak in July. 34.8% meteorological stations have incline rainfall erosivity. For most of the western region, there is a downward trend; while for more than half of the stations in eastern region is on the rise, among with Changbai station was significant. The changes of average rainfall erosivity of each decade in different regions are not consistent, and have fluctuation. The spatial distribution of rainfall erosivity contour in different decade differ slightly. [Conclusion] Rainfall erosivity in Jilin Province is consistent with the precipitation in temporal and spatial variation, and the trend of rainfall erosivity in different regions is different, and almost dose not pass the significance test.

rainfallerosivity;temporalandspatialvariations;meteorology;correlationcoefficient;climatetendencyrate

A

： 1000-288X(2017)04-0295-07

： S157.1

2016-06-02

：2016-06-16

吉林省任務(wù)下達(dá)式課題“吉林省水稻氣候品質(zhì)評(píng)價(jià)指標(biāo)及應(yīng)用”(2015020); 氣象行業(yè)專項(xiàng)(GYHY201206018); 中國(guó)氣象局沈陽大氣環(huán)境研究所2016年開放基金課題(2016SYIAE02)

邱美娟(1987—),女(漢族),遼寧省大連市人,碩士,工程師,主要從事農(nóng)業(yè)氣象災(zāi)害研究。E-mail:qmjcams@163.com。