Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection-diffusion equation. In the last years, special attention has been devoted to the task of searching analytical solutions for the advection-diffusion equation in order to simulate the pollutant dispersion in the planetary boundary layer (PBL). Presently, analytical solutions of the advection-diffusion equation are usually obtained by making strong assumptions about the eddy diffusivity coefficients and wind speed profiles. In this work, we present a general solution (i.e., for any wind and eddy diffusivity vertical profiles) of the two-dimensional steady-state advection-diffusion equation using the general integral Laplace transform technique (GILTT).