Aquifer characterization is an important problem in Geohidrology. The problem consits in estimating phenomenological parameters given piezometric data. In general eld data is obtained in observation wells, which are few and scattered in the aquifer. In this work we propose a solution to the problem of parameter estimation in two parts. First we use an interpolation technique by means of radial basis functions to generate suffcient piezometric data, then we apply the Differential System method for parameter estimation. To illustrate the method, we show a characterization of the Silao-Romita aquifer located in Guanajuato Mexico.