南方红壤区是我国水土流失的严重区域之一，其中面积在2.00×10⁵ km²以上的花岗岩区尤为突出。南方高温湿润的气候特点有助于母岩风化，形成深厚的风化壳^[1]。日益紧张的人地关系、滥砍乱发、无序开发建设等人为因素诱导、加速了低山丘陵区的土壤侵蚀。花岗岩结构松散，富含石英砂砾，具有不同于其他母质发育的土壤的特殊型。一旦花岗岩红壤坡面发生土壤侵蚀，坡面细沟侵蚀极易演变成沟蚀，形成危害巨大的“生态溃疡”——崩岗侵蚀。这类侵蚀治理难度大，会毁坏房屋道路，沙压农田，淤积水库塘堰，严重影响当地社会、生态、经济的发展和安全^[1]。从崩岗的侵蚀类型演变来看，崩岗大多数由地面径流形成的细沟侵蚀发育而来，是坡面侵蚀沟谷发育的高级阶段^[1-2]。降雨达到地表后，雨滴和薄层径流呈点状、片状或面状地剥蚀表面土壤。当坡面形成径流后，径流对花岗岩风化残积土产生下切侵蚀，在切穿表土后，坡面形态由无定型的细沟发育到永久性、固定型的浅沟、切沟和冲沟^[1]。尤其是在植被被破坏后，裸坡演化为沟蚀的速度更迅速。因此，花岗岩地区红壤坡面侵蚀的发生机制和防治工作近些年来一直是学者们的研究热点^[3]。如蔡崇法等^[4]研究了花岗岩红壤表土特征及对坡面侵蚀的影响，认为花岗岩红壤稳定的表土是一个上部多砾石、下部较紧密的紧实层，紧实层的存在影响地表入渗和坡面产流产沙量。张德谦等^[5]研究了花岗岩红壤各层次土体的坡面侵蚀特性和各粒径泥沙的分选规律，指出淋溶层具有较强的抗侵蚀能力，在侵蚀过程中坡面细颗粒流失严重。对于花岗岩风化岩壳上粘下砂、上紧下松的特点，表土对于表下层土壤具有重要的保护作用，一旦表土被切穿、剥离，表下层极易发生更严重的侵蚀^[5]。近年来，生物工程防治、坡面削坡开级、坡面覆盖物的生态治理等新兴措施逐渐兴起，越来越发挥重要的作用。因此，部分学者将红壤坡面侵蚀防治的研究重心集中在坡面覆盖因子的影响上，合适的覆盖物能够有效地减少坡面的土壤侵蚀和径流^[6-8]。目前，已有的研究中使用的坡面覆盖物主要为秸秆覆盖^[6-7]、砾石覆盖^[8]等措施。黄麻土工布是用低等麻或麻纺下脚料织成的网状物，是国内一种常见的农业工程材料。近些年来，土工布逐渐被引入了坡面水土保持技术中，在防治工程边坡侵蚀方面有一定应用前景^[9]。它含有良好的天然纤维，具有良好的透水性、隔离性、防渗防护，工程效益见效快，能迅速建立植被，在改善土壤结构，拦蓄径流泥沙等方面有巨大优势^[9-10]。张勇等^[11]报告了花岗岩红壤在黄麻土工布覆盖影响下的坡面侵蚀特性，包括坡面产流产沙、泥沙粒径分选等。然而，该试验只考虑了有黄麻覆盖与无黄麻土工布覆盖条件及陡坡条件，在不同坡度梯度和网格大小的黄麻土工布覆盖条件下的坡面侵蚀特性尚未见报道，尤其是花岗岩红壤表土的坡面侵蚀特性还需进一步的研究。基于以上研究背景，本试验以花岗岩红壤表土层为研究对象，通过室内模拟降雨试验的方法，探究黄麻土工布控制条件下花岗岩红壤的坡面侵蚀特性，以期为花岗岩红壤区坡面侵蚀的防治工作提供一定的科学参考。

1 材料与方法 1.1 供试土壤

