One of the difficulties in accurate characterization of unknown groundwater pollution sources is the uncertainty regarding the number and the location of such sources. Only when the number of source locations is estimated with some degree of certainty that the characterization of the sources in terms of location, magnitude, and activity duration can be meaningful. A fairly good knowledge of source locations can substantially decrease the degree of nonuniqueness in the set of possible aquifer responses to subjected geochemical stresses. A methodology is developed to use a sequence of dedicated monitoring network design and implementation and to screen and identify the possible source locations. The proposed methodology utilizes a combination of spatial interpolation of concentration measurements and simulated annealing as optimization algorithm for optimal design of the monitoring network. These monitoring networks are to be designed and implemented sequentially. The sequential design is based on iterative pollutant concentration measurement information from the sequentially designed monitoring networks. The optimal monitoring network design utilizes concentration gradient information from the monitoring network at previous iteration to define the objective function. The capability of the feedback information based iterative methodology is shown to be effective in estimating the source locations when no such information is initially available. This unknown pollution source locations identification methodology should be very useful as a screening model for subsequent accurate estimation of the unknown pollution sources in terms of location, magnitude, and activity duration.