Two-dimensional advection-diffusion equation with variable coefficients is solved by the explicit finite-difference method for the transport of solutes through a homogeneous, finite, porous, two-dimensional, domain. Retardation by adsorption, periodic seepage velocity, and a dispersion coefficient proportional to this velocity are permitted. The transport is from a pulse-type point source (that ceases after a period of activity). Included are the first-order decay and zero-order production parameters proportional to the seepage velocity, periodic boundary conditions at the origin and the end of the domain. Results are compared to analytical solutions reported in the literature for special cases and a good agreement was found. The solute concentration profile is greatly influenced by the periodic velocity fluctuations. Solutions for a variety of combinations of unsteadiness of the coefficients in the advection-diffusion equation are obtainable as particular cases of the one demonstrated here. This further attests to the effectiveness of the explicit finite difference method for solving two-dimensional advection-diffusion equation with variable coefficients in a finite media, which is especially important when arbitrary initial and boundary conditions are required.