An efficient numerical scheme for fully coupled flow and reactive transport in variably saturated porous media including dynamic capillary effects.