龙训建， 翁薛柔，2， 叶琰， 张逸轩，曾必顺， 叶勇， 徐廷兵

1.西南大学 资源环境学院，重庆 400715；2.中山大学 地理科学与规划学院，广州 510006；3.重庆市人工影响天气办公室，重庆 401120；4.西南大学 后勤保障部，重庆 400715

土壤侵蚀是农业生产中的重大生态环境威胁，任其发展不利于经济社会的可持续[1-3]．受全球气候变化和人类活动的共同影响，土壤侵蚀问题变得更加复杂，也受到更多关注，尤其是降雨引起的侵蚀问题(降雨侵蚀力)成为热点研究内容之一[4-6]．已有研究表明，在不利的地形地质条件、森林乱伐、过度耕作和放牧等因素影响下，土壤更易产生侵蚀与流失[2，7-10]．目前，水文循环速率加剧导致极端降雨事件频发，这预示着降雨侵蚀力和土壤侵蚀率的增加[11]．因此，准确评估降雨侵蚀力有利于采取针对性措施制止或减缓水力侵蚀．

降雨侵蚀力可描述不同降雨特征(雨量、历时、强度等)对土壤侵蚀的影响[12]，其计算模型有纯雨量模型和综合指标模型两类．纯雨量模型主要是以日雨量[13-14]、月雨量[15-16]、年雨量[17-18]为基础数据的模型，如Richardson等[19]、Renard等[20]、Lee等[21]分别基于日雨量、月雨量、年雨量资料，建立了估算美国、韩国等国家和地区降雨侵蚀力的简易模型；国内主要以章文波等[13]推荐的利用日雨量估算半月降雨侵蚀力，进而统计月、年降雨侵蚀力的算法应用最为广泛．综合指标模型不仅包括雨量因子，更考虑了雨强等指标的作用，如卜兆宏等[22]发现汛期雨量、年最大30 min雨强(I30)指标与土壤流失量间存在较高的相关性，建立了汛期雨量和年I30值乘积的年降雨侵蚀力模型．

以上两类模型对数据分辨率要求较低，其结果通常属于估算结果，不适宜应用于对精确性有要求的工作领域．在具备高时间分辨率和连续的降雨数据前提下，应用Wischmeier等[23]提出以次降雨总动能(E)与最大30 min雨强(I30)的乘积(EI30)作为降雨侵蚀力指标，并将之与通用土壤流失方程(USLE)相结合，降雨侵蚀力相关研究的准确性就能得到巨大提高．同时，王万忠等[24]通过比较EI30与EI10，EI60等计算方法，认为EI30符合我国绝大多数地区的降雨特性，建议采用EI30作为中国降雨侵蚀力的计算公式．基于此，本文以重庆市2009-2018年32个国家级基本气象台站的逐分钟降雨数据为基础，应用降雨侵蚀力EI30计算公式，利用MK秩相关检验、空间插值法和系统聚类法，分析探讨研究区降雨和降雨侵蚀力在时空尺度上的趋势变化特征，以期为区域生态环境保护提供可靠的数据基础和理论依据．

1 资料与方法

1.1 研究区概况

重庆市位于中国西南，是长江上游地区的经济、金融、科创、航运和商贸物流中心．重庆总面积8.24万km2，辖38个区县(自治县)；常住人口3 101.79万人，境内有长江、嘉陵江、乌江穿过，流域面积大于50 km2的河流374条，过境水资源丰富，降雨量相对充沛但分布不均，多年平均降雨量1 100～1 200 mm[25]．

1.2 数据来源

本次基于国家气象数据共享服务网，收集了重庆市32个国家级气象站点2009-2018年的逐分钟降雨量数据，各气象站点位置见图1．

1.3 研究方法

1.3.1 降雨侵蚀力计算

1) 次降雨

