基于数值方法的农田土壤零通量面的模拟

零通量面能够反映土壤水分运动情况,定量计算土壤水分的入渗和蒸散量。采用Richards方程描述了降雨和蒸发边界条件下的一维土壤水分运移,通过土壤传递函数估算得到Van-Genuchten水分特征曲线参数,借助参数化Hydrus-1D模型计算得到一维土壤空间不同层的土壤含水量和基质势的数值解,并采用田间实测数据对数值解进行了验证,根据不同层的水势梯度和水流通量的值判断零通量面出现位置以及零通量面类型,最后分别设置5个不同的降雨和蒸发强度,模拟分析了土壤零通量面的空间位置和类型。研究结果表明:(1)采用Hydrus-1D可以准确模拟降雨和蒸发条件下的一维土壤水分运移;(2)在降雨和入渗条件下,零通量面位置随之降低;(3)蒸发强度增加时,零通量面上部土水势降低,土壤吸力增大,使其位置上升。结果证明该方法可行,能够模拟零通量面的动态变化。