When the non-aqueous phase liquid (NAPL) contaminant is entrapped in soil pores, their release to the subsurface environment is mainly limited by interphase mass transfer such as dissolution and sorption. Sorption onto the aquifer material plays an important role as a retardation mechanism in subsurface contaminant transport processes. Considering the heterogeneity associated with aquifer properties as well as contaminant source distribution, long-term contamination is an inevitable consequence of mass transfer limitation of residual mass from soil micropores. A comprehensive numerical model is presented in this study for understanding the kinetic nature of dissolution and sorption of hydrophobic hydrocarbons and their interactions in order to estimate the rate and extent of mass removal from different phases. The results showed that sorption by soils and sediments having different physicochemical properties generally follow non-linear behavior. It is observed from the study that sorption non-linearity can be effectively incorporated by the two-site kinetic model rather than combination of linear and non-linear sorption isotherms. Even though initial phase of mass transfer is dissolution-dominated, the extended tailing of concentration at higher pore volumes is controlled by non-equilibrium sorption, which is better explained by the two-site kinetic model. © 2014, Saudi Society for Geosciences.