[1] This is a computational study of gravity-driven fingering instabilities in unsaturated porous media. The governing equations and corresponding numerical scheme are based on the work of Nieber et al. (2003) in which nonmonotonic saturation profiles are obtained by supplementing the Richards equation with a nonequilibrium capillary pressure-saturation relationship, as well as including hysteretic effects. The first part of the study takes an extensive look at the sensitivity of the finger solutions to certain key parameters in the model such as capillary shape parameter, initial saturation, and capillary relaxation coefficient. The second part is a comparison to published experimental results that demonstrates the ability of the model to capture realistic fingering behavior.Citation: Chapwanya, M., and J. M. Stockie (2010), Numerical simulations of gravity-driven fingering in unsaturated porous media using a nonequilibrium model, Water Resour. Res., 46, W09534,