Summary: We present a transient 2D coupled vertically averaged flow/transport model. The model deals with all bed geometries and guarantees global conservation and positive values of both water level and solute concentration in the transient solution. The model is based on an upwind finite volume method, using Roe’s approximate Riemann solver. A specific modification of Riemann solver is proposed to overcome the generation of negative values of depth and concentration, that can appear as a consequence of existing wetting/drying and solute advance fronts over variable bed levels, or by the generation of new ones when dry areas appear. The numerical stability constraints of the explicit model are stated incorporating the influence of the flow velocity, the bed variations and the possible appearance of dry cells. Faced to the important restriction that this new stability condition can impose on the time step size, we are able to present a different strategy to allow stability using a maximum time step, and in consequence a minimum computational cost.