Heavy metals in European soils: A geostatistical analysis of the FOREGS geochemical database

This paper presents the results of modeling the distribution of eight critical heavy metals (arsenic, cadmium, chromium, copper, mercury, nickel, lead and zinc) in topsoils using 1588 georeferenced samples from the Forum of European Geological Surveys Geochemical database (26 European countries). The concentrations were mapped using regression-kriging (RK) and accuracy of predictions evaluated using the leave-one-out cross validation method. A large number of auxiliary raster maps (topographic indexes, land cover, geology, vegetation indexes, night lights images and earth quake magnitudes) were used to improve the predictions. These were first converted to 36 principal components and then used to explain spatial distribution of heavy metals. The study revealed that this database is suitable for geostatistical analyses: the predictors explained from 21% (Cr) to 35% (Pb) of variability; the residuals showed spatial autocorrelation. The Principal Component Analysis of the mapped heavy metals revealed that the administrative units (NUTS level3) with highest overall concentrations are: (1) Liege (Arrondissement) (BE), Attiki (GR), Darlington (UK), Coventry (UK), Sunderland (UK), Kozani (GR), Grevena (GR), Hartlepool & Stockton (UK), Huy (BE), Aachen (DE) (As, Cd, Hg and Pb) and (2) central Greece and Liguria region in Italy (Cr, Cu and Ni). The evaluation of the mapping accuracy showed that the RK models for As, Ni and Pb can be considered satisfactory (prediction accuracy 45-52% of total variance), marginally satisfactory for Cr, Cu, Hg and Zn (36-41%), while the model for Cd is unsatisfactorily accurate (30%). The critical elements limiting the mapping accuracy are: (a) the problem of sporadic high values (hot-spots); and (b) relatively coarse resolution of the input maps. Automation of the geostatistical mapping and use of auxiliary spatial layers opens a possibility to develop mapping systems that can automatically update outputs by including new field observations and higher quality auxiliary maps. This approach also demonstrates the benefits of organizing standardized joint European monitoring projects, in comparison to the merging of several national monitoring projects.

[1]  P. de Caritat,et al.  Distinguishing between natural and anthropogenic sources for elements in the environment: regional geochemical surveys versus enrichment factors. , 2005, The Science of the total environment.

[2]  D. J. Walker,et al.  Heavy metal bioavailability in a soil affected by mineral sulphides contamination following the mine spillage at Aznalcóllar (Spain) , 2003, Biodegradation.

[3]  Heavy Metals in Waste , 2002 .

[4]  G. Heuvelink,et al.  Optimization of sample patterns for universal kriging of environmental variables , 2007 .

[5]  J. Strobl,et al.  SAGA-analysis and modelling applications , 2006 .

[6]  Walter H. F. Smith,et al.  A global, self‐consistent, hierarchical, high‐resolution shoreline database , 1996 .

[7]  Andreas Papritz,et al.  Modelling the spatial distribution of copper in the soils around a metal smelter in northwestern Switzerland , 2005 .

[8]  Philippe Renard,et al.  Geostatistics for environmental applications. Proceedings of the Fifth European Conference on Geostatistics for Environmental Applications. , 2005 .

[9]  R. Geary,et al.  Testing for Normality , 2003 .

[10]  G. Koptsik,et al.  Heavy metals in soils near the nickel smelter: chemistry, spatial variation, and impacts on plant diversity. , 2003, Journal of environmental monitoring : JEM.

[11]  P. Römkens,et al.  Calculation of critical loads for cadmium, lead and mercury; background document to a mapping manual on critical loads of cadmium, lead and mercury , 2005 .

[12]  A. Roth,et al.  The shuttle radar topography mission—a new class of digital elevation models acquired by spaceborne radar , 2003 .

[13]  P. Haastrup,et al.  Chapter 16 – A Water Reference for Europe , 2006 .

[14]  A. Facchinelli,et al.  Multivariate statistical and GIS-based approach to identify heavy metal sources in soils. , 2001, Environmental pollution.

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

[16]  Palle Haastrup,et al.  Environmental Data Exchange Network for Inland Water , 2007 .

[17]  M. Bidovec,et al.  Geochemical Atlas of Europe, Part 1, Background Information, Methodology and Maps , 2005 .

