This study deals with review of different improvements done in the formulation of the governing equations to simulate accurate solute transport in saturated porous media over the years. The traditional advection-dispersion equation (ADE) model is the simplest lumped model founded on the assumptions of Fick’s law of diffusion. But it typically underestimates the breakthrough concentration in leading and/or tailing region due to non-fickian transport. It is modified into mobile-immobile model (MIM) considering the medium having micropores with stagnant water pockets but allowing solute exchange by diffusion between mobile and immobile zone which is quantified by mass transfer coefficient. Multi-process non-equilibrium (MPNE) model further simulates for a system with both physical and chemical nonequilibrium by assuming instantaneous and rate-limited sorption in advective and non-advective domains. Using the concept of dual permeability, slow fast transport (SFT) model divides the liquid phase in the domain into three zones i.e. fast, slow and immobile. Here chemical interaction between the fluid and soil matrix takes place only in slow and immobile zones. Non-fickian solute transport does not follow Brownian motion rules so a random variable is required to explain it. Hence continuous time random walk (CTRW) model is used where solute transport is characterized by joint probability variable. Special case of CTRW with solute having considerable probability of moving long distances and follow power law gives Fractional advection-dispersion equation (FADE) model. These models varying from relatively simple to more complex formulations and assumptions are discussed here highlighting the merits and demerits of each.