原状土和填装土的水分特征与土壤孔隙分布.docx

原状土和填装土的水分特征与土壤孔隙分布 土壤含水量特征曲线是土壤含水量与基本势之间的关系曲线，它取决于土壤结构和土壤侵蚀的分布，直接或间接地与土壤通气、导水性和水分对作物的有效性有关。这是土壤水分和有机质运动、植物根系分布和土壤通气的最重要指标之一。沙箱法可以测定低吸力条件下(0~100hPa)土壤基质势与含水量的关系;压力板法可以测定高吸力条件下(100~15 000hPa)土壤基质势与含水量的关系。为此,同时利用沙箱和压力板法来测定土壤脱湿曲线,以明确原状土和填装土水分特征曲线及土壤孔隙分布的差异,为合理进行农业耕作提供依据。 1 材料和方法 1.1 原环刀的检测 试验所用土壤样品采自河南省农业科学院原阳试验田。用1cm高小环刀从田间采取原状土,底部用滤纸包被,并用橡皮筋固定用于测定原状土的水分特征曲线;待原状土测定完毕,把原状土磨碎并过2mm筛,重新填装到原环刀中,用于测定填装土的水分特征曲线。样品共有3个重复。 1.2 测定参数函数 采用沙箱法测定以下低吸力点:0、10、20、40、60、80hPa;压力板法测定以下高吸力点100、300、500、1 000、15 000hPa。采用van Genuchten建立的θ(h)函数式(1)拟合实际测定的土壤含水量和基质势数据,同时利用Excel中的规划求解命令,求得参数的最优值。 式中,θ为含水量;θsat为饱和含水量(基质吸力等于0时的土壤含水量);θres为残余含水量(基质吸力在15 000hPa时的土壤含水量),θres被设定为常数;h为基质吸力;θsat、m、n、α是需要拟合的参数[α0,m=1-1/n(n1,0m1)]。 1.3 土壤中的空隙大小和空隙度的计算 土壤孔隙当量孔径(D)计算公式为D=300/h。与当量孔径对应的孔隙度由土壤水分特征曲线上相应的含水量差计算求得。 2 结果与分析 2.1 不同原状土和填装土的拟合曲线 依据van Genuchten模型得到原状土和填装土的拟合参数(表1)和拟合曲线(图1)。由图1可以看出,无论是测定值还是拟合值,原状土和填装土水分特征曲线均存在明显差异。原状土的测定曲线和拟合曲线基本重合,拟合复相关系数(R2)高达0.999,而填装土的R2为0.981(表1)。 从图1可以看出,原状土和填装土的测定曲线相交于pF(土壤基质势绝对值的常用对数)为2.7处,当0≤pF2.7时,填装土含水量大于原状土;而当2.7pF≤4.2时,原状土含水量大于填装土。原状土和填装土的拟合曲线大致相交于pF为3.0处,当0≤pF3.0时,填装土含水量大于原状土,而当3.0pF≤4.2时,原状土含水量大于填装土。由此可见,当土壤颗粒大小分布和结构排列发生变化时(即原状土遭到机械破碎变成填装土),其在同一吸力下的持水能力产生差异。另外,原状土曲线的平滑性和自然弧度明显比填装土流畅,说明填装过程破坏了土壤孔隙结构的连续性。这种情况在频繁犁耕时易发生,因此,应合理安排犁耕强度与耕作频率,防止土壤孔隙连续性变差而影响作物有效吸水。 土壤水分特征曲线反映了土壤水和土壤固体的相互作用关系,是分析土壤水保持和运动最基本的资料之一,主要受土壤颗粒组成和孔隙结构等因素影响。本研究中,原状土和填装土的土壤颗粒组成相同,因此,导致两者水分特征曲线差异的原因主要是土壤结构的变化,即由土壤颗粒大小和排列变化引起的土壤孔隙含量分布和连续性的变化。 2.2 孔隙度的变化 由表2可以看出,土壤由原状土变成填装土后,孔隙大小、分布也随之变化,主要表现在50μm孔隙的孔隙度主要呈减少趋势,由原来的19.6%减少到17.8%,减少了9.2%,而37.5μm孔隙的孔隙度呈增加趋势,由原来的17.0%增加到29.2%,增加了71.7%。其中,150~300μm、75~150μm和50~75μm孔隙的孔隙度分别降低7.3%、32.4%和29.6%,而10~37.5μm、6~10μm和3~6μm孔隙的孔隙度分别增加了0.4倍、4.8倍和2.0倍。 前人定义孔径大于150μm的孔隙为在重力作用下可以自由排水的孔隙,30~150μm的孔隙为导水孔隙,0.2~30μm的孔隙为储水孔隙。导水孔隙和储水孔隙对改善土壤-作物-水分的关系、维持良好的土壤结构、保持作物和微生物利用的有效水分均起着非常重要的作用。由表2可以看出,土壤由原状变成填装后,重力排水孔隙的孔隙度基本没有变化,由原来的13.5%变为13.6%;导水孔隙(37.5~150μm)的孔隙度由7.6%减少到5.8%,减少了23.7%;储水孔隙(2~37.5μm)的孔隙度由17.0%增加到29.2%,增幅高达70.8%。总之,填装土导水和储水孔隙的总孔隙度较原状土增加(从24.6%增加到35.0%),增幅达42.3%,看似土壤导水和储水整体能