The paper presents an approach that extends the flexibility of the standard lattice Boltzmann single relaxation time scheme in terms of spatial variation of dissipative terms (e.g., diffusion coefficient) and stability for high Péclet mass transfer problems. Spatial variability of diffusion coefficient in SRT is typically accommodated through the variation of relaxation time during the collision step. This method is effective but cannot deal with large diffusion coefficient variations, which can span over several orders of magnitude in some natural systems. The approach explores an alternative way of dealing with large diffusion coefficient variations in advection-diffusion transport systems by introducing so-called diffusion velocity. The diffusion velocity is essentially an additional convective term that replaces variations in diffusion coefficients vis-à-vis a chosen reference diffusion coefficient which defines the simulation time step. Special attention is paid to the main idea behind the diffusion velocity formulation and its implementation into the lattice Boltzmann framework. Finally, the performance, stability, and accuracy of the diffusion velocity formulation are discussed via several advection-diffusion transport benchmark examples. These examples demonstrate improved stability and flexibility of the proposed scheme with marginal consequences on the numerical performance.