土壤生物结皮是由生长在土壤表面及其以下的细菌、真菌、苔藓和地衣类等个体微小的生物成分与土壤相互作用而形成的复杂复合土层^[1]，具有重要的生态功能。生物结皮一方面可改变表层土壤黏结力，促进土壤团聚体形成，提高土壤的抗侵蚀能力^[2-3]，另一方面可改善坡面入渗，影响水动力特性进而对坡面侵蚀产生影响^[4-5]。如生物结皮可增加表层斥水性，形成不透水层，影响产流过程^[6-7]，改变坡面流流态，降低坡面流流速^[8]，进而减少土壤侵蚀量。黄土高原实施退耕还林(草)工程后，以藻和苔藓为主要组成物的生物结皮得以广泛发育，其盖度可达70%左右^[9]。不同类型生物结皮组成和表面微形态具有差异，藻结皮黏连作用形成无机层使表面平滑，苔藓结皮植物体密集丛生^{[2, 10]}，不同类型生物结皮形成的下垫面不同，其水动力特性具有差异。此外，踩踏、冻融交替等均是自然界常见干扰，适度干扰及冻融交替可缓解区域水分状况，但也使原有结构在一定程度受到破坏，边缘翻起^[11-12]。研究^[13]表明生物结皮在11 Pa水流剪切力下会出现失效点，成层揭开，而结皮破坏后的表土由于无物理结皮保护，土壤侵蚀可能反而加剧。目前，有关生物结皮直接作用于土壤侵蚀的研究较多^[14-16]，而对于生物结皮覆盖及破坏状态下坡面流水动力特性研究较少。基于此，本研究采用人工模拟降雨的试验方法，探究不同类型生物结皮破坏前后坡面流水力学特征，以期为黄土高原丘陵区生物结皮坡面流侵蚀水动力研究提供理论依据。

1 材料与方法 1.1 试验方法 1.1.1 试验土槽及接种方法

试供土样采自黄土丘陵沟壑区安塞县纸坊沟流域，其土壤质地为粉质壤土，机械组成为黏粒15.74%，粉粒63.45%，砂粒20.81%(土壤颗粒组成采用马尔文MS2000激光粒度分析仪进行测定)。将试供土样过2 mm筛备用，填土前，将土槽调至水平，在土槽底部均匀装填0.1 m厚的天然沙，并在沙子表面铺上透水纱布，保证土壤能够有良好的透水性。按照容重1.2 g/cm³分4层每层10 cm厚填土，填装下层之前将表土打毛，消除土壤之间的分层现象。接种材料及方法：藻结皮种源为硅藻，接种量为50 g/m²，将藻粉溶于1 L水中并均匀喷洒于试验土槽表面；苔藓结皮种源采自纸坊沟流域铁杆蒿样地，选取具有成熟孢蒴的藓结皮，铲取藓结皮同1 cm的土层，剔出肉眼可辨的植物残渣、土块、石子并装入洁净塑料袋中，运回实验室自然风干，用粉碎机粉碎后制成种子土，接种量为1 000 g/m²，均匀洒于试验土槽的土壤表面，定期浇水保证结皮正常生长。

1.1.2 试验处理

本试验共设置5个处理(表 1)，其中每个处理2个重复，共10个试验土槽(2.0 m×0.5 m×0.5 m)。T₀为静置2个月的裸地对照处理，T₁为种植两个月的藻结皮处理(盖度为68.3%，厚度变化范围为1.18~1.87 mm)，T₂为种植两个月后铲除藻结皮处理，T₃为种植两个月的苔藓结皮处理(盖度为59.2%，厚度变化范围为4.00~4.18 mm)，T₄为种植两个月后铲除苔藓结皮处理。破坏处理均为轻铲掉生物结皮，使其与土壤分离，用小刷子扫去表面浮土。

1.1.3 试验降雨过程

