Assessing the diffuse pollutant loadings at watershed scale has become increasingly important when formulating effective watershed water management strategies, but the process was seldom achieved for heavy metals. In this study, the overall temporal-spatial variability of particulate Pb, Cu, Cr and Ni losses within an agricultural watershed was quantitatively evaluated by combining SWAT with sediment geochemistry. Results showed that the watershed particulate heavy metal loadings displayed strong variability in the simulation period 1981-2010, with an obvious increasing trend in recent years. The simulated annual average loadings were 20.21 g/ha, 21.75 g/ha, 47.35 g/ha and 21.27 g/ha for Pb, Cu, Cr and Ni, respectively. By comparison, these annual average values generally matched the estimated particulate heavy metal loadings at field scale. With spatial interpolation of field loadings, it was found that the diffuse heavy metal pollution mainly came from the sub-basins dominated with cultivated lands, accounting for over 70% of total watershed loadings. The watershed distribution of particulate heavy metal losses was very similar to that of soil loss but contrary to that of heavy metal concentrations in soil, highlighting the important role of sediment yield in controlling the diffuse heavy metal loadings. (C) 2014 Elsevier B.V. All rights reserved.