陈 鹏, 张铁钢, 董 智, 郭建英, 李红丽, 仇苏倩

(1.山东农业大学 林学院, 泰山森林生态系统定位研究站, 山东 泰安 271018;2.水利部 牧区水利科学研究所, 呼和浩特 010020)

砒砂岩区土壤侵蚀严重，水土流失现象剧烈，砒砂岩成岩程度较低，且遇水时极易膨胀，土壤颗粒组分中细颗粒含量较为匮乏，砂粒间胶结程度差，细小颗粒对大颗粒孔隙的阻塞、填充作用使得土壤孔隙结构遭到破坏，渗透性能严重降低[1]。这些特性导致砒砂岩分布区成为黄土高原治理难度最大、侵蚀强度最为剧烈的区域。合同沟小流域属于典型的砒砂岩强烈侵蚀区，生态环境较为恶劣，坡面土壤流失现象十分严重，是达拉特旗乃至鄂尔多斯多沙粗沙的代表区域[1]。虽然以往研究表明[2]，植被措施可以显著提升土壤抗蚀性，降低土壤的侵蚀程度，但对于水分承载力较弱的干旱、半干旱地区而言，植被的控流阻沙功能受到了严重挑战[3]。因此探究有限水资源支撑下的高效植被布设措施便显得尤为迫切。明晰不同条件下土壤侵蚀特征是研究土壤抗蚀性能、规避水土流失风险的重要前提。其中土壤剥蚀率是探究土壤侵蚀特征重要指标[4]，研究径流过程中土壤的剥蚀特征对表述水力侵蚀作用下坡面侵蚀过程以及水土流失评价具有重要意义。

流速、流量、剪切力和水流功率等水力参数是坡面土壤侵蚀的驱动力，对土壤剥离和泥沙输移过程影响显著[5-7]，其中水流功率、单位水流功率和水流剪切力作为探讨径流能量流和动力流的重要参数，在坡面土壤侵蚀过程和水蚀动力机理的研究中应用最为广泛[6-9]。揭示土壤剥蚀率与水力参数特征之间的关系对探究土壤侵蚀机制、合理改进土壤侵蚀预报模型具有重要意义[9]。但目前关于坡面水蚀的研究中，其侧重点较多的集中于坡度、雨强及冲刷流量等方面[2,9-10]，而利用水力参数，解析地表水文过程对不同植被格局下土壤侵蚀特征的研究则略显薄弱[11-12]。

基于此，本文选择砒砂岩区最为常见的灌草坡面为对象，选取无灌丛草地(GL)、灌丛位于坡上部(SU)、坡中部(SM)和坡下部(SL)4种灌草格局，以裸坡(BL)为对照，分析不同灌草格局下土壤剥蚀率和水动力参数的特征及其相互关系，以期揭示砒砂岩区不同灌草格局下土壤侵蚀特征的差异性及其作用机理，为合理完善水蚀预测模型提供科学理论与技术依据。

1 材料与方法

1.1 研究区概况

试验地点位于黄河一级支流罕台川的合同沟小流域，行政区划属于内蒙古达拉特旗，地理位置为39°59′58″—40°13′18″N，109°53′36″—110°06′53″E，受温带半干旱大陆性气候的影响，该区多年平均降雨介于200～400 mm，71.2%的降水集中于夏季(7—9月)，降水过程表现为“短历时、大雨强”的特征。监测最大降水量达223 mm/h，且土壤类型以栗钙土为主，质地为砂土，坡面坡度介于5°～15°。区域内主要有柠条(CaraganaKorshinskii)、油松(Pinustabuliformis)、小叶杨(Populussimonii)等造林树种，以及本氏针茅(Stipabungeana)、达乌里胡枝子(Lespedezadavurica)、羊草(Leymuschinensis)等建群植物，其中坡面以柠条种植面积最大。

1.2 试验装置