2019年6月在黄土高原土壤侵蚀与旱地农业国家重点实验室人工模拟降雨大厅进行模拟降雨试验。降雨试验前将土槽缓慢喷水至饱和，减少前期含水量差异对试验的影响，试验使用QYJY-503侧喷式人工降雨设备，降雨高度18 m，可调节雨强范围为30~300 mm/h，降雨均匀度大于80%，雨滴终点速度近似天然降雨。降雨前率定降雨强度，以确保最终调试的降雨强度的误差在5%内。本次试验降雨坡度为25°，试验土槽坡长为2 m，考虑到黄土高原极端降雨情况^[17]，平均雨强设置为120 mm/h，为减小试验过程中人为干扰的影响，从坡面产流4 min后开始计时，降雨历时为60 min，每隔2 min接一次径流泥沙样称重，并使用高锰酸钾染色示踪法测定2次水流流经坡面中段0.5~1 m的时间以计算流速，同时记录水流温度以计算运动黏性系数。

1.2 水动力参数计算

为了描述坡面流水动力特征的变化，本试验选取雷诺数(R_e)、傅汝德数(F_r)、流速(V，m/s)、水流剪切力(τ，Pa)、Darcy-Weisbach阻力系数(f)、水流功率〔ω，N/(m·s)〕参数指标进行水动力特征描述，各参数的计算公式为：

$ {R_e} = \frac{{Vh}}{\gamma } $

(1)

$ {F_r} = \frac{V}{{\sqrt {gh} }} $

(2)

$ V = k{V_m} $

(3)

$ h = \frac{q}{V} $

(4)

$ \tau = \rho gRJ $

(5)

$ f = \frac{{8gRJ}}{{{V^2}}} $

(6)

$ W = \tau V $

(7)

式中：V为径流平均流速(m/s)；h为坡面水深(m)；γ为运动黏性系数(m²/s)，是水流温度的函数；g为重力加速度，取9.8 m/s²；V_m为坡面流表层流速(m/s)；k为修正系数，取0.67^[18]；q为流经土槽的单宽流量〔m³/(s·m)〕；ρ为水的容重(kg/m³)；R为水力半径，用水深h代替水力半径(m)；J为水力坡度，取坡度的正切值。

2 结果与分析 2.1 流态

本试验中各处理坡面流平均雷诺数由大到小依次为：T₄＞T₂＞T₃＞T₁＞T₀，生物结皮覆盖及破坏处理均不同程度增大了雷诺数，且生物结皮破坏处理(T₂，T₄)的雷诺数高于生物结皮处理(T₁，T₃)，其中生物结皮处理(T₁，T₃)的雷诺数分布范围为22.98~46.90，生物结皮破坏处理(T₂，T₄)雷诺数分布范围为44.57~49.67。除藻结皮处理(T₁)雷诺数变化较大外，其余各处理分布都较为集中(图 1)，且均远小于500，表明各处理(T₁，T₂，T₃，T₄)在降雨过程中黏滞力占主导地位，坡面流均为层流。傅汝德数(Fr)常用来表示水流惯性力和重力的大小，其值与1的大小用来确定水流的流型，各处理坡面流平均傅汝德数由大到小依次为：T₀＞T₄＞T₂＞T₁＞T₃＞1，各处理(T₁，T₂，T₃，T₄)坡面流均为急流。生物结皮处理(T₁，T₃)的傅汝德数变化范围为1.21~3.06显著低于裸地(p＜0.05)，其中苔藓结皮处理(T₃)傅汝德数接近于1，显著区别于藻结皮处理(T₁)(p＜0.05)。生物结皮破坏处理(T₂，T₄)的傅汝德数变化范围为5.93~7.75小于裸地。研究表明生物结皮覆盖可显著降低坡面流傅汝德数(p＜0.05)，结皮破坏后坡面水流流态与裸地无显著差别。

2.2 流速和水流剪切力

