This study deals with transport of solutes through a saturated sub-surface rock formation with well-defined horizontal parallel fractures. For this purpose, a simplified conceptual model consisting of a single fracture and its associated rock-matrix is considered in the presence of a fracture-skin in order to study the mobility and mixing of solutes along the fracture. In this paper, a coupled fracture-skin-matrix system is modeled numerically using finite difference method in a pseudo two-dimensional domain with a constant continuous source at fracture inlet. Flow and transport processes are considered parallel to the fracture axis, while the transport processes in fracture-skin as well as in rock-matrix are considered perpendicular to the fracture axis. Having obtained the concentration distribution along the fracture, method of spatial moments is employed to study the mobility and spreading of solutes. Sensitivity analyses have been done to understand the effect of various fracture-skin parameters like porosity, thickness, and diffusion coefficient. Further, the influence of non-linear sorption and radioactive decaying of solutes are carried out for different sorption intensities and decay constants. Results suggest that the presence of fracture-skin significantly influences the mobility and spreading of solutes along the fracture in comparison with a coupled fracture-matrix system without fracture-skin. © 2012 Springer Science+Business Media B.V.