我国长江以南山地丘陵区地处热带亚热带，雨量充沛，由于自然环境、花岗岩母质特征及人为活动破坏，土壤侵蚀日益加剧，现已成为我国仅次于黄土高原的第二大侵蚀区，特别是花岗岩红壤区已成为我国水土流失严重的区域^[1-3]。土壤侵蚀是一个复杂的过程，受降雨和下垫面等因素的综合影响。在花岗岩深厚的风化壳中，各个风化层土壤性质不同，不同强度的水力侵蚀造成了不同层次土壤不同程度地暴露在外，其理化性质与抗侵蚀能力差异非常明显^[4-7]，造成坡面侵蚀产沙过程与泥沙颗粒特性存在较大差异^[8]。降雨是引起土壤侵蚀的驱动因子^[9]，大量研究表明，降雨强度越大土壤侵蚀越严重，特别是在细沟发育的活跃阶段^[10]。而短历时暴雨是造成该区坡面侵蚀的主要雨型^[11]，加之缺乏植物和有机质，裸露的花岗岩红壤坡面极易受到降雨和径流冲刷的影响^[12]。随着极端降水情况的频发，暴雨对土壤侵蚀的影响应引起人们的重视。

土壤颗粒随径流的输移和分选特征是坡面侵蚀的研究重点，一直被国内外学者所关注。当降雨动能较低时，土壤中细颗粒优先被分散^[13-14]；更高的降雨动能会增加径流搬运能力，使得泥沙颗粒中粗颗粒含量增加，黏粒含量相对减少^[15]。当坡面径流水流功率超过临界值时，径流分选作用降低，泥沙的颗粒组成与原生土壤逐渐趋于相似^[16]。侵蚀过程对泥沙颗粒的分选也存在一定的作用^[17]，当细沟发育后，降雨侵蚀颗粒粗化现象更为明显^[18]。泥沙的颗粒分布也受土壤母质和降雨侵蚀力影响，土壤母质相同的情况下，雨强和坡面条件会影响泥沙颗粒的比例，黏粒含量会随着坡长和雨强的增加而增加^[19]。Berger等^[20]在研究细沟发展与土壤侵蚀的关系中提出雨强和坡度对细沟泥沙颗粒的分布均有影响，且雨强的影响大于坡度。当土壤母质不同时，泥沙颗粒分布也会产生明显差异^[21]。

目前，关于花岗岩红壤的侵蚀特征研究中，有限的探索主要集中在红壤表土与崩积堆的泥沙颗粒分选特征方面^{[12, 22]}，而不同侵蚀程度土壤间的坡面侵蚀规律与泥沙分选特征研究不足。鉴于以上研究背景，本试验以花岗岩发育不同侵蚀程度红壤为对象，研究了人工模拟降雨和径流冲刷条件下的坡面侵蚀产流产沙特征，揭示了降雨强度和上方来水流量双因素与坡面泥沙颗粒特征的关系，以期为南方红壤区水土保持提供科学依据。

1 材料与方法 1.1 供试材料

供试土样采集于福建省长汀县和田镇（25°37′58.69″N，116°27′29.41″E）。该区属于亚热带季风湿润气候，年平均气温18.5 ℃，年平均降水量1 650 mm。采样时参照Wu等^[23]判断土壤侵蚀程度依据及立地崩岗侵蚀发生发育程度，在尚未发生崩岗侵蚀的花岗岩发育残积土坡面采集轻微侵蚀程度土体（E1），采样深度为10 cm，土壤层次为表土层，有丰富植被生长（主要为铁芒萁和马尾松），结构较好。在有崩岗侵蚀发生的集水坡面靠近坡顶位置采集中度侵蚀土体（E2），采样深度为10 cm，土壤层次为红土层，地表草本植被稀少（以马尾松为主），土体黏重，土粒细腻紧实，土壤结构严重破坏，透水性差。严重侵蚀土壤（E3）采集自崩岗崩壁，采样深度为200~230 cm，土壤层次为过渡层，呈红色夹杂零星黄色，土壤极为贫瘠，保肥持水能力差，土体结构疏松，用手搓会有沙质感。采样点及剖面详情见图 1。将土样运回试验室后，自然风干并过10 mm筛，剔除石块、根系等杂质后备用。土壤有机质采用重铬酸钾容量法测定，土粒密度采用比重瓶法测定，pH通过pH计测定；土壤粒径分布使用超声波振荡仪分散土壤颗粒，并通过湿筛法结合吸管法测定，采用美国制粒径分级，得到砾石（> 2 mm）、粗砂粒和细砂粒（0.05~2 mm）、粉粒（0.002~0.05 mm）、黏粒（< 0.002 mm）质量分数。土壤容重采用野外原位采集的环刀测定（100 cm³）。土样的基本理化性质见表 1。