流速作为坡面流观测的关键变量之一是水动力参数计算的基础^[19]，生物结皮处理(T₁，T₃)的流速变化较为平稳，生物结皮破坏处理(T₂，T₄)的流速波动变化(图 2)。各处理坡面流平均流速由大到小依次为：T₄＞T₂＞T₀＞T₁＞T₃，藻结皮处理(T₁)及苔藓结皮处理(T₃)较裸地(T₀)流速分别减少了48.74%和62.24%，藻结皮破坏处理(T₂)和苔藓结皮破坏处理(T₄)较裸地(T₀)流速分别增加了10.84%，16.46%。研究表明结皮覆盖较裸地可以显著降低流速(p＜0.05)，结皮破坏时流速显著增加(p＜0.05)。

水流剪切力作为一种基本的水动力学参数，对了解细沟侵蚀机理和建立侵蚀模型有重要的意义。各处理平均水流剪切力由大到小依次为：T₃＞T₁＞T₂＞T₄＞T₀，水流剪切力随着降雨历时的延长均呈上升的变化趋势，其中藻结皮处理(T₁)在前25 min变化平稳，在30 min时呈快速增加的变化趋势(图 3)。藻结皮处理(T₁)及苔藓结皮处理(T₃)较裸地(T₀)水流剪切力分别增加了130.01%和308.07%，藻结皮破坏处理(T₂)和苔藓结皮破坏处理(T₄)较裸地(T₀)水流剪切力分别增加了36.83%，42.44%。研究结果表明生物结皮覆盖下坡面流水流剪切力显著增加(p＜0.05)，结皮去除后水流剪切力较结皮覆盖减小但仍显著高于裸地(p＜0.05)。

2.3 Darcy-Weisbach阻力系数

Darcy-Weisbach阻力系数(f)的大小表征下垫面对水流的水力阻力状况。本试验中，阻力系数随着降雨历时的延长呈波动增加(图 4)。各处理坡面流Darcy-Weisbach阻力系数平均值由大到小依次为：T₃＞T₁＞T₂＞T₀＞T₄，藻结皮处理(T₁)及苔藓结皮处理(T₃)较裸地(T₀)的Darcy-Weisbach阻力系数分别增大了703.35%和2 532.01%，藻结皮破坏处理(T₂)及苔藓结皮破坏处理(T₄)与生物结皮处理(T₁，T₃)相比Darcy-Weisbach阻力系数分别减小了86.52%和96.44%，且与裸地(T₀)的Darcy-Weisbach阻力系数无显著差异，其中苔藓结皮破坏处理(T₄)的Darcy-Weisbach阻力系数甚至较裸地降低了6.16%。试验表明生物结皮覆盖可显著增大Darcy-Weisbach阻力系数(p＜0.05)，影响坡面阻力状况，结皮破坏后Darcy-Weisbach阻力系数较生物结皮覆盖显著降低且与裸地无显著差别。

2.4 水流功率

水流功率表征一定高度水流在顺坡流动时所具有的势能，反映剥蚀一定量土壤所需功率。在整个降雨过程中水流功率随降雨历时呈波动增加，其中藻结皮处理(T₁)水流功率波动较大，与水流剪切力变化趋势相同，在前25 min变化平稳，在30 min时快速增加(图 5)。水流功率平均值由大到小依次为：T₄＞T₂＞T₃＞T₁＞T₀。各处理(T₁，T₂，T₃，T₄)水流功率较裸地T₀均不同程度显著增加(p＜0.05)，其中藻结皮处理(T₁)及苔藓结皮处理(T₃)较裸地(T₀)水流功率分别增大21.16%和57.43%，藻结皮破坏处理(T₂)和苔藓结皮破坏处理(T₄)较裸地(T₀)水流功率分别增大60.90%和62.14%。试验表明生物结皮覆盖下水流功率显著增加(p＜0.05)，结皮破坏后水流功率进一步增大。

3 讨论

