The lattice Boltzmann method (LBM) has established itself as a valid numerical method in computational fluid dynamics. Recently, multiple-relaxation-time LBM has been proposed to simulate anisotropic advection-diffusion processes. The governing differential equations of advective-diffusive systems are known to satisfy maximum principles, comparison principles, the non-negative constraint, and t...