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.

[1]  V.J.G. Houba,et al.  Wageningen evaluating programmes for analytical laboratories (WEPAL), organization and purpose , 1996 .

[2]  D. Tibljaš,et al.  Discovery of Blueschists in the Medvednica Mountain (Northern Croatia) and Their Significance for the Interpretation of the Geotectonic Evolution of the Area , 1998 .

[3]  Pierre Goovaerts,et al.  A geostatistical approach to estimate probability of occurrence of Tuber melanosporum in relation to some soil properties , 2000 .

[4]  A. G. Fabbri,et al.  Quantitative analysis of mineral and energy resources , 1987 .

[5]  P. Goovaerts.  Factorial kriging analysis: a useful tool for exploring the structure of multivariate spatial soil information , 1992 .

[6]  A. Castrignanò,et al.  Study of spatial relationships among some soil physico-chemical properties of a field in central Italy using multivariate geostatistics , 2000 .

[7]  B. J. Alloway,et al.  An inventory of heavy metals inputs to agricultural soils in England and Wales. , 2003, The Science of the total environment.

[8]  Stefano Bocchi,et al.  Application of factorial kriging for mapping soil variation at field scale. , 2000 .

[9]  K. Korkka-niemi,et al.  Influence of parent sediments on the concentration of heavy metals in urban and suburban soils in Turku, Finland , 2007 .

[10]  Hans Wackernagel,et al.  Multivariate Geostatistics: An Introduction with Applications , 1996 .

[11]  Tomislav Hengl,et al.  Representing soil pollution by heavy metals using continuous limitation scores , 2007, Comput. Geosci..

[12]  B. J. Alloway,et al.  Heavy metals in soils , 1990 .

[13]  Dominique King,et al.  Scale of spatial dependence between chemical properties of topsoil and subsoil over a geologically contrasted area (Massif central, France) , 2003 .

[14]  D. Sparks.  Environmental Soil Chemistry , 1995 .

[15]  C. Reimann,et al.  Chemical elements in the environment , 1998 .

[16]  M. Romić,et al.  Heavy metals distribution in agricultural topsoils in urban area , 2003 .

[17]  H. Wackernagel.  Cokriging versus kriging in regionalized multivariate data analysis , 1994 .

[18]  A. Bellanca,et al.  Influence of inheritance and pedogenesis on heavy metal distribution in soils of Sicily, Italy , 2000 .

[19]  M. Romić,et al.  Heavy Metals Accumulation in Topsoils from the Wine-growing Regions Part 1. Factors which Control Retention , 2004 .

[20]  K. H. Wedepohl.  The Composition of the Continental Crust , 1995 .

[21]  Pierre Goovaerts,et al.  Scale-dependent correlation between topsoil copper and cobalt concentrations in Scotland , 1994 .

[22]  R. Webster,et al.  Coregionalization of trace metals in the soil in the Swiss Jura , 1994 .

[23]  A. Castrignanò,et al.  Accounting for extensive topographic and pedologic secondary information to improve soil mapping. , 2009 .

[24]  R. Webster,et al.  Geostatistical analysis of soil contamination in the Swiss Jura. , 1994, Environmental pollution.

[25]  M. Goulard,et al.  Linear coregionalization model: Tools for estimation and choice of cross-variogram matrix , 1992 .

[26]  B. Nicoullaud,et al.  Exploring the Spatial Relationships Between Some Soil Properties and Wheat Yields in Two Soil Types , 2004, Precision Agriculture.

[27]  Zlatan Bajraktarević,et al.  Stratigraphy and Palaeogeography of Miocene Deposits from the Marginal Area of Žumberak Mt. and the Samoborsko Gorje Mts. (Northwestern Croatia) , 2010, Geologia Croatica.

[28]  Hans Wackernagel,et al.  Geostatistical Techniques for Interpreting Multivariate Spatial Information , 1988 .