Thermal springs are the sites where the water temperature lies above ambient temperature. They are widely used for power generation, hot water spas, balneotherapy, agriculture, laundering, and aquaculture. In Nepal, many thermal springs are reported but scientific understanding on water quality and hydrogeochemistry of the springs is very limited. In this study, a total of 28 physico-chemical parameters were measured in water samples collected from 12 thermal springs from Gandaki Province, Nepal. Correlation matrix and multivariate statistical analysis such as principal component analysis (PCA) and cluster analysis were used to understand the water quality and hydrogeochemistry of the hot water springs. The pH, temperature, electrical conductivity, total dissolved solids, and turbidity in the hot water springs ranged from 7.3 to 8.8, 31.6-64.3 °C, 206-16270 μS/cm, 115-6637 mg/L, and 0.21-63.7 NTU; respectively. The dominance order of major anions and cations were: Cl > HCO > SO > NO > F and Na > Ca > K > Mg; respectively. Comparison of the water quality parameters with the WHO and National Water Quality Standards suggested that the majority of the parameters were within the safe limit. Out of 9 heavy metals and trace elements analyzed Zn, Ni, Cr, Cd, Hg (total), and Pb were found below the safe limit but Fe (total), As (total) and Cu (total) were found higher than the WHO safe limit in total of 3, 5 and 1 sampling sites; respectively. The water quality index (WQI), sodium absorption ratio (SAR), and magnesium hazard (MH) ratio in the hot water springs ranged from 40.9 to 573, 2.2-49.3, 7.1-70.8; respectively. The result of PCA analysis showed that four principal components are required to explain hydrogeochemistry. Cluster analysis suggested that the sampling sites can be grouped into three distinct clusters based on total dissolved solids. Interestingly, the classification of hydrochemical facies using a Piper diagram suggested that 7 out of 12 thermal springs have Na-Cl type water. Finally, a perspective on the suitability of the hot springs for hot water spas and balneotherapy and policy recommendation is provided.