土壤水分作为土壤的重要组成部分,是矿质养分迁移、化学物质溶解和热量传递的核心载体
[1-3]。土壤含水量是反映土壤水分动态变化关键指标,主要受降水量、降水强度、降水频率等指标的影响,同时也受土壤质地、土层深度的影响
[4],一定量、一定强度的降水发生后,土壤含水量会呈快速增大趋势,之后由于植物蒸腾、土壤蒸发等因素,土壤含水量呈缓慢降低趋势,这一过程无其他降水干扰情况称为土壤水分对降水的脉动响应
[5-6]。
目前,针对土壤水分与降水格局开展了大量相关研究,包括土壤初始状况对水平入渗的影响
[7-8],降水入渗和蒸发特征
[9-10],植物群落土壤水分动态特征
[11-14],降水入渗与再分配规律特征等
[15-16],上述研究主要围绕不同地区植被类型、土地利用类型及林分类型等展开了广泛的讨论,但已有研究中对不同沙化土地类型对降水的响应特征研究相对较少,在沙漠、沙地区域,土壤水分的主要补给方式为降水,土壤水分的时空变异性是植被生长发育的制约性因子
[17-18],因此研究沙地土壤水分对降水的响应对揭示区域水分动态变化和植被恢复均具有重要意义。
基于此,本文选取浑善达克沙地的流动沙地、半固定沙地、固定沙地为研究对象,基于降水和土壤含水量的连续动态监测,分析土壤含水量时间动态、空间分布以及降水入渗等特征,探讨不同沙地类型土壤含水量对降水的响应特征,以期为区域生态系统建设提供理论参考。
1 材料与方法
1.1 研究区概况
本研究在内蒙古多伦浑善达克沙地生态系统国家定位观测研究站内选取3种典型沙地类型作为研究对象(
图1),3种典型沙地均位于浑善达克沙地南缘,其中流动沙地(116°29′E,42°12′N)和固定沙地(116°29′E,42°8′N)行政区划隶属多伦县,固定沙地(115°9′E,42°44′N)行政区划隶属正镶白旗。研究区属中温带半干旱大陆性季风气候;降水量少且分布不均,正镶白旗年均降水量326.00 mm,多伦县年均降水量385.00 mm,7—9月降水量占全年降水量的65%以上。主要地带性土壤以栗钙土为主,非地带性土壤主要为风沙土,东部为草甸栗钙土或暗栗钙土,向西逐渐演变为淡栗钙土,到西北部则过渡为棕钙土,沙地内广泛分布的非地带性土壤风沙土大多呈坨(沙丘)、甸(丘间低地)相间分布或沙丘链与甸子地交错排列。植被类型为典型草原、草甸草原和沙地植被
[19-21],乔木植物种为白榆(
Ulmus pumila)、樟子松(
Pinus sylvestris),灌木植物种小叶锦鸡儿(
Caragana microphylla)、杨柴(
Hedysarum leave)、黄柳(
Salix gordejevii),草本植物种为羊草(
Leymus chinensis)、大针茅(
Stipa grandis)、冰草(
Leynus secalinus)、冷蒿(
Artemisia frigida)、菊叶委陵菜(
Potentilla tanace-
tifolia)、糙隐子草(
Cleistogenes squarrosa)。
1.2 试验设计
于研究区内选择流动沙地、半固定沙地和固定沙地设为试验观测场,观测时间2016—2021年,试验观测场分别设小型自动气象站和土壤含水量监测系统,其中流动沙地和固定沙地直线距离为7.4 km,两试验观测场为同一降水区域,共用一套小型自动气象站。小型自动气象站观测仪器为TYD-ZS2型自计雨量筒,雨量传感器为TR525M型,测量精度0.1 mm,记录仪时间30 min/次。土壤含水量监测系统观测仪器为HOBOH21型土壤水分自动记录仪,水分传感器型号为EC-5土壤水分传感器,土壤水分分辨率0.1%,传感器安置深度分别为20,40,60 cm土层的水平位置,记录仪时间30 min/次。
1.3 研究方法
(1) 降水事件选取:降水间隔时间大于24 h,作为两次独立的降水事件,降水历时为降水事件开始至结束的时间。并将降水事件按降水量分为Ⅰ(0.1~2 mm)、Ⅱ(2.1~5 mm)、Ⅲ(5.1~10 mm)、Ⅳ(10.1~18 mm)、Ⅴ(>18 mm)5个等级
[22-23]。
(2) 统计与分析。
式中:W为土壤储水量(cm);VWC为土壤体积含水量(%);h为土层厚度(cm)。ΔW为每次降水事件后土壤储水量增量(cm);W2 为降水事件发生后的最大储水量(cm);W1为降水事件发生前的初始储水量(cm)。
式中:V为土壤入渗速率(cm/h);ΔW为每次降水事件后土壤储水量增量(cm);T为降水发生的时间与土壤含水量峰值出现时间的差值(h)。
2 结果与分析
2.1 降水特征
分析试验观测场2016—2021年降水资料,并参考石明明
[22]、常昌明
[23]等将降水事件按降水量分为5个等级,详情见
表1。观测期内半固定沙地累计降水量1 049.40 mm,降水事件230次;Ⅰ,Ⅱ级降水事件占总降水事件67.83%,降水量占总降水量的23.84%。流动沙地、固定沙地累计降水量2 135.13 mm,降水事件418次;Ⅰ,Ⅱ级降水事件占总降水事件63.40%,降水量占总降水量的11.37%。半固定沙地、流动沙地和固定沙地以无降水天气为主,无降水日数分别占总日数的86.36%,80.74%,均呈现出为Ⅰ,Ⅱ级降水为主,Ⅲ,Ⅳ,Ⅴ级降水的频率较低,但Ⅲ,Ⅳ,Ⅴ级降水事件对总降水量的贡献较大。
2.2 土壤水分动态变化
浑善达克沙地年际降水量与土壤含水量动态变化见
图2,流动沙地和固定沙地属于同一降水区域,枯水年2020年降水量为276.52 mm;半固定沙地枯水年2020年降水量为130.60 mm。
流动沙地(
图2A)、半固定沙地(
图2B)、固定沙地(
图2C)土壤含水量年际变化均表现为先降低后增大的趋势,含水量分别为1.76%~3.12%,3.57%~6.37%,1.92%~8.04%,最低值均出现在2019年;年际间土壤平均含水量均表现为流动沙地<半固定沙地<固定沙地,不同土层土壤含水量差异表现为流动沙地土壤含水量受年内降水量和蒸发量的影响,0—60 cm土层深度土壤含水量高于0—20 cm,0—40 cm;半固定沙地、固定沙地土壤含水量受年内降水量、蒸发量、植物根系生长等共同作用,表现为0—20 cm深度土壤含水量高于0—40 cm,0—60 cm。
浑善达克沙地年内降水量与土壤含水量动态变化见
图3,流动沙地(
图3A)、半固定沙地(
图3B)、固定沙地(
图3C)观测期间降水量均表现为单峰形态,流动沙地、固定沙地降水量峰值时间为6月(180.61 mm),半固定沙地降水量峰值时间为7月(54.00 mm)。流动沙地土壤含水量年内变化表现为先增大后降低趋势,土壤含水量峰值时间均为7月;半固定沙地5—9月由于植被耗水作用土壤含水量呈降低趋势;固定沙地5—9月生长季内降水量增加及植被生长耗水共同作用,土壤含水量无明显规律性变化。
2.3 降水下渗特征
下渗速率表征土壤含水量对降水的响应时滞,沙地各降水强度下渗速率见
表2,Ⅰ,Ⅱ级降水事件因降水量、降水强度较低,降水对土壤补给作用微弱,未引起土壤下渗速率变化。Ⅲ,Ⅳ,Ⅴ级降水事件土壤下渗速率均随着降水级别的增加呈增加趋势,随着土层深度的增加呈增加趋势;流动沙地土壤下渗速率为0.4~7.8 cm/h,半固定沙地土壤下渗速率为0.30~8.35 cm/h,固定沙地土壤下渗速率为0.26~0.64 cm/h;20 cm,60 cm土层相同降水级别的土壤下渗速率表现为半固定沙地>固定沙地>流动沙地,40 cm土层表现为固定沙地>流动定沙地>半固定沙地。
2.4 土壤水分对降水的响应过程
选取观测期内典型独立降水事件分析土壤含水量对降水的响应过程,该降水事件前97~142 h无降水发生,降水事件后259~310.5 h无降水发生,因此忽略其他降水事件的影响。
图4为沙地土壤含水量对降水的响应过程,不同沙地类型土壤含水量对降水的响应均表现为先急剧增加后缓慢降低趋势,差异表现为不同沙地类型土壤含水量对降水的响应时间不同。
流动沙地(
图4A)20,40,60 cm土层土壤初始含水量分别为2.80%,2.90%,2.90%,降水事件开始14,42,86 h后土壤含水量增加至峰值,较土壤初始含水量增加96.43%,48.00%,10.34%。半固定沙地(
图4B)20,40,60 cm土层壤初始含水量分别为12.81%,10.27%,10.02%,降水事件开始8,12,120 h后土壤含水量增加至峰值,较土壤初始含水量增加57.30%,33.78%,3.19%。固定沙地(
图4C)20,40,60 cm土层土壤初始含水量分别为9.04%,9.25%,6.78%,降水事件开始22,28,88 h后土壤含水量增加至峰值,较土壤初始含水量增加65.04%,33.84%,6.49%。3种沙地土壤水分对降水事件的响应均呈现为迅速增加至峰值后缓慢降低趋势,由于流动沙地土壤初始含水量较低、固定沙地地表枯落物等原因,其差异表现为半固定沙地20,40 cm土层土壤峰值含水量的响应时间均低于流动沙地和固定沙地,流动沙地土壤峰值含水量增量均高于半固定和固定沙地。
3 讨 论
3.1 降水特征
沙地土壤水分的主要来源是大气降水,降水格局、植被特征、土壤类型等因素直接影响土壤水分状况
[24]。通过年降水资料分析表明,流动沙地、固定沙地的试验观测场2016—2021年总降水量为2 135.13 mm,降水事件418,无降水日数占总日数的80.74%;半固定沙地试验观测场总降水量为1 049.40 mm,降水事件230次,无降水日数占总日数的80.74%,流动沙地、固定沙地试验观测场降水特征优于半固定试验观测场,但均呈现出以小降水为主,发生大降水的频率较低,但大降水对总降水量的贡献较大,这与特日格勒等
[25]的研究一致。
浑善达克沙地2020年为枯水年,土壤水分随着降水量年际变化表现为先降低再增大的趋势。半固定沙地试验观测场,流动沙地、固定沙地试验观测场年内降水量呈单峰形态,峰值时间为6月,峰值为分别为32.80,180.61 mm,5—9月降水量占年降水量的92%以上;年内土壤水分表现为恢复期和波动期,在恢复期1—4月随着气温的升高,早春积雪融化,土壤解冻,土壤水分呈增加趋势,土壤水分维持在0.91%~11.35%,不同类型沙地土壤含水量表现为半固定沙地>固定沙地>流动沙地;波动期内由于降水量增加、土壤蒸发、植物蒸腾等因子共同作用,不同深度的土壤含水量波动存在差异,土壤水分表现为骤升骤降的脉冲式特点
[13]。
3.2 土壤水分对降水的动态响应
下渗速率反映土壤水分对降水响应时滞的重要指标,受降水量、降水历时、降水强度的直接影响,本研究结果显示,0~5 mm降水量对土壤水分补给作用微弱,未引起土壤含水量增加,魏雅芬等
[26]对降雨格局对库布齐沙漠土壤水分的补充效应研究指出,5.0 mm以下的降雨一般有增加空气湿度、降温的作用,一定程度上可以缓解旱情;5.0 mm以上的降雨才能有效补充土壤水分
[27]。随着土层深度的增加,下渗速率降幅明显,即滞后性的表现更加明显
[28]。
本研究结果显示,土壤水分对降水的响应过程表现为先快速增加后缓慢降低趋势,陈敏玲等
[29]研究指出表层土壤含水量的变化过程可用指数方程来模拟,本研究与其一致。差异表现为不同沙地类型、不同土层深度对降水的响应时滞性不同,流动沙地、半固定沙地、固定沙地20 cm土层深度土壤含水量对降水的响应峰值时间14,8,22 h,土壤峰值含水量较初始含水量增加96.43%,57.30%,65.04%;固定沙地由于植被冠层截留和凋落物层的影响,降水能被土壤吸收地较少,因此土壤含水量对降水的响应峰值时间高于其他沙地,土壤峰值水分较初始水分增加率低于其他沙地。
4 结 论
(1) 2016—2021年观测期内半固定沙地累计降水量1 049.40 mm,Ⅰ,Ⅱ级降水事件占总降水事件67.83%,降水量占总降水量的23.84%;流动沙地、固定沙地累计降水量2 135.13 mm,Ⅰ,Ⅱ级降水事件占总降水事件63.40%,降水量占总降水量的11.37%。均呈现出以Ⅰ,Ⅱ级降水为主,Ⅲ,Ⅳ,Ⅴ级降水的频率较低,但对总降水量的贡献较大。
(2) 观测期内枯水年为2020年,其流动沙地和固定沙地降水量为276.52 mm,半固定沙地降水量为130.60 mm。年内降水量均表现为单峰形态,流动沙地、固定沙地降水量峰值时间为6月(180.61 mm),半固定沙地降水量峰值时间为7月(54.00 mm)。年际间土壤平均含水量均表现为流动沙地<半固定沙地<固定沙地;流动沙地土壤含水量年内变化表现为先增大后降低趋势,土壤含水量峰值时间均为7月,半固定沙地5—9月由于植被耗水作用土壤含水量呈降低趋势;固定沙地5—9月生长季内降水量增加及植被生长耗水共同作用,土壤含水量无明显规律性变化。
(3) Ⅰ,Ⅱ级降水事件未引起土壤下渗速率变化;Ⅲ,Ⅳ,Ⅴ级降水事件土壤下渗速率均随着降水级别的增加呈增加趋势,随着土层深度的增加呈增加趋势,流动沙地土壤下渗速率为0.14~7.8 cm/h,半固定沙地土壤下渗速率为0.30~18.35 cm/h,固定沙地土壤下渗速率为0.26~11.64 cm/h。不同沙地类型土壤含水量对降水的响应均表现为先急剧增加后缓慢降低趋势,流动沙地20,40,60 cm土层土壤峰值含水量响应时间分别为14,42,86 h,较土壤初始含水量增加96.43%,48.00%,10.34%;半固定沙地20,40,60 cm土层土壤峰值含水量响应时间分别为8,12,120 h,较土壤初始含水量增加57.30%,33.78%,3.19%;固定沙地20,40,60 cm土层土壤峰值含水量响应时间分别为22,28,88 h,较土壤初始含水量增加65.04%,33.84%,6.49%。
内蒙古自治区林业科学研究院科研能力提升“揭榜挂帅”项目“‘三北’工程区林草生态产业模式研究与示范”(2024NLTS05)