1961—2019年江西省降雨侵蝕力時空變化分析

鄭瀟　白淑英　高吉喜

摘要 降雨侵蝕力（R）與降雨量、降雨歷時、降雨強度、降雨動能有關，反映了降雨特性對土壤侵蝕的影響，是影響土壤侵蝕的重要因素之一．本文利用1961—2019年江西省25個氣象站點的逐日降雨量數據，基于降雨侵蝕力模型，通過Mann-Kendall相關檢驗、小波分析和Kriging空間插值方法，分析江西省降雨侵蝕力空間分布及變化趨勢．結果表明：江西省年均降雨量和降雨侵蝕力從贛南到贛北逐漸增加，同時降雨空間分布和降雨侵蝕力空間分布具有相似性;春季和夏季降雨侵蝕力在3 000～6 000 MJ·mm·hm-2·h-1·a-1，且春季降雨侵蝕力最大值高于夏季，而秋冬季降雨侵蝕力明顯小于春夏季;江西省年降雨侵蝕力平均最大值分布在贛北區域，其次為贛中，贛南最小，侵蝕力由北向南逐漸遞減．研究表明江西省因降雨導致侵蝕力不斷增加，尤其以春夏季更為劇烈，其空間分布存在明顯差異．關鍵詞 降雨侵蝕力;降雨量;Mann-Kendall相關檢驗;小波分析

中圖分類號S157.1

文獻標志碼A

0 引言

降雨是引起土壤侵蝕的動力因子之一，降雨侵蝕力常作為衡量降雨導致土壤侵蝕的潛在能力［1］，其時空特征研究對于區域水土保持規劃、災害控制［2］以及生態敏感性指標評價具有重要指導意義［3］．當前全球氣候變化大背景下，極端氣候事件發生的概率明顯增多［4］，在某種程度上增加了土壤侵蝕的風險［5］，帶來的風險包括洪澇災害［6］等造成重大損失．降雨侵蝕力R值作為表征降雨對土壤潛在侵蝕能力的特征值，自20世紀50年代中后期提出并應用于通用土壤流失方程（USLE）以來，逐漸成為目前研究土壤水蝕對降雨變化響應的最佳選擇［7］．分析降雨侵蝕力時空變化規律，以及對后續土壤侵蝕力評估和水土整治等過程有利于采取針對性措施制止或減緩降雨侵蝕．

降雨侵蝕力的計算是由Wischmeier等［8］提出的，它利用次降雨動能和降雨強度的乘積來計算，其優勢在于計算精度高，但數據要求也高，需要長序列連續性降雨，處理過程繁瑣、耗時．國外學者將降雨侵蝕與農作物、徑流等結合并分析降雨侵蝕造成的影響．Johannsen 等［9］通過降雨侵蝕力研究奧地利不同農業生產區對土壤侵蝕的影響;Anderson等［10］研究降雨徑流對溝壑的侵蝕過程，并得出溝壑、降雨和植被之間的相互關系;Azari等［11］針對氣候變化對河流土壤資源和泥沙負荷的關系，利用回歸方程得出降雨量增加會導致泥沙負荷提高的結論．國內學者利用日、月、年降雨量計算降雨侵蝕力模型［12-14］，且計算方式得到廣泛應用［15］．溫天福等［16］利用駐日降雨資料采用TFPW-MK法分析流域降雨侵蝕變化規律，指出侵蝕力增加會引發水土流失風險加劇;周琳等［17］基于日降雨侵蝕模型分析了降雨侵蝕與山地、平原之間的相關性，以及降雨侵蝕與四季變化之間大小關系;索笑穎等［18］根據時間和空間侵蝕突變點分析了侵蝕季節變化和地形變化;曾瑜等［19］基于樹立統計及空間插值技術研究了降雨侵蝕力分配特征，給出了降雨侵蝕力整體分布格局．

江西屬于典型的中部平原、四周環山區域，贛南地區多以山地和丘陵為主［20］，并且全省分布較廣的第四紀紅壤，土壤質地黏重、透水性差極易造成水力侵蝕［21］，是全國水土流失最嚴重的省份之一．本文基于日降雨氣象資料，利用降雨侵蝕模型，結合克里金（Kriging）插值、Mann-Kendall突變性檢驗、小波分析等方法研究1961—2019年江西省降雨侵蝕力時空演替特征及變化規律，以期為后期水土流失監測與防治、環境保護等提供技術支撐．