1.2 研究方法

模拟降雨装置采用垂直旋转下喷式自动模拟降雨系统，该降雨器通过小、中、大三种规格的喷头组合实现15~220 mm·h^–1降雨强度范围，降雨均匀度大于90%。降雨高度为10 m，降雨雨滴可以达到终点速度^[24]。试验土槽规格为3.0 m（长）×0.8 m（宽）×0.45 m（高），坡度可在0°~30°自由调节。本试验基于开放小区，在降雨-径流耦合条件下进行模拟实验。径流通过恒压水泵和槽子后端的稳流槽提供，流量由供水管道上的两个控制阀调节。按照野外采集的环刀测定的容重，将3个层次土样分别按其容重作为试验土壤填入土槽。填装土样时，采用分层填装，每层5 cm，并在上层土壤填入前对下层土壤刮毛，逐层填至30 cm以保持土壤颗粒分布均匀贴合紧密。完成土壤装填后在土槽顶部覆盖一个15 cm长的纱布，以尽量减少径流冲刷和边缘效应的影响，试验重复2次，试验装置见图 2。

根据长汀暴雨经验公式，针对我国亚热带气候区典型的强风暴，设置了在90 mm·h^–1+2.00 L·min^–1和120 mm·h^–1+2.66 L·min^–1两种试验条件下（假定坡面上游有2 m²的集水区）对三个侵蚀程度土壤共进行了6组降雨-冲刷试验。丁光敏^[25]在对福建崩岗侵蚀成因的调查研究中发现10°~25°中坡为最易产生侵蚀的坡度，故本试验设置试验坡度为15°，也是崩岗侵蚀严重区集水坡面常见坡度。试验降雨强度根据当地设计暴雨强度公式计算了一年一遇、五年一遇、十年一遇下I_30min、I_45min、I_60min选定，其中90 mm·h^–1为多年平均暴雨降雨强度，120 mm·h^–1为极端暴雨条件下降雨强度。流量设置依据严冬春等^[26]不同坡度下细沟侵蚀发生的临界坡长，等比例添加上方来水使对应降雨强度下坡面汇水坡长达到5 m。试验前以小雨强（30 mm·h^–1）对坡面进行预湿润处理，使用纱布覆盖表土，减少雨滴溅蚀影响，直至产流，并放置24 h后开始试验。降雨历时1 h，每分钟采集3组泥沙样品，用于计算土壤侵蚀速率、泥沙浓度，并利用每三分钟时间段中第三分钟的三个样品（即第3、6、9…60分钟的三个样品）通过湿筛和吸管法测定泥沙颗粒粒径：砾石（> 2 mm）、砂粒（0.05~2 mm）、粉粒（0.002~0.05 mm）、黏粒（< 0.002 mm）。

1.3 参数获取

综合有关研究选取下列相关参数。土壤侵蚀速率^[27]（土壤剥蚀率）为单位时间、单位侵蚀面积上土壤在水流侵蚀动力的作用下被剥蚀的土壤颗粒质量，表征径流对坡面土壤的分离能力，其求解方法如下：

${D_r} = \frac{{{m_t}}}{{BLT}}$

(1)

式中，D_r为土壤侵蚀速率，kg·m^–2·min^–1；m_t为累积产沙量，kg；B为水宽，cm；L为沟长，m；T为时间，min。

为了研究泥沙粒径分选机制，计算了泥沙中各粒径含量和原土样对比的富集率（Enrichment ratio，ER）^[22]，具体计算公式如下：

${\rm{ER}} = \frac{{{P_s}}}{{{P_o}}}$

(2)

式中，ER为富集率；P_s为泥沙某一粒径的百分比，%；P_o为土壤某一粒径的百分比，%。ER > 1表示该颗粒富集，ER < 1表示该颗粒损耗。

