文献参数: 王逸凡, 李龙, 弥宏卓, 等.砒砂岩区不同林分土壤水分动态变化及其对降雨的响应[J].水土保持通报,2025,45(5):187-198. Citation:Wang Yifan, Li Long, Mi Hongzhuo, et al. Soil moisture dynamics in different forest species and their responses to rainfall in Pisha sandstone area [J]. Bulletin of Soil and Water Conservation,2025,45(5):187-198.
土壤水分是连接大气、土壤、植物和水体的关键纽带,在维持生态平衡方面具有重要作用
[1]。它也作为植被水分的主要供给源,影响植被生长、发育
[2]。因此土壤水分在物质循环与能量传递中占据着重要地位
[3]。土壤水分动态变化受到植被类型、温度、湿度、地形、土壤质地等影响
[4]。目前,国内外学者
[5-6]研究大多聚集在不同坡位、不同林分密度等单一影响因素对水分动态变化的影响。
降雨是影响土壤水分变化的一个重要因子,也是土壤水分重要来源,土壤含水量与降雨量、降雨强度关系密切
[7]。在植被生长季节,降雨量会显著影响土壤水分变化以及土壤水分分布
[8]。研究
[9]表明降雨量与降雨强度是影响土壤水分的主要因素。降雨强度越大,影响土层深度越深,进而影响土壤水分的蓄积。在半干旱地区,降雨可以直接影响到5—60 cm土层土壤水分,但是是在80 cm以下的土层受到降雨补给较为有限
[10-11],石建周等
[12]发现只有降雨量≥7 mm时,才可以有效补充土壤水分。成龙等
[13]研究发现,在毛乌素沙地降雨格局是影响土壤水分的重要因素,随着降雨量的增加,土壤水分入渗量也会随之增加。砒砂岩区处于半干旱区域,该地区土壤干层随着降雨量的增加出现消长状态
[14]。陈慧等
[15]研究发现,降雨强度越大,土壤下渗越快。
砒砂岩地处中国北方半干旱区,由于其独特的成岩过程,土壤-基岩水分运移规律复杂,董晓宇
[16]研究发现覆土砒砂岩入渗能力较差。目前,砒砂岩区关于土壤水分研究相对较少,主要集中在土壤水分影响因素以及土壤水分入渗特征等方面
[17-18]。砒砂岩地表土壤贫瘠,植被覆盖率低,生态环境脆弱,降雨量多集中在6—9月,降雨是补充土壤水分的重要来源。以往针对砒砂岩土壤水分的研究多集中在单独样地的水分状态研究层面
[19],难以揭示降雨输入-水分运移的链式机制,针对土壤水分的动态变化及其对多降雨的降雨关系尚未见报道。因此,本研究对砒砂岩区的雨季降雨特征、不同林分土壤水分动态变化以及土壤水分对降雨的响应展开分析,明确不同林分降雨规模下土壤水分变化规律,以期为砒砂岩区选择造林树种提供积极借鉴思路。
1 材料与方法
1.1 研究区概况
研究区位于内蒙古自治区鄂尔多斯市准格尔旗暖水乡,地理坐标为39°46′51″N,110°36′23″E,其平均海拔1 195.57 m。该地区属于温带大陆性气候,表现为冬季漫长且气候干燥,而夏季则相对短暂并带有温暖气息。多年平均下来,气温维持在7.2 ℃左右。年降水量为400 mm,其中降雨活动主要集中于夏季时段。具体而言,夏季的平均降雨量为256.4 mm,占该区全年降雨总量的64.1%,显示出显著的季节性分布特点。土壤基岩为砒砂岩,土壤类型主要以黄绵土为主,并伴随着以栗钙土和风沙土为主的砒砂岩。研究区内植被多以人工植被为主,主要有油松(Pinus tabuliformis)、樟子松(Pinus sylvestris)、侧柏(Platycladus orientalis)等;灌木主要有柠条(Caragana korshinskii)、沙棘(Hippophae rhamnoides)等。
1.2 样地设置
于2024年7—9月在砒砂岩区选择柠条纯林、油松纯林和柠条油松混交林作为试验样地,这3种林分的植被均为2013年人工种植,土壤类型以黄绵土为主,选取坡面,海拔,坡向等特征相近的样地,以减少影响,样地基本情况详见
表1。
1.3 数据获取
1.3.1 土壤含水率的动态监测
2024年7—9月在所选的3种林分样地安装土壤水分设备,3种林分均使用HOBO数据采集器(传感器为S-SMD-M005)对土壤进行连续的土壤水分检测,土壤水分探头通过监测土壤介电常数来测量土壤体积含水量。将土壤水分探头通过挖土壤剖面埋置坑中,分为10个梯度(0—10,10—20,20—30,30—40,40—50,50—60,60—70,70—80,80—90,90—100 cm)埋于样地内。将水分探头水平插入每个梯度土壤中,并设置监测频率为30 min记录1次,获取土壤水分数据。为验证土壤水分自动监测设备数据的准确性,同时采用铝盒烘干法计算每层的土壤含水率,所得数据与水分自动监测设备比较,以降低设备误差。样地内设置翻转式雨量筒进行监测降雨情况。
(1) 降雨强度。根据中国气象局对24 h内降雨量划分为标准,按照降雨量划分为小雨(0.1~9.9 mm)、中雨(10.0~24.9 mm)、大雨(25.0~49.9 mm)、暴雨(50.0~99.9 mm)。
式中:I为降雨强度(mm/h); P为降雨量(mm); T为降雨历时(h)。
(2) 变异系数(Cv )。采用垂直变异系数来分析不同植物群落样地土壤水分的垂直变化特征。计算公式为
式中:δ为标准差; μ为均值。
(3) 土壤水分对降雨的响应。土壤水分的增加量表示降雨对土壤的补给量,补给速率可以用土壤水分增量、与土壤水分达到峰值的时间之比表示。计算公式为
式中:SWmax为降雨后的土壤水分峰值(%); SW0为雨前土壤水分; ∆SW r 为降雨后的土壤水分增量(%)。
(4) 土壤储水量。土壤储水量是指在一定深度范围内土壤体积含水量。计算公式为
式中:W为各层土壤储水量(mm); θi 为第i层土壤体积含水量(m3/m3); di 为土层厚度(mm)。
1.3.2 土壤理化性质的计算
测定土壤的容重、含水率、田间持水量、总孔隙度、有机质含量。
(1) 土壤容重采用环刀法测定。将带土环刀放入烘箱中在105 ℃下烘至恒重,待冷却后称烘干土重。计算公式为
式中:为土壤容重(g/cm³); V为环刀容积(cm3); m1为环刀质量(g); m2为环刀及烘干土总质量(g)。
(2) 土壤含水率采用烘干法测定。将带土的铝盒开盖放入烘箱中在105 ℃下烘至恒重,待冷却后称烘干土重。计算公式为
式中:为土壤含水率(%); w1 为原土重(g); w2 为烘干土重(g)。
(3) 田间持水量采用环刀法测定。将取回的土壤样品放在1个大托盘中,将清水缓缓倒入托盘中,当水面低于环刀口1~2 mm时注水停止,在水中浸泡24 h,称取饱和土质量,将饱和土样置于沙盘静置12 h使重力水完全排除,放入烘箱中在105 ℃下烘至恒重,待冷却后称取烘干土质量。计算公式为
式中:F为土壤田间持水量(%); a1 为饱和土质量(g); a2 为烘干土质量(g)。
(4) 总孔隙度根据土壤容重进行计算。计算公式为
式中:Z为土壤总孔隙度(%); 为土壤容重(g/cm³); 为土壤密度,取值2.65 g/cm³。
(5) 土壤粒径通过马尔文3000激光粒度分析仪测定粒径质量分数。
(6) 土壤有机质测定采用重铬酸钾外加热法测定。
1.4 数据分析
利用Excel 2020进行数据整理,利用Origin 2024对整理后的数据进行绘图。
2 结果与分析
2.1 观测期降雨量变化特征
由
表2可知,研究区观测期内共计降雨26场,总降雨量为254.8 mm。降雨等级为小雨共发生17场,共计降雨量为26.4 mm;中雨发生3场,总降雨量为37.8 mm;大雨发生4场,总降雨量为107.4 mm;暴雨发生1场,总降雨量为78.2 mm。其中降雨等级为小雨发生频率最高,占总降雨场次的65.38%。降雨等级为大雨的降雨量最多,占总降雨量的42.1%。研究区内7—9月最高气温为39.8 ℃,最低气温为3.86 ℃,平均温度为20.76 ℃。
2.2 不同林分土壤含水率时空变化
图1为7—9月柠条纯林、油松纯林、柠条×油松混交林3种不同林分在0—100 cm土层的水分特征描述箱线图。由
图1可知,柠条×油松混交林平均土壤含水率最高(26.15%)、油松纯林次之(24.15%)、柠纯林平均土壤含水率最低(16.43%),3种林分在0—10 cm土层处,土壤含水率均呈现:油松纯林>柠条纯林>柠条×油松混交林,表明油松土壤表层保水效果优于其他2种林分。变异系数的大小可以描述土壤水分动态变化剧烈程度,变异系数分为3个等级,弱变异(变异系数<10%)、中等变异(10%>变异系数<100%)和强变异(变异系数>100%)。3种林分均属于中等变异,柠条纯林变异系数最大(33%),油松纯林变异系数次之(32.8%),柠条×油松混交林变异系数最低(28.2%)。表明柠条×油松混交林土壤水分较其他两种林分更稳定。3种林分在不同土层深度中,土壤含水率差异较大,大致呈现出随着土层深度增加,土壤含水率表现为先增加后降低再增加的波动趋势。
如
图2所示,7月3种林分土壤含水率最低,随着降雨事件的发生,3种林分在8月土壤含水率发生大幅度提升,9月土壤含水率达到峰值,在10月开始衰退,12月底,3种林分土壤含水率下降到最低。
2.3 不同林分对降雨量的响应
由
图3可知,不同林分土壤含水率对降雨等级响应程度表现为:暴雨>大雨>中雨>小雨。其中降雨等级为小雨(2.4 mm)时,只有在0—10 cm的土层深度作出响应,土壤水分增量表现为:柠条纯林(1.4%)>油松纯林(1%)>柠条×油松混交林(0.6%)。降雨等级为中雨(13.4 mm)时,在0—60 cm的土层深度作出响应,柠条纯林0—60 cm土壤水分增量分别为4.8%,1.8%,1.8%,0.9%,0.6%,1.2%。油松纯林在0—60 cm土壤水分增量分别为2.7%,3.3%,0.6%,0.8%,0.2%,0.6%。柠条×油松混交林在0—60 cm土壤水分增量分别1.9%,0.7%,0.6%,0.3%,1.4%,0.7%。降雨等级为大雨(30.6 mm)与暴雨(78.2 mm)时,3种林分均对0—100 cm的土层深度作出响应。混交林在降雨等级为暴雨时,80—100 cm的土层响应程度较低。随着降雨量的增加,土壤水分增量越大。在降雨等级为小雨和中雨时,随着土层深度的增加,降雨对土壤水分的影响大致呈现减少的趋势。
本研究以0—20 cm土层为土壤浅层,20—100 cm岩石层为深层。由
图4可知,3种林分岩石层土壤储水量高于0—20 cm土层的储水量,不同林分土壤储水量大致表现为:柠条×油松混交林>油松纯林>柠条纯林。在降雨等级为小雨时,3种林分只有土壤浅层的储水量增加。在降雨等级为中雨时,油松样地浅层储水量增加了最多(9.3 mm),柠条样地岩石层储水量增加了最多(24.6 mm)。在降雨等级为大雨时,混交林样地平均储水量增加最多(115.8 mm)。在降雨等级为暴雨时,柠条纯林、油松纯林和柠条×油松混交林浅层储水量分别增加了32.1,45.0和28.6 mm,岩石层储水量分别增加了408.9,421.4,397.4 mm。
2.4 不同林分水分动态变化特征
如
图5所示,降雨等级为暴雨(降雨量为78.2 mm)时,3种林分在0—50 cm土层中,均在雨后立刻对降雨作出响应,土壤含水率大幅度升高,在雨后第2天呈现出平稳下降趋势。在50—70 cm土层中,柠条纯林与油松纯林均在雨后第1天出现响应,但柠条×油松混交林只有在50—60 cm的土层中,对降雨做出响应。在60—70 cm土层中,柠条×油松混交林土壤水分的增加出现滞后性,在雨后第3天才出现响应。相比其他2种林分,柠条×油松混交林在70—100 cm的土层中,土壤含水率出现极小幅度的增加。在70—100 cm的土层中,油松纯林的响应时间比其他两种林分快,其中油松纯林70—90 cm的土层中,含水率在雨后第2天作出响应,90—100 cm处土壤含水率在雨后第4天作出响应。柠条纯林70—100 cm的土壤含水率依次在雨后第3,4,5 d作出响应,出现明显的滞后效应。
如
图6所示,在8月25日降雨等级为暴雨条件下,油松对降雨的响应程度优于柠条纯林,柠条×油松混交林对降雨响应程度最差。在柠条纯林样地10—60 cm土层深度响应最强烈,70—100 cm土层深度随着时间的推移,土壤水分出现滞后性。油松纯林0—80 cm土层深度在雨后土壤含水率出现大幅度升高,80—100 cm土层深度的土壤含水率只有小幅度提升。在混交林样地,0—10 cm,70—100 cm土层深度在雨后变化不明显,只有小幅度提升,20—60 cm土层深度变化强烈。
2.5 降雨量与土壤水分增量关系
在降雨发生后,土壤水分含量会相应增加,在不同降雨量下,土壤水分增量也会不同。由
图7可知,3种样地降雨量与土壤水分增量之间呈现出高度拟合,土壤水分增量与降雨量大小呈现显著正相关关系,随着降雨量的增加,土壤水分增量也相应增加。由线性回归拟合方程可知,降雨量每增加1 mm,柠条纯林土壤水分增量增加0.873%,油松纯林土壤水分增量增加0.607%,柠条×油松混交林土壤水分增量增加0.539%。
2.6 水分运移的影响因素
由图
8—
10与表3—5可知,3种林分的土壤体积含水率均与土壤有机质呈现出显著正相关关系(
p<0.05)。柠条纯林、油松纯林、柠条×油松混交林的土壤有机质含量范围分别为5.75~8.74 g/kg,3.69~10.24 g/kg和3.62~10.52 g/kg,其中10—40 cm的土层中,有机质含量相较于其他土层偏高,导致砒砂岩区的3种林分在20—40 cm土层处,含水率较高。相比于其他2种林分,柠条纯林的土壤含水率与土壤黏粒也呈现显著正相关关系(
p<0.05)。其中柠条纯林在20—30 cm的土层中,土壤含水率随着土壤黏粒含量的增加而增加。
3 讨 论
3.1 降雨等级对水分动态变化的影响
不同的降雨量、降雨强度对土壤含水率具有显著影响。
Heisler-White J L等人
[20]对半干旱地区长时序土壤水分动态研究结果表明,该地区由于蒸发强烈,土壤渗透性差,较小强度降雨事件通常仅对表层土壤产生影响。只有降雨强度较大时,才可以对深层土壤水分进行补给。
本研究发现土壤含水率随着降雨量的增加也相应呈现出上升趋势,不同林分土壤含水率对降雨等级响应程度表现为:暴雨>大雨>中雨>小雨,然而,余榕等
[21]在石漠化地区发现不同林分土壤含水率对降雨等级响应程度为:中雨>小雨>大雨>暴雨。造成这一问题原因可能是,砒砂岩是松散岩层结构,暴雨的高强度水流快速渗入深层。当降雨等级较低时,雨水仅能湿润砒砂岩地区的表层土壤,3种林分样地的降雨等级为小雨时,只有0—10 cm的土层受到湿润,这与于晓燕等
[22]在黄土丘陵区研究结果一致。
由于砒砂岩区土壤基质不同,砒砂岩的土壤保水能力较弱
[23],雨水会很快的在太阳辐射和风力等作用下蒸发,使得土壤水分不能有效地补充。随着降雨量的增大,雨水开始下渗到较深层的土壤。在降雨等级为中雨时,3种林分在0—60 cm的土层深度受到影响。杨凯悦等
[24]发现18 mm以上降雨量能影响高寒沙地人工林0—60 cm土壤。本研究发现3种林分随着土层深度的加深,土壤含水率的变化率逐渐减少。由于砒砂岩的特殊地质条件,下渗的水分也会因为重力作用和土壤特性等因素,一定时间内继续在土壤孔隙中重新分布,深层的含水量也会逐渐升高,为植被根系吸收水分提供更好的条件。在降雨等级为大雨和暴雨时,3种林分0—100 cm土层深度均受到影响,土壤水分均受到有效补充。在毛乌素沙地,降雨量≥50 mm时,湿润深度可以达到110 cm,即较高的降水会使土壤补充更多的水分
[25]。
3.2 水分运移的影响因素
水分运移与其他土壤理化性质具有密切的联系。本研究发现,在3种林分样地中,水分运移与土壤理化性质具有相关性。3种林分的土壤含水率均与土壤有机质呈现显著性正相关关系(
p<0.05),土壤有机质的增加可以使土壤水分增加
[26]。研究区柠条×油松混交林、油松纯林和柠条纯林在土层深度为20—40 cm平均土壤有机质分别为8.98,9.43和7.99 g/kg,相比于其他土层较高,因此3种林分土壤含水率在20—40 cm处高于其他土层土壤含水率。有机质含量和黏粒还关系到土壤入渗情况,有机质含量越高,土壤入渗越好
[27]。在降雨等级为暴雨时,3种林分在0—100 cm土层中,水分运移响应情况呈现出油松纯林最好,柠条纯林次之,柠条×油松混交林稍差。尤其是在70—100 cm处,混交土壤含水率增加极少,是由于油松纯林0—100 cm土层平均土壤有机质为(7.6 g/kg)>柠条纯林平均土壤有机质(7.2 g/kg)>柠条×油松混交林平均土壤有机质(6.7 g/kg)。油松纯林在0—100 cm的土层中,平均土壤黏粒含量最高(2.66%),柠条纯林平均土壤黏粒次之(2.09%),柠条×油松混交林平均土壤黏粒含量(1.54%)最少。良好的黏粒与有机质结合形成的团聚体结构,增加了土保水能力,土壤黏粒与有机质含量越高,土壤持水能力越强,透水能力越弱,进而土壤含水率越高。
本研究结果表明,砒砂岩区雨季在0—100 cm的土层中,柠条×油松混交林(26.15%)的土壤含水率最高,油松纯林(24.15%)次之,柠条纯林(16.43%)最低,出现这一现象的原因可能是受到土层厚度的影响。在雨季,土层厚度越厚,土壤含水率呈现出越高的趋势
[28]。研究表明
[29],厚土层更有利于植被生长发育,提高了土壤保水蓄水能力。柠条×油松混交林、油松纯林和柠条纯林土层厚度分别为27 cm,23 cm,21 cm,越厚的土层,土壤水分状况越好,与冯喆等
[30]研究一致。3种林分在不同土层深度中,土壤含水率差异较大,大致呈现出随着土层深度增加,土壤含水率先增加后降低再增加的波动趋势。这与李龙等
[31]在太行低山区探究荆条土壤变化含水率时得出结论相一致。植被覆盖度也是影响土壤水分的关键因子。Guo Jiahuan等
[32]发现混交林土壤肥力和酶活性均高于纯林,从而导致混交林土壤保水性要高于纯林。由于地表凋落物以及植物根系可以影响土壤水分含量
[33-34]。研究表明
[35]混交林凋落物种类数量更加丰富,有利于水分储存和渗透,从而混交林土壤含水率要高于纯林土壤含水率。油松凋落物多于柠条凋落物,且油松根系较柠条发达,保水能力更强。因此油松纯林土壤含水率高于柠条土壤含水率。
4 结 论
(1) 不同林分对土壤含水率具有影响,本研究表砒砂岩区雨季土层深度0—100 cm的土壤含水率柠条×油松混交林最高(26.15%),油松纯林次之(24.15%)、柠条纯林最低(16.43%)。3种林分土壤含水率随着土层深度的增加呈现出先增加后减少再增加的波动趋势。
(2) 降雨规模为小雨和中雨时,土壤水分补充深度分别为10 cm,60 cm,降雨规模为大雨和暴雨时,0—100 cm土壤水分均补充。3种林分土壤含水率对降雨等级响应程度表现为:暴雨>大雨>中雨>小雨。在不同降雨规模下,柠条纯林土壤水分增量最高,油松次之,柠条×油松混交林的土壤水分增量最低。
(3) 在暴雨降雨规模下,油松纯林与柠条×油松混交林0—50 cm土壤水分当天响应,而柠条纯林0—60 cm土壤水分当天响应,在雨后第2天呈现出小幅度平稳下降趋势,柠条纯林80—100 cm土壤水分在雨后6 d依旧呈现上升趋势。
(4) 在降雨等级为暴雨时,3种林分在0—100 cm土层中,水分运移响应情况呈现出油松纯林最快,柠条纯林次之,柠条×油松混交林稍慢。这是由于相比于其他2种林分,油松纯林平均有机质含量与黏粒含量最高(7.6 g/kg,2.66%),导致土壤持水能力最强,透水能力最弱,土壤运移状况最好。
内蒙古自治区科技计划项目(内蒙古自治区风沙物理与防沙治沙工程重点实验室)“砒砂岩沟坡侵蚀过程与防治机制研究”(2025KYPT0043)
内蒙古自治区科技计划项目“砒砂岩沟坡崩塌机制与沟-坡-岸综合防治关键技术研究”(2025YFSH0001)
内蒙古黄河流域多沙粗沙区林草植被恢复技术体系研究采购项目(ZB2023021185)