采样点位于湖北省咸宁市通城县五里镇(29°12′39″N，113°46′26″E)，该地区气候属于北亚热带季风气候，气候温和，年平均气温为17.1 ℃，无霜期为260 d左右，平均年降雨1 604 mm。通城县广泛出露着花岗岩母质发育的松散土体，花岗岩占全县地层岩石面积的80%以上。采样点的土地利用方式为林地，主要植物有马尾松、金樱子、铁芒萁等。在前期野外调查的基础上，选取了典型土壤剖面，并将其上部0—15 cm厚度的表层土壤采集并运回试验室。同时，3个100 cm³环刀原状土样被采集用于测定自然状态下土壤容重。将土样自然风干并过10 mm筛网，剔除石块、根系等杂质。土壤理化性质采用常规方法测定^[12]，土壤容重和土壤含水率通过烘干环刀样获得；pH值通过1:2.5土水比混合浸提后，电极法测定；土壤的机械组成(国际制)采用湿筛法和吸管法测定；有机质为重铬酸钾外加热法。上述试验每组重复3次进行，结果见表 1。

1.2 试验方法 1.2.1 试验设备

模拟降雨通过自动模拟降雨系统提供，该降雨器通过控制供水压力和不同的喷头组合可以实现15~220 mm/h降雨强度范围，且降雨均匀度＞90%，降雨高度为10 m，雨滴到达坡面表面时可以达到终点速度。试验土槽规格为2.0 m(长)×1.2 m(宽)×0.4 m(高)，坡度可在0°~30°自由调节(见图 1)。

1.2.2 土槽制备

在填装土壤前，将土槽底部填装15 cm鹅卵石与细沙以便坡面自由排水。根据野外采集的环刀测定的容重，将试验土壤填入土槽。土壤采用分层填装的方法，以每层5 cm的增量进行填装，逐层填至40 cm。为保证土壤颗粒分布均匀、贴合紧密，在填入上层土壤前将下层土壤刮毛处理。每组试验重复进行2次。土槽完成填装后，在土壤表面覆盖一层纱布以减少预湿润降雨对坡面的影响。试验开始前以30 mm/h的雨强对坡面进行预湿润处理，直至坡面开始产流，然后将其放置24 h自由排掉重力水。当通过TDR(time domain reflectometry)测定的坡面土壤含水率将至20%时，开始进行降雨试验。

1.2.3 试验设置

试验于2018年3月至2019年3月在华中农业大学降雨大厅进行，每个处理试验重复2次。针对我国亚热带气候区典型的强风暴^[13]，试验设置了90 mm/h的降雨强度。丁光敏等^[14]对福建崩岗侵蚀成因的调查研究发现10°~25°是最易产生侵蚀的坡度，同时也是崩岗侵蚀区集水坡面的常见坡度，因此，试验设置了缓坡5°和陡坡15°的2个坡度。本试验中使用的黄麻土工布与张勇等^[11]的材料一样，黄麻直径为5 mm。依据坡面覆盖度和黄麻土工布的网格大小，试验依次设置了对照无覆盖组(CK)、6 cm×6 cm组(T₁)和3 cm×3 cm组(T₂)3个处理(图 1)。

1.2.4 观测指标和数据处理

在降雨试验开始前，将土槽调整至目标坡度，并将降雨强度校准到目标强度(误差＜5%)。当坡面开始产流后，试验开始计时，每场降雨持续60 min。试验过程中，通过放置在土槽四周的雨量筒实时监测降雨强度。降雨过程中通过铝盒与塑料桶收集泥沙样，采集间隔为2 min。将铝盒中的泥沙样称重、烘干后，结合塑料桶收集的泥沙样，获得坡面产流量、含沙量等数据。径流系数为坡面产流量与坡面实时降雨量之比。在本研究中，考虑到花岗岩红壤土壤团聚体的弱稳定性，选择最终粒径分布(分散处理)来表征泥沙的粒径分布，其通过湿筛法和吸管法得到。由于＞0.1 mm泥沙粒径在沉降筒中沉降太快，因此将上述烘干铝盒(每4 min)中的土壤转移进烧杯，移除有机质并分散后，过0.1 mm的筛网。粒径＞0.1 mm的泥沙过2 mm的筛网进一步筛分，粒径＜0.1 mm的泥沙通过吸管法获取泥沙的粒径分布。泥沙粒径被划分为黏粒(＜0.002 mm)、粉粒(0.002~0.05 mm)、砂粒(0.05~2 mm)及砾石(2~10 mm)。

