Development of hot spring sites in Malaysia has been targeted by authorities to enhance its touristic and recreational assets. Previous experiences showed that some hot springs were damaged, and the flow of water had stopped, or the water temperature had decreased. To solve this problem, the knowledge of the subsurface structural model may be needed. The modelling should be carried out using non-invasive methods according to the present land-use conditions. Hence, the geophysical methods are adopted as they have two advantages; first, they are non-invasive; second, they are cost effective. The present work was conducted to aide in both modelling the subsurface structures and guiding the development plans of the recreation complex. Besides, the study focusses on the shallow part of the geothermal system that pave the way for more deeper study to weigh the potential of establishing a geothermal power plant near the site. Geophysical methods used in this study are the 2D ERT and potential fields (magnetic and gravity). Pseudo 3-D ERT block was created from the 2-D profiles to compare with the of the potential fields. The final model of the area deduced that the area is affected by a shear zone. The present-day hot spring is located near the boundary of the shear zones. The model also shows that there are several sites in the complex with stable subsurface that can be used for development without affecting the hot spring system.