1 研究區域概況

江西省位于中國東南部，在長江中下游南岸，地處中亞熱帶，季風氣候顯著，四季化分明．江西省的主要河流為贛江、撫河、信江、饒河和修水，流域面積達到16.22萬km²，這5大水系又以鄱陽湖為中心，經鄱陽湖集納后流向長江．鄱陽湖是中國最大的淡水湖，贛江流域面積79 173 km²，是長江最重要的支流之一．江西境內水熱條件差異較大，多年平均氣溫自北向南依次增高，南北溫差約3 ℃;年降水量1 300～2 100 mm，地區分布上是北多南少，東多西少，且降雨的年際變化很大，空間分布差異十分明顯，多雨與少雨年份相差也較大，大約50％的蒸發發生在7月中旬到11月，氣溫與降雨的變化是導致江西旱澇頻發的原因之一．

2 數據與研究方法

2.1 數據來源

降雨量是表征氣候特征的最基本要素．江西省1961—2019年25個站點的降雨量數據來源于中國氣象數據網（http：／／data.cam.cn），站點分布均勻（圖1），數據可靠．利用Kriging空間插值方法可以推算出該區域氣候要素情況，該方法在我國降雨量研究上有很好的擬合優度［22］．本文采用Kriging空間插值方法生成年均降雨量和降雨侵蝕力柵格數據．

2.2 降雨侵蝕力（R）

2.3 Mann-Kendall相關檢驗

2.4 小波分析

3 結果分析

3.1 降雨量與降雨侵蝕力時空分析

3.1.1 江西省降雨量和降雨侵蝕力年際分析

江西省的降雨侵蝕力和降雨量空間分布變化與該地區的地形地貌有很大關系，高程對氣象要素的影響是主要的，復雜山地環境則會影響氣象要素的地形空間尺度和分布格局［28］．贛南多數為丘陵地帶，贛中、北平原地帶相較多［29］，廬山海拔較高，加上贛江和鄱陽湖分布廬山四周，濕潤氣流受到地形阻擋，易形成長時間降雨．因此，廬山降雨侵蝕力比其他地區大，從贛北向贛南呈依次遞減的空間分布特征．

3.1.2 江西省降雨侵蝕力季節性分析

降雨分布較為均勻，且春季降雨侵蝕力最大值高于夏季，但夏季降雨侵蝕力最小值高于春季降雨侵蝕力最小值．秋季（圖9）和冬季（圖10）降雨侵蝕力明顯小于春夏季降雨侵蝕力，秋季降雨主要分布在江西西南部（即井岡山和遂川一帶），大部分地區降雨侵蝕力較為平穩，冬季降雨分布在贛北地區，且由北向南逐漸遞減．贛北地區降雨量較大，主要原因是受到北方冷空氣向南方偏移的影響．

3.2 降雨侵蝕力檢驗分析

本文運用Mann-Kendall突變檢驗方法對江西省58 a來降雨侵蝕力年際變化趨勢進行突變性檢驗，設上下顯著性水平為0.05，紅色實線代表UF統計量，藍色虛線代表UB統計量．通過圖11可知：1970年后UF大于0，表明降雨隨年際變化有增加趨勢，且在1996—2006年超過0.05的顯著性水平檢驗;1970—1990年存在多次震蕩相交，1990年之后逐漸平緩，表明江西省的降雨侵蝕力在1970—1990年發生顯著突變．根據Mann-Kendall變異系數和Z值空間插值分析降雨侵蝕力分布情況（圖12），各站點變異范圍在0.25～0.38之間，各地區變化差異比較明顯，以廬山為代表的贛北地區，以德興、上饒和玉山為代表的贛東北地區和以寧都為代表的贛南地區變異系數相對較大，而贛西南和西北地區變異系數相對較小，主要原因是贛北地區降雨變化不穩定，而贛南和贛西南地區降雨相對穩定，且變異系數基本小于0.3．由降雨侵蝕力統計量Z值得出，除贛州、廣昌、波陽和寧都的降雨侵蝕力表現出降低的趨勢外，其余各站點都表現出增長趨勢．但是，無論降低還是增加的變化，統計量|Z|<1.96，均未通過顯著性檢驗，侵蝕力均表現為不顯著的變化趨勢．

