Abstarct: The need for groundwater resources to have sustainable development in a country is undoubted. In this research work, it is attempted to recognize and explore fresh groundwater resources and also, to investigate and analyze the quality of groundwater in Mashah –Chenaran plain for drinking and agricultural consumptions. For this purpose, the method of continuous vertical electrical sounding (CVES) has been used in which 132 sounding points on 7 lines have been surveyed in the area, and then, one-dimensional (1-D) and two-dimensional (2-D) modeling and interpretation of the results have been made. The electrical data have been interpreted using standard and auxiliary curves and also VES, IX1D and Res2dinv software packages. As a result, resistivity, thickness and depth of water-bearing formation or aquifer have been determined in each sounding location and also, the aquifer extent has been determined along each survey line. The interpretation results indicate that the depth of water table and also water quality (with respect to its total dissolved solids (TDS), electrical conductivity (EC) and clay content) considering variation of the resistivity of water-bearing formation are varied in the study area. In northern parts of the area, water table is located at a low depth so that it is as low as 5 meters in some parts. Towards the southern parts of the area, the depth of water table increases and reaches to a maximum of 130 meters. Also, the interpretation results indicate that the water quality is low in most of northern parts of the area is low and the water possesses remarkable clay content and salinity. In the southern parts of the area, considering higher resistivity values of the water-bearing formation, the water quality improves comparatively. The bedrock in northern parts of the area is mainly formed of shale type and is located at a low depth. However, the bedrock in southern parts of the area includes coarse-grain sediments located at a relatively high depth. Since, nowadays, most of spatial datasets are used inefficiently without employing suitable systems of data management, and thus, are wasted, an attempt has been made in this research work to digitize, organize, model and integrate in GIS environment all the data or maps (i.e. resistivity, depth and thickness of water-bearing formation) obtained from the interpretation of sounding data. For this, to model and integrate the data in GIS environment, first, each of the data or quantities comprising of resistivity, depth and thickness of water-bearing formation is classified after a preliminary processing, and then, each class is weighted separately and as a result, relevant evidence maps are obtained. In the next stage, considering different importance of each quantity in modeling and integration of the data or results to obtain the potential map of groundwater in the area, different weights are attributed to the quantities. Finally, using the capabilities of GIS and the index overlay integration method, the zones with high potential of groundwater for consumption are introduced. Furthermore, using trapezoidal membership functions corresponding to the study target (i.e. groundwater), the maps of all three quantities are fuzzified in GIS environment, and then, for integration of these fuzzy maps, considering their importance in the groundwater exploration and investigation, various fuzzy operators are used. The zones with high potential of groundwater in the area are introduces as a result of investigation of the final integration maps. Hence, fuzzy values more than 0.6 are seen in most of central and southern parts, and also in small portion of northern parts of the study area. This generally indicates a high groundwater reservoir and quality (considering high resistivity values of aquifer) and also a low depth for the water table, or an easy access to the groundwater reservoir, in these locations.