In this study, Ordinary Kriging (ok), and Adaptive Neuro Fuzzy based Inference System (anfis) are evaluated for assessing hydraulic head distribution in an aquifer unit covering 40 km2. Cartesian coordinates of the samples were used as inputs of anfis. Calibrated models are used to interpolate the hydraulic head distribution on a 50 m square - grid. Both simulations have realistic pattern (R2 > 0.97) even if ok performs slightly better than anfis at sampling location. The two methods capture different patterns. The Comparison of the two distributions allow for identifying area of estimate uncertainty, what can be used to improve the sampling network.