川中低山丘陵区不同方法测量土壤的饱和导水率.docxVIP

川中低山丘陵区不同方法测量土壤的饱和导水率 0 调查方法和仪器 土壤饱和导水率是土壤水动力学研究的重要参数，影响着土壤水的入渗、径向流和蒸发之间的分配关系。土壤饱和导水率易受土壤容重、质地、土壤结构、有机质含量、土地耕作利用方式等诸多因素的影响,且空间变异性强。因此,如何可靠测定土壤饱和导水率一直是土壤水动力学的一大热点问题。目前,研究土壤渗透性的方法较多,室内方法主要有定水头法和降水头法,田间测定方法则主要有单环法、双环法、圆盘渗透仪法、Guelph入渗仪法、Philip-Dunne入渗仪法等。单环刀是省去双环的外环,该法测定过程快速简单,可应用于坡耕地导水率的测定。但由于尺寸小,插入土中对土壤的扰动大。双环刀耗水量大,实际操作费力繁琐,测量精度不高,且很难运用于坡耕地上。其中单环法运用较多,入渗仪的使用相对较少。 Guelph入渗仪是利用恒定水头三维流原理快速测量田间饱和导水率的设备,其方法研究在近几年得到了快速发展,但未见该法在中国西南紫色土低山丘陵区的应用报道。因此,本文采用Guelph入渗仪测定不同土地利用方式(林地、坡耕地)、不同土壤深度(0~20、20~40 cm)的稳定下渗率,对比分析单水头法(LA法、USH法)和双水头法(TH法、BH法)所测得饱和导水率(Ks)的差异,同时比较了田间原位Guelph法与室内降水头法测定结果的差异,以期为进一步研究该地区土壤水分运动特征提供可靠的饱和导水率实测值。 1 材料和方法 1.1 中丘—研究区概况 中国科学院盐亭紫色土农业生态试验站位于四川盆地中北部的盐亭县林山乡,105°27′E,31°16′N,海拔400~600 m,属中丘—深丘地貌,植被为桤柏混交林散生乔木及黄茅草坡等。土壤以钙质紫色土为主,中亚热带气候,年均气温17.3℃,年均降雨量826 mm。 1.2 学习方法 1.2.1 压力水头法 土壤饱和导水率是土壤中的孔隙全部被水充满,土壤水分连续性最好,基质势为0时的导水率。如果不考虑溶质势和温度势,饱和土壤具有2种势的作用,即压力势与重力势,且这种势的梯度为正,在土壤结构和质地不变情况下饱和导水率是一个定值。根据对饱和条件下水分运动的驱动力的不同考虑,许多学者提出了不同的计算公式。目前,主要分为2大类:单水头和双水头法。 1)单水头法 (1)拉普拉斯法(laplace analysis,LA) Reynolds等基于拉普拉斯方程推导出饱和导水率的计算公式。他同时考虑了压力势和重力势对导水率的影响,忽略了毛细管力的作用。数学表达式如下 式中,Ks为饱和导水率,cm/s;C为形状系数,无量纲;H为压力水头,cm;Qs为水流通量,cm3/s,Qs=x R(x为储水管横截面积,cm2;R为土壤稳定下渗率,cm/s);r为洞口半径,cm。 (2)新单水头法(updated single head analysis,USH) Reynolds等利用Richards方程的近似解来处理土壤水分运动中的稳定流问题,他们同时考虑压力势、重力势和毛细管力对土壤饱和流的作用,更全面地反映了实际水流在土壤中的运移情况。其通量公式表达式为 令 则 式中,φm为土壤基质势,cm。 上述2方法的参数C取决于井的半径r(cm)和压力水头H(即井中水深,cm)。由此,Zhang等提出了3种不同土壤质地(砂土,壤土和黏土,无结构的黏土)的经验计算公式,其数学表达式如下 该方法计算初期需要根据土壤质地对α*进行人为赋值(α*≥0.12 cm-1,砂土;α*=0.04 cm-1,壤土和黏土;α*=0.01 cm-1,无结构的黏土),进而根据α*对应的经验公式计算参数C。Reynolds等指出若α*值选择得不合适,计算误差可高达25%。α*越大,C越大,Ks也因而越高。 2)双水头测量法 (1)双水头法(two head analysis,TH) 双水头法饱和导水率的计算公式如下 基质势计算公式 其中,Qs1,Qs2分别为压力水头为H1(cm)和H2(cm)对应的稳定水流通量,cm3/s。 由式(9)、式(11)有可能获得负的Ks或φm估算值。这可能是由于土壤水力性质的空间异质性所导致。当φm0,Ks0时,说明上层土壤透水性大于下层土壤;若φm0,Ks0,情况相反。若两者同时为正且Ks大于由单水头LA法得到的计算值,φm大于公式φm=CQs/2πH的计算值,则由TH法得到的Ks、φm估算值更接近真实值。 当双水头TH法中计算得到的Ks,φm有一个为负值时,不应采用双水头法。建议采用USH法(新单水法)分别求两个定水头的饱和导水率,并取2者平均值。即 其中,Ks1、Ks2分别为单水头USH法中定水头分别为H1(cm)和H2(cm)的饱和导水率,cm/s。 1.2.2 其他部件设计 Guelph入渗仪(20