圖13為1961—2019年江西地區降雨侵蝕力小波分析［30］結果，圖中紅色代表降雨強度變化顯著，越接近藍色，則表示降雨強度變化穩定．由圖13可見，年際降雨量在時間尺度上存在復雜的周期結構變化，在降雨演變的過程中，降雨強度也存在5～10 a、3～6 a的周期變化規律，在1970年后生成6 a左右的震蕩周期，其中在5～10 a中存在4次震蕩，3～6 a中存在5次震蕩．小波系數為正值時，表明降雨強度提高;小波系數為負值時，表明降雨強度明顯減弱．1970—1975年小波系數逐漸增加，降雨強度增大;1975—1980年小波系數逐漸減小，降雨強度逐漸下降．

4 討論

降雨對于降雨侵蝕力而言是主要的因素，江西省降雨量與降雨侵蝕力呈現正相關，和夏正兵［31］對黃河流域降雨侵蝕力研究、伊力哈木·伊馬木［32］對新疆維吾爾自治區降雨侵蝕力研究結論一致．因此，針對不同降雨資料選擇不同的降雨侵蝕力模型，可以得出大致相同的降雨侵蝕力結果．

江西的土壤類型以紅壤為主，再加上江西降雨分布不均，超過70％的降雨量發生在4—9月之間［33］，但降雨侵蝕力較強地區分布在贛北一帶，其次是贛中，最后是贛南．產生的原因有：1）不同地區所處的季風帶不相同;2）地區高程的大小影響季風的走勢，進而造成降雨量的不同．贛北地區降雨侵蝕力在春夏季出現峰值，原因在于地處亞熱帶季風氣候，受大陸季風影響，導致春夏季降雨侵蝕力增強．另外，并非所有的降雨都會帶來侵蝕影響，尤其是生態脆弱性區域和氣候影響顯著地區．

5 結論

本文以江西省為研究對象，從時間和空間上分析了降雨侵蝕力的變化趨勢，得到以下結論：

1） 江西省年降雨侵蝕力空間分布趨勢為由贛北向贛南逐漸遞減，年平均降雨量與年平均降雨侵蝕力分布基本一致．隨著降雨量的增加，降雨侵蝕力逐漸增大;隨著地貌地形的變化，降雨侵蝕力也有相應改變．

2） 春季和夏季（3—8月）降雨侵蝕力波動較大，侵蝕力較高地區主要分布在贛北及贛東，且春季降雨侵蝕力的最大值高于夏季降雨侵蝕力的最大值，夏季降雨侵蝕力波動較小，降雨侵蝕分布比夏季降雨侵蝕力均勻;秋冬季降雨侵蝕力明顯小于春夏季降雨侵蝕力，秋季降雨主要分布在江西西南部，大部分地區降雨侵蝕力較為平穩，而冬季降雨分布在贛北地區，主要原因是受到北方冷空氣的影響．

3） 根據Mann-Kendall突變檢驗、變異系數和Z值以及小波分析表明，1970年后降雨雖然有增加趨勢，但是變化不顯著，且變化范圍在0.25～0.38之間．通過小波分析可得在1990—1995年UF與UB變化方向一致，且UF在1995年超過0.05顯著水平．空間分布上，吉安、樟樹站點Z值平均超過1，其余站點小于1，說明吉安和樟樹區域年降雨變化較顯著，且降雨侵蝕分布贛北>贛中>贛南．因此，江西省應該進一步加強對贛北區域進行水土防治工作，減少降雨引起的水土流失危害．

參考文獻 References