分形维数通过表征土壤颗粒粒径分布、土壤质地均一程度、土壤团粒结构团聚体、水稳性团聚体和粒径组成等^[28]，可用于反映坡面土壤质地及肥力对侵蚀过程的响应^[29]及坡面土壤结构对侵蚀的影响^[30]。Mandelbrot ^[31]建立了二维空间颗粒大小分形特征模型。此后，Tyler和Wheatcraft^[32]在此基础上建立了三维空间的体积分形维数模型：

$\frac{{W(\delta > {d_i})}}{{{W_{\rm{o}}}}} = 1 - {\left( {\frac{{{d_i}}}{{{d_{\max }}}}} \right)^{3 - D}}$

(3)

对分形维数模型进行推导，两边取对数，以左右两边为因变量和自变量进行线性拟合，拟合的直线斜率即为3$ - $D，从而可求出分形维数D。

$\lg \frac{{W(\delta > {d_i})}}{{{W_{\rm{o}}}}} = (3 - D)\lg \frac{{{d_i}}}{{{d_{\max }}}}$

(4)

式中，d_i=（d_j+d_j+1）/2，W（$\delta $ > d_i）为粒径大于d_i累积土粒质量；W_o为各粒级土粒质量之和；d_i为筛分粒径范围[d_j，d_j+1]平均值；d_max为某个粒径级[d_j，d_j+1]最大粒级平均直径。

1.4 数据处理

数据处理与分析在Excel 2016、Origin 2017与SPSS 19.0软件中完成，相关分析采用多元逐步回归分析。

2 结果与讨论 2.1 不同降雨强度和上方来水流量对产流产沙的影响

不同降雨-冲刷强度下坡面产流率变化如图 3。产流率随产流时间表现为先增加、后相对稳定或略有波动。其原因是降雨初期，土壤含水率较低，未达到饱和，雨水下渗作用延缓了径流的产生。随着降雨历时的延长，土壤含水率逐渐饱和，加之雨滴击溅导致地表粗颗粒破碎分离，细颗粒逐渐堵塞土壤孔隙，促进径流的产生，当坡面逐渐达到稳定下渗率时，产流率也趋于稳定。两种降雨-冲刷强度下产流率表现为E2 > E1 > E3，差异显著（表 2）。其原因是坡面侵蚀过程取决于降雨特征、土壤性质和地表条件等因素^{[23, 33]}，不同侵蚀程度土体的容重存在显著差异，容重越大产流率也会越大（r=0.91，P < 0.01），而降雨-冲刷强度的增加会减少到达稳定阶段需要的时间，高降雨-冲刷强度下产流率也会更高。

不同降雨-冲刷强度下各侵蚀程度土壤侵蚀速率变化如图 4所示。除E1在90 mm·h^–1+2.00 L·min^–1 条件下，其余各试验条件下土壤侵蚀速率随着产流历时呈增加趋势。各土体土壤侵蚀速率差异显著，表现为E3 > E2 > E1（表 2），E3在90 mm·h^–1+2.00 L·min^–1条件下土壤侵蚀速率为0.324 kg·m^–2·min^–1，为E2的1.57倍和E1的9.95倍；在高雨强下E3土壤侵蚀速率为0.740 9 kg·m^–2·min^–1，为E2的1.76倍和E1的2.40倍。随着降雨-冲刷强度的增加，各侵蚀强度土体土壤侵蚀速率也显著增加，E1、E2及E3分别增加了9.48倍、2.03倍、2.28倍。究其原因为E1在90 mm·h^–1+2.00 L·min^–1条件下，侵蚀过程主要为雨滴击溅与层状侵蚀，未发育出细沟。随着土壤表面可被搬运的细颗粒被输移耗尽，产流率增加加剧了坡面粗化程度，使得表土层中丰富的砾石裸露在地表对下层土壤起到了保护作用。当降雨-冲刷强度增加后，坡面水流功率增加，超过了地表对粗颗粒的阻力，使其推移搬运。随着细沟开始发育，坡面集中水流侵蚀能力增大，土壤侵蚀速率随之增加^[34]，并伴有细沟沟壁崩塌，使得土壤侵蚀速率存在波动。由于水流功率较高，E2、E3土体与E1之间基本性质存在差异，粗颗粒含量也远少于E1，导致E2、E3坡面细沟发育过程剧烈，土壤侵蚀量也高于E1。

