第42卷第2期2022年4月
水土保持通报
B u l l e t i no f S o i l a n d W a t e r
C o n s e r v a t i o n
V o l .42,N o .2
A p
r .,2022
收稿日期:2021-08-14 修回日期:2021-10-26
资助项目：中国科学院 西部之光 西部青年学者A 类项目 黄土高原退耕草地近地表特征变化及其对土壤侵蚀过程的影响 (X A B 2019A 04)；国家自然科学基金重点项目 黄土高原刺槐人工林结构与生态功能相互作用机制及调控 (42130717
) 第一作者：张慧妮(1996 )，女(汉族)，陕西省西安市人，硕士研究生，研究方向为土壤侵蚀㊂E m a i l :
Z h a n g h u i n i @n w a f u .e d u .c n ㊂ 通讯作者：王兵(1982 )，男(汉族)，陕西省武功县人，博士，研究员，博士生导师，主要从事植被恢复与土壤侵蚀相关研究㊂E m a i l :b w a n g @m s .i s w c .a c .c n
㊂黄土丘陵区生物结皮对坡面流水动力特征的影响
张慧妮1，王兵1,
2
，张宝琦3，韩剑桥1，许欢欢1(1.
西北农林科技大学水土保持研究所黄土高原土壤侵蚀与旱地农业国家重点实验室，陕西杨凌712100;2.中国科学院大学，北京100049;3.黄河水土保持西峰治理监督局，甘肃庆阳745000)摘 要:[目的]探究生物结皮对坡面流水动力特征的影响，为黄土高原丘陵区生物结皮防蚀机理研究提供理论指导㊂[方法]设置生物结皮覆盖及破坏处理，进行人工模拟降雨试验㊂[结果]①生物结皮覆盖及破坏处理下坡面流均属于层流急流㊂②生物结皮可显著降低坡面流流速，增加水流剪切力，阻力系数和水流功率(p <0.05)㊂其中，藻结皮处理和苔藓结皮处理较裸地流速分别减少了48.74%和62.24%，水流剪切力㊁阻力系数㊁水流功率较裸地分别增加了130.01%和308.07%;703.35%和2532.01%;21.16%和
57.43%㊂③生物结皮破坏后其阻力系数与裸地无显著变化，流速㊁水流剪切力㊁水流功率均显著高于裸地(p <0.
05)㊂其中，藻结皮破坏和苔藓结皮破坏处理较裸地流速㊁水流剪切力㊁水流功率分别增加了10.84%和16.46%;36.83%和42.44%;60.90%和62.14%㊂[结论]生物结皮覆盖会增大坡面流阻力系数，降低流速㊂生物结皮破坏后坡面流阻力系数与裸地无差异，流速增加㊂此外，生物结皮覆盖及破坏处理均会增大坡面流水流剪切力及水流功率，进而增加径流侵蚀能量㊂关键词：生物结皮；藻结皮；苔藓结皮；水动力学；结皮破坏
文献标识码:A 文章编号:1000-288X (2022)02-0016-06
中图分类号:S 157.1
文献参数：张慧妮，王兵，张宝琦，等.黄土丘陵区生物结皮对坡面流水动力特征的影响[J ].
水土保持通报,2022,42(2):16-21.D O I :10.13961/j .c n k i .s t b c t b .2022.02.003;Z h a n g H u i n i ,W a n g B i n g ,Z h a n g B a o q
i ,e t a l .E f f e c t s o f b i o c r u s t s o n h y d r o d y n a m i c c h a r a c t e r i s t i c s o f s l o p e f l o w i n h i l l y r e g i o n s o f L o e s s P l a t e a u [J ].B u l l e t i no f S o i l a n d W a t e rC o n s e r v a
t i o n ,2022,42(2):16-21.
E f f e c t s o fB i o c r u s t s o nH y d r o d y
n a m i cC h a r a c t e r i s t i c s o f S l o p eF l o w i nH i l l y R e g
i o n s o fL o e s sP l a t e a u Z h a n g H u i n i 1,W a n g B i n g 1,
2,Z h a n g B a o q i 3,H a n J i a n q
i a o 1,X uH u a n h u a n 1
(1.S t a t eK e y L a b o r a t o r y o f S o i lE r o s i o na n dD r y l a n dA g
r i c u l t u r e o n t h e L o e s sP l a t e a u ,I n s t i t u t e o f S o i l a n d W a t e rC o n s e r v a t i o n ,N o r t h w e s tA&F U n i v e r s i t y ,
Y a n g l i n g ,S h a a n x i 712100,C h i n a ;2.U n i v e r s i t y o f C h i n e s eA c a d e m y o f S c i e n c e ,B e i j i n g 1
00049,C h i n a ;3.S u p e r v i s i o nB u r e a uo f S o i l a n d W a t e rC o n s e r v a t i o n ,X i f e n g
G o v e r n a n c e ,Q i n g y a n g ,
G a n s u 745000,C h i n a )A b s t r a c t :[O b j e c t i v e ]T h e e f f e c t s o f b i o l o g i c a l s o i l c r u s t s (b i o c r u s t s )o n t h e h y d r o d y
n a m i c c h a r a c t e r i s t i c s o f s l o p e f l o w w e r e i n v e s t i g a t e d i no r d e r t o p r o v i d e t h e o r e t i c a l g u i d a n c e f o r t h es t u d y o
fb i o c r u s t sa n t i -e r o s i v e m e c h a n i s mi nh i l l y r e g i o n s o f t h eL o e s s P l a t e a u .[M e t h o d s ]T h e t r e a t m e n t s o f b i o c r u s t s c o v e r a g e a n dd e s t r u c t i o n w e r e s e t ,a n dt h ea r t i f i c i a l r a i n f a l l s i m u l a t i o ne x p e r i m e n t sw e r ec o n d u c t e d i nt h i s s t u d y
.[R e s u l t s ]①T h e s l o p ef l o w o fb i o c r u s t sc o v e r a g
ea n dd e s t r u c t i o nt r e a t m e n t s w a st h el a m i n a r j e ts t r e a m.②T h eb i o c r u s t s c o v e r a g es i g n i f i c a n t l y r e d u c e ds l o p ev e l o c i t y ,w h e r e a si n c r e a s e ds h e a rs t r e s s ,r e s i s t a n c ec o e f f i c i e n t ,a n d s t r e a m p o w e r (p <0.
05).C o m p a r e dw i t h t h e b a r e l a n d ,t h e s l o p e v e l o c i t y o f a l g a l c r u s t a n dm o s s
c r u s t t r e a t m e n t d e c r e a s e db y 48.74%a n d62.24%.T h es h e a r f o r c e ,r e s i s t a n c ec o e f f i c i e n t ,a n ds t r e a m p o w e r i n c r e a s e db y
130.01%a n d 308.07%,703.35%a n d2532.01%,21.16%a n d57.43%,r e s p e c t i v e l y .③C o m p
a r e dw i t ht h e
b a r e l a n d ,t h e r e s i s t a n
c e c o e f f i c i e n tw a s n o s i g n i f i c a n t
d i f f
e r e n c e a
f t e r t h e b i o c r u s t sw e r e d e s t r o y
e d ,w h e r e a s
t h e f l o wv e l o c i t y,s h e a r s t r e s s,a n ds t r e a m p o w e rw e r es i g n i f i c a n t l y h i g h e r t h a nt h eb a r e l a n d(p<0.05). T h e s l o p e v e l o c i t y,s h e a r s t r e s s,a n d s t r e a m p o w e
r o f t h e a l g a e c r u s t a n dm o s s c r u s t t r e a t m e n t s i n c r e a s e d b y
10.84%a n d16.46%,36.8%a n d42.44%,60.8a n d62.14%,r e s p e c t i v e l y.[C o n c l u s i o n]T h eb i o c r u s t s c o v e r a g e
i n c r e a s e d s l o p e f l o wr e s i s t a n c ec o e f f i c i e n t a n dd e c r e a s e df l o wv e l o c i t y.A f t e r t h eb i o c r u s t sw e r ed e s t r o y e d, t h e s l o p e f l o wr e s i s t a n c ec o e f f i c i e n tw a sn od i f f e r e n t f r o mt h eb a r e l a n d,b u t t h ef l o w v e l o c i t y i n c r e a s e d.
F u r t h e r m o r e,t h e b i o c r u s t s c o v e r a g e a n d d e s t r u c t i o n t r e a t m e n t i n c r e a s e d s h e a r f o r c e a n d s t r e a m p o w e r o f t h e s l o p e f l o w,t h e r e b y i n c r e a s i n g t h e r u n o f f e r o s i o ne n e r g y.
K e y w o r d s:b i o l o g i c a l s o i l c r u s t s;a l g a e c r u s t;m o s s c r u s t;h y d r o d y n a m i c s;b i o c r u s t s d e s t r u c t i o n
土壤生物结皮是由生长在土壤表面及其以下的细菌㊁真菌㊁苔藓和地衣类等个体微小的生物成分与土壤相互作用而形成的复杂复合土层[1]，具有重要的生态功能㊂生物结皮一方面可改变表层土壤黏结力，促进土壤团聚体形成，提高土壤的抗侵蚀能力[2-3]，另一方面可改善坡面入渗，影响水动力特性进而对坡
面侵蚀产生影响[4-5]㊂如生物结皮可增加表层斥水性，形成不透水层，影响产流过程[6-7]，改变坡面流流态，降低坡面流流速[8]，进而减少土壤侵蚀量㊂黄土高原实施退耕还林(草)工程后，以藻和苔藓为主要组成物的生物结皮得以广泛发育，其盖度可达70%左右[9]㊂不同类型生物结皮组成和表面微形态具有差异，藻结皮黏连作用形成无机层使表面平滑，苔藓结皮植物体密集丛生[2,10]，不同类型生物结皮形成的下垫面不同，其水动力特性具有差异㊂此外，踩踏㊁冻融交替等均是自然界常见干扰，适度干扰及冻融交替可缓解区域水分状况，但也使原有结构在一定程度受到破坏，边缘翻起[11-12]㊂研究[13]表明生物结皮在11P a水流剪切力下会出现失效点，成层揭开，而结皮破坏后的表土由于无物理结皮保护，土壤侵蚀可能反而加剧㊂目前，有关生物结皮直接作用于土壤侵蚀的研究较多[14-16]，而对于生物结皮覆盖及破坏状态下坡面流水动力特性研究较少㊂基于此，本研究采用人工模拟降雨的试验方法，探究不同类型生物结皮破坏前后坡面流水力学特征，以期为黄土高原丘陵区生物结皮坡面流侵蚀水动力研究提供理论依据㊂
1材料与方法
1.1试验方法
1.1.1试验土槽及接种方法试供土样采自黄土丘陵沟壑区安塞县纸坊沟流域，其土壤质地为粉质壤土，机械组成为黏粒15.74%，粉粒63.45%，砂粒20.81%(土壤颗粒组成采用马尔文M S2000激光粒度分析仪
进行测定)㊂将试供土样过2mm筛备用，填土前，将土槽调至水平，在土槽底部均匀装填0.1m 厚的天然沙，并在沙子表面铺上透水纱布，保证土壤能够有良好的透水性㊂按照容重1.2g/c m3分4层每层10c m厚填土，填装下层之前将表土打毛，消除土壤之间的分层现象㊂接种材料及方法：藻结皮种源为硅藻，接种量为50g/m2，将藻粉溶于1L水中并均匀喷洒于试验土槽表面；苔藓结皮种源采自纸坊沟流域铁杆蒿样地，选取具有成熟孢蒴的藓结皮，铲取藓结皮同1c m的土层，剔出肉眼可辨的植物残渣㊁土块㊁石子并装入洁净塑料袋中，运回实验室自然风干，用粉碎机粉碎后制成种子土，接种量为1000g/m2，均匀洒于试验土槽的土壤表面，定期浇水保证结皮正常生长㊂
1.1.2试验处理本试验共设置5个处理(表1)，其中每个处理2个重复，共10个试验土槽(
2.0mˑ0.5mˑ0.5m)㊂T0为静置2个月的裸地对照处理, T1为种植两个月的藻结皮处理(盖度为68.3%，厚度变化范围为1.18~1.87mm),T2为种植两个月后铲除藻结皮处理,T3为种植两个月的苔藓结皮处理(盖度为59.2%，厚度变化范围为4.00~4.18mm),T4为种植两个月后铲除苔藓结皮处理㊂破坏处理均为轻铲掉生物结皮，使其与土壤分离，用小刷子扫去表面浮土㊂
1.1.3试验降雨过程2019年6月在黄土高原土壤侵蚀与旱地农业国家重点实验室人工模拟降雨大厅进行模拟降雨试验㊂降雨试验前将土槽缓慢喷水至饱和，减少前期含水量差异对试验的影响，试验使用Q Y J Y-503侧喷式人工降雨设备，降雨高度18m，可调节雨强范围为30~300mm/h，降雨均匀度大于80%，雨
滴终点速度近似天然降雨㊂降雨前率定降雨强度，以确保最终调试的降雨强度的误差在5%内㊂本次试验降雨坡度为25ʎ，试验土槽坡长为2m，考虑到黄土高原极端降雨情况[17]，平均雨强设置为120mm/h，为减小试验过程中人为干扰的影响，从坡面产流4m i n后开始计时，降雨历时为60m i n，每隔2m i n接一次径流泥沙样称重，并使用高锰酸钾染示踪法测定2次水流流经坡面中段0.5~1m的时间以计算流速，同时记录水流温度以计算运动黏性系数㊂
71
第2期张慧妮等：黄土丘陵区生物结皮对坡面流水动力特征的影响
表1生物结皮各处理试验设置及说明
处理编号处理名称试验处理处理表征
T0裸地静置2月的裸地通用对照
T1藻结皮培育2月藻结皮藻结皮及物理结皮的影响T2藻结皮破坏培育2月后铲除藻结皮藻结皮失效后的影响
T3苔藓结皮培育2月苔藓结皮苔藓结皮及物理结皮的影响T4苔藓结皮破坏培育2月后铲除苔藓结皮苔藓结皮失效后的影响
1.2水动力参数计算
为了描述坡面流水动力特征的变化，本试验选取雷诺数(R e)㊁傅汝德数(F r)㊁流速(V,m/s)㊁水流剪切力(τ,P a)㊁D a r c y-W e i s b a c h阻力系数(f)㊁水流功率 ω,N/(m㊃s) 参数指标进行水动力特征描述，各参数的计算公式为:
R e=V hγ(1)
F r=V g h(2)
V=k V m(3)
h=q V(4)
τ=ρg R J(5)
f=8g R J V2(6)
W=τV(7)式中:V为径流平均流速(m/s);h为坡面水深(m);γ为运动黏性系数(m2/s)，是水流温度的函数;g为重力加速度，取9.8m/s2;V m为坡面流表层流速(m/s);k为修正系数，取0.67[18];q为流经土槽的单宽流量
m3/(s㊃m) ;ρ为水的容重(k g/m3);R 为水力半径，用水深h代替水力半径(m);J为水力坡度，取坡度的正切值㊂
2结果与分析
2.1流态
本试验中各处理坡面流平均雷诺数由大到小依次为:T4>T2>T3>T1>T0，生物结皮覆盖及破坏处理均不同程度增大了雷诺数，且生物结皮破坏处理(T2,T4)的雷诺数高于生物结皮处理(T1,T3)，其中生物结皮处理(T1,T3)的雷诺数分布范围为22.98~ 46.90，生物结皮破坏处理(T2,T4)雷诺数分布范围为44.57~49.67㊂除藻结皮处理(T1)雷诺数变化较大外，其余各处理分布都较为集中(图1)，且均远小于500，表明各处理(T1,T2,T3,T4)在降雨过程中黏滞力占主导地位，坡面流均为层流㊂傅汝德数(F r)常用来表示水流惯性力和重力的大小，其值与1的大小用来确定水流的流型，各处理坡面流平均傅汝德数由大到小依次为:T0>T4>T2>T1>T3>1，各处理(T1, T2,T3,T4)坡面流均为急流㊂生物结皮处理(T1,T3)
的傅汝德数变化范围为1.21~3.06显著低于裸地(p<0.05)，其中苔藓结皮处理(T3)傅汝德数接近于1，显著区别于藻结皮处理(T1)(p<0.05)㊂生物结皮破坏处理(T2,T4)的傅汝德数变化范围为5.93~7.75小于裸地㊂研究表明生物结皮覆盖可显著降低坡面流傅汝德数(p<0.05)，结皮破坏后坡面水流流态与裸地无显著差别㊂
图1不同处理坡面水流流型、流态特征
2.2流速和水流剪切力
流速作为坡面流观测的关键变量之一是水动力参数计算的基础[19]，生物结皮处理(T1,T3)的流速变化较为平稳，生物结皮破坏处理(T2,T4)的流速波动变化(图2)㊂各处理坡面流平均流速由大到小依次为:T4>T2>T0>T1>T3，藻结皮处理(T1)及苔藓结皮处理(T3)较裸地(T0)流速分别减少了48.74%和62.24%，藻结皮破坏处理(T2)和苔藓结皮破坏处理(T4)较裸地(T0)流速分别增加了10.84%,16.46%㊂研究表明结皮覆盖较裸地可以显著降低流速(p<0.05)，结皮破坏时流速显著增加(p<0.05)㊂
水流剪切力作为一种基本的水动力学参数，对了解细沟侵蚀机理和建立侵蚀模型有重要的意义㊂各处理平均水流剪切力由大到小依次为:T3>T1>T2
81水土保持通报第42卷
>T4>T0，水流剪切力随着降雨历时的延长均呈上升的变化趋势，其中藻结皮处理(T1)在前25m i n变化平稳，在30m i n时呈快速增加的变化趋势(图3)㊂藻结皮处理(T1)及苔藓结皮处理(T3)较裸地(T0)水流剪切力分别增加了130.01%和308.07%，藻结皮破坏处理(T2)和苔藓结皮破坏处理(T4)较裸地(T0)水流剪切力分别增加了36.83%,42.44%㊂研究结果表明生物结皮覆盖下坡面流水流剪切力显著增加(p<0.05)，结皮去除后水流剪切力较结皮覆盖减小但仍显著高于裸地(p<0.05)㊂
图2生物结皮各处理流速过程及平均值图3生物结皮各处理水流剪切力过程及平均值
2.3D a r c y-W e i s b a c h阻力系数
D a r c y-W e i s b a c h阻力系数(f)的大小表征下垫面对水流的水力阻力状况㊂本试验中，阻力系数随着降雨历时的延长呈波动增加(图4)㊂各处理坡面流D a r c y-W e i s b a c h阻力系数平均值由大到小依次为: T3>T1>T2>T0>T4，藻结皮处理(T1)及苔藓结皮处理(T3)较裸地(T0)的D a r c y-W e i s b a c h阻力系数分别增大了703.35%和2532.01%，藻结皮破坏处理(T2)及苔藓结皮破坏处理(T4)与生物结皮处理(T1,T3)相比D a r c y-W e i s b a c h阻力系数分别减小了86.52%和96.44%，且与裸地(T0)的D a r c y-W e i s b a c h 阻力系数无显著差异，其中苔藓结皮破坏处理(T4)的D a r c y-W e i s b a c h阻力系数甚至较裸地降低了6.16%㊂试验表明生物结皮覆盖可显著增大D a r c y-W e i s b a c h 阻力系数(p<0.05)，影响坡面阻力状况，结皮破坏后D a r c y-W e i s b a c h阻力系数较生物结皮覆盖显著降低且与裸地无显著差别㊂
2.4水流功率
黄土高原的范围水流功率表征一定高度水流在顺坡流动时所具有的势能，反映剥蚀一定量土壤所需功率㊂在整个降雨过程中水流功率随降雨历时呈波动增加，其中藻结皮处理(T1)水流功率波动较大，与水流剪切力变化趋势相同，在前25m i n变化平稳，在30m i n时快速增
加(图5)㊂水流功率平均值由大到小依次为:T4>T2 >T3>T1>T0㊂各处理(T1,T2,T3,T4)水流功率较裸地T0均不同程度显著增加(p<0.05)，其中藻结皮处理(T1)及苔藓结皮处理(T3)较裸地(T0)水流功率分别增大21.16%和57.43%，藻结皮破坏处理(T2)和苔藓结皮破坏处理(T4)较裸地(T0)水流功率分别增大60.90%和62.14%㊂试验表明生物结皮覆盖下水流功率显著增加(p<0.05)，结皮破坏后水流功率进一步增大㊂
91
第2期张慧妮等：黄土丘陵区生物结皮对坡面流水动力特征的影响
图4生物结皮各处理D a r c y-W e i s b a c h 阻力系数过程及平均值图5生物结皮各处理水流功率过程及平均值
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]，导致坡面径流量仍多余裸地,
02水土保持通报第42卷
发布评论