根据达拉特旗1960—2016年气象资料记录的实际降雨情况，以及该区“短历时，大雨强”的降水特征，选择1.5，2.0，3.0 mm/min雨强，根据径流小区面积(5 m×2 m)换算后确定试验冲刷流量为15，20,30 L/min。以该区平均坡度8°作为本研究的试验坡度。根据柠条灌丛在坡面分布的差异性，分别设置无灌丛(即单一草地坡面，GL)、灌丛位于坡上部(SU)、坡中部(SM)、坡下部(SL)和裸坡(BL)共5种试验径流小区，小区建造后进行围封，经3个月自然沉降后进行试验。各小区基本概况见表1。

表1 径流小区基本概况

本研究中径流小区的规格为5 m×2 m(按水平投影计)，试验之前使用流量计对流量进行率定，由水泵将自来水泵入供水稳流装置中，消除水流扰动，小区与溢流槽的交界处放置一个2 m×0.2 m的钢制平板，以保证水流在较低的动能水平下均匀而平稳的溢入小区。预试验发现产流经过60 min径流基本趋于稳定，因此设定每场试验时间为60 min，径流流出小区时开始试验计时。试验前期每分钟收集一次径流样品，10 min后每5 min收集一次径流样品，同步测定流速、水宽和水深，其中流速采用染料示踪法测定(本研究中染色剂为高锰酸钾)，染料前缘速度乘以修正系数0.8，得到平均流速[4]，水宽和水深采用精度为1 mm的直尺测量，将泥沙样品风干去除根系等杂质后测定重量，用以土壤剥蚀率的计算。不同处理间进行3次重复，试验数据取其平均值，共计45场。每次试验之前测定各径流小区的土壤含水量，消除土壤水分条件的差异，保证试验结果的有效性。

1.3 数据分析

1.3.1 土壤剥蚀率计算 土壤剥蚀率(Dr)采用下式计算[13]：

(1)

式中:Dr为土壤剥蚀率[g/(m2·s)];M为径流过程中的产沙量(g);b为径流平均水宽(m);L为径流长度(m),本研究中L为5 m;T为径流时间(s)。

1.3.2 水动力参数计算 水流功率(ω)、单位水流功率(P)和水流剪切力(τ)采用下式计算[13]：

ω=τV

(2)

P=VJ

(3)

τ=ρghJ

(4)

式中:ω为水流功率(W/m2);P为单位水流功率(m/s);τ为水流剪切力(Pa);V为径流平均流速(m/s);J为水力坡度;ρ为水的密度(kg/m3);g为重力加速度(g/m);h为径流深度(m),本研究使用坡面坡度的正弦值近似代替J[13]。

本研究使用Excel 2016进行试验数据的计算处理，SPSS 24.0对进行各指标的统计分析。采用Origin进行图件绘制。

2 结果与分析

2.1 不同灌草格局下土壤剥蚀特征

本试验条件下径流平均宽度整体介于0.995～0.424 m(图1)，并呈现出随冲刷流量增大而逐渐增大的趋势。灌草格局对径流宽度具有一定影响，但不同冲刷流量下，随灌草格局的改变径流宽度的大小并未呈现出一致的变化趋势。

图1 不同灌草格局下径流平均宽度特征

灌草格局对土壤剥蚀特征具有显著影响(p<0.05)。不同格局土壤剥蚀率整体介于0.055～0.469 g/(m2·s)(图2)，并且随冲刷流量的提高，土壤剥蚀率均表现出增大的趋势。3种冲刷流量下土壤剥蚀率均以BL格局最大，其中BL格局土壤剥蚀率介于0.277～0.469 g/(m2·s)，均值为0.386 g/(m2·s)。其他格局土壤剥蚀率的范围介于0.055～0.359 g/(m2·s)，土壤剥蚀率均值仅为BL格局的44.63%，表明不同灌草格局均能够显著降低土壤剥蚀率。不同灌草格局间土壤剥蚀特征具有一定差异，GL格局土壤剥蚀率显著高于SU，SM和SL格局，3种冲刷流量下，不同灌草格局表现出一致的趋势，土壤剥蚀率的大小关系均为GL>SU>SM>SL，其中GL格局土壤剥蚀率均值为分别是SU，SM，SL格局的1.55，1.77,2.73倍。SL土壤剥蚀率均值仅为0.098 g/(m2·s)，是降低土壤剥蚀率的最佳灌草布设措施。

图2 不同灌草格局下土壤剥蚀特征

2.2 不同灌草格局下水动力参数特征

