Water quality experiments are difficult, costly, and time-consuming. Therefore, different modeling methods can be used as an alternative for these experiments. To achieve the research objective, geospatial artificial intelligence approaches such as the self-organizing map (SOM), artificial neural network (ANN), and co-active neuro-fuzzy inference system (CANFIS) were used to simulate groundwater quality in the Mazandaran plain in the north of Iran. Geographical information system (GIS) techniques were used as a pre-processer and post-processer. Data from 85 drinking water wells was used as secondary data and were separated into two splits of (a) 70 percent for training (60% for training and 10% for crossvalidation), and (b) 30 percent for the test stage. The groundwater quality index (GWQI) and the effective water quality factors (distance from industries, groundwater depth, and transmissivity of aquifer formations) were implemented as output and input variables, respectively. Statistical indices (i.e., R squared (R-sqr) and the mean squared error (MSE)) were utilized to compare the performance of three methods. The results demonstrate the high performance of the three methods in groundwater quality simulation. However, in the test stage, CANFIS (R-sqr = 0.89) had a higher performance than the SOM (R-sqr = 0.8) and ANN (R-sqr = 0.73) methods. The tested CANFIS model was used to estimate GWQI values on the area of the plain. Finally, the groundwater quality was mapped in a GIS environment associated with CANFIS simulation. The results can be used to manage groundwater quality as well as support and contribute to the sustainable development goal (SDG)-6,