包括次降雨过程的界定和侵蚀性雨量的界定．次降雨过程的界定参照USLE方程推荐的雨量历时间隔标准，即将降雨间歇≤6 h的降雨视为一个次降雨过程，否则视为下一次降雨．侵蚀性雨量采用国内通用的侵蚀性雨量分级标准和史冬梅等[18]对重庆市降雨侵蚀的多年研究经验，即研究区侵蚀性降雨标准为次降雨雨量≥12 mm．

2) 次降雨侵蚀力

次降雨侵蚀力计算公式如下[26]：

R=EI30=E·I30

(1)

式中，I30是单场次降雨最大30 min雨强，单位mm/h；E为单场次降雨的动能，是降雨-历时函数，单位MJ/ha．

(2)

式中，D是次降雨历时，单位为h；v是次降雨量，单位mm；i是雨强值，单位mm/h．

1.3.2 趋势变化分析

离差系数(Cv)和趋势系数(r)两项指标通常被用来描述降雨侵蚀力的年际变化特征，计算公式为：

(3)

(4)

1.3.3 空间插值分析

空间分布是空间对象群体在一定空间区域的散布特性，是空间分布研究的主要内容之一．考虑研究区内各气象站点为离散分布，本次选用精确性插值法——反距离权重法(Inverse Distance Weighted，IDW)来实现降雨侵蚀的空间分布．反距离权重法是加权移动平均方法，以插值点与样本点间的距离为权重进行加权平均，距离插值点越近的样本点赋予的权重越大[27-28]，计算公式为：

(5)

式中，υ0是未知点的估计值，υi是第i(i=1，2，…，n)个样本点的值，di是采样点与未知点之间的距离，k是距离的幂值，默认值为2，一般取值范围在0.5～3之间可获得最合理结果．

1.3.4 Ward系统聚类法

Ward系统聚类法是根据选定的特征来识别相对均一的个案(变量)组，基于方差分析思想，以各类别之间的欧式距离最小作为分类标准，从单独聚类中的每个个案(或变量)开始对各聚类进行组合，直至剩下一个类别．Ward系统聚类法可以判断研究区内各站点的降雨及其侵蚀力在同一时间尺度内的相似性，以验证空间分布特征分析的准确性．

1.3.5 重庆市四季和分区划分

四季划分一般按照气候划分标准，根据当地30年以上的气温平均值，划出四季交换的日期．根据重庆市气候中心统计数据，重庆为四季分明区，入春的平均时间为3月2日，平均在6月1日入夏，入秋时间为10月2日，12月22日是平均入冬时间．据此，将重庆市的3-5月划为春季，6-9月划为夏季，10-12月划为秋季，1-2月划为冬季．

重庆市一般行政分区涉及4个区域，主城片区、渝西片区、渝东南片区和渝东北片区；从地形地貌特征上一般分为3个区域．相比于行政区域，降雨侵蚀力与地形地貌更为直接相关，本次将研究区整体划为渝西区域、渝东南区域和渝东北区域3部分，各分区对应行政区范围见表1所示．

表1 重庆市分区范围划分

2 结果与分析

2.1 降雨量时空变化特征分析

2.1.1 降雨的年际变化特征

重庆市2009-2018年多年平均降雨量为1 082.98 mm，各分区年平均面雨量最小值发生在渝西片区的2011年，仅758.24 mm，最大发生在渝东北片区的2017年，为1 476.02 mm，两者倍差为1.95倍；整体上，渝东北多年平均值最大，为1 130.55 mm，渝西最小，为1 047.02 mm，两者倍差为1.08倍(图2)．为判断各站点间多年平均降雨量间的差异，表2给出了重庆市2009-2018年各区域面降雨量的年际变化特征参数值．可以看出，渝东北区域标准差、Cv值均最大，说明该区域降雨年际差异最大；研究时段内，各区域年平均面降雨量均呈增加趋势，其中渝西地区增加幅度最为明显．

