Dissolution from residual source zones of benzene poses serious threat to the groundwater quality. Proper understanding of fate and migration of dissolved benzene is a prerequisite for planning remediation strategies to reduce groundwater contamination. In the present study, an attempt has been made to numerically model the dissolution of benzene and to investigate the transport of aqueous phase benzene in a saturated fracture-matrix system under steady-state flow condition. In addition to dissolution mass transfer, advection, dispersion and matrix diffusion of aqueous benzene has been considered along the fracture. In the present numerical model, residual phase benzene is considered to be present along the entire length of the fracture and residual phase benzene is assumed to be immobile for the flow conditions considered in the analysis. Transport equations for fracture and rock-matrix are solved using implicit finite difference method. Transport equation for aqueous benzene within the fracture has been solved in a one-dimensional domain and transport equation for aqueous benzene within rock-matrix has been solved in a pseudo-two-dimensional domain. Sensitivity studies have been conducted to investigate the impact of variation of flow velocity, dispersivity, fracture aperture, inlet benzene concentration, rock-matrix diffusion coefficient, and half fracture spacing on transport of aqueous benzene concentration within the fracture. From the present study, it can be concluded that the addition of slow liquid benzene dissolution into aqueous phase influences the benzene breakthrough curves/profiles at different velocity, fracture aperture, initial benzene concentration, mass transfer rate, rock dispersivity and half fracture spacing. © 2016, Springer International Publishing Switzerland.