[18]  Yonglong Lu,et al.  Landscape ecology of the Guanting Reservoir, Beijing, China: multivariate and geostatistical analyses of metals in soils. , 2007, Environmental pollution.

[19]  A.J.H. Visschedijk,et al.  Incremental cost and remaining emission in 2010 of Heavy Metals (HM) resulting from the implementation of the draft HM Protocol under the UN/ECE Convention on Long Range Transboundary Air Pollution , 1998 .

[20]  P. Sutton Modeling population density with night-time satellite imagery and GIS , 1997 .

[21]  P. Döll,et al.  Development and validation of a global database of lakes, reservoirs and wetlands , 2004 .

[22]  Donald L. Gautier,et al.  Map showing geology, oil and gas fields, and geologic provinces of Europe including Turkey , 1997 .

[23]  Dar A. Roberts,et al.  A Comparison of Nighttime Satellite Imagery and Population Density for the Continental United States , 1997 .

[24]  M. Kane Correlation of major eastern earthquake centers with mafic/ultramafic basement masses , 1977 .

[25]  C. Reimann,et al.  Geochemical background--concept and reality. , 2005, The Science of the total environment.

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

[27]  H. Kern,et al.  Earthquakes and strength in the laminated lower crust —Can they be explained by the “corset model”? , 2008 .

[28]  J. L. Parra,et al.  Very high resolution interpolated climate surfaces for global land areas , 2005 .

[29]  R. Swennen,et al.  Assessment of acid neutralizing capacity and potential mobilisation of trace metals from land-disposed dredged sediments. , 2004, The Science of the total environment.

[30]  R. Salminen Geochemical atlas of the Eastern Barents region , 2004 .

[31]  Gerard B. M. Heuvelink,et al.  About regression-kriging: From equations to case studies , 2007, Comput. Geosci..

[32]  B. Gawlik,et al.  Background values in European soils and sewage sludges , 2006 .

[33]  Edzer J. Pebesma The Role of External Variables and GIS Databases in Geostatistical Analysis , 2006, Trans. GIS.

[34]  C. Micó,et al.  Assessing heavy metal sources in agricultural soils of an European Mediterranean area by multivariate analysis. , 2006, Chemosphere.

[35]  T. Sterckeman,et al.  Of the necessity of knowledge of the natural pedo-geochemical background content in the evaluation of the contamination of soils by trace elements. , 2001, The Science of the total environment.

[36]  S. Pirc,et al.  Geochemical atlas of Europe. Part 2, Interpretation of geochemical maps, additional tables, figures, maps, and related publications , 2006 .

[37]  I. Vanha-Majamaa,et al.  Understorey vegetation along a heavy-metal pollution gradient in SW Finland. , 2001, Environmental pollution.

[38]  Chaosheng Zhang,et al.  Using multivariate analyses and GIS to identify pollutants and their spatial patterns in urban soils in Galway, Ireland. , 2006, Environmental pollution.

[39]  Edzer Pebesma,et al.  GSTAT: a program for geostatistical modelling, prediction and simulation , 1998 .

[40]  G. Heuvelink,et al.  A generic framework for spatial prediction of soil variables based on regression-kriging , 2004 .

[41]  Hettelingh Jp,et al.  Critical Loads and Dynamic Modelling results, CCE Progress Report 2004 , 2004 .

[42]  Tomislav Hengl,et al.  A practical guide to geostatistical mapping of environmental variables , 2007 .

[43]  J. Eriksson,et al.  Agricultural soils in Northern Europe: a geochemical atlas. , 2003 .

[44]  N. Skarpelis Lateritization processes of ultramafic rocks in Cretaceous times: The fossil weathering crusts of mainland Greece , 2006 .

[45]  I. Thornton,et al.  Application of factorial kriging analysis to the FOREGS European topsoil geochemistry database. , 2008, The Science of the total environment.

[46]  Yu-Pin Lin,et al.  Multivariate analysis of soil heavy metal pollution and landscape pattern in Changhua county in Taiwan , 2002 .

[47]  Tomislav Hengl,et al.  Finding the right pixel size , 2006, Comput. Geosci..

[48]  B. Henderson,et al.  Australia-wide predictions of soil properties using decision trees , 2005 .

[49]  C. W. Thornthwaite An Approach Toward a Rational Classification of Climate , 1948 .