Abstract In this present work, a well-structured and limit-based derivative definition of fractional derivative term, known as conformable derivative, is employed to develop a local generalized form of time variable fractional convection-dispersion equation (FCDE). The fractional models of convection-dispersion equation have been widely established as more authentic mode to characterize the pollutant transport in geological structures. In this work, models are formulated corresponding to spatial dependency of velocity and dispersion coefficient and also for the temporally varying decay rate, that dealt with more realistic phenomenon of the pollutant transport in groundwater reservoir. Due to the non-linearity of the problem, homotopy analysis method (HAM) is adopted to investigate the complex solutions of FCDE. As the FCDE comprises the classical convection-dispersion equation (CDE) as a special case corresponding to the fractional order α = 1 , so the obtained solutions are validated by the corresponding numerical solution and exact analytical solution for α = 1 and solutions are also verified for different fractional values of α . The effect of conformable derivative order α is properly visible over concentration strength. Hence, this investigation helps to interpret the accurate description of time dependent behaviour of contaminant transport in porous structure.