Advection-diffusion processes can be simulated by the Lattice Boltzmann method. Two formulations have been proposed in the literature. We show that they are not fully correct (only first order accurate). A new formulation is proposed, which is shown to produce better results, both from the point of view of the Chapman-Enskog expansion or when comparing simulations with an exact time-dependent solution of the advection-diffusion equation.