为研究降雨-冲刷强度对各侵蚀程度土体侵蚀特征的影响，进行了多元逐步回归分析。结果表明，容重、降雨-冲刷强度和土壤粉粒含量对产流率起正向作用（R²=0.800，P < 0.01），容重起主要作用。降雨-冲刷强度和土壤粉粒含量对泥沙浓度起正向作用，土壤黏粒含量对泥沙浓度起负向作用且起主要作用（R²=0.689，P < 0.01）。土壤侵蚀速率受降雨-冲刷强度、砂粒含量和土壤粉粒含量的影响（R²=0.715，P < 0.01），土壤砂粒含量和粉粒含量起主要作用。降雨强度在土壤侵蚀中起着双重作用，不仅使团聚体消散破碎，而且形成径流促进泥沙搬运。花岗岩红壤随着侵蚀程度的加深，土壤性质随之发生变化，表现为粉粒和砂粒含量的增加与黏粒含量的减少。

2.2 不同降雨强度和上方来水流量下泥沙分选差异

图 5是各侵蚀程度土体泥沙颗粒分布随产流时间的变化。E1各粒级泥沙分布随产流历时变化均较大，黏粒和粉粒含量逐渐减少，而粒径较大的砂粒和砾石含量逐渐增加。其原因是试验初期坡面侵蚀为片蚀阶段，悬移-跃移是主要的搬运机制^[35]，径流会优先选择易侵蚀的细颗粒，而坡面形成的薄层水流在雨水冲击作用下，径流雷诺数的增加加剧了细颗粒的侵蚀。随着细沟侵蚀阶段开始，径流功率增加，加强了对大颗粒泥沙或团聚体的搬运能力，使大粒径的砾石发生输移^[36]。集中水流作用下细沟内泥沙分选作用下降，细沟间侵蚀仍以雨滴击溅为主，但产流初期细颗粒已被大量侵蚀，导致了细颗粒随降雨时间搬运量下降，致使坡面发生粗化。而泥沙中粒径较大的砾石在初始阶段较低也是由于其质量较大，水流产生的拖拽力小于泥沙颗粒与坡面的摩擦力，无法使土粒在坡面上发生滑动或滚动而发生输移。而中度侵蚀土体和严重侵蚀土体在降雨前期较短的时间变化较大，此后有波动但整体变化趋势并不大。这是由于随着土体侵蚀程度的加深，其抗侵蚀能力相应下降，细沟在极短的时间产生，集中水流功率迅速增加，分选作用大大下降。

随着降雨-冲刷强度增加砾石和黏粒含量均增加，砂粒和粉粒含量减少（表 3）。这主要是因为更高的降雨-冲刷能量具有更大动能带动大颗粒的运移，同时也会增加细沟间区域的侵蚀强度，引起黏粒和粉粒含量的改变。随着降雨-冲刷强度增加，不同粒径颗粒增加比例也存在差异。高降雨-冲刷条件下E1、E2、E3砾石含量分别增加了13.61%、7.54%、13.60%，而黏粒分别增加了6.73%、0.39%、1.48%。其差异表明高降雨-冲刷条件对细沟侵蚀作用的影响作用明显高于细沟间侵蚀。砂粒含量减少较大，通过对砂粒粒径细化可以看出，其改变的主要原因是细砂粒（0.05~0.25 mm）含量的减少和粗砂粒（0.5~2 mm）含量的增多，使其分布更趋近于原土。高降雨-冲刷条件下细沟会更早的形成，细沟水流功率增加，较大粒径的粗砂粒侵蚀量增加，而水层厚度的增加削弱了雨滴击溅的影响，导致细砂粒搬运量减少。泥沙中各侵蚀程度砂粒含量均很高，这也与原土砂粒含量较高有关。

对比两个降雨-冲刷强度下各侵蚀程度土体泥沙的富集率可以发现：在较低降雨-冲刷强度下，泥沙中黏粒和粉粒的富集率较高、砂粒和砾石富集率均较低（表 3），这与杨伟等^[37]的研究结果一致。细颗粒物质的富集和土壤表面团聚体的破碎程度有关^[38]，在较低降雨-冲刷强度下，受到雨滴动能的影响，大粒径颗粒被剥蚀为较小粒径颗粒，小颗粒以悬移的方式流失，而大颗粒较难被搬运。这间接表明粗颗粒在坡面地表富集，这一现象与在野外花岗岩红壤崩岗侵蚀区集水坡面观察到的现象一致。在更高的降雨强度下砾石富集率增幅较大，E1、E2、E3分别增加了2.52倍、3.96倍、31倍，其中E1砾石富集率接近1，而E2、E3则远大于1。这一现象可能归因于泥沙输移机制的差异（悬移、跃移和推移过程）^[39]，粒径较大的颗粒通过推移被搬运，粒径较小的颗粒通过悬移-跃移被搬运^[40]。当降雨-冲刷强度增加后，细沟更快的发生发育，水流功率及挟沙能力的增加，明显增大了推移搬运的比重。

