资助项目:国家自然科学基金(51679115); 兰州大学西部环境教育部重点实验室开放基金及兰州大学中央高校基本科研业务费专项基金(lzujbky-2018-kb01)
第一作者:魏宁(1981—),女,陕西扶风人,硕士,讲师,主要从事应用数学与统计学方面的研究。E-mail:weining@nwafu.edu.cn 通信作者:魏霞(1980—),女,陕西扶风人,博士,副教授,主要从事土壤侵蚀与水土保持、水文学及水资源学等方面的研究。E-mail:weix@lzu.edu.cn
(1.西北农林科技大学 理学院, 陕西 杨陵 712100; 2.兰州大学 资源环境学院, 兰州 730000; 3.Agronomy Dep., Purdue Univ., and USDA-ARS National Soil Erosion Research Lab., 275S Russell St, West Lafayette, IN 47907-2077,USA)
(1.College of Science, Northwest A&F University, Yangling, Shaanxi 712100, China; 2.College of Earth and Environmental Sciences, Lanzhou University, Lanzhou 730000, China; 3. Agronomy Department, Purdue University, USDA-ARS National Soil Erosion Research Laboratory, West Lafayette, IN 47907-2077, USA )
effects of freeze-thaw; variable flow scouring; soil critical shear stress; soil erodibility
冻融是指土层由于温度降到0℃以下和升至0℃以上而产生冻结和融化的一种物理地质作用和现象[1-2]。土壤的冻结和融化实质是土壤水的冻结和融化,即土壤中水分相态的变化过程。交替冻融指的是由于季节或昼夜热量变化在冻土表层及以下一定深度形成的反复冻结——融化的土壤过程,它作为一种自然现象普遍存在于高纬度、高海拔地区以及大多中纬度地区。交替冻融对土壤物理性质的影响,主要表现在水与冰转化过程中体积胀缩对土壤团聚结构的破坏。交替冻融对土壤物理性质的影响主要取决于交替冻融的速率、交替冻融温度、土壤含水率、土壤容重和交替冻融循环的次数等[3]。研究表明,交替冻融作用下土壤理化性质、结构和质地等均会发生改变,使土壤更容易成为侵蚀的物质来源,进而加剧土壤侵蚀程度[3-8]。既然交替冻融作用破坏土壤团聚体结构,降低土壤的团聚体稳定性[8-15],那么必然会影响到土壤的临界剪切力和土壤细沟可蚀性,进而影响土壤侵蚀量。众所周知,土壤临界剪切力和土壤细沟可蚀性是表征土壤抗侵蚀能力的两个重要参数,反映土壤对侵蚀外营力剥蚀和搬运的易损性和敏感性,是影响土壤流失量的内在因素和定量研究土壤侵蚀的基础[16-23]。然而,目前有关冻融的研究主要集中在交替冻融作用对土壤团聚体稳定性[8,10-12]、土壤容重[13-14]、土壤抗冲性等[15]方面,对冻融作用影响下的土壤临界剪切力和土壤细沟可蚀性这两个参数的研究较少涉及。鉴于此,本文利用室内冻融循环试验和变流量放水冲刷试验,研究冻融作用影响下土壤临界剪切力和土壤细沟可蚀性的变化特征,旨在揭示交替冻融作用影响下土壤自身抗侵蚀能力的两个重要参数——土壤临界剪切力和土壤细沟可蚀性的变化趋势与规律,为冻融作用影响下土壤侵蚀的预报提供重要参考。
本试验是在美国农业部土壤侵蚀国家重点实验室进行。供试土壤取自美国印第安纳州西拉法叶普渡大学动物科学研究和教育中心(87°00'54″W,40°29'24″N)的表层土壤(0—10 cm),供试土壤>0.05,0.05~0.002,<0.002 mm的粒径质量百分比为14%,66%,20%。将采集的土壤剔除植物根系和杂物后自然风干,然后过5 mm筛,配置成试验设计的土壤质量含水率(10%和20%),放置在冻融循环监测系统中进行交替冻融循环,将满足冻融循环周期设计的土壤再次风干,作为试验用土备用。冻融循环试验的最高温度为4℃,最低温度为-12℃,交替冻融的温度变化速率分别是-12℃/24 h和4℃/24 h,冻融循环的周期为1,3,6。该冻融循环格局(温度、温度的变化速率、冻融循环周期等)的设计与该土壤所在地区的主要交替冻融方式相符合。
试验土槽[22-23]尺寸为3个0.5 m×0.045 m×0.13 m底部开孔的土槽,每个土槽的两端有一个“V”型通道。装土前,在土槽下部铺填一层0.08 m厚的天然砂砾,以保持试验土的透水状况接近天然坡面,再在天然沙层上铺设一层粗纤维,然后装填试验用土,试验用土的厚度为0.05 m,试验设计坡度3%。具体试验处理见表1,除此之外,为了保证降雨前期土壤初始含水率一致,以不经受任何冻融作用的风干土为参照试验(CK)进行对照。
小流量的情况下,径流剪切力小于土壤抗侵蚀的临界剪切力,几乎不会发生土壤侵蚀,当土壤表面径流的强度增加时,径流剪切力增加,当超过土壤临界剪切力时就会产生侵蚀。因此依据土壤临界剪切力出现的实际,本试验设计从小流量向大流量依次梯度变化。冲刷试验采用逐渐增大流量的方法进行,在稳流水箱与冲刷试验槽之间安装小型流量计进行流量控制。试验过程中,总共有100,150,200,250,300,350,400,450,500,550,600,700 ml/min 12个梯级流量。试验过程中每分钟接取一个径流样品,每种流量下接取3个径流样品。冲刷方式为连续冲刷。
试验结束后,将试验径流样品称重,加入饱和明矾溶液沉淀24 h,然后倒掉上清液,将倒掉清液的径流泥沙样放进烘箱,105℃烘24 h,等彻底干燥后,取出称重,计算径流泥沙样品。根据土壤侵蚀率和水流剪切力的关系确定土壤临界剪切力和土壤细沟可蚀性。
土壤临界剪切力计算公式如下[24-26]:
Dr=Kr(τ-τc)(1)
式中:Dr为土壤分离能力[kg/(m2·s)]; Kr为土壤细沟可蚀性(s/m); τ为水流剪切力(Pa); τc为土壤临界剪切力(Pa)。
水流剪切力计算公式如下[27-28]:
τ=ρgHS(2)
式中:τ为水流剪切力(Pa); ρ为水的密度(kg/m3); g为重力加速度(m/s2); H为水深(m); S为水槽坡度(%)。
土壤分离能力计算公式如下[29-30]:
Dr=(Ms)/(tA)(3)
式中:Dr为土壤分离能力[kg/(m2·s)]; Ms为时间t内被分离的土壤干重(kg); t为坡面径流冲刷时间(s); A为采样器面积(m2)。
采用Microsoft Excel 2010和IBM SPSS Statistics 20统计软件对试验数据进行统计分析。采用单因素方差分析(one-way ANOVA)和最小显著差异法(least-significant difference,LSD)比较相同含水率不同冻融循环周期间的差异。
众所周知,土壤侵蚀的发生必须具备两个条件,即外部营力——侵蚀动力和内在物质——侵蚀对象。对于水力侵蚀而言,侵蚀的外营力即是降雨及其形成的径流,侵蚀对象即为地表土壤。径流在沿坡面梯度方向运动过程中,在其运动方向上有径流剪切力产生,径流剪切力是径流在坡面位置和时间的函数[24-25]。只有当土壤表面径流的强度增加,径流剪切力增加,超过临界剪切力τc时就会产生侵蚀。与之相伴随的临界剪力τc是表征土壤抗侵蚀能力的一个参数,其值相对比较固定,与土壤的比重、含水率、结构等因素有关,是土壤侵蚀计算式中的一个非常重要的参数。土壤细沟可蚀性Kr也是表征土壤自身抗侵蚀能力大小的一个重要因子,是评价土壤对侵蚀敏感程度的重要指标,表示了土壤被侵蚀的难易程度,反映土壤对侵蚀外营力剥蚀和搬运的易损性和敏感性。但如何确定τc和Kr,至今仍没有一个广为接受的理论或试验的方法。目前,通常的做法是将侵蚀率与径流剪应力进行回归分析,取回归直线在剪应力轴上的截距为τc值,取回归直线的斜率为Kr值。
因此,本文也采用目前比较通用的做法来确定土壤的临界剪切力τc和土壤细沟可蚀性Kr。土壤侵蚀率采用“V”型水槽的放水冲刷试验,观测得到土壤侵蚀量,进而转换得之。径流剪切力由公式(2)确定。将计算的每种试验处理,不同交替冻融循环周期和不同土壤质量含水率组合下的土壤的侵蚀率与径流剪切力点绘在直接坐标系中,见图1。可以看出,土壤侵蚀率随水流剪切力的增大而增大,土壤侵蚀率与水流剪切力呈线性关系。其函数的一般形式如下:
y=Ax-B(4)
式中:y为土壤侵蚀率[g/(m2·s)]; x为水流剪切力(Pa); A,B为系数。
为了进一步分析不同交替冻融循环周期和不同的土壤质量含水率对土壤临界剪切力的影响,将图1各种处理下的侵蚀率与径流剪切力的关系式(1)进行变形,变形为与公式Dr=Kr(τ-τc)[25]形式一致的关系式,如下所示:
y=A(x-B/A)(5)
式中:y为土壤侵蚀率[g/(m2·s)]; x为水流剪切力(Pa); A为土壤细沟可蚀性(s/m); B/A为水流临界剪切力(Pa)。也就是说对于每场试验,只有当水流剪切力大于相应的B/A时,才会有侵蚀发生,该值也为每种处理下的土壤临界剪切力值。
表2列出了根据图1中的各个拟合关系方程计算的土壤临界剪切力和土壤细沟可蚀性值。由表2可知,参照试验CK的土壤临界剪切力为0.300 Pa。当土壤质量含水率为10%,冻融循环周期为1,3,6时,对应的土壤临界剪切力的值分别为0.380,0.259,0.267 Pa。当土壤质量含水率为20%,冻融循环周期为1,3,6时,对应的土壤临界剪切力的值分别为0.338,0.238,0.242 Pa。分析发现,两种土壤质量含水率下,土壤经过1个周期的冻融循环后,土壤临界剪切力与参照试验的相应值相比有所增大,当土壤经过3个周期的冻融循环后,土壤临界剪切力显著降低,随着土壤的冻融循环周期增至6时,土壤临界剪切力又有所增加,但其值仍小于参照试验CK的土壤临界剪切力。
由表2可知,参照试验CK的土壤细沟可蚀性为0.003 0 s/m,当土壤质量含水率为10%,冻融循环周期为1,3,6时,对应的土壤细沟可蚀性的值分别为0.007 9,0.002 7,0.001 5 s/m,即当土壤经历1个周期的冻融循环后,土壤细沟可蚀性相对参照试验显著增大,当土壤经历3个和6个周期的冻融循环后,土壤临界剪切力降低,且均小于参照试验CK的土壤细沟可蚀性值。同理,当土壤质量含水率为20%,冻融循环周期为1,3,6时,对应的土壤细沟可蚀性的值分别为0.006 8,0.002 1,0.003 3 s/m,即当土壤经历1个周期的冻融循环后,土壤细沟可蚀性增大,当土壤继续经历冻融循环作用,其周期增大至3时,土壤细沟可蚀性出现明显的减小趋势,当冻融循环周期增大至6时,土壤细沟可蚀性又有所增加,但冻融循环周期为3,6时,土壤的细沟可蚀性均小于参照试验。
由以上分析可知,交替冻融循环周期对土壤临界剪切力和土壤细沟可蚀性的影响显著,且随着交替冻融循环周期的不同,土壤临界剪切力和土壤细沟可蚀性的影响趋势与程度不同。两种土壤含水率水平下,当交替冻融循环周期为1时,交替冻融作用使土壤临界剪切力和土壤细沟可蚀性显著增大,当交替冻融循环周期为3,6时,交替冻融作用使土壤临界剪切力和土壤细沟可蚀性减小。
由表2可知,当交替冻融循环周期为1时,10%和20%两种不同土壤含水率下的土壤临界剪切力的值分别为0.380,0.338 Pa,均大于参照试验的0.300 Pa,即在两种不同含水率下,经过1个周期的冻融循环后,土壤临界剪切力增大,且低含水率时土壤临界剪切力增大的幅度较高含水率时土壤临界剪切力增加的幅度大。当交替冻融循环周期为3时,10%和20%两种不同土壤含水率下的土壤临界剪切力分别为0.259,0.238 Pa,均小于参照试验的0.300 Pa,且低含水率下的土壤临界剪切力较高含水率下土壤临界剪切力的值大。当交替冻融循环周期为6时,10%和20%两种不同土壤含水率下的土壤临界剪切力分别为0.267,0.242 Pa,均小于参照试验的0.300 Pa,且10%含水率下的土壤临界剪切力较20%含水率下的土壤临界剪切力大。
同理,当交替冻融循环周期为1时,10%和20%两种不同土壤含水率下的土壤细沟可蚀性分别为0.007 9,0.006 8 s/m,均大于参照试验的0.003 0 s/m,即在两种不同含水率处理下,土壤经过1个周期的冻融循环后,土壤细沟可蚀性增大,且低含水率时土壤细沟可蚀性增大的幅度较高含水率时土壤细沟可蚀性增大的幅度大。当交替冻融循环周期为3时,10%和20%两种不同土壤含水率下的土壤细沟可蚀性分别为0.002 7,0.002 1 s/m,均小于参照试验的0.300 Pa,且低含水率下的土壤细沟可蚀性较高含水率下的土壤细沟可蚀性大。当交替冻融循环周期为6时,10%和20%两种不同土壤含水率下土壤细沟可蚀性分别为0.001 5,0.003 3 s/m,均小于参照试验的0.300 s/m,且10%含水率下的土壤细沟可蚀性较20%含水率下的土壤细沟可蚀性大。
综上可知,不同土壤含水率对土壤临界剪切力和土壤细沟可蚀性的影响程度不同,土壤质量含水率水平为10%时的土壤临界剪切力和土壤细沟可蚀性的变化幅度均大于土壤质量含水率水平为20%时的土壤临界剪切力和土壤细沟可蚀性的变化幅度。
本研究是探索交替冻融作用对土壤临界剪切力和细沟可蚀性的影响,结果表明,交替冻融能够改变土壤的临界剪切力和细沟可蚀性。当冻融循环周期为1时,10%和20%含水量处理下的土壤临界剪切力和细沟可蚀性相比对照组都显著增大。随着冻融循环周期继续增大至3时,两种含水量处理下的土壤临界剪切力和细沟可蚀性显著降低,且明显低于对照组。当冻融循环周期继续增大至6时,10%含水量的土壤临界剪切力反而显著增大,细沟可蚀性显著减小,但其值小于控制组; 20%含水量的土壤临界剪切力和细沟可蚀性显著增大。并且随着冻融循环周期的增加,土壤临界剪切力和细沟可蚀性的变化速率减小,即随着冻融循环周期的增加,土壤临界剪切力和细沟可蚀性趋于稳定。已有研究结果也表明,历经一次冻融循环,土壤的团聚体稳定性显著增大,随着冻融循环周期的增加,土壤团聚体的稳定性降低,最终趋于稳定。
此外,土壤临界剪切力和土壤细沟可蚀性是表征土壤抗侵蚀能力的两个重要参数,土壤细沟可蚀性越大,土壤越易侵蚀,土壤临界剪切力越小,土壤越易侵蚀。从理论上来讲,土壤临界剪切力和细沟可蚀性对交替冻融的响应应该相反。但本研究的结果表明,在10%含水量时,随着冻融循环周期的增加,土壤临界剪切力和细沟可蚀性的变化趋势相反,但在20%含水量时,随着冻融循环周期的增加,土壤临界剪切力和细沟可蚀性的变化趋势相同。这可能是由于初始冻融含水量相比冻融循环周期对土壤临界剪切力和细沟可蚀性影响较大缘故,已有研究也表明,对于同一种土壤,团聚体稳定性可能随着F-T循环的增加而增加或降低,具体取决于冻融初始土壤含水量[10-12]。
(1)交替冻融循环周期不同,其对土壤临界剪切力和土壤细沟可蚀性的影响程度和趋势也不同,当冻融循环周期为1时,交替冻融作用使土壤临界剪切力和土壤细沟可蚀性增大,当冻融循环周期为3,6时,交替冻融作用使土壤临界剪切力和土壤细沟可蚀性减小。
(2)土壤含水率不同,其对土壤临界剪切力和土壤细沟可蚀性的影响程度不同,本研究所涉及的两种土壤质量含水率下,当交替冻融循环周期相同时,低含水率(10%)下的土壤临界剪切力与土壤细沟可蚀性均大于高含水率(20%)下的土壤临界剪切力与土壤细沟可蚀性。
(3)各种试验处理下的土壤侵蚀率与径流剪切力之间的关系仍然呈良好的线性关系。不同处理下,交替冻融作用对土壤临界剪切力与土壤细沟可蚀性的影响趋势相同。
致谢:感谢同在美国农业部土壤侵蚀国家重点实验室吴秋菊博士和杨韶洋硕士在试验过程中给予的帮助!