图2 重庆市及各分区年降雨量折线

表2 各分区平均面降雨量特征参数

2.1.2 降雨的年内分布特征

重庆市及各分区降雨量季节分布及比例情况见表3．从季节比例角度来看，研究区内冬季降雨量占全年降雨量的比例最小，仅为全年降雨量的5%左右，这与冬季气温低，大气层结稳定，绝对水汽含量低，产生降雨的概率也较低有关；其次为秋季，仅占全年降雨量的12.39%～15.29%；夏季降雨量占全年降雨量的比例最大，达49.37%～55.15%，这一时期的7-8月，夏季西太平洋副热带高压西伸北移，控制重庆大部分地区，导致区域内盛行下沉气流，产生降雨的概率较小，但是，受地形地貌影响，大气层结不稳定，绝对水汽含量高，易产生对流性降雨，从而表现出高的强降雨比例．同时，在此气候作用下，重庆市夏季频繁发生暴雨、洪涝、山体滑坡，秋季多伏旱，冬季则降雨频次低、雨量少．

在区域分布上，各季节年平均降雨量存在一定的差异，其中渝东北最大，渝东南次之，渝西最小．在季节分布上，各分区分布特征与全市相似，夏季最多，冬季最少．渝西片区夏季降雨量比例在3个片区中最高，达到了55.15%，秋季降雨量比例最低，这与渝西区域夏季多暴雨，秋季伏旱突出的特征相一致．

表3 各分区年平均降雨量季节分布及比例

审图号：GS(2019)3333号(自然资源部监制)．图3 年平均降雨量空间分布

2.1.3 降雨量的空间变化

应用反距离权重法得到重庆市多年平均降雨量的空间分布．结果显示，降雨量总体上呈现西低东高的分布趋势．高值区主要分布在渝东南的酉阳、渝东北的城口-开州-万州-云阳一带，其多年平均降雨量均超过了1 200 mm；低值区主要分布在渝西的潼南-大足-荣昌-永川一带，介于950～970 mm之间；此外，长江沿线的河谷区域降雨也大多低于两侧山区．这与方德贤等[29]、马丽珲等[30]对重庆市近年来降雨量空间分布特征的研究结果一致(图3)．

从重庆市多年平均季节降雨量的空间分布可知，不同季节降雨量的空间分布存在显著差异．在春季，高值区主要分布在渝东南酉阳县，达379.11 mm，低值区主要分布在渝西片区的潼南-大足-荣昌-永川一带，降雨量在206.99～224.70 mm．在夏季，西部降雨明显增强，整个区域呈现两边高、中间低的分布，渝东北区域降雨最多，渝西次之，区域中部较低，这与夏季西南涡活动较为频繁以及西太平洋副热带高压的西伸、北抬有关[31-32]．秋季降雨较夏季整体偏弱，呈西低东高的空间分布，高值区位于渝东北云阳区，低值区主要分布在渝西片区的潼南-大足-荣昌-永川一带．冬季降雨呈西高东低分布，高值中心在渝东南的酉阳，低值中心为渝东北的奉节．季节分布特征表明渝西地区的降雨量显著少于其他区域，这也解释了渝西地区易发生季节性干旱的特征(图4)．

为进一步研判降雨量的空间分布相似特征，应用Ward系统聚类法对研究区32个站点的月平均降雨量进行分类．由结果可知，研究区降雨量可分为3类，其中类Ⅰ包括潼南、大足、荣昌、永川、铜梁、北碚、合川、渝北、璧山和江津，主要分布在渝西片区；类Ⅱ包括万盛、沙坪坝、巴南、南川、綦江和酉阳，多聚集于重庆中部；其余区县归为类Ⅲ．研究区月降雨的分类特征具有空间相似性，与对应时间尺度降雨空间分布特征一致(图5)．