［1］殷水清，薛筱嬋，岳天雨，等.中國降雨侵蝕力的時空分布及重現期研究［J］.農業工程學報，2019，35（9）：105-113YIN Shuiqing，XUE Xiaochan，YUE Tianyu，et al.Spatiotemporal distribution and return period of rainfall erosivity in China［J］.Transactions of the Chinese Society of Agricultural Engineering，2019，35（9）：105-113

［2］ 辜倩，馬旭林，何寧，等.湖南1961—2010年降雨侵蝕力時空變化特征［J］.科學技術與工程，2018，18（6）：184-189GU Qian，MA Xulin，HE Ning，et al.Temporal and spatial variation characteristics of rainfall erosion force in Hunan province from 1961 to 2010［J］.Science Technology and Engineering，2018，18（6）：184-189

［3］ 劉軍會，鄒長新，高吉喜，等.中國生態環境脆弱區范圍界定［J］.生物多樣性，2015，23（6）：725-732LIU Junhui，ZOU Changxin，GAO Jixi，et al.Location determination of ecologically vulnerable regions in China［J］.Biodiversity Science，2015，23（6）：725-732

［4］ IPCC.Climate change 2013：the physical science basis［R］.Cambridge：Cambridge University Press，2014

［5］ 朱大運，楊倩，陳海，等.1960—2017年貴州省不同水系降雨侵蝕力時空變異特征［J］.水土保持通報，2021，41（1）：6-14ZHU Dayun，YANG Qian，CHEN Hai，et al.Spatiotemporal variations in rainfall erosivity of different river systems in Guizhou province during 1960-2017［J］.Bulletin of Soil and Water Conservation，2021，41（1）：6-14

［6］ 歐陽千林，占承德.近30年江西省洪災損失等級評估及對降雨響應研究［J］.江西水利科技，2020，46（5）：325-331OUYANG Qianlin，ZHAN Chengde.Research on flood disaster losses grade and response to rainfall in Jiangxi province over 30 years［J］.Jiangxi Hydraulic Science ＆ Technology，2020，46（5）：325-331

［7］ 馬良，姜廣輝，左長清，等.江西省50余年來降雨侵蝕力變化的時空分布特征［J］.農業工程學報，2009，25（10）：61-68MA Liang，JIANG Guanghui，ZUO Changqing，et al.Spatial and temporal distribution characteristics of rainfall erosivity changes in Jiangxi province over more than 50 years［J］.Transactions of the Chinese Society of Agricultural Engineering，2009，25（10）：61-68

［8］ Wischmeier W H，Smith D D.Predicting rainfall erosion losses：a guide to conservation planning［M］／／Agriculture Handbook No.537.Washington，DC：U.S. Department of Agriculture，1978：285-291

［9］ Johannsen L L，Schmaltz E M，Mitrovits O，et al.An update of the spatial and temporal variability of rainfall erosivity （R-factor） for the main agricultural production zones of Austria［J］.Catena，2022，215：106305

［10］ Anderson R L，Rowntree K M，le Roux J J.An interrogation of research on the influence of rainfall on gully erosion［J］.Catena，2021，206：105482

［11］ Azari M，Oliaye A，Nearing M A.Expected climate change impacts on rainfall erosivity over Iran based on CMIP5 climate models［J］.Journal of Hydrology，2021，593：125826

［12］ 邵祎婷，何毅，穆興民，等.秦巴山區降雨侵蝕力時空變化特征［J］.長江流域資源與環境，2019，28（2）：416-425SHAO Yiting，HE Yi，MU Xingmin，et al.Spatiotemporal variation of rainfall erosivity in Qin-Ba Mountains region［J］.Resources and Environment in the Yangtze Basin，2019，28（2）：416-425

［13］ 卜兆宏，宮世俊，阮伏水，等.降雨侵蝕力因子的算法及其在土壤流失量監測中的選用［J］.遙感技術與應用，1992，7（3）：1-10BU Zhaohong，GONG Shijun，RUAN Fushui，et al.Preliminary study on the algorithms of erosivity factor of rainfall and their application in monitoring soil loss［J］.Remote Sensing Technology and Application，1992，7（3）：1-10

［14］ 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

［15］ Duan X W，Bai Z W，Rong L，et al.Investigation method for regional soil erosion based on the Chinese soil loss equation and high-resolution spatial data：case study on the mountainous Yunnan province，China［J］.Catena，2020，184：104237

