  We introduce in this paper two time discretization schemes tailored for a
range of Wasserstein gradient flows. These schemes are designed to preserve
mass, positivity and to be uniquely solvable. In addition, they also ensure
energy dissipation in many typical scenarios. Through extensive numerical
experiments, we demonstrate the schemes' robustness, accuracy and efficiency.
