The present work solves two-dimensional Advection-Dispersion Equation (ADE) in a semi-infinite domain. A variable source concentration is regarded as the monotonic decreasing function at the source boundary (x=0). Depth-dependent variables are considered to incorporate real life situations in this modeling study, with zero flux condition assumed to occur at the exit boundary of the domain, i.e....