不同灌草格局条件下水流剪切力介于15.458～22.750 Pa，与BL格局相比(21.503 Pa)，各灌草格局均能显著降低水流剪切力，其中SL格局水流剪切力最低，仅为BL格局的78.99%，但随植被格局的改变水流剪切力的变化趋势并不明显。不同格局下水流功率介于1.727～5.830 W/m2，随冲刷流量的增大，不同格局间的水流功率显著增加。15，20，30 L/min冲刷流量时，BL格局水流功率分别为3.967，4.483，5.830 W/m2，是其余灌草格局的1.80，1.56，1.71倍，不同灌草格局可以发现水流功率呈现BL>GL>SU>SM>SL的分布趋势。单位水流功率介于0.014～0.055 m/s，随冲刷流量的提高，单位水流功率呈增大趋势。BL格局单位水流功率均值为0.047 m/s，是其余灌草格局的2.29倍，不同灌草格局单位水流功率分布特征与水流功率相同(图3)。

图3 不同灌草格局下水动力参数特征

2.3 土壤剥蚀率与水动力学参数特征的关系

水流剪切力、水流功率和单位水流功率与土壤剥蚀率均呈正相关关系，其中水流功率与土壤剥蚀率的线性函数Dr=0.1275ω-0.1881及幂函数Dr=0.0167ω2.176的决定系数均达到了0.92以上，即两者均可以对该区土壤剥蚀特征进行较好的预测，其次则是单位水流功率的线性函数以及水流剪切力的幂函数(R2=0.912，0.661)。由回归分析结果可知，水流功率是预测土壤剥蚀率的最佳水动力参数，水流剪切力的R2值分别只有0.572，0.661，对土壤剥蚀率的描述效果不佳(图4)。

图4 土壤剥蚀率与水动力参数的关系

通径分析是研究自变量对因变量的直接及间接重要性多元统计技术[13]，可以为统计决策提供可靠依据。通径分析结果表明，水流功率和剪切力对土壤剥蚀率的直接、间接通径系数均为正数，贡献度分别为86.85%，8.33%，单位水流功率对土壤剥蚀的贡献度最低(4.81%)，水流剪切力的间接通径系数达到了0.817。说明水流功率和水流剪切力均对土壤剥蚀率起增进作用，水流功率是土壤剥蚀率的主要影响因素，水流剪切力的间接作用在土壤剥蚀过程中起重要影响作用(表2)。

表2 不同灌草格局下水动力参数与土壤剥蚀率的通径分析

3 讨 论

3.1 灌草格局对土壤剥蚀特征的影响

雨滴或水流驱动下引起的土壤剥离过程是外界侵蚀营力与土壤抗蚀能力间相互作用的外在表现，当侵蚀应力大于土壤的抗蚀能力时就会产生土壤侵蚀[14]。以往研究表明，植被影响是坡面土壤剥蚀特征的重要因素之一[8]，Zhang等[8]研究认为灌丛及草本植被均能有效的降低径流过程中坡面土壤的剥蚀程度，坡耕地的土壤剥蚀率是不同植被条件下的6.9～47.8倍。这在本研究中也得到了证实。本研究发现植被类型的差异也是影响土壤剥蚀特征的重要因素，其中SU，SM，SL格局土壤剥蚀率显著低于GL格局，这是由于较草地而言灌丛的存在增大了坡面的根系密度，根系的缠绕及其分泌物对土壤颗粒的胶结作用，可以固持土壤、提高表土及附近土壤的稳定性[15]；根系提升土壤有机质的作用，可以促进水稳性团聚体的形成[6]，提升表土粗糙度，降低径流对土壤的剥蚀；其次，在试验过程中观测发现，在研究区风水复合的特殊的侵蚀环境下，灌丛所在区域形成了凸起的丘体[16]，径流通过灌丛所在位置时均出现绕过丘体向坡下汇集的现象。丘体对地表径流的拦蓄作用可以在一定程度上降低径流的流态，削弱径流对坡面的剥蚀[17]。从而在对土壤抗蚀性能的提升效果方面灌草坡面优于草地。