坡面径流作为三维、非恒定、非均匀的沿程变量流，其水动力参数非定值常数^[20]，由于试验技术限制，当前研究坡面流水动力多采用借鉴明渠水流研究理论，本试验中也均借鉴明渠水流研究理论讨论坡面平均状况下的水动力特性。经典的明渠水流研究理论中雷诺数小于500时坡面流为层流^[21]，本试验中各处理平均雷诺数均远小于500坡面流均为层流，薄层径流受到降雨影响虽属层流但其流体紊动强烈，当前对坡面流流态归属尚未有定论，但均认为坡面薄层水流的层流与明渠水流层流存在区别^[22]。此外，当前研究表明生物结皮可显著改变坡面流流型，降低坡面径流能量^[5]，本研究中生物结皮覆盖下流型虽未改变，但其坡面流傅汝德数仍显著降低(p＜0.05)，其径流流型接近于较为简单的缓流，降低坡面沟蚀发生的可能，结皮破坏后坡面流傅汝德数与裸地无显著差异。

生物结皮生长使土壤表层形成了一个渗透性和抗侵蚀性不同于下层的土层^[3]，从而影响水动力过程。当前研究表明藻结皮生长发育会使表面糙度增加^[23]，而苔藓结皮作为刚性植物的一种也可通过吸水体积膨胀，高度增高造成随机糙度的增大^{[12, 24]}，生物结皮通过增大表面糙度，影响坡面阻力状况，使流体动能减小，降低坡面流流速。然而，不同于藻结皮，苔藓结皮为具有一定高度的簇状生长，为径流提供通道的同时，对坡面径流形成阻水作用，使坡面形成壅水，增加径流能量消耗，因此，苔藓结皮坡面流流速较藻结皮显著降低(p＜0.05)。此外，当前研究表明生物结皮覆盖可降低土壤饱和导水率减少入渗，增加坡面径流量^[25]，进而使坡面流势能增加，本试验中苔藓促进产流作用强于藻结皮，径流量的增加及流速的减小使坡面流水深增加，坡面流势能增大，表征侵蚀动力的水流剪切力及水流功率增加，生物结皮覆盖下径流潜在侵蚀能力提高。这与当前研究生物结皮可减少侵蚀的研究结果不符^{[14, 26]}，主要是由于生物结皮增加径流侵蚀力的同时，亦会影响土壤黏结力与团聚体含量，提高土壤抗冲性减少土壤侵蚀发生概率^[16]。

生物结皮对环境改变、人为干扰较为敏感，生物结皮破坏后土壤可蚀性增大^[15]，坡面水蚀过程取决于径流侵蚀力^[27]。本试验中结皮破坏后阻力系数与裸地无差异，流速、水流剪切力、水流功率均显著高于裸地(p＜0.05)，其主要原因是结皮破坏后表面粗糙度降低，阻力状况与裸地无差异，坡面流消耗的能量减少造成流速的增加，且试验土壤为不同于沙土的粉质壤土，生物结皮破坏后土壤对水分渗透性的抑制作用不能立刻解除^{[3, 28]}，导致坡面径流量仍多余裸地，造成了结皮破坏后流体能量的增加，水流剪切力及水流功率增高，坡面径流侵蚀力增强。

4 结论

(1) 生物结皮覆盖及破坏下坡面流均属于层流急流，生物结皮可显著降低坡面流傅汝德数(p＜0.05)，使径流运动状态接近于较简单的缓流，生物结皮破坏后坡面流流态与裸地无显著区别。

(2) 生物结皮一方面可改变地表粗糙度，增加坡面流阻力系数，增大径流消耗，降低坡面流流速。另一方面可通过增加坡面径流，使流体势能增加，增大坡面流水流剪切力及水流功率，增加径流侵蚀风险。此外，由于苔藓结皮与藻结皮结构的差异，两种类型生物结皮下坡面流水动力特征也有显著区别(p＜0.05)。

(3) 生物结皮破坏后坡面阻力状况与裸地相近，流体能量消耗减少，坡面水蚀过程取决于水流侵蚀力，坡面流流速、水流剪切力、水流功率均显著高于裸地(p＜0.05)，径流侵蚀能力增强。