本研究选取下列相关参数。土壤侵蚀速率〔g/(m²·min)〕为单位时间、单位侵蚀面积上土壤在坡面径流侵蚀力的作用下被剥蚀的土壤颗粒质量^[15]：

$ E = \frac{{{m_1}}}{{{\rm{BLT}}}} $

(1)

式中：E为土壤侵蚀速率〔g/(m²·min)〕；m_t为T时间内的累计产沙量(kg)；B为水面宽度(这里用土槽宽度代替)(cm)；L为坡长(m)；T为时间(min)。

泥沙颗粒的富集率，用于研究土壤侵蚀过程中泥沙中各粒级泥沙颗粒的富集现象^[6]，具体计算公式为：

$ {E_R} = \frac{{{P_s}}}{{{P_o}}} $

(2)

式中：E_R为富集率；P_s为各粒径泥沙颗粒的质量占泥沙总质量百分比(%)；P_o为相应粒径的土壤颗粒占原始土壤的百分比(%)。富集率＞1表示泥沙中该粒径颗粒比原始土壤中相应颗粒所占的比例更大，富集率＜1则表示泥沙中该粒径颗粒小于原始土壤中相应颗粒所占的比例。平均重量(MWD)直径(mm)为用以表征泥沙和土壤粒径分布的指标：

$ {\rm{MWD}} = \sum\limits_1^{n + 1} {\frac{{{r_{i + 1}} + {r_i}}}{2}} \times {m_i} $

(3)

式中：r_i为第i个筛子孔径(mm), r₀=r_i，r_n=r_n+1；m_i为第i个筛网中颗粒的质量分数；n为筛子的数量。

本试验中，数据处理通过Excel与Origin 2020软件完成。

2 结果与分析 2.1 坡面产流产沙过程

在本试验中，坡面侵蚀属于溅蚀—片蚀(细沟间侵蚀)过程的范畴，泥沙的分离与输移是由雨滴的击溅作用和降雨引起的坡面薄层径流所驱动。图 2表示不同试验条件下坡面径流系数随降雨历时的变化。由图 2可知，坡面径流系数随降雨历时呈逐渐增加的趋势。在降雨初期，雨水下渗延缓了径流的产生^[6]。随着降雨的持续，土壤含水率逐渐饱和，同时大团聚体破碎形成的细颗粒逐渐堵塞土壤孔隙，从而促进了径流的产生与增加，当坡面入渗率趋于稳定时，径流率也趋于稳定^[5]。另外，坡面径流量随着坡度的增加呈现增加的趋势，而随黄麻土工布覆盖密度(无覆盖CK至6 cm×6 cm网格T₁，3 cm×3 cm网格覆盖T₂)的增加呈现减小的趋势。如在T₁覆盖条件下缓坡与陡坡的平均径流系数分别为0.56和0.46，而当坡面覆盖由无覆盖(CK)增加到T₁，T₂时，缓坡条件下坡面的平均径流系数分别为0.64，0.56及0.49。

如图 3所示，本试验采用了土壤侵蚀速率对坡面侵蚀产沙进行了定量化表达。在各试验条件下，土壤侵蚀速率随降雨历时呈现逐渐减小的趋势。结合上述坡面径流随降雨历时逐渐增加的规律，表明坡面侵蚀过程是一个侵蚀分离受限的过程。和径流系数的规律相似，土壤侵蚀速率随着坡度的增加而增加。坡度的增加一方面使坡面土体稳定性降低，另一方面增加了径流在坡面方面的重力分力，使其侵蚀力和能量增加，更易侵蚀及搬运土壤^{[6, 16]}。同时，当坡面覆盖由无覆盖(CK)增加到T₁，T₂时，土壤侵蚀速率逐渐减小，如在陡坡条件下3种覆盖处理的坡面土壤侵蚀速率分别为55.05，38.20及22.92 g/(m²·min)，相比之下，2种黄麻土工布处理使土壤侵蚀分别减少了28.0%和56.8%。本试验中尽管与其他研究所使用的土壤类型和覆盖物不同，但同样达到了相应的减流减沙效果。

