The final publication is available at Springer via http://dx.doi.org/10.1007/s12665-013-2704-5 ; Characterizing the spatial patterns of variability is a fundamental aspect when investigating what could be the causes behind the spatial spreading of a set of variables. In this paper, a large multivariate dataset from the southeast of Belgium has been analyzed using factorial kriging. The purpose of the study is to explore and retrieve possible scales of spatial variability of heavy metals. This is achieved by decomposing the variance-covariance matrix of the multivariate sample into coregionalization matrices, which are, in turn, decomposed into transformation matrices, which serve to decompose each regionalized variable as a sum of independent factors. Then, factorial cokriging is used to produce maps of the factors explaining most of the variance, which can be compared with maps of the underlying lithology. For the dataset analyzes, this comparison identifies a few point scale concentrations that may reflect anthropogenic contamination, and it also identifies local and regional scale anomalies clearly correlated to the underlying geology and to known mineralizations. The results from this analysis could serve to guide the authorities in identifying those areas which need remediation. ; Benamghar, A.; Gómez-Hernández, JJ. (2014). Factorial kriging of a geochemical dataset for the heavy-metal spatial-pattern characterization The Wallonian Region. Environmental Earth Sciences. 71(7):3161-3170. doi:10.1007/s12665-013-2704-5 ; S ; 3161 ; 3170 ; 71 ; 7 ; Bartholomé P et al (1977) Métallogénie de la Belgique, des Pays-Bas et du Luxembourg, Rapport nr 1, Belgium, pp 38 ; Candeias C, Ferreira da Silva E, Salgueiro AR, Pereira HG, Reis AP, Patinha C, Matos JX, Avila PH (2011) The use of multivariate statistical analysis of geochemical data for assessing the spatial distribution of soil contamination by potentially toxic elements in the Aljustrel mining area (Iberian Pyrite Belt, Portugal) ; da Silva EF, Avila PF, ...