对聚类后的各类月降雨量的中心距离进行计算，得到各分类的降雨量年内变化曲线．类Ⅰ的年内变化曲线较尖瘦，在7月达到了峰值，属于类1的区域在7月需尤为注意暴雨事件的频发；类Ⅱ的曲线较扁平，6-8月雨量较多；类Ⅲ在6月达到降雨量峰值，8月次之．这说明重庆西部降雨集中发生在7月，重庆中部降雨在年内的集中度较西部弱，重庆东部降雨的年内分布与中、西部有明显差异(图6)．

2.2 降雨侵蚀力时空变化特征分析

2.2.1 降雨侵蚀力年际分布特征

应用公式(1)和(2)计算得到各站点次降雨侵蚀力计算结果，根据区域分区原则，采用面积加权法求得各分区的降雨侵蚀力．由公式(3)和(4)计算得到研究区降雨侵蚀力的年际变化统计参数．结果表明，重庆市2009-2018年多年平均降雨侵蚀力为1 816.99 MJ·mm/(hm2·h)，分区降雨侵蚀力最大值出现在渝东南片区的2016年，达3 266.88 MJ·mm/(hm2·h)，最小值则出现在该区域的2009年，仅949.75 MJ·mm/(hm2·h)，两者倍差3.44倍．不同分区多年平均降雨侵蚀力存在一定差异，渝东北区域年均降雨侵蚀力最大，为1 857.68 MJ·mm/(hm2·h)，侵蚀潜在危险性较高，渝西次之，为1 808.82 MJ·mm/(hm2·h)，渝东南最小，为1 763.37 MJ·mm/(hm2·h)．标准差变化范围为382.72～734.12，最大值出现在渝东南，最小值出现在渝西．渝东南地区Cv值是渝东北和渝西地区的2倍，其降雨侵蚀力的年际波动显著高于其他区域．趋势系数变化范围为0.27～0.70，各分区降雨侵蚀力均呈增加趋势(表4)．

审图号：GS(2019)3333号(自然资源部监制)．图4 各季节降雨量空间分布

表4 各分区降雨侵蚀力演变趋势参数

图5 系统聚类结果

图6 各分类降雨量的年内分布曲线

由重庆市2009-2018年各分区降雨侵蚀力的年际变化曲线图可知，各地区降雨侵蚀力变化幅度大且存在较大差异，重庆市降雨侵蚀力在2009-2016年整体呈上升趋势，在2016年达最大值后，近两年有所下降；渝东南地区在2016年和2017年降雨侵蚀力反差强烈，这主要与渝东南的降雨有关(图7)．

2.2.2 降雨侵蚀力年内分布特征

根据研究区降雨侵蚀力在不同季节的均值及比例可以看出，统计时段内，研究区降雨侵蚀力的季节分布表现出不均匀性．其中，夏季比例均超过70%，渝西地区更是达到了79.39%；春季降雨侵蚀力比例仅次于夏季，为17.25%～24.89%；秋季较少，降雨侵蚀力比例为3.05%～4.66%；冬季降雨侵蚀力最低，比例不足0.5%．根据统计结果，秋冬季侵蚀力普遍较低，这与该时段内降雨整体比例较低，能引起土壤侵蚀的有效降雨量较少相关．无论是全市整体还是分区比较，降雨侵蚀力的高值期均出现在夏季，应加强该时段内水土流失的监测、防治与监督，以最大程度降低水土流失量，减小水土流失发生的可能性(表5)．

图7 重庆市及各分区年降雨侵蚀力折线

表5 研究区降雨侵蚀力季节分布特征

审图号：GS(2019)3333号(自然资源部监制)．图8 年降雨侵蚀力空间分布

2.2.3 降雨侵蚀力空间分布特征

采用反距离权重法对研究区多年平均降雨侵蚀力进行插值，得到空间分布结果．可以看出，年降雨侵蚀力呈中间低、两端高的空间分布，高值区主要分布在渝东北的梁平、开州一带，低值区主要位于研究区中部(图8)．

