Vol. 4, No. 3, 2015, pp. 127-133. doi: 10.11648/j.ajep.20150403.13 Abstract: Sophisticated prediction of aquifer performance requires numerical simulation. To date, no comprehensive simulation has been reported on groundwater modeling. Most available simulators are not applicable for fractured aquifer, and do not account for contaminant leaching and degradation, particularly in the vadose zone. Consequently, studying contaminant transport in a fractured or vuggy formation offers a forthidable challenge. This paper addresses the problem of contaminant transport in carbonate aquifer, in the presence of fractures. Most of aquifers in UAE are of limestone or carbonate origins. A series of experiments was conducted using potassium nitrate as the contaminant. Dynamic adsorption and desorption tests were carried out using both homogeneous and fractured formation models. Initial modeling and experiments were carried out for a range of initial concentration values. The concentration at the outlet was measured with the Flame Ionization technique. A numerical model was developed using the surface excess theory, combined with a non-Fickian dispersion coefficient. Numerical results agreed favorably with experimental results. It was found that the non-Fickian model was necessary for modeling fracture flow results and with this version, there was no need to use the dual porosity/dual permeability formulation. Strong dependence of adsorption on initial concentration was observed and was justified with the numerical model.