2.2 泥沙分选过程

整体上，在缓坡和陡坡条件下各粒径泥沙颗粒的变化呈现了相似的规律(见图 4)。泥沙中黏粒和粉粒颗粒的含量随黄麻土工布覆盖密度的增加而增加，而砂粒和砾石颗粒的含量则相反。各粒径颗粒随降雨历时的变化则表现出相反的规律，黏粒和粉粒颗粒的百分比均随降雨历时而逐渐减少，粒径较大的砂粒和砾石颗粒的百分比均随降雨历时而逐渐增加。其原因是在坡面片蚀阶段，坡面的土壤被薄层径流层状剥蚀，径流搬运能力不足，悬移是主要的搬运机制，细颗粒倾向于被径流优先选择性搬运^[7]，石英粗颗粒由于其质量较大而难以搬运。然而，随着降雨的持续，一方面土壤表面细颗粒逐渐耗尽，另一方面土壤表面形成结皮，侵蚀受到限制，径流能量逐渐更多地消耗在粗颗粒的搬运上^[16]。虽然，随着降雨的持续粗颗粒含量在逐渐增加，但泥沙中仍以细颗粒为主。

各试验条件下侵蚀泥沙颗粒的富集率变化如图 5所示。从图 5可以看出，泥沙中黏粒和粉粒颗粒呈现富集的规律，而砂粒和砾石颗粒呈现亏损的规律。这一结果间接表明随着降雨的持续，坡面土壤中的石英粗颗粒逐渐倾向于被留在坡面，并不断在坡面富集形成一层粗粒覆盖层，这与我们在野外观察到的现象一致。另外，在本试验的短时强降雨条件下，缓坡更有利于泥沙中细颗粒的富集，坡度的增加使泥沙趋于粗化(见图 5)。如T₁处理下，当坡度由5°增加到15°时，泥沙中黏粒的富集率由1.27减少到了1.17，而砂粒的富集率由0.60增加到了0.74。坡度的增加使难以搬运的粗颗粒开始倾向于被搬运，加之径流能量增加，导致泥沙中粗颗粒所占的比例增加，进而造成细颗粒的相对含量下降^[17]。

黄麻土工布的存在明显有利于泥沙中细颗粒的富集(见图 5)。随着黄麻土工布覆盖密度的增加，泥沙中细颗粒的富集率呈现增加的趋势，粗颗粒的富集率则相反。在缓坡条件下，当黄麻土工布覆盖密度增加时，黏粒的富集率由1.19增加到了1.27，1.38。和无覆盖试验相比，即使在陡坡条件下，黄麻土工布对径流搬运的大颗粒也具有明显的拦截效应，其结果导致石英粗颗粒更多地富集在坡面表面。泥沙颗粒的平均重量直径(MWD)提供了更直观的证据(见图 6)。泥沙颗粒的MWD越大，说明泥沙越粗糙。对MWD分析可知，缓坡条件下，3个处理的泥沙MWD平均值分别为0.37，0.30，0.22 mm，陡坡条件下泥沙MWD平均值分别为0.44，0.37和0.31 mm，表明泥沙随坡度的增加呈现越加粗糙的规律，随着黄麻土工布覆盖密度的增加呈相反的规律。和无覆盖处理相比，黄麻土工布覆盖处理进一步使土壤中的粗粒富集在坡面，对坡面的侵蚀过程具有重要影响。

3 讨论

水力侵蚀涉及了由降雨溅蚀和地表径流引起的土壤分离、泥沙输移、泥沙沉积3个过程^[7]。雨滴和薄层径流的能量用于克服土壤颗粒之间的黏结力，并将剥离的土壤颗粒运输到其他位置^[18]。在本试验中，坡面径流量和土壤侵蚀速率随着黄麻土工布覆盖密度(CK至T₁，T₂)的增加而减小。试验中使用的网格状黄麻土工布与Shi等^[7]的研究中使用秸秆覆盖措施的效果相似，他指出当坡面覆盖率分别增加12.7%和86.6%时，坡面稳定产流率分别减少了8.2%和75.5%，同时坡面覆盖极大地减少了土壤侵蚀，尤其是在高密度覆盖条件下坡面含沙量减少了20多倍。Lin等^[6]报告了在高密度覆盖条件下由于入渗的增加导致土壤侵蚀量出现了转折。究其原因，一方面黄麻土工布能够增加土壤表面糙度，增加入渗能力，减少径流量及流速，削减径流能量，从而减少径流的侵蚀能力及挟沙能力；另一方面黄麻土工布能够削减、分散降雨能量，阻止雨滴直接打击土壤表面，减少土壤颗粒的侵蚀和破碎，抑制土壤结皮的发生与发展^[8]。

