Assessing heavy metal contamination in soils of the Zagreb region (Northwest Croatia) using multivariate geostatistics

The assessment of soil contamination and location of pollution sources represent a crucial issue in soil remediation. Topsoil samples were collected in the Zagreb area (Northwest Croatia) and the total contents of trace and major elements were determined. A multivariate geostatistical analysis was used to estimate soil chemical composition variability. Factorial Kriging Analysis (FKA) was used to investigate the scale-dependent correlation structure of some variables by modelling co-regionalization of ten chemical variables, co-kriging specific factors and mapping them. The FKA provided two regionalized factors at different spatial scales of variability: the first factor at shorter range for Zn, Pb, Cd, Cu and Ni indicated different sources of anthropogenic contamination, whereas Ca (mainly loading on the longer range factor) was related to the lithology and parent material composition. The methodology used has proved to be a useful tool to separate geological and anthropogenic causes of variation in soil heavy metal content and to identify common pollution sources.