This study presents a novel mathematical model for describing the transport of the remedial reagent in a vertical circulation flow field in an anisotropic aquifer. To develop the mathematical model, the radial and vertical components of the pore water velocity are calculated first by using an analytical solution for steady-state drawdown distribution near a vertical circulation well. Next, the obtained radial and vertical components of the pore water velocity are then incorporated into a three-dimensional axisymmetrical advection-dispersion equation in cylindrical coordinates from which to build the reagent transport equation. The Laplace transform finite difference technique is applied to solve the three-dimensional axisymmetrical advection-dispersion equation with spatial variable-dependent coefficients. The developed mathematical model is used to investigate the effects of various parameters such as hydraulic conductivity anisotropy, longitudinal and transverse dispersivities, the placement of the extraction and injection screened intervals of the vertical circulation well and the injection modes on the transport regime of the remedial reagent. Results show that those parameters have different degrees of impacts on the distribution of the remedial reagent. The mathematical model provides an effective tool for designing and operating an enhanced groundwater remediation in an anisotropic aquifer using the vertical circulation well technology. (C) 2010 Elsevier B.V. All rights reserved.