［16］ 溫天福，張范平，胡建民，等.1958—2014年贛江流域降雨侵蝕力時空變化特征［J］.水土保持研究，2020，27（1）：7-14，20WEN Tianfu，ZHANG Fanping，HU Jianmin，et al.Spatiotemporal variation of rainfall erosivity during the period 1958-2014 in Ganjiang River Basin［J］.Research of Soil and Water Conservation，2020，27（1）：7-14，20

［17］ 周琳，許武成，尹義星.近61年四川省降雨侵蝕力的時空變化趨勢［J］.水土保持研究，2018，25（4）：8-14ZHOU Lin，XU Wucheng，YIN Yixing.Temporal and spatial variations of rainfall erosivity in Sichuan province in recent 61 years［J］.Research of Soil and Water Conservation，2018，25（4）：8-14

［18］ 索笑穎，劉玉春，趙光耀，等.河北省山區降雨侵蝕力的時空變化特征［J］.水土保持通報，2020，40（6）：259-265SUO Xiaoying，LIU Yuchun，ZHAO Guangyao，et al.Temporal and spatial variation characteristics of rainfall erosion in mountainous areas of Hebei province［J］.Bulletin of Soil and Water Conservation，2020，40（6）：259-265

［19］ 曾瑜，厲莎，胡煜彬.1961—2014年鄱陽湖流域降雨侵蝕力時空變化特征［J］.生態與農村環境學報，2019，35（1）：106-114ZENG Yu，LI Sha，HU Yubin.Spatial-temporal variation of distribution of rainfall erosivity in the Poyang Lake Basin during 1961-2014［J］.Journal of Ecology and Rural Environment，2019，35（1）：106-114

［20］ Shao C J，Liu Y，Lan H X，et al.Spatiotemporal distribution characteristics，causes，and prevention advice of fatal geohazards in Jiangxi province，China［J］.Science of the Total Environment，2022，834：155337

［21］ 齊述華，蔣梅鑫，于秀波.基于遙感和ULSE模型評價1995—2005年江西土壤侵蝕［J］.中國環境科學，2011，31（7）：1197-1203QI Shuhua，JIANG Meixin，YU Xiubo.Evaluating soil erosion in Jiangxi province with USLE model and remote sensing technology during 1995-2005［J］.China Environmental Science，2011，31（7）：1197-1203

［22］ 馬軒龍，李春娥，陳全功.基于GIS的氣象要素空間插值方法研究［J］.草業科學，2008，25（11）：13-19MA Xuanlong，LI Chune，CHEN Quangong.Study on the method of GIS based spatial interpolation of climate factors in China［J］.Pratacultural Science，2008，25（11）：13-19

［23］ 于延勝，陳興偉.基于Mann-Kendall法的水文序列趨勢成分比重研究［J］.自然資源學報，2011，26（9）：1585-1591YU Yansheng，CHEN Xingwei.Study on the percentage of trend component in a hydrological time series based on Mann-Kendall method［J］.Journal of Natural Resources，2011，26（9）：1585-1591

［24］ Gülü Y S.Improved visualization for trend analysis by comparing with classical Mann-Kendall test and ITA［J］.Journal of Hydrology，2020，584：124674

［25］ 陳曉瑜，孫曉航，黃奇曉，等.1970—2017年福建省降雨侵蝕力的時空變化特征［J］.福建農林大學學報（自然科學版），2021，50（1）：125-133CHEN Xiaoyu，SUN Xiaohang，HUANG Qixiao，et al.Spatial and temporal variation of rainfall erosivity in Fujian province during 1970-2017［J］.Journal of Fujian Agriculture and Forestry University （Natural Science Edition），2021，50（1）：125-133

［26］ 韓知明，賈克力，楊芳，等.基于Morlet小波的呼倫湖流域降水多時間尺度分析［J］.水土保持研究，2018，25（1）：160-166HAN Zhiming，JIA Keli，YANG Fang，et al.Analysis on multi time scale yearly precipitation of Hulun Lake watershed based on wavelet Morlet analysis［J］.Research of Soil and Water Conservation，2018，25（1）：160-166

