A numerical study of a negatively buoyant plume intruding into a neutrally stratified porous medium has been undertaken using finite different methods. Of particular interest has been to ascertain whether the experimentally observed gravitational instabilities that form along the lower edge of the plume are reproduced in the numerical model. The model has been found to faithfully reproduce the mean flow as well as the gravitational instabilities in the intruding plume. A linear stability analysis has confirmed the fact that the negatively buoyant plume is in fact gravitationally unstable and that the stability depends on two parameters: a concentration Rayleigh number and a characteristic length scale which is dependent on the transverse dispersivity. ¿ American Geophysical Union 1994