Application of Leak 2D to Describe Preferential Water Flow in a Soil Containing Artificial Macropores

Leak 2D is a new two-dimensional dual permeability mathematical model for the simulation of the preferential flow in the vadose zone. In this model, water flow in the soil matrix domain is described by the two-dimensional h-based Richards’ equation. Water flow in the fracture domain is estimated using the kinematic wave approach. Richards’ equation is solved by a combination of the alternating direction implicit (A.D.I.) method and the Douglas and Jones predictor−corrector method. The kinematic wave equation is solved explicitly. In the present paper, Leak 2D is calibrated and validated with data obtained in a Hele–Shaw apparatus filled with sand. Preferential flow is achieved by inserting four artificial macropores of various sizes into the soil. Six irrigations of various intensities and durations were used for the calibration and validation process. The water content at various depths was recorded by five sensors that were inserted into the soil. A comparison of the simulated water content with the measured profiles shows that Leak 2D can sufficiently describe preferential flow into the unsaturated zone of the soil, even under extreme irrigation conditions.