Changes in groundwater quality, is caused destruction of other resources, either directly or indirectly, due to bad management of groundwater extractions. Considering the importance of groundwater resources, especially in arid and semiarid regions, this research was done to modeling the spatial distribution of some groundwater qualitative factors with emphasis on drinking using geostatistical methods. Qualitative variables selected are: total hardness, chloride, sulfate and nitrate ions. First of all, normality of data investigated based on Kolmogrov-Smirnov test and non-normal variables were transformed to normal by using logarithmic transform. Then variography and cross-variogram analysis of variables was done. Accuracy assessment of results was investigated based on cross validation and root mean square error (RMSE). Results showed that Cokriging Method (using auxiliary variables) compared with other Geostatistical methods is a more accurate estimator, so interpolating and zoning of variables were done by using Cokriging Method and drinking water standards.