Experimental and Numerical Studies of Saturation Overshoot during Infiltration into a Dry Soil. Issue 1 (3rd May 2019)

• Record Type: Journal Article
• Authors: Zhuang, Luwen; Hassanizadeh, S. Majid; van Duijn, C.J.; Zimmermann, Susanne; Zizina, Irina; Helmig, Rainer
• Is Part Of: Vadose zone journal. Volume 18:Issue 1(2019)
• Page Start: 1
• Page End: 13
• Publication Date: 2019-05-03
• Subjects: Soil science -- Periodicals; Zone of aeration -- Periodicals; Groundwater flow -- Periodicals; Groundwater flow; Zone of aeration; Periodicals; Electronic journals; 631.4
• DOI: 10.2136/vzj2018.09.0167
• Languages: English
• ISSNs: 1539-1663

Abstract

Abstract : Core Ideas: Experiments measured saturation overshoot during infiltration of water into dry soil. Two models based on Richards' equation with a dynamic capillarity term simulated the data. The IFA model is an alternative approach to simulate saturation overshoot. The IFA model and dynamic capillarity equation were combined to simulate experiments. Downward infiltration of water into almost dry soil, when there is no ponding at the soil surface, often occurs in the form of fingers, with saturation overshoot at the finger tips. While this is well known, there is still uncertainty about the exact saturation pattern within fingers. We performed a series of one‐dimensional water infiltration experiments into a dry soil to study the non‐monotonicity of the saturation. We observed that saturation showed a non‐monotonic behavior as a function of time. The overshoot was somewhat plateau shaped at relatively low flow rates but was quite sharp at higher flow rates. Two mathematical models, referred to as the extended standard (ESD) model and the interfacial area (IFA) model, were used to simulate the experimental results. Both models were based on extended forms of the Richards equation by including a dynamic capillary term. In the ESD model, standard equations for hysteresis were used. In the IFA model, the specific interfacial area was introduced to simulate hysteresis. Parameter values for both models were obtained from preliminary experiments or using empirical formulas. Only one parameter, the dynamic capillarity coefficient τ, was optimized to model saturation overshoot. While the ESD model did not reproduce the form of saturation overshoot for any combination of parameter values, the IFA model could provide good agreement with the data. To our knowledge, this is the first time where a combination of the IFA model and the dynamic capillarity equation has been used to simulate a set of experiments.