紫色土是三峡库区农业用地的主要土壤类型,坡耕地是三峡库区水土流失的发源地
[1],土壤侵蚀不仅导致表土的养分流失和库区河沙淤积,而且会造成水质恶化并危害生态健康
[2]。三峡库区土壤侵蚀剥离富含氮、磷、有机质的表层土壤,导致坡耕地耕层厚度年均减少3~5 mm
[3]。三峡水库蓄水后,回水变动区每年淤积泥沙约1.2亿t,减少有效库容,降低防洪与发电效益
[4]。三峡库区也是长江上游面源污染的高负荷区,降雨径流监测表明,一次30 mm的暴雨可造成坡耕地氮、磷流失负荷分别达到4.5~7.8 kg/hm
2和0.6~1.2 kg/hm
2[5-6];库区19条一级支流的入江口在汛期TN、TP浓度普遍超过国家地表水Ⅲ类标准1~2倍
[7]。土壤侵蚀过程中土壤氮、磷等发生迁移与富集,成为农业面源污染问题的重要源头,对三峡库区水环境生态安全构成严重威胁
[8]。氮、磷经库湾缓流区富集后,促进蓝藻、硅藻异常增殖,形成100~300 km²的水华带,威胁重庆、宜昌等沿江城市1 000余万人的饮用水源安全。三峡库区悬沙以≤0.062 mm粉砂为主,其磷通量占82%,是颗粒磷的主要载体;非汛期细颗粒含磷更高,库区对黏粒磷滞留率>60%
[9]。更严峻的是,侵蚀泥沙携带的颗粒态氮、磷可在库区干支流中再次悬浮,释放溶解态营养盐,使水体总氮、总磷浓度在汛期升高1.5~3.0倍,诱发水华,降低溶解氧,破坏鱼类产卵场和浮游生物群落结构,形成“侵蚀-淤积-富营养化”的级联生态风险
[10]。因此,在氮、磷进入受纳水体前,通过适当的技术手段阻断氮、磷等营养物质是控制农业原料污染的主要策略之一
[11]。
植物篱是将植物沿等高线行列种植形成的植物篱笆墙或植物条带,枝叶缓冲雨滴,茎秆阻滞径流,深根固土,使表土不被剥蚀,促雨水入渗,减少泥沙与氮磷流失
[12],因此,植物篱技术是目前水土保持领域治理农业面源污染较为有效且可持续性的植物措施。植物篱可通过减少坡耕地养分物质输出量,从而达到控制非点源污染的效果。有研究
[13]发现,等高植物篱通过控制<2 µm团聚体流失,有效减少氮、磷养分迁移,因其为N、P富集主要载体,是防控面源污染的关键。生态草沟是一种沿坡面或汇水线开挖并密植耐湿草本或灌木的浅沟系统
[14]。其根系能够固持水土,截留泥沙,吸附养分,净化水质。相比于传统水泥排水沟渠,草沟沟渠减少径流速度及产流量,同步实现排洪、减蚀、控污与景观提升,是一种低维护、可持续的坡面生态措施
[15-16]。生态草沟同样对于坡面径流中的各形态N素有较好的拦截作用,有研究
[17]在绿狐尾藻草沟的试验中发现,草沟对总氮的去除率可达48.7%。农业沟渠种植水芹等水生植物,不仅能通过植物自身吸收富集土壤团粒中的氮、磷,还可有效拦截携带大量氮、磷的泥沙,从而增强养分滞留效果
[18]。
目前,针对植物篱-草沟系统在坡面径流过程中对泥沙颗粒的分选机制及氮、磷等养分拦截效应的研究仍较为薄弱,尤其缺乏对其协同作用过程的系统解析。深入探究该系统对泥沙团聚体粒径分布的调控规律及养分迁移路径,不仅有助于揭示其水土保持内在机理,也为优化生态治理模式、提升坡地资源可持续利用水平提供科学依据,具有重要的理论价值与实践意义。因此,本研究在室内构建不同茎秆直径的植物篱-草沟系统,设计5种植物篱茎秆直径与2种草沟直径的组合条件,开展室内人工模拟降雨试验,探究不同茎秆直径植物篱-草沟系统坡面产流产沙及氮磷迁移特征,研究结果为防治三峡库区坡地水土流失提供有效理论依据。
1 材料与方法
1.1 试验区概况
试验用土取于重庆市西南大学水土保持试验基地的紫色土坡耕地(29°48′42″N,106°24′20″E)。该地属于亚热带季风气候,海拔266.3 m,年平均降水量为1 173.6 mm
[19],无霜期359 d,年平均气温为18.6 ℃,年日照时间达1 006.2 h。供试土壤采自长江三峡库区典型坡耕地,其母质为中生代侏罗纪沙溪庙组紫色砂泥页岩,属典型弱酸性紫色土,是该区域分布最广、耕作最为普遍的农业土壤类型,约占库区耕地总面积的78%。样品采集于表层0~20 cm耕作层。野外取样后,土壤经自然风干处理,以避免高温烘干对土壤结构和理化性质造成破坏,随后手工拣除石块、根系等杂物,并通过10 mm孔径筛进行均质化处理,确保颗粒分布相对均匀,满足模拟降雨试验对土壤样本的标准化要求。试验土壤体积质量为(1.25±0.05)g/cm
3,TN质量分数为(0.60±0.33)g/kg,TP质量分数为(0.78±0.12)g/kg,砂粒质量分数为(50.87±1.29)%,粉粒质量分数为(38.05±1.41)%,黏粒质量分数为(11.08±0.12)%。
1.2 试验设计
按照野外实地测量不同生长阶段植物篱带茎秆直径,成比例构建植物篱-草沟系统,利用柔韧仿真软胶模拟草改变单株草径围来模拟植物篱带的5种茎秆直径,单株草茎秆直径为0(对照,CK)、0.1、0.2、0.4、0.8 cm,草沟设置有草沟(0.5 cm)和无草沟(0 cm)2种处理方式。植物篱-草沟系统茎秆直径布设见
图1。
软胶模拟草可重复使用,且比自然植被更容易控制茎秆直径,能够在外力(如水流冲击)作用下产生类似真实植被的弯曲与摆动行为,从而更真实地模拟植物在降雨-径流条件下的动态阻水固土过程。每株草单独插入土壤中,插入深度为5 cm,植物篱带布置在上方土槽距离末端0.5 m处,面积为0.25 m2(0.5 m×0.5 m)。
在西南大学水土保持试验基地内进行室内人工模拟降雨试验。采用长4.8 m、宽1 m的钢制变坡土槽模拟布设植物篱模式的坡面,采用长2 m、宽0.5 m的钢制变坡土槽模拟布设生态草沟的排水沟道(
图2)。根据三峡库区降雨资料设置雨强为90 mm/h,通过实地调查紫色土区农业坡耕地和草沟坡度,设置植物篱坡面15°,生态草沟坡面16°,按照完全组合的方法进行模拟降雨试验,各组合条件重复2次,共20场降雨。
1.3 样品采集与测定
将试验土壤按体积质量分层装填到钢制土槽内,填土深度为20 cm。植物篱坡面和生态草沟填土体积质量分别为1.25、1.35 g/cm3。装填结束后以30 mm/h小雨强湿润土壤直至水分饱和,待沉降24 h后用土壤湿度计测量土壤含水量,试验开始前保持土壤含水量为标准田间持水量。降雨试验采用侧喷式人工模拟降雨装置,该装置由中国科学院水利部水土保持研究所研发制造,专用于坡面水土流失过程的模拟研究。试验中,土槽置于降雨器中央,两侧喷嘴水平间距为1.5 m,喷嘴与土壤表面之间的垂直距离约为4.0 m,可形成8 m×3 m的有效降雨面积,确保降雨覆盖整个试验小区并减少边缘效应。该侧喷式模拟降雨器通过调节供水压力和更换不同直径的喷嘴,实现30~230 mm/h内连续可调的降雨强度,能够覆盖三峡库区典型暴雨事件的强度范围。每次试验前测定雨强为90 mm/h,降雨均匀度可达到80%以上。准备工作结束后,开始正式降雨试验,生态草沟集流槽产流后开始用秒表计时,并用20 L刻度桶收集径流,降雨时长设置为60 min,前10 min内每2 min用刻度桶采集1次径流泥沙样,后50 min每5 min采集1次径流泥沙样,单场降雨共收集15次径流泥沙样品,待每场降雨结束后,每桶径流泥沙混匀后用590 mL塑料水瓶采集1瓶用于泥沙团粒粒径分布测定,每场共收集15瓶径流泥沙样。桶内泥沙沉降完全后倒掉上清液,将每桶底部的泥沙分别采集带回实验室,放入烘箱进行低温40 ℃烘干后,研磨过1 mm筛进行泥沙N、P质量分数测定。每场降雨试验结束后更换表层10 cm深土壤,并进行之前相同操作。
模拟降雨试验中所收集到的径流泥沙样品带回实验室后用Rise系列激光粒度仪分析泥沙团粒的粒径分布,将泥沙团粒共分为黏粒<2 µm、粉粒2~50 µm、砂粒50~1 000 µm 3个粒级
[20],然后采用湿筛法配合吸管法进行不同粒级泥沙团粒的分离,首先将采集的土壤样品经自然风干后,挑除根系和碎石等杂质,随后将其轻轻研磨至通过2 mm筛,获取用于分析的细土。在测定前,取适量过筛样品置于套筛中,在水中缓慢振荡进行湿筛,以避免机械破坏团聚体结构。湿筛完成后,将各筛级中的团聚体收集并烘干称重,计算各级团聚体的含量。将收集的径流泥沙样按照3个粒级进行烘干、称重,测定产沙量、TN和TP质量分数。样品测定时,先摇匀水样,采用碱性过硫酸钾消解-紫外分光光度法测定TN质量分数;将摇匀水样静置60 min取上清液,采用碱性过硫酸钾-钼酸铵分光光度法测量TP质量分数。
1.4 数据分析
采用SPSS 26.0软件进行数据整理、计算和统计分析。采用单因素方差(ANOVA)分析不同茎秆直径变化对坡面侵蚀泥沙粒径分布、TN、TP的影响;采用Pearson相关分析揭示不同茎秆直径下粒级泥沙团粒与TN、TP的关系,使用Microsoft Excel 2021和Origin 23.0软件绘制图表。
2 结果与分析
2.1 迁移泥沙粒径平均值分布
在有、无草沟条件下,CK处理坡面迁移泥沙的砂粒质量分数均为最高,分别为24.21%和31.33%,比0.1、0.2、0.4、0.8 cm茎秆直径坡面迁移泥沙的砂粒质量分数分别显著提高162.65%、44.55%、101.89%、93.13%和188.95%、160.94%、78.08%、108.76%(
图3a)。在有、无草沟条件下,CK处理坡面迁移泥沙的粉粒质量分数均为最低,分别为66.67%和59.62%,比0.1、0.2、0.4、0.8 cm茎秆直径的粉粒质量分数分别显著降低22.81%、12.87%、18.88%、20.45%和32.80%、32.10%、23.51%、24.80%(
图3b)。在有、无草沟条件下,0.8 cm茎秆直径坡面迁移泥沙的黏粒质量分数均显著最低,分别为7.16%和7.11%,比CK、0.1、0.4 cm茎秆直径的黏粒质量分数分别显著降低27.35%、24.36%、22.23%和27.32%、40.35%、23.39%(
图3c)。
2.2 迁移泥沙中TN、TP分布
2.2.1 迁移泥沙中TN质量分数平均值分布
在有草沟条件下,0.2 cm茎秆直径下坡面迁移泥沙中砂粒的TN质量分数最高,为1.02 g/kg,比CK和0.1 cm茎秆直径下的TN质量分数分别显著提高107.73%、144.59%(
图4a);在无草沟条件下,CK、0.1、0.2、0.4、0.8 cm茎秆直径下坡面迁移泥沙中砂粒的TN质量分数无显著性差异(
图4a)。
在有草沟条件下,0.2 cm茎秆直径下坡面迁移泥沙中粉粒的TN质量分数最高,为1.26 g/kg,比CK、0.1、0.4、0.8 cm茎秆直径下的粉粒的TN质量分数本部显著提高71.18%、23.42%、12.74%、51.53%(
图4b);在无草沟条件下,CK处理坡面迁移泥沙中粉粒的TN质量分数最低,为0.58 g/kg,比0.1、0.2、0.4 cm茎秆直径的粉粒的TN质量分数分别显著降低71.69%、70.43%、43.77%(
图4b)。
在有草沟条件下,0.2 cm茎秆直径下坡面迁移泥沙中黏粒的TN质量分数最高,为1.25 g/kg,比CK、0.4、0.8 cm茎秆直径黏粒的TN质量分数分别显著提高68.03%、12.41%、35.55%(
图4c);在无草沟条件下,0.1 cm茎秆直径下坡面迁移泥沙中黏粒的TN质量分数最高,为0.68 g/kg,比CK和0.8 cm茎秆直径下的TN质量分数分别显著提高58.78%、37.24%(
图4c)。
2.2.2 迁移泥沙中TP质量分数平均值分布
在有草沟条件下,0.8 cm茎秆直径下坡面迁移泥沙中砂粒的TP质量分数最低,为0.33 g/kg,比CK、0.1、0.2、0.4 cm茎秆直径砂粒的TP质量分数分别显著降低58.82%、87.55%、55.47%、43.06%(
图5a);在无草沟条件下,0.8 cm茎秆直径下坡面迁移泥沙中砂粒的TP质量分数同样最低,为0.43 g/kg,各茎秆直径下坡面迁移泥沙中砂粒的TP质量分数无显著性差异(
图5a)。
在有草沟和无草沟条件下,0.8 cm茎秆直径下坡面迁移泥沙中粉粒的TP质量分数均为最低,其中有草沟条件下粉粒的TP质量分数为0.33 g/kg,比CK、0.1、0.2、0.4 cm茎秆直径分别显著降低50.14%、56.21%、50.43%、54.93%(
图5b);在无草沟条件下粉粒的TP质量分数为0.49 g/kg,比CK、0.1、0.2、0.4 cm茎秆直径分别显著降低31.47%、32.65%、31.50%、30.48%(
图5b)。
在有草沟和无草沟条件下,0.8 cm茎秆直径下坡面迁移泥沙中黏粒的TP质量分数均最低,为0.44 g/kg,比0.1、0.2、0.4 cm茎秆直径黏粒TP质量分数分别显著降低71.89%、63.30%、62.01%(
图5c);在无草沟条件下,0.8 cm茎秆直径下坡面迁移泥沙中黏粒TP质量分数最低,为0.54 g/kg,比CK、0.1、0.4 cm茎秆直径黏粒的TP质量分数分别显著降低22.06%、28.37%、19.20%(
图5c)。
2.3 茎秆直径与泥沙粒径、TN、TP的关系
在有、无草沟时,茎秆直径与泥沙粒径、TN、TP的关系见
图6。有草沟条件下,黏粒TN质量分数与粉粒TN质量分数呈极显著正相关(
p<0.01),砂粒TP质量分数和黏粒TP质量分数、粉粒TP质量分数呈显著正相关(
p<0.05),粉粒TP质量分数和黏粒TP质量分数呈显著正相关(
p<0.05);砂粒质量分数和粉粒质量分数呈极显著负相关(
p<0.01),茎秆直径和砂粒TP质量分数呈显著负相关(
p<0.05)。
无草沟条件下,黏粒TN质量分数和粉粒TN质量分数呈极显著正相关(p<0.01),黏粒TP质量分数和粉粒TP质量分数呈显著正相关(p<0.05),黏粒质量分数和粉粒TP质量分数、黏粒TP质量分数呈显著正相关(p<0.05);茎秆直径和黏粒质量分数、粉粒TP质量分数、砂粒TP质量分数呈显著负相关(p<0.05),粉粒质量分数和砂粒质量分数呈极显著负相关(p<0.01)。
3 讨 论
本试验雨强为90 mm/h,属于暴雨,此时雨滴的溅蚀效应能显著使较大团聚体发生崩裂,产生较多细团粒,且交替的植被格局也会引起团聚体的崩裂,这种双重作用使得沟床内不易发生细沟侵蚀。从迁移泥沙粒径平均值分布来看,植物篱-草沟系统下侵蚀泥沙主要以粒径2~50 µm的团粒为主,平均占比为75.25%。在有草沟条件下,篱带茎秆直径在增大的过程中,相较0 cm茎秆直径条件下,其他茎秆直径下的侵蚀泥沙中50~1 000 µm团粒质量分数均显著减少(44.55%~162.65%),而2~50 µm团粒质量分数为12.87%~22.81%和<2 µm团粒质量分数为2.40%~27.35%,均呈增加趋势,表明植物篱-草沟系统具有泥沙分选的作用。其原因为植物篱带通过增加地表糙度显著降低径流剪切力,使50~1 000 µm大团聚体因重力沉降优先在篱前沉积;而草沟内缓流环境配合根系-微生物,促进2~50 µm微团聚体与<2 µm黏粒通过絮凝作用形成稳定悬浮流,随水流向下游迁移,从而表现出明显的粒径分选效应
[21]。当迁移泥沙通过坡面篱带的分选后汇入生态草沟,草沟沟床内更高密度的植被格局继续对迁移泥沙进行分选
[22],且对径流动能产生二次削弱,使得沟床内很难形成细沟侵蚀,故被搬运的泥沙团粒主要以2~50 µm团粒为主。
在有草沟条件下CK处理径流坡面迁移泥沙的砂粒质量分数最高为24.21%,说明相较于无草沟条件下,植物篱带能有效地拦截更多径流中粒径较大的团粒,与肖海等
[20]的研究结果相似,其在90 mm/h雨强和15°坡面条件下发现,植物篱能显著增强坡面的泥沙分选作用,且相较于裸坡地,侵蚀泥沙中50~1 000 µm团粒质量分数减少30.75%~47.28%。除0.2 cm茎秆直径下的植物篱,CK、0.1、0.4、0.8 cm茎秆直径的侵蚀泥沙中砂粒质量分数均表现为无草沟显著高于有草沟条件(17.61%~46.69%),说明草沟对径流中砂粒的拦挡作用,减小土壤表面砂粒的迁移速度,促进颗粒物的沉积,与LEES等
[23]的结果类似。与传统水泥沟渠相比,草沟由草沟壁和沟内植物组成,能有效降低径流中泥沙等颗粒物质的沉积,并且土壤孔隙被草沟堵塞,降低土壤的入渗能力。有研究
[24-25]表明,在侵蚀运行过程中,径流能量和草沟土壤基质的抗冲刷力共同决定着土体大团粒的运动,当径流能量小于草沟土壤基质的抗冲刷力时,草沟使得土壤大团粒发生沉积,草沟的沟床也发生演变,最终使泥沙悬浮与沉积达到动态平衡;吴凤至等
[26]研究降雨动能对泥沙分选的试验结果中指出,侵蚀泥沙中黏粒含量随降雨历时的增大而减小,可能因为其试验雨强超过90 mm/h,且坡面不设置水土保持措施,降雨动能的显著增大使径流搬运的泥沙团粒以大团聚体为主,导致侵蚀泥沙中细团粒质量分数降低,并且泥沙在搬运的过程中会发生絮凝效应,即细团粒重新组合成较大团粒。彭怡等
[27]研究得出,随着降雨的进行,砂粒质量分数逐渐增加,粉粒质量分数基本保持不变,黏粒质量分数在逐渐减小,与本研究结果存在差异,可能是因为植物篱和生态草沟的泥沙分选作用改变不同粒级泥沙团粒搬运-沉积的过程,也说明土壤侵蚀中泥沙团粒的迁移过程存在复杂性
[28]。
从迁移泥沙中TN质量分数结果来看,在有草沟条件下,CK处理下坡面迁移泥沙中3个不同粒径的TN质量分数均最低,随着茎秆直径增大,粒径的TN质量分数均呈现先上升后下降趋势,且0.2 cm茎秆直径下的3个粒径TN质量分数均达到最高,表明不同茎秆直径下的植物篱-草沟系统能增加土壤TN、有机质等养分,与蒲玉琳等
[29]研究结果相似。从迁移泥沙中TP质量分数结果来看,在有、无草沟的条件下,3个粒径下的TP均呈现0.1 cm茎秆直径下的侵蚀泥沙TP质量分数最高,0.8 cm茎秆直径下的TP质量分数最低。本研究发现,随着茎秆直径增大,植物篱-草沟系统显著提高径流TP的拦截效果,其原因是篱-沟复合系统的协同效应显著提升土壤有机质积累,改善土壤团粒结构,使得>250 µm团聚体质量分数提高15.00%,增强土壤对磷的吸附固定能力,从而促进土壤的熟化进程。当生态草沟茎秆直径为0.5 cm时,虽然各粒级泥沙团粒平均质量分数相较土质沟时没有显著变化,但随降雨历时的各粒级泥沙质量分数的变化相较土质沟时更稳定,波动性显著减小,说明生态草沟在生长前后对于土壤侵蚀的控制能力存在差异,在实际农业生产过程中,可针对性地对草沟生长初期采取额外水土保持措施
[30-32]。土壤团粒可表征土壤的粒径分布,反映坡面侵蚀过程对土壤质地的影响
[33]。植物篱-草沟系统中,侵蚀泥沙随篱带茎秆直径增大呈增大趋势。茎秆直径和黏粒、粉粒TP、砂粒TP质量分数呈显著负相关,与50~1 000 µm团粒质量分数呈显著负相关,与兰龙焱等
[34]在赣南草地发现的随植被盖度增大,土壤粗颗粒减少、细颗粒增加现象一致,表明根系-茎秆对水动力的削弱作用可降低大团聚体的搬运距离。本试验基于人工土槽开展,试验设计在尽可能还原自然坡面水文与侵蚀过程的基础上,设置可控的降雨强度、坡度和植被覆盖条件,以探究紫色土坡耕地在不同情景下的水土流失机制。研究结果同时表明谢贤健
[35]在泥石流频发区的分形维数与砂粒质量分数呈负相关,而与黏粒态磷(TP)呈正相关,进一步证明植物篱-草沟系统通过“截砂-保磷”作用,促进侵蚀泥沙向细颗粒富磷组分的方向分选,具有显著的养分拦截与水土保持协同效应。然而,受限于实验室空间与土槽尺寸,本研究尺度较小,难以完全再现自然流域中复杂的地形异质性、深层土壤结构及长期水文响应过程。同时,土槽的刚性边界可能引发边缘效应,干扰水分侧向运移与产流产沙的空间分布,影响观测结果的自然性;人工模拟降雨在雨滴谱、冲击动能等方面亦与天然降雨存在差异,可能导致侵蚀动力过程的模拟偏差。因此,尽管本试验为揭示植物篱-草沟系统的泥沙与养分调控机理提供一定的数据支撑,其结论在推广至更大空间尺度或不同生态区域时仍需谨慎。未来研究应结合野外大尺度径流小区观测与长期定位监测,开展多尺度对比验证,进一步检验分形特征与颗粒组成、养分分布之间的关系,优化系统参数,为三峡库区乃至南方丘陵山区的水土保持与生态修复提供更为科学、坚实的实践依据。
4 结 论
1)在植物篱-草沟系统中,侵蚀泥沙以黏粒2~50 µm团粒为主。相较裸坡,茎秆直径0.8 cm时显著增加黏粒、粉粒质量分数,显著减小砂粒质量分数。
2)在有、无草沟条件下,CK处理坡面迁移泥沙的砂粒质量分数均显著最高(p<0.05);有草沟时,0.8 cm处理下的各粒径泥沙的TP质量分数均最低(p<0.05),对TP的拦截效果最好。
3)有草沟时,CK处理坡面迁移泥沙中各粒径的TN质量分数均最低,表明植物篱-草沟系统的设置可降低氮流失,其中0.1、0.2 cm的削减效果显著(p<0.05)。
4)有草沟时,茎秆直径和砂粒TP质量分数呈显著负相关(p<0.05)。无草沟时,茎秆直径和黏粒及粉粒、砂粒TP质量分数呈显著负相关(p<0.05)。
本试验可为三峡库区建设合理且科学的植物篱-草沟系统种植模式提供理论支撑,对于促进三峡库区水土流失和面源污染生态治理有重要现实意义。
国家重点研发计划项目(2023YFF1305204)
重庆市水利科技项目(CQSLK-2023022)
国家自然科学基金项目(42107347)
重庆市硕士研究生科研创新项目(CYS240153)