按照季节划分，研究区各季节降雨侵蚀力的空间分布差异也存在统计学意义．在春季，高值区主要分布在渝东南片区的酉阳，低值区位于渝西的潼南-大足-荣昌-永川一带；夏季降雨侵蚀力呈中部低、两端高的分布特征；在秋季，渝西地区降雨侵蚀力较低，渝东北的垫江、梁平和渝东南的降雨侵蚀力较高；在冬季，仅有酉阳-秀山一个高值中心，其余地区降雨侵蚀力均较低．其中，渝东北的高值区主要为山地丘陵区，土壤类型多属于黄壤、石灰岩土、黄棕壤，易于发生侵蚀；渝东南高值区主要分布在海拔高差较大、易发生土壤流失的酉阳、秀山一带；渝西片区的高值区主要分布在土地复种指数较高的区县；这些区域的降雨侵蚀潜在危险性都较高，加之地形地貌的影响，降雨导致的水土流失、土壤侵蚀等问题都亟需注意(图9)．

审图号：GS(2019)3333号(自然资源部监制)．图9 各季节侵蚀力分布

为进一步探讨研究区内侵蚀力的分布类型，应用Ward联结法对研究区32个站点的月降雨侵蚀力进行系统聚类，以多年平均月侵蚀力为初始聚类，32个站点为聚类成员进行聚类分析，形成距离谱系聚类图．以15为划分距离，研究区降雨侵蚀力可分为3类，其中类Ⅰ包括城口、潼南、大足、荣昌、永川和江津，此类在空间上主要归属于渝西片区，类Ⅱ包括开州区、云阳、巫溪、梁平、垫江、彭水和酉阳，此类多数聚集于重庆东部，其余区县归属于类Ⅲ(图10)．

图10 各站点月平均降雨侵蚀力聚类树

图11 各分类降雨侵蚀力的年内变化曲线

各分类降雨侵蚀力的年内变化曲线如图11所示，类Ⅰ与类Ⅱ在年内的分布特征相似，在6月达到峰值，7月曲线下凹，8月回升，之后持续下降，但类Ⅱ在年内的变化趋势较类Ⅰ更加平缓．类Ⅲ在5月达到降雨侵蚀力最大值，6月与5月基本持平，7-8月明显下降，9月有小幅回升，与其他两类变化趋势不同．这说明重庆西部和东部降雨侵蚀力在年内变化趋势相似，但渝西变化更加剧烈．

对比分析研究区降雨量和降雨侵蚀力的聚类结果发现，以15为划分距离，降雨量和降雨侵蚀力均可分为3大类，但分类结果存在差异．降雨量分类结果存在较强的区域性，同一类型的地理位置明显聚集，而同一类型降雨侵蚀力的站点分布相对分散．

各分类降雨和降雨侵蚀力的年内变化曲线分别呈单峰、双峰特征，变化趋势有显著差别．重庆各地区降雨量多在6-7月达到峰值，而降雨侵蚀力的最高值为5-6月．重庆西部和中部的降雨量年内变化较为相似，而重庆西部和东部的降雨侵蚀力年内变化相近．这说明尽管降雨量是引发水力侵蚀的重要因素，但降雨侵蚀力受其他降雨特征参数的影响明显，仅通过降雨量研究降雨侵蚀力分布特征缺乏精确性．

3 讨论

3.1 降雨与降雨侵蚀力时间变化特征对应关系

通过对研究区2009-2018年的降雨量和降雨侵蚀力的变化特征进行分析，降雨量和降雨侵蚀力均存在年际和年内间的波动变化，且空间分布不均的特征．

