The present study analyses spatial variation and covariation for seven heavy metal concentrations in the agricultural topsoil of the Duero basin (N Spain) in an attempt to identify their sources. Factorial kriging is used in conjunction with a principal components analysis(PCA).