For various kinds of waters including surface water, shallow groundwater (＜70 m deep) and deep groundwater (500∼810 m deep) from the Pungki area, an integrated study based on hydrochemical, multivariate statistical, thermodynamic, environmental isotopic (tritium, oxygen-hydrogen, carbon and sulfur), and mass-balance approaches was attempted to elucidate the hydrogeochemical and hydrologic characteristics of the groundwater system in the gneiss area. Shallow groundwaters are typified as the 'Ca-HCO$_3$'type with higher concentrations of Ca, Mg, SO$_4$and NO$_3$, whereas deep groundwaters are the 'Na-HCO$_3$'type with elevated concentrations of Na, Ba, Li, H$_2$S, F and Cl and are supersaturated with respect to calcite. The waters in the area are largely classified into two groups: 1) surface waters and most of shallow groundwaters, and 2) deep groundwaters and one sample of shallow groundwater. Seasonal compositional variations are recognized for the former. Multivariate statistical analysis indicates that three factors may explain about 86% of the compositional variations observed in deep groundwaters. These are: 1) plagioclase dissolution and calcite precipitation, 2) sulfate reduction, and 3) acid hydrolysis of hydroxyl-bearing minerals(mainly mica). By combining with results of thermodynamic calculation, four appropriate models of water/ rock interaction, each showing the dissolution of plagioclase, kaolinite and micas and the precipitation of calcite, illite, laumontite, chlorite and smectite, are proposed by mass balance modelling in order to explain the water quality of deep groundwaters. Oxygen-hydrogen isotope data indicate that deep groundwaters were originated from a local meteoric water recharged from distant, topograpically high mountainous region and underwent larger degrees of water/rock interaction during the regional deep circulation, whereas the shallow groundwaters were recharged from nearby, topograpically low region. Tritium data show that the recharge time was the pre-thermonuclear age for deep groundwaters (＜0.2 TU) but the post-thermonuclear age for shallow groundwaters (5.66∼7.79 TU). The $\delta$$\^$34/S values of dissolved sulfate indicate that high amounts of dissolved H$_2$S (up to 3.9 mg/1), a characteristic of deep groundwaters in this area, might be derived from the reduction of sulfate. The $\delta$$\^$13/C values of dissolved carbonates are controlled by not only the dissolution of carbonate minerals by dissolved soil CO$_2$(for shallow groundwaters) but also the reprecipitation of calcite (for deep groundwaters). An integrated model of the origin, flow and chemical evolution for the groundwaters in this area is proposed in this study.