Diffuse pollution is usually temporally and spatially uncertain, and thus hard to analyze. In many cases, discretizing a diffuse source of pollution into individual point sources can ease diffuse pollution modelling and analysis, and therefore reduce high uncertainty especially in the spatial distribution of pollution loads. This is however a difficult task, since quite a number of sub-drainage areas, wit h complex structures and land-use properties, has to be delineated. Watershed models can be used to delineate the sub-drainage areas in a watershed with high accuracy and locate the related outlets which connect the sub-drainage areas to the main waterbody in a watershed. In this study, such an approach has been used on a case study to model the diffuse nutrient loads carried to streams that reach to a medium-sized lake in Turkey. The annual nutrient loads, which were calculated by using mathematical models, were then converted to a load-map with the help of a geographical information system.