A hydraulically and hydrologically based estimation method of pollutant delivery load for water quality modeling is proposed. The proposed method works on grid basis and routes overland flows from one cell to the next following the maximum downslope directions. The method is able to consider spatially-varied data of source pollutant, topography, land slopes, soil characteristics, land use and aspects, which can be extracted from geographic information systems (GIS) and from digital elevation models (DEMs). Because of this feature, the proposed method can be expected to be used for evaluating the impacts of various practices on watershed management for water quality.