进一步分析发现，3种冲刷流量下土壤剥蚀率均呈现SU格局>SM格局>SL格局的趋势。对此Zhang等[2]认为，径流对坡面的侵蚀主要集中于上坡下部到下坡下部的广泛区域，径流过程中坡面侵蚀主要发生在中下部位置，因此相较于灌丛位于坡上和坡中部而言，灌丛位于坡面下部时更有利于坡面重点侵蚀区域的保护；并且在干旱半干旱区，由于气候的原因有机质分解速率较慢，坡面存在一定量的凋落物组织，并在径流作用下形成一种类似“坝体”的结构[18]，此类屏障的形成对径流路径及运动状态会产生显著影响。本研究中灌丛所在位置的凋落物明显高于坡面其他部位，未分解的茎秆与小的凋落物组织及植物根茎部位结合形成足以拦截径流的“自然屏障”，极大的降低了径流速率，消耗水流能量，削弱了灌丛附近土壤剥蚀的强度[14]。因而在多种因素的共同作用下，土壤剥蚀率表现SU格局>SM格局>SL格局的分布特征。

3.2 水动力参数对土壤剥蚀率影响

水流功率、单位水流功率和水流剪切力是土壤侵蚀模型中描述土壤剥离过程常用的水动力参数[8,19-20]。运用水动力参数对不同灌草格局下土壤剥蚀率进行量化，有利于建立适用于本研究区土壤剥蚀率的预测模型。Wang等[7]研究发现，土壤剥蚀率与水流功率、单位水流功率和水流剪切力均呈正相关关系，这与本研究结果相符，不同灌草格局下水流功率和单位水流功率均表现为SL格局

土壤剥蚀率与水动力参数的变化趋势基本相同，说明不同植被措施均能够通过降低水流动能来削弱坡面的剥蚀程度[7]，但在土壤剥蚀特征的预测研究中，最佳水动力参数选择仍具有一定争论，Zhang等[21]认为水流功率是估算土壤剥蚀率的最有效水动力参数，而一些研究[19-20]却发现，水流剪切力、单位水流功率是土壤剥蚀率的最佳预测参数，对此Wang等[12]认为，研究条件的变化可能会导致最适宜水力参数的选取出现差异。针对本研究而言，简化式中水流剪切力是水深和坡度的函数[13]，本研究试验地点为野外自然坡面，与室内试验相比，自然坡面在多种因素干扰下平整度较差，径流填洼的现象会使水深对灌草格局及冲刷流量的反馈降低，导致不同灌草格局下水流剪切力与土壤剥蚀率拟合精度较差；此外本研究条件下，不同格局冲蚀细沟的发育也是造成差异的重要原因，Stefanovic等[22]研究认为侵蚀泥沙的来源并非全是由水流剪切作用对坡面土壤的剥离作用产生的，细沟发育过程中的沟岸坍塌、沟头溯源等活动作为细沟发育的重要组成部分[22]，也是泥沙的重要贡献来源[23]，而水流剪切力的剪切效应主要作用于细沟底部，这在一定程度上导致了水流剪切力对土壤剥蚀现象响应的敏感性降低；对于水流功率、单位径流功率等能量指标而言，沟岸崩塌、沟头溯源等现象可以影响水流剥离搬运的能量分配，因而对土壤剥蚀特征的响应更为灵敏。但对水动力参数和土壤剥蚀率进行通径分析后发现，水流剪切力的间接通径系数达到了0.817，且对土壤剥蚀率贡献率仅次于水流功率，说明在研究坡面径流对土壤的侵蚀作用时，水流剪切力对土壤剥蚀率的间接作用也不应被忽视。

4 结 论

(1) 植被可以通过降低水流能量来削弱径流对坡面的剥蚀程度，3种冲刷流量下，砒砂岩区坡面土壤剥蚀率介于0.055～0.469 g/(m2·s)，并且呈现随流量的增大逐渐上升的趋势。其中SL格局是降低坡面土壤剥蚀率的最佳植被布设措施。

(2) 土壤剥蚀率、水动力参数均值整体呈BL格局>GL格局>SU格局>SM格局>SL格局的趋势。水流功率是描述土壤剥蚀率的最佳水动力参数，水流剪切力对土壤剥蚀率的间接作用也不应被忽视。