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 rain gauges 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 method has slightly a higher correlation coefﬁcient and lower mean relative error than the other methods for the two seasons that make the method the most optimal one for mapping of rainfall distribution in the Faria catchment. Moreover, the spatial distribution of raingauges between the two seasons affects the performance of different interpolation methods in the catchment.