The aim of this study is to develop analytical solutions for one-dimensional advection-dispersion equation in a semi-infinite heterogeneous porous medium. The geological formation is initially not solute free. The nature of pollutants and porous medium are considered non-reactive. Dispersion coefficient is considered squarely proportional to the seepage velocity where as seepage velocity is considered linearly spatially dependent. Varying type input condition for multiple point sources of arbitrary time-dependent emission rate pattern is considered at origin. Concentration gradient is considered zero at infinity. A new space variable is introduced by a transformation to reduce the variable coefficients of the advection-dispersion equation into constant coefficients. Laplace Transform Technique is applied to obtain the analytical solutions of governing transport equation. Obtain results are shown graphically for various parameter and value on the dispersion coefficient and seepage velocity. The developed analytical solutions may help as a useful tool for evaluating the aquifer concentration at any position and time.