Abstract. In this paper, results from three different interpolation techniques based on Geostatistics (ordinary kriging, kriging with external drift and conditional simulation) and one deterministic method (inverse distances) for mapping total monthly rainfall are compared. The study data set comprised total monthly rainfall from 1998 till 2001 corresponding to a maximum of 121 meteorological stations irregularly distributed in the region of Galicia (NW Spain). Furthermore, a raster Geographic Information System (GIS) was used for spatial interpolation with a 500×500 m grid digital elevation model. Inverse distance technique was appropriate for a rapid estimation of the rainfall at the studied scale. In order to apply geostatistical interpolation techniques, a spatial dependence analysis was performed; rainfall spatial dependence was observed in 33 out of 48 months analysed, the rest of the rainfall data sets presented a random behaviour. Different values of the semivariogram parameters caused the smoothing in the maps obtained by ordinary kriging. Kriging with external drift results were according to former studies which showed the influence of topography. Conditional simulation is considered to give more realistic results; however, this consideration must be confirmed with new data.