In this study, seven types of first-order and one-variable grey differential equation model (abbreviated as GM (1, 1) model) were used to predict the hardness of four groundwater monitoring stations in Kaohsiung City of Taiwan. The mean absolute percentage error (MAPE) was used to evaluate the predicting performance. The results indicated the minimum MAPE of 4.71 %, 3.15 %, 2.66 %, and 16.63 % could be achieved when predicting hardness of Fonsi, Datung, Shaukang, and Chihsien stations, respectively. According to the results, it revealed that GM (1, 1) was an efficiently early warning tool for providing groundwater quality information to the competent authority.