黄土丘陵区作为宁夏脆弱生态的典型区域,植被恢复、发展林草是该区生态建设工作的重要内容。植被能够通过多重生态过程对水分下渗过程产生调控作用。草本作为干旱半干旱地区生态恢复的主要植被,适应性强,生长快,耗水少,根系复杂多样,大部分根集中分布于土壤浅层(0—30 cm),这些根与土壤紧密结合并常常以复合体的形式发挥生态功能。草本根土复合体不仅能够抵抗地表水流冲刷、分散,还具有渗水特性,在降雨特别是短时强降雨过程中减少地表径流的同时影响水分下渗。土壤入渗是描述水分通过地表向下的流动过程,土壤入渗率是指单位时间、单位面积土壤表面入渗的水量
[1],土壤入渗性能对调节径流、涵养水源有重要作用。因此明确根土复合体的渗透性是了解不同类型植被水文功能的重要方面,研究不同草本根土复合体入渗性能对新时期改造及恢复植被、防治水土流失具有重要意义。
在林草生态系统中,根系与土壤相互作用,密不可分,根系影响土壤的理化性质,土壤也会通过物理、化学和生物等作用影响根系的生长、功能及分布。根土复合体作为一个“组合体”,其渗透性能就与根、土壤具有一定的联系。不同植被由于独特的生物学特性,其根系对土壤理化性质、水文生态功能的影响存在差异
[2]。关于土壤渗透性及影响因素的研究成果很多,在土壤和渗透性方面:云慧雅和李平等
[3-4]研究得出,黄土区不同林分类型土壤入渗性能与孔隙状况或非毛管孔隙度、>0.25 mm水稳性团聚体含量以及有机质含量具有极显著正相关性,与容重、土壤含水量呈负相关性;张轩铭等
[5]研究表明,三峡库区紫色砂岩不同发育程度土壤入渗速率与总孔隙度、毛管孔隙度、通气孔隙度、砂粒含量呈正相关关系,与容重、黏粒含量、有机质含量呈负相关关系。苏世平等
[6]研究得出,兰州市4种灌木类型土壤渗透性与总孔隙度、土壤有机质含量、土壤全氮、土壤无机氮、土壤速效磷、土壤速效钾分别呈极显著正相关,与土壤初始含水率、土壤非毛管孔隙度分别呈显著或极显著正相关,与土壤容重呈极显著负相关;还有研究
[7]认为,土壤的机械组成和非毛管孔隙度是影响荒漠草原土壤入渗的主要因素。丁海晶等
[8]研究渗透性的区域变化规律及影响因子认为,容重、有机质含量、非毛管孔隙度等是影响土壤渗透性的主导因子,但各因子对土壤渗透性的影响存在一个阈值(分别为30%,1.0 g/cm
3,30%),一旦超过这个值其对渗透性的影响减弱。可见针对不同研究区域,不同研究对象取得的研究结果有一定差异。在根系和渗透性方面:一些研究
[1,9]指出,根系的根长密度、根表面积密度,及根体积密度和根生物量与土壤入渗特征指标均存在显著或极显著相关关系。根系构型也会影响土壤渗透性能,王鑫皓等
[10]研究表明,在根系体积相同的情况下,各根系构型下的土壤稳渗速率为V型>M型>R型>W型>VH型>H型。根系提高了土壤的入渗性能,根长密度、根体积分数、根体积分形维数与土壤稳渗速率均呈正相关关系。有研究者
[11]提出,植物根系有利于改良土壤物理性质,但也可能会堵塞、挤压孔隙进而阻碍水的下渗。Ghestem等
[12]认为水平根系能够阻止水分的下渗,倾斜根和垂直根能够促进土壤水分流向深层土壤。而Hinz等
[13]认为,根的粗细程度以及它在孔隙中所占的比例是决定根系如何影响水分运移的关键因素。还有研究者
[14]指出,根系特征对土壤入渗能力的影响较复杂,受土壤质地以及根系直径的影响。闫东锋等
[1]研究认为根系对土壤渗透能力的增强作用主要由0.5~5 mm径级的根系表现出来,娄义宝等
[15]研究表明,细根(
d≤1 mm)是植物群落改善土壤孔隙状况、提高土壤入渗和增强土壤抗冲性的有效根系,桑凯新
[16]研究表明,直径介于0.5~2 mm的根系对土壤渗透性的改善作用最好。然而,根系、土壤往往会共同影响土体的渗透性能。一些研究
[13,17-19]就表明根系特征(RLD, RSD, RVD, RB)、团聚体稳定性、土壤孔隙度/土壤大孔隙、土壤各粒径比例、饱和持水率是决定土壤入渗能力的主要因子。我们发现,很多相关研究是单独研究根系或者土壤对渗透性的影响,将根土作为一个整体,系统研究其渗透性及影响因素的很少
[1,3-6,8-9],因此本文以宁夏南部黄土丘陵区3种典型草本的根土复合体为研究对象,探究不同根土复合体渗透性的差异、规律,弄清影响渗透性能的因素,为黄土丘陵区植被恢复建设中草种的选择配置(选择利于降雨入渗、增强土壤蓄水能力的特定植物种类)及区域水土保持功能评价提供一定参考。
1 材料与方法
1.1 试验区概况
试验区在彭阳县中庄流域,位于宁南山区彭阳县草庙乡和白阳镇境内,流域总面积128.24 km²,海拔高度在1 470~1 848 m。该区是典型的黄土丘陵区,地貌以梁状丘陵为主,部分地域梁峁交错,侵蚀沟道遍布,地形破碎,沟道间偶有小型平缓地带,土质疏松,是水土流失重点区域。气候特征为典型的温带大陆性季风气候,年均气温7.6 ℃,无霜期140~160 d。年均降水量420~500 mm,年平均潜在蒸发量1 360 mm,降雨年内和年间分布不均,主要集中在7—9月。主要土壤类型为黄绵土,土壤容重1.1~1.3 g/cm3,有机质3.0~35.2 g/kg。试验区草本植物主要有百里香(Thymus mongolicus Ronn.)、星毛委陵菜(Potentilla acaulis L.)、长芒草(Stipa bungeana Trin.)、赖草(Leymus secalinus (Georgi) Tzvel.)、二裂委陵菜(Potentilla bifurca Linn.)、猪毛蒿(Artemisia scoparia)、达乌里胡枝子(Lespedeza davurica)、糙隐子草(Cleistogenes squarrosa (Trin.) Keng)、冷蒿(Artemisia frigida Willd.)、铁杆蒿(Artemisia gmelinii)等。
1.2 试验方法
1.2.1 样品采集
于2023年6—8月,在彭阳中庄小流域选择以百里香、星毛委陵菜、长芒草等为优势种的自然坡面(
表1)进行随机取样,每种草本类型选取7个点。取样时剪去植株地上部分,露出长1 m,宽0.3 m的地面作为取样面,同时挖取30 cm深度的土壤剖面。每个剖面按照0—5 cm,5—10 cm,10—20 cm分三层分别使用100 cm
3环刀、直径 7 cm的特制取样器进行取样,每层3个重复,每种草本共计取126个样品,其中42个环刀样用于渗透性测定,另42个环刀样用于容重、孔隙度及水分测定,42个铝盒样用于团聚体测定,所有样品确保原状带回实验室统一测定。同时采集分析土样,将取样点的每层土壤进行多点取样并混合,每种草本共取散土样21份,用于土壤有机质、机械组成等理化性质的测定。
1.2.2 根土复合体渗透性及根、土特征测定方法
(1) 渗透性。采用环刀法测定。将取回的环刀样品在水槽内浸泡8 h,之后取出静置30 min,去掉上盖,在上面套一个空环刀,接口处用胶布密封,防止漏水。将双环刀放于漏斗中,漏斗下承接有烧杯,向空环刀中缓缓加水,从漏斗滴下第一滴水时开始计时,以后每间隔1 min更换漏斗下的烧杯,分别记录渗出水量Q1,Q2,Q3,…,Qn,直到渗出水量相等为止,同时记录水温。为便于比较,渗透总量统一取前55 min内的渗透量。测定结束后计算:初渗率f0(mm/min),稳渗率fc (mm/min),平均入渗速率fv (mm/min)、渗透系数(Kθ )。
(2) 土壤团聚体及理化性质。团聚体,采用Yoder湿筛法测定。有机质,采用重铬酸钾外加热法测定;土壤机械组成,使用Mastersizer 3000激光衍射粒度分析仪测定,土壤粒径按照美国制(United States Department of Agriculture, USDA)标准划分为砂砾(0.05~2.0 mm)、粉粒(0.002~0.05 mm)和黏粒(<0.002 mm)3级;容重、孔隙度、含水量,采用环刀烘干称重法测定。
(3) 根系。渗透性试验结束后,将环刀样置于网筛中用清水反复冲洗,直至洗出所有根系,洗净的根用滤纸吸干水分,展开平铺放入根盘,使用浙江万深 LA-S系列根系扫描分析系统测定根系长度、表面积、体积等指标,最后将根装入纸袋后放入烘箱,在 60 ℃下烘干至恒重,称取生物量,计算根长密度、表面积密度、体积密度、比根长。
1.2.3 土壤入渗模型
本文在借鉴前人研究成果
[3-4,20]的基础上,选取常用的4个渗透性模型,分别是Philip模型,其公式为
,通用经验模型公式为
,Kostiakov模型表达式为
, Horton模型表达式为
。式中:
S(mm/min
-0.5)为吸渗率,表征土壤吸水能力的强弱;
A为稳渗率(mm/min);
t为入渗时间(min);
a、
a1、
f0为初始入渗速率;
b、
fc 为稳定入渗速率;
b1、
k为衰减系数和衰减指数;
n为拟合参数。
1.3 数据分析
通过Excel 2016,DPS17.0软件对土壤理化性质和土壤入渗指标进行统计分析,在DPS17.0中使用单因素方差分析法和Duncan新复极差法分析检验处理间的差异及显著性,利用上述4种模型在Origin 2021中对不同草本类型根土复合体的入渗过程进行拟合,利用相关性分析研究土壤、根与复合体入渗的关系。采用Origin 2021软件进行作图。
2 结果与分析
2.1 不同草本样地土壤理化性质
土壤有机质含量、容重、孔隙度、团聚体组成等是影响土壤渗透性的主要因素
[4,7-8,14]。对不同草本复合体的土壤理化性质进行测定,结果显示(
表2),3种草本0—20 cm土层根土复合体中土壤的各理化指标存在显著差异(
p<0.05)。其中百里香的有机质、含水量最大,星毛委陵菜,长芒草本较之分别低10.2%,70.1%和8.6%,25.9%,草本间差异达到显著水平(
p<0.05)。百里香、星毛委陵菜土壤总孔隙度、非毛管孔隙度显著高于长芒草1.77%~3.44%,13.75%~20.11%。草本间土壤容重、毛管孔隙度排序与总孔隙度、非毛管孔隙度正好相反。土壤机械组成方面,长芒草黏粉粒最多(65.02%),星毛委陵菜砂粒最多(38.37%),草本间差异显著(
p<0.05)。团聚体方面,星毛委陵菜>0.25 mm水稳性团聚体含量达70.05%,百里香、长芒草较之低7.24%,38.53%,草本间差异也显著(
p<0.05)。可见,百里香相比星毛委陵菜土壤有机质含量更高、土粒更细,但星毛委陵菜土壤团聚体发育更好,长芒草土壤容重最大,有机质、孔隙度及水稳性大团聚体含量最低。
2.2 不同草本根系特征
3种草本0—20 cm土层各径级根的根长密度为0~10.79 cm/cm
3(
图1),草本间大小为百里香>星毛委陵菜>长芒草。对不同草本相同径级根的根长密度进行方差分析显示,0—5 cm土层,百里香各径级根的根长密度显著大于其他两种草本(
p<0.05),星毛委陵菜与长芒草大部分根差异不显著。5—10 cm和10—20 cm土层,除0~0.5 mm根外,百里香大部分根的根长密度与星毛委陵菜差异不显著,但显著大于长芒草。可见,百里香根数量多且分布密集,长芒草根数量最少,分布稀疏。
另外,各径级根的根长密度随根径增粗在减小,0~0.5 mm根最大,占总根长密度的78.0%~90.9%,具有显著优势。草本各径级根的根长密度垂直分布有差异,百里香各径级根的根长密度随土层深度增加逐渐降低,变幅为52.9%~92.9%。星毛委陵菜、长芒草0~0.5 mm根的根长密度随土层加深先升后降(倒“V”型),其他径级根则随土层深度增加逐渐减小,变幅为34.8%~89.6%。
3种草本0—20 cm土层各径级根的表面积密度为0~78.71 mm
2/cm
3(见
图1),草本间排序与根长密度相同。不同草本间差异显著性结果显示,0—10 cm土层表面积密度差异与根长密度相同,但到10—20 cm土层时,草本间差异均不再显著。
总体来看,百里香的根与土体接触最为密切,这从侧面也反映出百里香根系最为丰富,之后依次是星毛委陵菜、长芒草。各草本不同直径根的表面积密度随根径、土层的变化与根长密度相同,其中百里香变幅53.9%~93.7%,其他草本37.6%~90.9%。各径级根中0~0.5 mm根表面积密度占比达50.5%~76.5%,具有明显优势。
结果显示(
图1),3种草本0—20 cm土层各径级根的体积密度为0~7.57 mm
3/cm
3,百里香最大,其次是星毛委陵菜,草本间根体积密度差异同表面积一致。随土层深度增加,草本间差异在变小,尤其是百里香与星毛委陵菜,但总体仍然是百里香最大。各草本不同直径根的体积密度随根径变化规律与表面积密度等明显不同,百里香在0—5 cm土层时,1.0~2.0 mm根的体积密度最大,其次是≥2.0 mm根,直径0~0.5 mm根最小,粗根(≥2.0 mm)占比达27.6%。5—10 cm土层和10—20 cm土层时,体积密度随根径增粗逐渐减小,粗根仅占6.8%~7.4%。星毛委陵菜在0—5 cm土层和5—10 cm土层时,各径级根体积密度占比18.4%~31.3%,差异不大。10—20 cm土层时,体积密度随根径加粗逐渐减小,≥2.0 mm根只占4.2%。长芒草在0—5 cm土层时,0.5~1.0 mm根的体积密度最大,其次是0~0.5 mm根,≥2.0 mm根最小,其他土层根体积密度随根径的变化与百里香一致。另外,草本各径级根体积密度垂直分布有差异,百里香各径级根体积密度随土层深度增加逐渐减小,与根长密度等相同。星毛委陵菜0~1.0 mm根、长芒草0~0.5 mm根体积密度随土层加深先增大后减小,其他径级根随土层的变化与百里香相同。
3种草本0—20 cm土层根系比根长为19.42~84.05 m/g,大小为长芒草>百里香>星毛委陵菜,说明长芒草对生境内资源的利用效率及空间拓展的能力最强。随土层深度增加比根长逐渐增大,与大部分根的根长密度、表面积密度正好相反。
综上可知,3种草本根系均以直径0~0.5 mm根占绝对优势。百里香根细数量多且分布密集,星毛委陵菜根粗数量较少,长芒草几乎无粗根,比根长最大。草本大部分根的根长密度、表面积密度随根径增粗逐渐减小。百里香根系主要集中在0—10 cm土层,随土层深度增加根系减少,根长密度等变幅最大。星毛委陵菜、长芒草根长密度等随土层增加变幅较小。
2.3 不同草本根土复合体渗透性能
2.3.1 根土复合体入渗过程
3种草本根土复合体的入渗过程如
图2所示,不同草本根土复合体入渗速率存在一定差异但呈现的规律基本一致,前0~6 min迅速递减,6~30 min缓慢波动,30~55 min趋于稳定。百里香根土复合体入渗率在0~6 min期间由3.70 mm/min降至3.02 mm/min,降低18.42%,降幅最大;星毛委陵菜根土复合体在0~6 min期间由4.17 mm/min降至3.68 mm/min,降低11.92%;而长芒草根土复合体由2.98 mm/min下降至2.67 mm/min,变幅为10.26%。
2.3.2 根土复合体渗透性
土壤初始入渗速率(
f0)、稳定入渗速率(
fc )和平均入渗速率(
fv )是反映土壤入渗性能的重要指标
[3],比较其大小关系可以了解不同草本类型根土复合体入渗能力差异。试验结果显示(
图3):3种草本0—20 cm土层根土复合体的初渗速率、稳渗速率、平均入渗速率及总渗透量大小排序均是星毛委陵菜(3.92 mm/min,3.52 mm/min,3.53 mm/min,197.69 mm)>百里香(3.38 mm/min,2.80 mm/min,2.93 mm/min,159.11 mm)>长芒草(2.80 mm/min,2.02 mm/min,2.47 mm/min,142.50 mm)。渗透性各指标变异系数为0.17~0.29。草本根土复合体的渗透性都表现为初渗速率>平均入渗速率>稳渗速率的特征。方差分析显示,各草本间的初渗速率、平均入渗速率差异显著(
p<0.05),稳渗速率、单位面积渗透总量百里香与长芒草差异不显著。由此可见,星毛委陵菜根土复合体渗透性能最好,其次是百里香,长芒草复合体渗透性能最差。
3种草本根土复合体渗透性各指标随土层的变化有所不同。对于星毛委陵菜,渗透性各指标随土层加深呈“V”型变化,5—10 cm最小,与其他土层差异显著(p<0.05),0—5 cm与10—20 cm差异不显著。对于百里香和长芒草,渗透性各指标随土层加深逐渐降低,百里香各指标3个土层差异均显著(p<0.05),长芒草初渗速率、稳渗速率及渗透总量0—5 cm与其他土层差异显著(p<0.05),5—10 cm与10—20 cm差异不显著,平均入渗速率3个土层差异均不显著。
2.3.3 不同草本根土复合体入渗过程拟合
将不同草本根土复合体入渗速率(
y)随时间(
t)的变化规律用不同模型进行拟合,结果如
表3所示。Phllip模型中
S表征土壤吸水能力的强弱,与土壤入渗能力正相关。
S值表现为百里香>星毛委陵菜>长芒草,
A值表现为星毛委陵菜>百里香>长芒草,模型拟合的入渗过程与实测数据呈现规律有所差异。通用经验模型中,
a值表现为星毛委陵菜>百里香>长芒草,拟合结果与实测结果呈现的规律一致。
b值拟合的结果百里香最大,其次是星毛委陵菜,而实际测量中星毛委陵菜大于百里香,与实际情况有一定差异。Horton模型的拟合结果与通用经验模型相似,初始入渗速率
f0与实测结果一致,但稳定入渗速率
fc 结果与实际情况不符。Kostiakov模型拟合的
a1值为2.900 4~4.046 6,且
a1值表现为星毛委陵菜>百里香>长芒草,与实测数据呈现的规律一致。
R2 表示不同草本根土复合体的入渗过程在不同模型下的拟合效果,
R2 越大,拟合效果越好
[3]。Philip模型的拟合精度为0.886~0.943,Kostiakov模型拟合精度为0.891~0.919,Horton模型的拟合精度为0.586~0.646,通用经验模型拟合精度均在0.930以上,综合以上认为,Kostiakov模型对3种草本根土复合体入渗过程的拟合效果最好,适合描述本研究区域草本复合体的入渗特征。
2.4 根土复合体渗透性能影响因子分析
2.4.1 土壤理化性质与渗透性的关系
选取与土壤特性密切相关的10个因子与渗透指标进行相关性分析,剔除影响不显著指标后结果显示(
表4、
表5):0—5 cm土层时,百里香土壤有机质(Soil Organic Matter, SOM)含量与根土复合体渗透性各指标(
f0,
fc,
fv )呈显著(
p<0.05)或极显著(
p<0.01)正相关,非毛管孔隙度(Non-Capillary Porosity, NCP)与稳渗速率
fc 、平均入渗速率
fv 呈显著正相关。10—20 cm土层时,土壤粉粒(Silt)、砂粒(Sand)含量与渗透性各指标(
f0,
fc,
fv )分别呈显著负相关、显著正相关,说明当表层土壤有机质含量越高,非毛管孔隙度越大,百里香根土复合体渗透性越强。到土壤下层时砂粒越多,粉粒越少,根土复合体渗透性越强。星毛委陵菜0—5 cm土层SOM与根土复合体渗透性各参数显著正相关,土壤容重(Bulk Density, BD)与
f0显著负相关。5—10 cm土层毛管孔隙度(Capillary Porosity, CP)与
fc 显著正相关,土壤含水量(Soil Water Content, SWC)与
fv 显著负相关。10—20 cm土层CP, NCP分别与复合体渗透性各指标极显著正相关、极显著负相关,说明星毛委陵菜表层土壤有机质含量越高,复合体渗透性越强,容重越小,初渗速率越大。中层毛管孔隙度越大,稳渗速率越大,含水量越低,平均入渗速率越快,下层毛管孔隙度越大,非毛管孔隙度越小,复合体渗透性越强。长芒草5—10 cm土层土壤黏粒(Clay)含量与渗透性各参数显著正相关,意味着它的黏粒越多,根土复合体渗透性越强。由此可见,百里香、星毛委陵菜根土复合体渗透性均受土壤孔隙状况、有机质含量显著影响,百里香、长芒草根土复合体渗透性均受土壤质地显著影响,草本种类、土层不同,影响根土复合体渗透性的土壤因素有差异。
2.4.2 根系与渗透性的关系
选取14个根系指标与渗透性参数进行相关性分析,剔除共线性强且影响不显著指标后的结果显示:0—5 cm土层时,百里香根土复合体中0~1.0 mm根的根长密度(Root Length Density, RLD)与渗透性各参数(
f0,
fc,
fv )均呈显著负相关(
p<0.05)。10—20 cm土层,根土复合体中0~0.5 mm根的RLD与渗透性各参数也显著负相关,5—10 cm土层,根土复合体的根生物量(Root Mass, RM)与
f0显著负相关,比根长(Specific Root Length,SRL)与
fc,
fv 显著正相关(
表6)。说明百里香的细根(<2.0 mm)会阻碍水分下渗,单位体积内根系长度越大,复合体渗透性越小。中层根系生物量越大,复合体初始入渗速率越小,比根长越大,稳渗速率和平均入渗速率也越大。星毛委陵菜0—5 cm土层根土复合体中≥2.0 mm根的RLD与
f0显著负相关,5—10 cm土层根土复合体的RM与渗透性各参数显著正相关(
表7),说明星毛委陵菜根土复合体特定根及根生物量显著影响水分入渗速率,表层单位体积≥2.0 mm根的长度越长,复合体初渗速率越小。长芒草0—5 cm土层根土复合体中直径≥2.0 mm根的RLD与渗透性各参数极显著正相关(
p<0.01)(
表8),表明表层单位体积内≥2.0 mm根的长度越长,其根土复合体的渗透性越强,这与星毛委陵菜表层根土复合体中≥2.0 mm根的作用正好相反。10—20 cm土层复合体中0~0.5 mm根的RLD与
f0,
fc 显著负相关,与百里香、星毛委陵菜的细根作用相似。可见,3种草本根土复合体中的根系均会显著影响复合体的渗透性能。
综合相关性分析结果认为,3种草本根土复合体渗透性能与二元系统组成部分土壤、植被根系密切相关。对百里香根土复合体而言,土壤孔隙度、有机质含量和粉砂粒含量,根系生物量及0~2.0 mm细根显著影响其渗透性。对于星毛委陵菜,土壤孔隙度、有机质含量,根系生物量及≥2.0 mm粗根和0.5~1.0 mm细根能显著影响复合体渗透性。对于长芒草,土壤黏粒含量、粗根及0~0.5 mm细根显著影响复合体渗透性。不同草本、不同土层影响渗透性的土壤因素有差异,不同草本,不同直径根系对渗透性的作用效果(正向或负向、影响程度)及土层也有差异。
3 讨 论
3.1 3种草本根土复合体的土壤、根系特征
土壤有机质等理化性质及团聚体等土壤结构是水土保持功能评价的重要指标,对防治水土流失具有重要意义
[3]。本研究中3种典型草本根土复合体土壤理化性质差异明显,百里香相比星毛委陵菜土壤有机质含量高10.2%、土粒更细,但星毛委陵菜>0.25 mm水稳性大团聚体含量较百里香高7.24%,团聚体发育更好,长芒草土壤容重最大,有机质、孔隙度及水稳性大团聚体含量显著最低,土壤养分及结构较差,推断3种草本复合体中长芒草土壤不利于水分入渗
[3,5,9],而草本间的土壤差异可能与植被生物学特性、根系构型差异等有关
[10]。
不同根系类型的草本植物具有不同的分布特征,会导致产生不同的力学效应和水文效应
[21]。根长密度、表面积密度和体积密度等是衡量植物根系分布和形态特征的重要指标。根长密度反映根系的密集程度,根长密度高的区域根系网络更密集
[14]。根表面积密度反映根与土体接触的紧密程度
[21]。根体积密度一定程度上反映根系的粗细。比根长是细根生长和吸收能力的重要指标,反映植被对特定生境有限资源的利用策略。明确根系特征对于了解不同类型草本植物固土保水能力具有至关重要的作用。研究发现,宁南山区3种典型草本根土复合体的根系形态及分布特征有一定差异,3种草本根系均以直径0~0.5 mm根占绝对优势,根长密度、表面积密度、体积密度均为百里香>星毛委陵菜>长芒草,而比根长正好相反。各草本根系特征为,百里香根细数量多且分布密集,主要集中在0—10 cm土层,其根系接触和影响的土壤多,有助于土壤孔隙的改善
[19]。星毛委陵菜根较粗,或对大孔隙及深层的土壤孔隙状况影响较大
[19,21]。长芒草几乎无粗根,但比根长(SRL)最大,对资源包括水分的利用效率最高。
3.2 3种草本根土复合体渗透性能
3种草本根土复合体渗透性具有显著差异,星毛委陵菜根土复合体渗透性能最好,其次是百里香,长芒草最差。星毛委陵菜复合体渗透性能最好,可能是它的土壤团聚体发育最好,加之根系相对少,总孔隙丰富,有利于水分的下渗
[14]。各草本根土复合体的渗透性都呈现出初渗速率>平均入渗速率>稳渗速率的特征,与李志等、白永会等
[20,22]研究结果一致。3种草本根土复合体渗透性各指标随土层的变化有所不同,对于百里香和长芒草,各指标随土层加深逐渐降低,表层渗透性最好,与许多研究结果一致
[1,6-7,20,22-23],可能是表层土壤的有机质含量高、土壤疏松,利于水分下渗。对于星毛委陵菜,渗透性各指标随土层加深呈“V”型变化(与各级根的根长密度、表面积密度变化相反),5—10 cm土层最小,可能是影响渗透性的因素除土壤外,还与其中的根系有关
[8,14]。4种常见的入渗模型中Kostiakov模型对3种草本根土复合体入渗过程的拟合效果最好(
R2≥0.890),可用来描述本研究区域典型草本根土复合体的入渗特征,与桑凯新
[16]、张志华等
[24]研究结果一致。不同根系构型植物适用的模型也会不同,王鑫皓等
[10]研究指出Horton公式适用于H型根系,Kostiakov公式适用于R型和W型根系,M型根系采用2种方式拟合均可。百里香、星毛委陵菜和长芒草根系构型分属于“W型”和“M型”,适用的模型与王鑫皓等的研究结果基本一致。
3.3 3种草本根土复合体渗透性影响因素
许多研究证实:土壤、根系都会影响渗透性,且土壤通常为直接作用,而根系并不直接影响土壤水力性质,而是通过引起土壤性质变化而产生间接影响
[14]。
本研究在土壤与渗透性关系方面发现,百里香0—5 cm土层根土复合体渗透性各指标(
f0,
fc,
fv )与土壤有机质含量显著或极显著正相关,与很多研究
[4,6,8,20,22-23]结果一致。稳渗速率、平均入渗速率与非毛管孔隙度显著正相关,与云慧雅等
[3,8,25]研究结果一致,丁海晶等
[8]认为当非毛管孔隙度在0~30%时,随着非毛管孔隙的增加,土壤稳渗速率显著提高,但超过阈值时其影响减弱。百里香10—20 cm土层根土复合体渗透性各指标与土壤粉粒、砂粒分别呈显著负相关、显著正相关关系,与许多研究
[5,23]结果一致。表层土壤有机质含量越高,非毛管孔隙度越大,百里香复合体渗透性越强,可能是表层土壤含有较多的有机质,有利于胶结土粒形成团粒结构,增加土壤孔隙,进而增强复合体渗透性能。到土壤下层时砂粒越多,粉粒越少,复合体渗透性越强,原因是较多的砂粒和较少的粉粒可以提高土壤粒间孔隙(非毛管孔隙度),利于复合体的透水透气。星毛委陵菜0—5 cm土层根土复合体初渗速率与土壤容重显著负相关,与很多研究
[5,20,23,25]结果一致。容重反映土壤的紧实程度和孔隙状况,容重越大土壤越紧实,孔隙度越小,土壤透气透水能力弱。有研究指出,容重与土壤渗透力呈现抛物线关系,土壤容重大于1.0 g/cm
3时,土壤渗透力随着土壤容重的升高而下降。容重在0.35~1.0 g/cm
3时,土壤渗透力受容重影响不大
[8]。也有研究在比较不同土壤类型入渗性能时发现,容重越大土壤入渗性能最好
[7]。5—10 cm土层根土复合体平均入渗速率与含水量显著负相关,与云慧雅等
[3]研究结果一致,但与杨浩等
[25]研究结果相反,可能是土壤含水量较大时,大量(非毛管)孔隙被水充满形成相对封闭的环境,阻碍了水分的快速渗透。10—20 cm土层根土复合体渗透性各指标与毛管孔隙度正相关,与李平等
[4-5]研究结果一致,毛管孔隙具有毛细管作用,可以持水从而使水分长久留在土壤中,因此在对饱和土体施加水头压力时,毛管孔隙度越大,复合体渗透性越强。根土复合体渗透性各指标与非毛管孔隙度呈极显著负相关,通常情况下非毛管孔隙越多土壤的渗透性越好,如百里香表层二者的正相关关系及灰钙土由低非毛管孔隙度造成的低入渗速率
[7]。非毛管孔隙能帮助降水下渗,影响土壤透气与渗水,本研究中二者呈负相关关系,可能是非毛管孔隙中还存在根系、动物和真菌等,在土体饱和情况下它们阻碍了水分的下渗
[17]。长芒草5—10 cm土层根土复合体渗透性各参数与土壤黏粒含量显著正相关,与张轩铭等
[5,14,23]研究结果一致,一般情况下,土壤细颗粒含量越高,渗透性就越低
[22,26],本结果的出现可能是研究区域土壤本身黏粒含量少,增加黏粒有助于土壤形成团聚结构,提高土壤孔隙度,进而增大根土复合体渗透性
[3-4]。
在根系与渗透性关系方面发现,百里香表层(0—5 cm)及下层(10—20 cm)根土复合体中<2 mm根参数与渗透性各指标显著负相关,即单位体积内根的长度越长,根土复合体渗透性越小,细根会阻碍水分的下渗,与闫东锋等
[1,16,23]研究结果相反。一些研究
[13,17]指出,根长密度增加不一定增加入渗能力,甚至降低土壤入渗能力,因为根系也可能会堵塞、挤压孔隙从而阻碍水的下渗。Vergani等
[27]研究认为根系发育的增加对应着土壤入渗能力的增加,但可能会存在一个阈值,根系发育一旦超出阈值,土壤入渗能力反而会随着根系的增加而降低。本研究结果可能一方面是百里香细根较多,随根系长度增加形成致密的根系网络,根网挤压土壤造成孔隙度减小,渗透性降低
[14,28],另一方面是土体在饱和情况下细根堵塞了孔隙。此外还发现,随根系增粗,百里香根影响渗透性的能力在减弱。闫东锋等
[1]研究认为,根系对土壤渗透能力的增强作用主要由0.5~5 mm径级的根系表现出来,相对径级过大或过小的根系对土壤渗透能力的增强作用都会减弱,与本研究结果相似。百里香5—10 cm土层复合体中根生物量与初渗速率显著负相关,与许多研究
[16,19]结果相反,与根长密度相似,它丰富的细根在饱和土体入渗初期堵塞孔隙导致复合体渗透性降低。星毛委陵菜根土复合体特定根及根生物量显著影响水分入渗速率,其根系特征指标与渗透性参数显著负相关,当表层(0—5 cm)单位体积≥2.0 mm根长度越长,根土复合体初渗速率越小,与张金武等
[29]研究结果一致,可能是星毛委陵菜根系为根蘖型
[21],表层为横走根分布的主要区域,这些粗根会阻碍水分下渗
[10,12]。细根降低渗透性与百里香表层细根作用相似,可能是在饱和土体中细根堵塞了孔隙。中层根生物量与渗透性各指标呈显著正相关,即根系生物量越大,水分入渗速度越快,与蒋忙舟
[30]、祁子寒等等
[19]研究结果一致,但与百里香中层根系生物量与渗透性的关系正好相反。有研究
[31]指出,根系系统功能的发挥在不同时间、不同环境条件下会有很大的差异。长芒草0—5 cm土层根土复合体中≥2.0 mm根参数与渗透性各指标极显著正相关,这与星毛委陵菜表层中≥2.0 mm根的作用正好相反,可能是长芒草根系为根茎型,表层分布的粗根为竖直向下的侧根而非横卧根,由其形成的大孔隙利于水分入渗
[13,21,28]。
虽然本研究明确了研究区3种典型草本根土复合体的特征、入渗性能及相关影响因素,但究竟哪些因素占主导 每个因素的贡献率是多少还不清楚,下一步应对此展开分析,确定关键性指标,以期更好地为宁夏南部山区生态质量提升提供参考。另外,除了明确复合体构成元素特征外,还可以利用CT扫描技术测定其内部结构,从微观视角直观分析渗透性差异原因,因此,未来渗透性研究还需要多手段、多尺度观测。
4 结 论
(1) 不同草本根土复合体土壤理化特征存在明显差异,百里香相比星毛委陵菜土壤有机质含量更高、土粒更细,但星毛委陵菜土壤团聚体发育更好,长芒草土壤容重最大,有机质含量等最低,土壤养分及结构较差。
(2) 不同草本的根系形态及分布有差异。各草本均以直径0~0.5 mm根占绝对优势,百里香根细数量多且分布密集,与土壤接触程度最高,根系主要集中在0—10 cm土层,根长密度等参数随土层增加逐渐减少,变幅最大。星毛委陵菜根粗数量相对较少,长芒草几乎无粗根但比根长最大,对资源利用效率最高,星毛委陵菜、长芒草根长密度等参数随土层增加变幅较小。
(3) 3种草本根土复合体渗透性表现为星毛委陵菜>百里香>长芒草。随土层加深,百里香、长芒草渗透性各指标逐渐减小,星毛委陵菜呈“V”型变化。Kostiakov模型能较好地拟合3种草本根土复合体的入渗过程,拟合精度在0.890以上。
(4) 3种草本根土复合体渗透性与土壤和根系密切相关。对百里香而言,土壤孔隙度、有机质含量、粉砂粒含量、根系生物量及0~2.0 mm细根显著影响复合体渗透性。对于星毛委陵菜,土壤孔隙度、有机质含量,根系生物量、≥2.0 mm粗根和0.5~1.0 mm细根显著影响复合体渗透性。对于长芒草,土壤黏粒含量、≥2.0 mm粗根及0~0.5 mm细根显著影响复合体渗透性。不同草本、不同土层影响渗透性的土壤指标有差异,不同草本,不同直径根系对渗透性的作用效果(正向或负向、影响程度)也有差异。
宁夏回族自治区农业科技自主创新资金项目“黄土丘陵区林草植被稳定性及退化人工林改造研究与示范”(NGSB-2021-14-01)
宁夏回族自治区自然科学基金项目“基于稳定同位素技术的宁南黄土区典型人工林水分利用策略研究”(2025AAC031075)