［27］ 紀小敏，劉煒偉，聶青，等.基于Morlet小波分析的長江大通站徑流年際變化周期研究［J］.江蘇水利，2016（6）：26-29，34JI Xiaomin，LIU Weiwei，NIE Qing，et al.Research on runoff annual variation cycle at Datong hydrological station in Yangtze River based on Morlet wavelet analysis［J］.Jiangsu Water Resources，2016（6）：26-29，34

［28］ 蔣育昊，劉鵬舉，夏智武，等.PRISM模型在復雜地形月降雨空間插值中的可行性研究［J］.水土保持研究，2018，25（1）：57-61，71JIANG Yuhao，LIU Pengju，XIA Zhiwu，et al.Study on the feasibility of PRISM in spatial interpolation of rainfall in complex terrain［J］.Research of Soil and Water Conservation，2018，25（1）：57-61，71

［29］ 何紹浪，何小武，李鳳英，等.近60年來江西省各等級侵蝕性降雨與降雨侵蝕力的關系［J］.水土保持研究，2018，25（2）：8-14HE Shaolang，HE Xiaowu，LI Fengying，et al.Relationship between erosive rainfall and rainfall erosivity of each rainfall grade in Jiangxi province in recent 60 years［J］.Research of Soil and Water Conservation，2018，25（2）：8-14

［30］ Torrence C，Compo G P.A practical guide to wavelet analysis［J］.Bulletin of the American Meteorological Society，1998，79（1）：61-78

［31］ 夏正兵.黃河流域極端氣候下降雨侵蝕力時空特征研究［J］.水電能源科學，2021，39（7）：16-19XIA Zhengbing.Study on temporal and spatial characteristics of rainfall erosion under extreme climate in the Yellow River Basin［J］.Water Resources and Power，2021，39（7）：16-19

［32］ 伊力哈木·伊馬木.新疆維吾爾自治區1981—2018年降雨侵蝕力的空間變化特征［J］.水土保持通報，2020，40（4）：1-9Yilihamu·Yimamu.Spatial variation of rainfall erosivity in Xinjiang Uygur Autonomous Region from 1981 to 2018［J］.Bulletin of Soil and Water Conservation，2020，40（4）：1-9

［33］ Duan J，Liu Y J，Tang C J，et al.Efficacy of orchard terrace measures to minimize water erosion caused by extreme rainfall in the hilly region of China：long-term continuous in situ observations［J］.Journal of Environmental Management，2021，278：111537

Spatial-temporal variation of rainfall erosivity in Jiangxi province from 1961 to 2019

ZHENG Xiao BAI Shuying GAO Jixi

1School of Remote Sensing ＆ Geomatics Engineering，Nanjing University of Information Science ＆ Technology，Nanjing 210044

2Nanjing Institute of Environmental Sciences，Ministry of Ecology and Environment，Nanjing 210042

Abstract As one of the important factors affecting soil erosion，rainfall erosivity （R） is related to rainfall characteristics of amount，duration，intensity，and kinetic energy.Using the daily rainfall data of 25 meteorological stations in Jiangxi province from 1961 to 2019，we analyzed the spatial distribution and variation trends of rainfall erosivity in Jiangxi province via rainfall erosivity model，as well as Mann-Kendall correlation test，wavelet analysis and Kriging spatial interpolation.The results showed that the annual average rainfall and rainfall erosivity varied similarly in spatial distribution，which were gradually increased from southern Jiangxi to northern Jiangxi;the rainfall erosivity in spring and summer ranged from 3 000 to 6 000 MJ·mm·hm-2·h-1·a-1，which was significantly higher than that in autumn and winter，with the maximum rainfall erosivity occurring in spring;the average annual rainfall erosivity in Jiangxi decreased from north to south，which maximized in northern Jiangxi and followed by central Jiangxi then southern Jiangxi.The research shows that the rainfall-induced soil erosion was increasing in Jiangxi province，especially in spring and summer，and there were obvious differences in their spatial distribution.

Key words rainfall erosivity;rainfall;Mann-Kendall correlation test;wavelet analysis