The electrical resistivity tomography method has been widely used in geophysics for many purposes such as determining geological structures, water movement, saltwater intrusion, and tectonic regime modeling. Karstic springs are important for water basin management since the karst systems are highly complex and vulnerable to exploitation and contamination. An accurate geophysical model of the subsurface is needed to reveal the spring structure. In this study, several karst springs in the Gökova Bay (SW, Turkey) were investigated to create a 3D subsurface model of the nearby karstic cavities utilizing electrical resistivity measurements. For this approach, 2D resistivity profiles were acquired and interpreted. Stratigraphically, colluvium, conglomerate, and dolomitic-limestone units were located in the field. The resistivity values of these formations were determined considering both the literature and field survey. Then, 2D profiles were interpolated to create a 3D resistivity model of the study area. Medium-large sized cavities were identified as well as their locations relative to the springs. The measured resistivities were also correlated with the corresponding geological units. The results were then used to construct a 3D model that aids to reveal the cavity geometry in the subsurface. Additionally, several faults are detected and their effect on the cavities is interpreted.