在研究区域内，从降雨量的年际变化看，各站点多年平均降雨量为950.28～1296.17 mm，统计特征参数中，渝东北变差系数最大，主要是由于时段内最大降雨量达最小降雨量的1.6倍，具有热带降雨特点．从年内降雨分布特征看，主要集中于春、夏两季，在西太平洋副热带高压作用下，降雨次数和降雨量都几乎集中于这一时段．这些特征与付友行等[33]、史东梅等[18]对重庆市内涉及的小流域和单一行政区域的研究成果基本一致．对于降雨侵蚀力，变差系数最大值则出现在渝东南，这主要与计算时段内，年平均降雨侵蚀力最大值、最小值均发生在渝东南有关．

由此，对比分析研究区内平均面降雨量和平均降雨侵蚀力的极值发现，年降雨量大的区域，并不一定会产生最大的降雨侵蚀力，可能是由于单次降雨侵蚀力受单次降雨量、降雨动能和降雨强度综合影响．

3.2 降雨与降雨侵蚀力空间特征对应关系

年平均降雨量和年降雨侵蚀力空间分布结果表明，两者最大值均位于渝东北城口-开州-万州-云阳一带和渝东南的酉阳，这与区域单次降雨量和单次降雨侵蚀力最大值出现在该区域有关，由此也可解释次降雨量越大越易于产生更多的降雨侵蚀．但在降雨侵蚀力的高值分布中，渝西片区也有零星分布，这主要是受区域极端暴雨事件影响．如铜梁区在2016年6月30日、7月13-14日和7月18-19日发生了3场暴雨事件，单次降雨量分别为96.7，78.3和54.5 mm，对应产生的降雨侵蚀力分别为264.32，353.05和136.15 MJ·mm/(hm2·h)．3场降雨量占当年总降雨量的18.89%，对应的累积降雨侵蚀力占当年总降雨侵蚀力的32.99%，这一现象也可从季节空间分布结果中得到印证．尤其是在夏季，尽管降雨量高值区几乎仅分布于渝东北的城口一带，但降雨侵蚀力在渝西片区的高值区表现十分明显．与此同时，基于月尺度的多年平均降雨量和降雨侵蚀力聚类结果表明，空间分布上，两者的类Ⅰ具有高度区域相似性．

另外，根据殷水清等[34]基于全国范围内774个气象站1961-2016年逐日降雨数据计算得到的多年平均年和次降雨侵蚀力空间分布结果，重庆区境内的多年平均降雨侵蚀力值为2 000～3 000 MJ·mm/(hm2·h)，与本文采用具有高时间分辨率的计算结果为1 400～2 300 MJ·mm/(hm2·h)不同，由此可说明，利用更小时间尺度的降雨数据计算降雨侵蚀力更有利于区域生态环境保护工作的开展．

4 结论

1) 重庆市2009-2018年平均降雨量为1 082.98 mm，区域内各站点多年平均降雨量为950.28～1 296.17 mm，最大值是最小值的1.36倍，差异有统计学意义．降雨量的季节分布不均，从大到小依次为夏季、春季、秋季、冬季，其中夏季占据全年的50%以上．空间分布上降雨量高值区分布在渝东北片区，低值区主要分布于渝西片区．

2) 研究区各分区降雨侵蚀力的演变趋势统计结果表明，降雨侵蚀力最大值出现在渝东南，最小值出现在渝西；渝东南降雨侵蚀力的年际波动高于渝东北和渝西．受地形地貌、土地复垦频率和其他人类活动等因素影响，降雨侵蚀力高值区主要分布在渝东北的梁平、开州一带，渝东南的酉阳、秀山一带，低值区位于重庆中部．

3) 通过Ward系统聚类法对研究区降雨量和降雨侵蚀力进行空间分类，结果显示，研究区降雨和降雨侵蚀力均可划分为3类，但分类结果不同，同一类型降雨具有明显的地域聚集特征，而降雨侵蚀力相对分散．重庆中部降雨在年内的集中度较西部弱．重庆西部和东部降雨侵蚀力在年内变化趋势相似，但渝西变化更加剧烈．