GIS is a powerful tools in analysing the multiple thematic maps at a time. GIS analysis is carried out in the present study to locate favorable shallow and deep groundwater bearing zones. To do this subsurface investigation and its interpretation were carried out through geophysical survey. Fifteen vertical electrical soundings (VES) survey were carried out in the study area. The field data were interpreted by RESIST software to determine the resistivity and thickness of the different layers. By using conventional GIS method, the spatial distribution maps for fracture zone resistivity and thickness zones were prepared. Integration of the said themes was done it in GIS. First layer fracture zone thickness and resistivity maps were overlaid and the polygon combinations were brought out. 9 combinations of 1LR-1LT, 1LR-1MT, 1LR-1HT, 1MR-1LT, 1MR-1MT, 1MR-1HT, 1HR-1LT, 1HR-1MT and 1HR-1LT were arrived, and designated as output map 1. This map shows that the shallow depth groundwater combination of (1LR-1HT and 1LR-1MT) covers an area of 50.39 Km² and 28.73 Km². Similarly, the second fracture zone thickness and resistivity maps were superposed. The result map designated as output map 2 gives a combination of 9 polygons of 2LR-2LT, 2LR-2MT, 2LR-2HT, 2MR-2LT, 2MR-2MT, 2MR-2HT, 2HR-2LT, 2HR-2MT and 2HR-2LT. This map shows that the deeper depth groundwater combinations of (2LR-2HT and 2LR-2MT) for deep groundwater zones were delineated and occupies an area of 3.16 Km2 and 13.10 Km².