The protection of groundwater resources is worldwide considered a task of high importance and imperative necessity due to their rapid deterioration and the subsequent risk of public health. A task that in its most common practical implementation is equivalent to the protection of water-supply wells through methods and techniques that result in the determination of site-specific protection areas around the individual wells. In the interior of these areas activities that might deteriorate the fresh water physical status, and particularly its chemical properties, are controlled and very often prohibited. In this study, we delineate water field protection zones for the Lusaka aquifer system. A numerical modelling process that consists of: a) the simulation of groundwater flow in the Lusaka Dolomite aquifer by applying the three-dimensional finite difference model MODFLOW and b) the delineation of protection zones for the domestic water-supply wells by applying the particle tracking post-processing package MODPATH. Regarding the flow model, a steady state model was developed and a calibration process was performed. The estimation of the hydraulic conductivity distribution was optimized using a trial-and-error technique. With regard to the particle-tracking process, two simulations were performed. The first simulation tracked from the water supply wells, while the second one tracked from point pollution sources. Finally, suggestions are given for the control and management of the identified potential point sources of pollution, especially if they are located within the protection zones of the examined water-supply wells.