A better understanding of the hydrogeochemical evolution of groundwater in vulnerable aquifers is 18 important for the protection of water resources. To assess groundwater chemistry, groundwater 19 sampling was performed from different representative aquifers in 2012-13. A Piper trilinear diagram 20 showed that the groundwater types can be classified into Na-SO4 and Na-Cl types. Only one 21 groundwater sample was Na-HCO3 type. The dominant cations for all samples were Na +. However, 22 2 the dominant anions varied from HCO3to SO4 2-, and as well Cl-. The mean total dissolved solid 23 (TDS) content of groundwater in the region was 1,889 mg/L. Thus, only 20% of groundwater 24 samples meet Chinese drinking water standards (<1,000 mg/L). Principal component analysis (PCA) combined with hierarchical cluster analysis (HCA) and self-organizing maps (SOM) were 26 applied for the classification of the groundwater geochemistry. The three first principal components 27 explained 58, 20, and 16% of the variance, respectively. The first component reflects sulfate 28 minerals (gypsum, anhydrite) and halite dissolution, and/or evaporation in the shallow aquifer. The 29 second and third components are interpreted as carbonate rock dissolution. The reason for two 30 factors is that the different aquifers give rise to different degree of hydrogeochemical evolution 31 (different travel distances and travel times). Identified clusters for evolution characteristic and 32 influencing factors were confirmed by the PCA-HCA methods. Using information from eight ion 33 components and SOM, formation mechanisms and influencing factors for the present groundwater 34 quality were determined.