In this paper, we present a new analytical method for solving problems in water pollution. The method provides a simulation using two differential constraints of first order, from which are found auto-Bäcklund transformations for the equation in two-dimensional advective-diffusive steady state. The main feature of the proposed formulation is the reduced processing time required to obtain the analytical solutions. Results of numerical simulations are presented in graphical format