土壤水力性质是决定土壤中水分保持和传导能力的关键,对于确保农业和环境的可持续发展具有至关重要的意义
[1]。一般而言,土壤饱和导水率直观体现土壤在饱和状态下传导水分的能力,其值越高说明水分的渗透速度较快,土壤透水性能较好;而土壤水分特征曲线(SWCC)则精准描绘出土壤水吸力与含水率之间的关系,曲线越平缓说明土壤持水能力越好;土壤水分常数(如田间持水量、凋萎系数、最大有效含水量)是土壤在特定条件下含水量保持相对稳定的数值,反映土壤持水能力
[2]。这些指标从不同角度揭示土壤水分的运移规律和特征差异,是土壤水力性质的核心内容。一般而言,土壤水力性质与土壤性质密切相关,随着粉粒、黏粒含量的增加,土壤饱和导水率下降、SWCC在高吸力阶段变得平缓,低吸力阶段变得陡峭,田间持水量则有所提高
[3];而土壤孔隙度特别是大孔隙度的增加,会导致土壤饱和导水率的增加
[4];土壤有机质含量的增加可提高田间持水量、有效含水量及饱和导水率
[1]。受植被恢复驱动的土壤性质变化,土壤水力学性质也存在显著差异。相关研究
[5]表明,植被林分结构越复杂,根系分布范围越广,越有利于提高土壤水分的保持和渗透。与单一树种的纯林相比,混交林的多树种组合形成浅根与深根的互补,促进根系在土壤中的层次分布,从而提高土壤孔隙度和土壤饱和导水率等
[6]。此外,固氮与非固氮树种的组合则提高土壤肥力,促进田间持水量等提升
[7],使得混交林在土壤水力性质方面呈现出更加复杂且多变的变化。然而,这些变化受到混交林的树种组合、混交模式等多重因素的影响。张甜等
[8]研究表明,朴树×榉树的土壤持水能力大于油松×朴树;张益等
[9]研究则发现,侧柏×灌木的土壤饱和导水率大于侧柏×五角枫。因此,混交林的不同树种组合对土壤水力性质的影响调节机制仍有待进一步研究。
黄土高原是我国“两屏三带”生态安全战略格局的关键区域,自1999年退耕还林(草)工程实施以来,人工林种植面积大幅增加,显著改善土壤理化性质。刺槐是黄土丘陵沟壑区分布最广、面积最大的造林树种,占该区人工林总面积的90%以上
[10]。吴起县地处黄土丘陵区的森林草原过渡带,低降雨、高蒸发导致植被生长受限,不利于刺槐的生长,而营造刺槐混交林促进植被的良好生长,表明刺槐混交林可有效提高水资源利用效率,对该地区的植被恢复产生积极的影响
[11]。然而,土壤水力学性质作为研究刺槐混交林土壤水分利用的基础,其对刺槐混交林的响应机制仍不清楚。因此,本研究以黄土高原丘陵沟壑区刺槐混交林为研究对象,探究刺槐混交林土壤水力性质差异,分析土壤性质对土壤水力学性质的影响,明确土壤水力性质对刺槐林混交的响应机制,以期为黄土丘陵沟壑区土壤水分高效利用和植被优化配置提供科学依据。
1 材料与方法
1.1 研究区概况
研究区位于陕西省延安市吴起县(36°33′33″~37°24′27″N,107°38′57″~108°32′49″E),处于水蚀风蚀交错区,总面积达到3 791.5 km2,地貌类型为黄土高原梁状丘陵沟壑区,海拔1 233~1 809 m,绝大多数土壤为黄绵土。吴起县属于温带大陆性气候,四季分明,年平均气温8.23 ℃,极端最高气温37.1 ℃,极端最低气温-25.1 ℃,年平均降水量466.36 mm,主要集中在7—9月,占全年降水量的60%以上。自1998年实施封山禁牧退耕还林(草)后,林草盖度显著增加,通过直接营造刺槐(Robinia pseudoacacia)与山杏(Prunus sibirica)混交林或在沙棘(Hippophae rhamnoides)纯林中补植刺槐、油松(Pinus tabuliformis)等方式以提高混交林的占比,以改善土壤理化性质,提高生态系统稳定性。
1.2 样点布设与样品采集
在吴起县选取2000年左右营造的刺槐混交林(刺槐×山杏、刺槐×山桃、刺槐×油松、刺槐×沙棘、刺槐×紫穗槐),以刺槐纯林为对照,每个样地设置3个乔木样方(10 m×10 m)并调查坡度、坡向、树种组成、密度、树高、胸径、林龄、林下植被等信息(
表1)。
土壤样品采集:在样方内沿对角线设置样点(土壤饱和导水率和田间持水量为5个样点,土壤质量含水量和SWCC为3个样点),将采样点表面的枯落物等清除后用内径5 cm、体积100 cm
3的环刀采样,每次采样在降雨结束3~5 d后进行。在室内进行土壤质量含水量、土壤饱和导水率、田间持水量、SWCC的测定。其中土壤含水量采用土钻法进行采样,取样深度为500 cm(0~100 cm,每隔10 cm取1个样品;100~500 cm,每隔20 cm取1个样品)(
图1)。
1.3 土壤性质分析
1.3.1 土壤水分测定
将样品中的砾石等杂质剔除掉后,用烘干法测定土壤质量含水量(W),计算公式为:
式中:W为土壤质量含水量,%;mw为烘干前的湿土重,g;ms为烘干后的干土重,g。土壤质量含水量为体积含水量与体积质量的比值。
变异系数(CV)通常用于衡量数据的离散程度,采用变异系数法来划分土壤水分变化的活跃程度层次,计算公式为:
式中:CV为变异系数;
S为标准差;
为平均值。当CV≥30%为速变层,20%≤CV<30%为活跃层,10%≤CV<20%为次活跃层,0<CV<10%为相对稳定层
[12]。
1.3.2 田间持水量(θf )测定
采用威尔科克斯法
[13]测定。土壤样品浸泡至饱和直至重量不再发生变化后将其标记为A。在相同土层采集的土壤风干后过1 mm的筛子后将其填入环刀内(体积质量1.27 g/cm
3),并标记为B。将环刀A底盖打开,连同滤纸一起放到环刀B上使其紧密接触并将结合处密封,减少水分在试验过程中蒸发和渗漏。8 h后从环刀A内取出20 g左右的原状土测量含水量,此值为田间持水量。所有样品的测定在室内温度和湿度等环境条件一致下进行。
1.3.3 土壤水分特征曲线(SWCC)测定
采用离心法测定SWCC。高速离心机(GR21G)在4 ℃的恒温下,在0、1、10、20、40、60、80、100、200、400、600、800 kPa的土壤基质势下测定土壤体积含水率。SWCC采用Van Genuchten(VG)模型进行拟合,拟合公式为:
式中:为土壤体积含水率,%;为土壤饱和体积含水率,%;为土壤残余体积含水率,%;当土壤基质势为1 500 kPa时为凋萎含水量();最大有效含水量=田间持水量-凋萎含水量;为土壤基质势,kPa;为尺度参数,值越小,土壤平均孔隙直径越大;n为土壤水分特征曲线形状参数,其值越小,曲线越平缓;。
1.3.4 饱和导水率(kt)测定
采用常水头法
[14]测定
kt。将样品浸泡至饱和后取掉上盖并放上滤纸,用另一空环刀接在其上并将连接处密封好,用马氏瓶向空环刀加水并保持5 cm常水头,从滴下第1滴水开始计时,每5 min记录读数直至达到稳定状态(最后5个读数没有显著差异),计算公式为:
式中:为饱和导水率,mm/min;为n个时段入渗的水量,mL;为土层厚度,cm;为第n时段的时长,min;h为水层厚度,cm。
1.4 数据处理
采用 Excel 2019、SPSS 25.0软件对数据进行整理和分析,用Origin 2021b软件进行曲线拟合与图表绘制。
2 结果与分析
2.1 刺槐混交林土壤含水量垂直分布特征
刺槐纯林及其混交林土壤质量含水量随土层深度逐渐降低且趋于一致,为3.59%~10.84%,0~100 cm土层土壤质量含水量较高且具有较大的变异性,主要为活跃层和次活跃层(
图2)。90~100 cm土层土壤质量含水量较表层0~10 cm降低1.98%~3.90%,而480~500 cm土层较90~100 cm土层也仅降低0.53%~2.33%。尽管如此,0~100 cm土层平均土壤质量含水量依然较100~200、200~300、300~400、400~500 cm土层分别高1.66%~3.49%、2.21%~3.98%、2.66%~4.18%、2.11%~4.20%。刺槐×油松土壤水分在0~500 cm土层整体最高,尤其是在0~100 cm,其土壤质量含水量较其他样地增加4.02%~31.53%;其次是刺槐×沙棘和刺槐×山桃,0~100 cm土层土壤质量含水量较刺槐纯林平均高出19.09%、25.56%;对于刺槐×山杏,在0~30 cm土层其土壤质量含水量高于刺槐纯林0.78%~1.38%,而后均不同程度低于刺槐纯林;刺槐×紫穗槐与刺槐×山杏则相反,在0~20 cm土层土壤质量含水量较刺槐纯林低0.14%~0.61%、而后略高于刺槐纯林。说明刺槐混交林对提高土壤水分利用效率的影响主要集中在浅层土壤,并随着土层深度的增加,其影响效果逐渐减小。
2.2 刺槐混交林土壤水分特征曲线变化特征
刺槐纯林及其混交林SWCC呈“L”形,土壤含水量随土壤基质势的增加先急剧下降而后缓慢下降(
图3)。在低吸力段(0~40 kPa),所有刺槐混交林土壤含水量较刺槐纯林高2.89%~12.09%,该阶段刺槐混交林的土壤持水能力高于刺槐纯林。在中高吸力段(>40 kPa),与刺槐纯林相比,刺槐×山桃、刺槐×油松、刺槐×紫穗槐的土壤体积含水量高3.08%~14.69%,刺槐×山杏与刺槐×沙棘低10.98%、5.27%。总体而言,刺槐×油松、刺槐×山桃和刺槐×紫穗槐土壤持水能力强于刺槐纯林和其他混交林。
VG模型拟合SWCC参数见
表2,各刺槐人工林
R2均为0.99,说明采用VG模型拟合SWCC效果较好;
θr为5.90%~7.95%,刺槐×山桃与刺槐×紫穗槐比刺槐纯林高0.42%、0.71%,其他混交林比刺槐纯林低0.46%~1.34%;
θs为45.14%~51.49%,所有刺槐混交林与刺槐纯林相比提高0.10%~6.35%;
α为32.21~71.65,其值越低说明在低吸力阶段土壤持水能力越高,刺槐×山杏和刺槐×沙棘分别比刺槐纯林低2.38、8.84;
n为1.48~1.65,其值越低,曲线越平缓,刺槐×油松、刺槐×山桃、刺槐×紫穗槐与刺槐纯林相比分别低0.15、0.05、0.06。说明刺槐混交林的土壤持水能力高于刺槐纯林。
2.3 刺槐混交林土壤水分常数及饱和导水率特征
刺槐纯林及其混交林的田间持水量、最大有效含水量、土壤饱和导水率存在差异,其中刺槐混交林多高于刺槐纯林(
图4)。与刺槐纯林相比,所有刺槐混交林田间持水量均显著提高11.97%~56.47%(
p<0.05);所有刺槐混交林(除刺槐×紫穗槐外)最大有效含水量显著提高67.46%~90.78%(
p<0.05),刺槐×紫穗槐提高4.78%但不显著(
p>0.05);刺槐×山桃和刺槐×紫穗槐的凋萎含水量显著提高19.53%、23.45%(
p<0.05),其他刺槐混交林与刺槐纯林的差异不显著(
p>0.05);刺槐×油松、刺槐×山桃土壤饱和导水率显著提高150.00%、83.33%(
p<0.05),刺槐×山杏、刺槐×沙棘、刺槐×紫穗槐提高16.67%~33.33%,但不显著(
p>0.05),其中刺槐×油松的田间持水量(22.36%)、最大有效含水量(16.77%)、土壤饱和导水率(0.15 mm/min)均最高,表明刺槐混交林对提高土壤持水性能和透水性能具有显著优势,特别是刺槐×油松混交林。
2.4 土壤水力性质的影响因素分析
将土壤理化性质(
D50、黏粒、粉粒、砂粒、有机质含量、体积质量)作为解释变量、将土壤水力性质(土壤饱和导水率、田间持水量、凋萎含水量、最大有效含水量)作为响应变量进行冗余分析(
图5)。刺槐混交林(除刺槐×紫穗槐外)与刺槐纯林存在显著差异,说明刺槐混交林对土壤水力性质有显著影响。RDA1 和RDA2分别解释了土壤水力性质指标变异的85.87%、12.55%,总解释度为98.42%。有机质含量与土壤饱和导水率呈极显著正相关(
p≤0.01),与田间持水量和最大有效含水量呈显著正相关(
p≤0.05);黏粒含量与土壤饱和导水率和最大有效含水量呈显著正相关(
p≤0.05),与凋萎含水量呈显著负相关(
p≤0.05);
D50与凋萎含水量呈极显著正相关(
p≤0.01),与土壤饱和导水率、最大有效含水量呈显著负相关(
p≤0.05)。
3 讨 论
3.1 土壤水分对刺槐混交林的响应
土壤含水量与植被类型和土壤剖面深度密切相关
[15]。刺槐纯林及其混交林的土壤含水量主要在100 cm土层呈现明显下降,而随土层深度的增加土壤含水量逐渐稳定且含量较低。主要是因为黄土高原土层深厚没有地下水补给,降雨作为土壤水分唯一的来源,其补给深度主要受植物根系的影响
[16]。本研究中,林下灌木和草本的根系多集中在浅层土壤
[17],其生长过程中穿透于土壤增加土壤孔隙度,有利于水分的入渗
[18]。也正因如此,林下植被对浅层土壤水分利用也相对较多,加上蒸散发作用,导致浅层土壤含水量虽然较高但变异较大
[19]。而深层土壤由于植物根系较少,上层土体压实导致土壤容重较大、孔隙度较小,水分难以入渗
[20],因此土壤含水量总体较低且变异小。刺槐混交林类型差异也可显著影响土壤水分含量。本研究发现,刺槐×油松的土壤含水量高于其他样地,是因为油松的耗水量小于刺槐,且根系主要集中在0~40 cm土层,而刺槐根系集中在40~80 cm土层,能够减少对同一土层中水分的竞争
[21],加上油松枝叶由于具有较高的油脂和木质素,凋落物分解速度慢,枯落物蓄积量高,林下草本少,有利于提高林地对降雨的截留,减少土壤水分蒸发和植物蒸腾
[22-23]。本研究还发现,刺槐×沙棘和刺槐×紫穗槐土壤含水量呈先降低后增加再降低的变化趋势,有可能是因为沙棘和紫穗槐均为灌木,根系分布较浅 对浅层土壤水分的消耗较高
[24],同时刺槐、紫穗槐、沙棘均有较强的固氮能力,能够有效地提高土壤肥力,有利于林下植物的生长
[25-26],而草本根系多集中在0~20 cm土层,进而增加对表层土壤水分的消耗
[27],因而表层土壤水分较低且一直降低。而刺槐×山杏的土壤含水量低于刺槐纯林,可能是因为刺槐×山杏的枯落物单位厚度蓄积量低于刺槐纯林和其他刺槐混交林
[23],减少水分截留的同时增加土壤水分蒸发量,因而土壤含水量低。
3.2 土壤水力性质对刺槐混交林的响应
土壤的水力性质受到土壤理化性质的显著影响,如土壤饱和导水率、田间持水量与土壤有机质含量和孔隙度呈正相关、与土壤体积质量呈负相关,而地上凋落物、地下根系则显著影响土壤理化性质
[28]。本研究结果表明,刺槐混交林的土壤持水性能(田间持水量、最大有效含水量、SWCC)和透水性能(土壤饱和导水率)较高,与ZHANG等
[29]研究结果一致。首先,是因为混交林不同树种凋落物的化学组成不同,如山杏、沙棘凋落物的碳氮比、木质素/氮低更易分解,能够快速释放养分,进而提高土壤养分含量,同时混交林中不同的凋落物分解速度有所差异,后能够为微生物提供长期的营养物质,更有利于微生物对凋落物的持续分解,促进土壤有机质的积累、团粒结构的形成和容重的降低等
[30]。其次,混交林根系分布互补性及其分泌物对改善土壤理化性起到关键作用,例如,紫穗槐为灌木,根系分布较刺槐乔木浅,二者混交后根系有利于形成复杂且互补的网络分布格局,提高土壤孔隙度和连通性,而不同树种的根系分泌物促进土壤颗粒胶结以形成团粒结构,提升土壤有机质含量,降低土壤体积质量等
[31],因此刺槐混交林具有较高的土壤饱和导水率、田间持水量、最大有效含水量,SWCC的斜率更低。
本研究还表明,刺槐×油松的土壤饱和导水率、田间持水量、最大有效含水量均高于刺槐纯林和其他刺槐混交林,可能是因为油松较刺槐、沙棘的生长速度慢,能够减少对土壤养分的消耗
[32],加上油松叶中较高的木质素和单宁等物质,导致分解速度慢,凋落物的蓄积量较高,能够长期为微生物提供营养物质,有利于养分的长期释放
[33],进而提高土壤有机质含量等,土壤持水能力和透水能力强于其他刺槐混交林。然而,深层土壤的透水性能和持水性能对土壤剖面水分平衡有重要影响,能够影响植物的生长发育
[34],而本研究仅进行浅层土壤水分和表层土壤水力性质的研究,无法全面了解土壤水分在剖面中的分布规律及影响因素,因此,在未来的研究中可进一步探讨土层深度、季节变化等对土壤水力性质的影响,从而为提高土壤水资源利用效率提供更为全面的科学依据。
4 结 论
对于5种刺槐混交林及其纯林对照,土壤含水量在0~500 cm土层范围内均表现为先降低后趋于平稳的变化趋势,水分补给深度主要集中在0~100 cm土层,刺槐×油松和刺槐×山桃的土壤含水量较刺槐纯林提升明显,具有较高的水分利用效率。同时,刺槐×油松和刺槐×山桃的持水性能和透水性能均优于刺槐纯林。土壤持水性能和透水性能变化的驱动因素主要由土壤质地、有机质含量和体积质量解释。总体而言,刺槐×油松和刺槐×山桃对土壤水力性质的改善效果较为明显,二者可作为当地混交林建设的推荐模式。