2.3 不同降雨强度和上方来水流量下泥沙颗粒分形维数变化

已有研究表明土壤颗粒分形维数越高，土壤细颗粒化越明显，土壤质地相对会更好^[41]。在泥沙颗粒分选研究中引入分形维数的概念，有助于更好地表征不同粒径颗粒的变化。泥沙颗粒质量分形维数与其黏粒、粉粒和砂粒质量分数之间的关系如图 6所示。由图可知，分形维数与黏粒和砂粒质量分数之间关系均可以采用线性模型进行较好拟合（P < 0.01），其中，分形维数与黏粒质量分数之间为极显著正相关关系，与粉粒和砂粒质量分数之间为负相关关系（其中与砂粒质量分数呈极显著负相关），泥沙颗粒分形维数主要由黏粒的质量分数决定，这与刘淼等^[42]、王国梁等^[43]的研究结果相符。因此，暴雨驱动下分形维数相对较高反映了泥沙中砂粒比重的下降与黏粒比重的升高，可以作为评价不同侵蚀程度坡面侵蚀泥沙的颗粒粒级分布特征的指标。

图 7为不同降雨-冲刷条件下泥沙颗粒分形维数随产流历时的变化。E1泥沙质量分形维数在较低降雨-冲刷强度下随着产流历时的延长逐渐降低，当降雨-冲刷强度增加后分形维数随产流历时表现为先下降后趋于平稳。在90 mm·h^–1降雨强度和2 L·min^–1冲刷流量耦合条件下，初期雨滴击溅对土粒的分散作用较强，坡面薄层水流的搬运和层状剥蚀作用使得侵蚀泥沙中细颗粒增多。随着降雨历时延长，坡面粗颗粒露出减少了雨滴击溅的影响，分形维数随之逐渐减小。在120 mm·h^–1降雨强度和2.66 L·min^–1冲刷流量耦合条件下，雨滴击溅分散土壤颗粒的作用更强，致使细颗粒增多，分形维数相较于较低降雨-冲刷下更大；同时坡面细沟开始发育后，集中水流对泥沙颗粒分选作用反而降低^[16]，分形维数随之稳定。E2、E3侵蚀颗粒分形维数在试验初期有下降趋势，其后变化不明显，随着降雨-冲刷强度变化分形维数增加不大。其原因是E2、E3相对于砾石含量更少，抗侵蚀能力更低，两个降雨-冲刷条件下均在较早时间内发育出细沟。分形维数存在波动，是因为沟岸坍塌造成了泥沙颗粒分布有较大变动。其中90 mm·h^–1降雨强度和2 L·min^–1冲刷流量耦合条件下E3侵蚀泥沙分形维数在12~18 min增加，原因是试验条件为降雨-冲刷，抗侵蚀能力较弱的E3不仅存在细沟的溯源侵蚀，坡面顶部在水流冲刷作用下也会出现跌坎，当沟头与跌坎汇合加剧水流对细颗粒的侵蚀，导致分形维数增加，随着细沟的稳定发育，分形维数逐渐降低。

3 结论

随着降雨-冲刷强度的增加，土壤侵蚀速率、泥沙浓度和产流率均增加，轻度侵蚀土壤受降雨-冲刷强度影响最大。不同侵蚀程度土体间，土壤侵蚀速率表现为E3 > E2 > E1，而产流率表现为E2 > E1 > E3，通过多元逐步回归分析发现降雨强度越大，黏粒含量越高、粉粒含量越低，土体侵蚀速率越小，而容重越大产流率越高。随着降雨-冲刷强度和降雨历时增加，泥沙中砾石和黏粒含量均增加，砂粒和粉粒含量减少。黏粒和粉粒的富集率较高、砂粒和砾石富集率较低，随着降雨-冲刷强度增加砾石富集率增幅较大。E1泥沙质量分形维数随产流历时降低，随降雨-冲刷强度增加而增大。E2、E3侵蚀颗粒分形维数下降集中在试验初期，细沟发生后其变化不明显，随着降雨-冲刷增加分形维数增加不明显。