Abstract. A simple numerical solution procedure – namely the method of lines combined with an off-the-shelf ordinary differential equation (ODE) solver was shown in previous work to provide efficient, mass-conservative solutions pressure-head form Richards' equation. We implement such a our model openRE. developed novel quantify boundary fluxes that reduce water balance errors without negative ...