We investigate effectiveness of an acceleration method applied to the modified Picard iteration for simulations of variably saturated flow. We solve nonlinear systems using both unaccelerated and accelerated modified Picard iteration as well as Newton’s method. Since Picard iterations can be slow to converge, the advantage of acceleration is to provide faster convergence while maintaining advan...