The availability of rainfall data is essential for the prediction of extreme rainfall events, hydrological modeling and design of water resources systems. In Palestine, the limited available raingauges and the available rainfall data are insufficient to precisely characterize the areal rainfall distribution over the entire area. Hence, it is essential to find proper methods to utilize the available rainfall data at the existing raingauges to predict rainfall values at certain locations where rainfall measurements are not available. This study aims at identifying the most suitable interpolation method that can be used to produce interpolated rainfall surfaces at the Faria catchment, Palestine. Faria catchment (320 km2), located in the northeastern part of the West Bank, Palestine, is characterized by high temporal and spatial rainfall variations. Six, GIS-based, spatial interpolation methods have been tested by visual observation and cross validation. These methods are: Inverse Distance Weighted, Spline completely Regularized, Spline with Tension, Ordinary Kriging, Universal Kriging, and Regression. Available seasonal rainfall data recorded at 16 locations for two rainy seasons; 2016/2017 and 2017/2018 were utilized. Results indicate that the Universal Kriging and Inverse Distance Weighted are the most suitable interpolation methods for the rainy seasons 2016/2017 and 2017/2018, respectively. Moreover, the spatial distribution of rain gauges highly affects the suitability of different interpolation methods in the catchment.