坡面的侵蚀过程是一个泥沙颗粒的分选过程，泥沙的粒径分布与富集率可以为坡面的原位和异位侵蚀提供更多的信息^[7]。泥沙中黏粒和粉粒颗粒的含量随黄麻土工布覆盖密度的增加而增加，而砂粒和砾石颗粒的含量则相反。黄麻土工布覆盖密度的增加不仅导致了坡面径流量的减少，也增加了土壤表面糙度，致使径流剥蚀、搬运泥沙时需要消耗更多的能量，使易蚀的细颗粒被优先剥蚀、搬运。然而，泥沙中黏粒和粉粒颗粒的含量通常高于原始土壤，砂粒和砾石颗粒的含量低于原始土壤，这个结果直接反映在各粒径颗粒的富集率上：黏粒和粉粒颗粒的富集率＞1，而砂粒和砾石颗粒的富集率＜1。这与Lin等^[6]采用崩积堆土壤获得的研究结果相似。尤其是，黄麻土工布覆盖密度对泥沙中细颗粒的富集有积极的作用。黄麻土工布覆盖密度越高，泥沙中黏粒和粉粒颗粒的富集率越高，相应地，石英粗颗粒随之更多地富集在坡面表面。Shi等^[7]指出在坡面秸秆覆盖率较低时(≤30%)坡面分离受限机制占主导地位。在本试验中，分离受限机制和黄麻土工布对降雨和径流能量的削减有关，黄麻土工布覆盖密度越高，径流能量越小，越加难以搬运质量较大、启动难度大的石英粗颗粒，这与前文的分析结果一致。研究^[8]表明，砾石覆盖对土壤的物理性质起重要作用。石英粗粒覆盖层的存在类似于坡面覆盖物的作用，不仅能够直接削减降雨和径流能量，减少雨滴的直接打击作用，还增加了坡面糙度，能够拦截被侵蚀的泥沙，尤其是泥沙中的大颗粒^[7]。对于花岗岩红壤这种特殊类型的土壤，石英粗颗粒的不断富集也是造成坡面侵蚀分离受限机制、土壤侵蚀速率随降雨的持续而逐渐降低的原因(见图 3)。相比裸坡，黄麻土工布覆盖导致坡面原位在后续侵蚀过程中有更高的侵蚀抗性，并形成一个正反馈过程。但是，原位产生的含沙量少的径流有可能会增加下坡异位的侵蚀。

4 结论

(1) 各试验条件下，坡面径流系数和土壤侵蚀速率随坡度的增加呈现增加的趋势，随黄麻土工布覆盖密度的增加呈现减少的趋势。土壤侵蚀速率随降雨历时的增加呈现减小的趋势，径流系数则相反。黄麻土工布覆盖的存在能够有效地减少坡面径流量与土壤侵蚀速率。

(2) 侵蚀泥沙中黏粒和粉粒颗粒的百分比均随降雨历时逐渐减少，砂粒和砾石颗粒的百分比则相反，同时泥沙中黏粒和粉粒颗粒的含量通常高于原始土壤，砂粒和砾石颗粒的含量低于原始土壤。黏粒和粉粒颗粒倾向于被优先输移。

(3) 不同粒径大小的泥沙颗粒的输移过程均具有选择性，泥沙中黏粒和粉粒颗粒呈现富集的规律，而砂粒和砾石颗粒呈现亏损的规律。这间接表明坡面发生了粗颗粒的富集，尤其出现在在缓坡与高密度黄麻土工布覆盖条件下。泥沙颗粒的MWD直观地证明了上述结果。