The spatial pattern of precipitation is a complex variable that strongly depends on other geographic and topographic factors. As precipitation is usually known only at certain locations, interpolation procedures are needed in order to predict this variable in other regions. The use of multivariate interpolation methods is usually preferred, as secondary variables generally derived using GIS tools correlated with precipitation can be included. In this paper, a comparative study on different univariate and multivariate interpolation methodologies is presented. Our study area is centred in the region of Valencia, located to the eastern Spanish Mediterranean coast. The followed methodology can be divided in three steps. First, secondary variables having significant correlations with the precipitation were derived, where the hillsides were used as influence areas of certain variables. Secondly, precipitation was interpolated with different methodologies. Finally, the derived models were compared in terms of predicted errors. Models were achieved for seasonal scales, considering a total of 179 raingauges; data of another 45 raingauges were also used to predict errors. Results prove that there is no ideal method for all the cases but it will depend on one hand, on the number of geographical factors that influence the rainfall and, on the other hand, on the major or minor spatial correlation within the rainfall.