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. Moreover, large eddies are able to mix scalar quantities in a manner that is counter to the local gradient. We present a general solution of a two-dimension steady state advection-diffusion equation, considering non-local turbulence closure using the General Integral Laplace Transform Technique. We show some examples of applications of the new solution with different vertical diffusion parameterisations.