模拟降雨条件下含碎石土壤入渗过程的研究.docxVIP

模拟降雨条件下含碎石土壤入渗过程的研究 土壤入渗和流量特征可以反映梯度过程和坡面输沙能力，这是代表降水侵蚀的重要参数。为深入研究坡面侵蚀过程,很多学者建立了坡度、坡长、坡形、地表覆盖以及土地利用方式与土壤侵蚀量的关系,并以此来预测试验区的土壤侵蚀量。但相关研究大多是基于均质土壤,较少考虑土壤颗粒组成的复杂性。由于成土过程和人类活动的影响,含砾石土壤在自然界分布广泛,如我国西南喀斯特地区的土石山区,黄土高原北部地区成土过程中出现的风化程度不同的钙质结核等。砾石能影响土壤结构和水力学特性,进而改变坡面水土流失过程。国外研究表明,地表土壤中砾石对入渗和径流的影响存在不确定性。国内研究大多是在室内用一维垂直定水头入渗法描述土石混合介质的入渗过程,着重分析砾石含量和直径对土壤水分、溶质运移的影响。还有研究认为在模拟降雨条件下土壤入渗率随砾石含量变化可能存在一个阀值。很少有学者分析含砾石土壤的坡面产流和入渗特征,相关研究不仅是量化非均质土壤坡面水土过程的基础,而且在土壤侵蚀、坡面水文模型研究和实际应用中具有一定的价值。 本研究在模拟降雨条件下,对不同坡度下含砾石土壤的坡面产流和入渗特征进行了分析,旨在探讨砾石含量和坡度对坡面径流和水分入渗的影响及其可能原因,为该类型土壤的侵蚀过程以及水土资源的合理利用提供科学依据。 1 材料和方法 1.1 土样与砾石的制备 试验在黄土高原土壤侵蚀与旱地农业国家重点实验室人工降雨大厅进行。试验土槽为铁制,规格3 m(长)×1 m(宽)×0.25 m(深)。供试土壤取自西北农林科技大学试验地表层塿土(钙积土垫旱耕人为土),土壤质地为粉壤土,有机质含量14 g kg-1,土壤颗粒组成见表1(国际土壤质地分类系统)。土样风干后测定其风干质量含水率(平均为3%),过5 mm筛,试验所用砾石取自渭河滩,为光滑椭圆形硬质鹅卵石,对水几乎不通透。在将土样和砾石混合前,去掉较大的石块,过20 mm筛,留取直径为20 mm~50 mm的砾石作为供试砾石,砾石平均直径为4.07 cm,砾石级配比例见表2。用水冲洗,除去砾石上附有的杂物,风干备用。 1.2 土石质量浓度和降雨时坡面流速 试验中,土壤容重控制为1.35 g cm-3,实测砾石密度为2.65 g cm-3;风干土与砾石按质量比设10∶0、9∶1、8∶2、7∶3和6∶4等5个处理,孔隙度分别为0.49、0.44、0.39、0.34和0.29。然后将土壤和砾石混合均匀(砾石质量含量则分别为0%、10%、20%、30%和40%)。将土石混合物分层装入土槽,为防止土体垂直分层,每两层之间打毛。 降雨过程中雨强控制在60 mm h-1左右(实际按率定雨强来计算),土槽设5°、10°、15°等3个坡度处理。降雨开始后记录土槽产流时间。从产流开始后每隔1 min收集一次径流,每5 min测定一次坡面水流流速(采用KMnO4染色剂法,记录水流流过固定坡长的时间,然后换算成流速)。降雨历时为产流后持续60 min。降雨结束后,用称重法和烘干法分别得到每分钟的径流量和产沙量。入渗率(i)用以下公式计算: 式中,i-入渗率(cm min-1);R-降雨量(cm);T-降雨历时(min);θ-土槽坡度(°);F-单位时间内的产流量(g);S-单位时间内的产沙量(g);k-换算系数(1g cm-3);A-土槽截面积(cm2)。 2 结果与分析 2.1 不同坡度下产流时间的变化 产流时间是坡面综合效应的反映,主要取决于土壤初始含水率。在降雨的最初阶段,降雨全部渗入土壤,坡面无径流产生。当雨强大于土壤入渗能力时,产生超渗雨。超渗雨形成坡面积水,然后沿坡面向低处流动,称为坡面漫流。影响产流时间的因素有土壤特性、地面坡度、地表覆盖等。图1为模拟降雨条件下坡面产流时间随砾石含量的变化关系。 图1显示,在坡度为5°时,随着砾石含量的增加,产流时间变化不大,当砾石含量较高时,产流时间急剧降低;而在10°和15°时,产流时间随砾石含量的增加而降低。3个坡度下,当砾石含量相同时,产流时间(T)随坡度变化大致呈T15°T10°T5°;当坡度相同时,产流时间随砾石含量变化大致呈T0%T10%T20%T30%T40%,这种变化趋势在坡度为10°和15°时特别明显。其原因在于:在降雨的最初阶段,主要是雨滴浸润土壤和砾石的过程。由于砾石的存在,一部分雨滴击溅砾石,使砾石表面湿润。随着时间的推移,砾石包裹一层水膜使土壤与砾石的周边接触成为入渗点,进而逐渐增大了水分湿润土壤的面积,入渗随之增加。当砾石被完全浸润后,入渗不再以入渗点的形式进行。这实际上是由于砾石的存在增大入渗孔隙的弯曲度,从而延缓了入渗,增加了坡面汇流,使产流时间减小。而且当土壤中砾石含量越大,坡面暴露的砾石也就越多,砾石的这种作用也就越明显。当砾石含量一