peer reviewed ; Accurate precipitation maps are essential for ecological, environmental, element cycle and hydrological models that have a spatial output component. It is well known that topography has a major influence on the spatial distribution of precipitation and that increasing topographical complexity is associated with increased spatial heterogeneity in precipitation. This means that when mapping precipitation using classical interpolation techniques (e.g. regression, kriging, spline, inverse distance weighting, etc.), a climate measuring network with higher spatial density is needed in mountainous areas in order to obtain the same level of accuracy as compared to flatter regions. In this study, we present a mean total annual precipitation mapping technique that combines topographical information (i.e. elevation and slope orientation) with average total annual rain gauge data in order to overcome this problem. A unique feature of this paper is the identification of the scale at which topography influences the precipitation pattern as well as the direction of the dominant weather circulation. This method was applied for Belgium and surroundings and shows that the identification of the appropriate scale at which topographical obstacles impact precipitation is crucial in order to obtain reliable mean total annual precipitation maps. The dominant weather circulation is determined at 260°. Hence, this approach allows accurate mapping of mean annual precipitation patterns in regions characterized by rather high topographical complexity using a climate data network with a relatively low density and/or when more advanced precipitation measurement techniques, such as radar, aren't available, for example in the case of historical data. © 2016 Elsevier B.V.