Frontpage
Articles
Statistics
Maps
Compare Maps
A contribution to IUGS/IAGC Global Geochemical Baselines

PDF version

DISTRIBUTION OF ELEMENTS IN SUBSOIL AND TOPSOIL

W. De Vos1, V. Gregorauskiene2, K. Marsina3, R. Salminen4, I. Salpeteur5, T. Tarvainen4, P.J. O’Connor6, A. Demetriades7, S. Pirc8, M. J. Batista9, M. Bidovec10 with contributions by A. Bel-lan11, M. Birke12, N. Breward13, B. De Vivo14, M. Duris15, J. Halamic16, P. Klein17, A. Lima14, J. Locutura11, J. Lis18, A. Mazreku19, R.T. Ottesen20, A. Pasieczna18, V. Petersell21, S. Reeder13, U. Siewers12, I. Slaninka3

1Geological Survey of Belgium, Brussels, Belgium; 2Geological Survey of Lithuania, Lithuania; 3Geological Survey of Slovak Republic, Slovak Republic; 4Geological Survey of Finland, Espoo, Finland; 5Geological Survey of France, France; 6Geological Survey of Ireland, Dublin, Ireland; 7Institute of Geology and Mineral Exploration, Athens, Greece; 8Geology Department, University of Ljubljana, Ljubljana, Slovenia; 9Geological Survey of Portugal, Portugal; 10Geological Survey of Slovenia, Ljubljana, Slovenia; 11Geological Survey of Spain, Madrid, Spain; 12Bundesanstalt für Geowissenschaften und Rohstoffe, Hannover, Germany; 13British Geological Survey, Keyworth, Nottingham, United Kingdom; 14Dipartimento di Scienze della Terra, Universita' di Napoli "Federico II", Naples, Italy; 15Czech Geological Survey, Prague, Czech Republic; 16Institute of Geology, Croatia, Zagreb, Croatia; 17Geological Survey of Austria, Wien, Austria; 18Polish Geological Institute, Warsaw, Poland; 19Centre of Civil Geology, Tirana Albania; 20Geological Survey of Norway, Trondheim, Norway; 21Geological Survey of Estonia, Tallinn, Estonia.

The FOREGS soil samples do not represent all soil types of Europe. According to the FOREGS Geochemical Mapping Field Manual (Salminen, Tarvainen et al. 1998), the sample should be taken from residual and sedentary (i.e., undisplaced) soil, which generally reflects the underlying lithology. Residual soil samples were collected from the small, second order drainage basins at suitable sites above their alluvial plain and base of slope, where alluvium and colluvium are respectively deposited. Residual soil may have been developed either on bedrock or on till, as is the case in glaciated terrains. Soil samples were mostly collected from forest soils. However, agriculture or pastureland sometimes constitutes the dominant land use type in the small catchments selected for sampling in southern Europe (Map 1). Photographs from each soil sampling site are available from the electronic version of the Atlas.

Map 1

The main purpose of the FOREGS geochemical survey is not to document the geochemistry of different soil types, but to discover the large pattern of geochemical signatures on the scale of the continent, and to investigate the different factors influencing this pattern, notably bedrock geology, climate and human influences. Detailed pedogenetic processes fall outside the scope of this investigation. From this point of view, soil geochemistry helps us to understand geology on the continental scale, and not the other way around.

The distribution of most elements in soil shows a pattern related to geology and/or mineralisation. A major geochemical dichotomy appears on the various maps between soil developed on crystalline and Palaeozoic basement rocks, and those on unfolded cover rocks that are naturally impoverished in most trace elements (except for areas of intense agricultural activity). Palaeocene and Eocene climates caused strong argillic and ferralitic weathering of some basement rocks, with a change in mineralogy, and a modification of vertical and lateral distribution patterns of most major elements.

On crystalline bedrock, there is also a systematic difference between the glacial till area encompassing Fennoscandia and Scotland, and the rest of Europe. Till is mostly quite fresh and chemically unweathered detritus, while in other areas the surficial deposits contain old and chemically strongly weathered material. Moreover, weathering of carbonates and silicates shows a different pattern. The anomaly pattern shown on the Na2O map (higher content over glacial till areas) can be found also for elements such as Ca, Sr, Ba, K and some others. Since these elements are also derived from other sources (carbonate in the case of Ca and Sr, shale in the case of K and Ba) high element values can be found in non-glacial southern Europe as well. The dilution of many element concentrations in the marginal area of the last glaciation is another geological process. The amount of weathered material in till increases towards the ice marginal area, where till contains a large amount of silica.

The European map of the weathering index (Map 2) shows that weathering in Fennoscandia is generally less advanced than elsewhere in Europe. The weathering index is defined as:

CIA =

100 x [Al2O3 / (Al2O3 + CaO* + Na2O + K2O)],

where CaO* = Ca in silicate fraction (McClennan et al. 1990).

Map 2

A high CIA value indicates a long weathering history of the soil. It is noted that the underlying bedrock also affects the CIA index value, e.g., Al-rich sedimentary rocks tend to have higher values, whereas Al-poor limestone always has low values (usually <40).

In the section on soil, for the sake of uniformity, and as a necessary simplification in descriptions, the following definitions were adopted with reference to the coloured maps and histograms in Part 1 of the Geochemical Atlas of Europe (Salminen et al. 2005)

Low values group the three lowest shades of blue in the colour scale, corresponding to the range from the minimum value up to the 25th percentile, defined as “very low” and “low background” concentrations in Part 1 (Tarvainen et al. 2005, p.97), and

High values group the three highest shades of red in the colour scale, corresponding to the range of values from the 75th percentile up to the maximum, defined as “high”, “very high” and “highly anomalous” concentrations in Part 1 (Tarvainen et al. 2005, p.97).

Correlation coefficients were calculated with Pearson’s product-moment linear correlation method (Table available in electronic format on website www.gtk.fi/publ/foregsatlas) after deletion of outliers and subsequent pairwise deletion of absent data. For a given element or determinand, outliers were defined here as values exceeding by a factor of 1.5 other nearby results, when all analytical results are ranked. They are generally visible on the histogram accompanying each map in Part 1 of the Geochemical Atlas. A maximum of four outliers were removed in this work for the calculation of linear correlation coefficients. A list of outliers is given separately for subsoil (Table 2) and topsoil (Table 3).

Table 2. Outliers of the subsoil data. Criterion: an outlier has a value exceeding by factor of 1.50 other nearby results when all analytical results are ranked. A maximum of four outliers were removed for the calculation of linear correlation coefficients.
Sample Country Element Value Unit Next value Factor
N36E03C1 Netherlands S 32 768 mg kg-1
N32E05C2 Switzerland S 21 069 mg kg-1 8 956 2.35
N37E05C2 Germany S 8 956 mg kg-1
N32E04C3 Switzerland S 8 657 mg kg-1 3 496 2.48
N26W03C1 Spain Hg 0.93 mg kg-1
N29W04C3 Spain Hg 0.90 mg kg-1 0.36 2.50
N36W02C1 UK TOC 48.52 %
N32E04C2 France TOC 41.3 % 23.1 1.79
N30E05C4 Italy Co 170 mg kg-1 77.7 2.19
N40E04C4 Norway Nb 133 mg kg-1
N31E06C4 Italy Nb 88.4 mg kg-1
N31E01C5 France Nb 70.4 mg kg-1 32.9 2.14
N30E05C4 Italy Ni 2 400 mg kg-1
N26E14C2 Greece Ni 1 790 mg kg-1
N27E12C1 Greece Ni 1 270 mg kg-1 513 2.48
N30E03C5 France Pb 938 mg kg-1
N26E14C4 Greece Pb 559 mg kg-1
N30E03C3 France Pb 531 mg kg-1 178 2.98
N28W05C1 Portugal As 593 mg kg-1 203 2.92
N26E12C4 Greece Cd 14.2 mg kg-1 6.25 2.27
N41E03C5 Norway Mo 17.2 mg kg-1
N31E06C3 Italy Mo 16.10 mg kg-1 9.62 1.67
N30E03C5 France Sb 30.3 mg kg-1
N26E14C4 Greece Sb 16.9 mg kg-1
N29W04C3 Spain Sb 16.6 mg kg-1 10.2 1.63
N32E06C3 Austria Te 1.63 mg kg-1 0.35 4.66
N30E03C5 France Tl 21.3 mg kg-1 6.47 3.29
N26E12C4 Greece P2O5 1.66 %
N35E07C2 Germany P2O5 1.56 % 0.91 1.71
N28E14C3 Greece Sn 106 mg kg-1
N28W05C4 Portugal Sn 103 mg kg-1
N28W04C4 Portugal Sn 86 mg kg-1 46 1.87
N32E05C2 Switzerland Sr 2 010 mg kg-1 1 167 1.72
N30E03C5 France Zn 3 062 mg kg-1
N30E03C3 France Zn 989 mg kg-1 624 1.58
Table 3. Outliers of the topsoil data. Criterion: an outlier has a value exceeding by factor of 1.5 other nearby results, when all analytical results are ranked. A maximum of four outliers were removed for the calculation of linear correlation coefficients.
Sample Country Element Value Unit Next value Factor
N28W01T4 Spain S 112 280 mg kg-1
N25E08T1 Italy S 6 518 mg kg-1
N36W02T1 UK S 5 052 mg kg-1
N32E04T3 Switzerland S 4 566 mg kg-1 2 657 1.72
N28E08T1 Italy Hg 1.35 mg kg-1
N33E11T2 Slovakia Hg 1.29 mg kg-1
N38W02T3 UK Hg 1.18 mg kg-1
N26W03T1 Spain Hg 1.04 mg kg-1 0.56 1.86
N30E05T4 Italy Co 249 mg kg-1 135 1.84
N37W03T4 UK Dy 44.9 mg kg-1
N27W04T1 Spain Dy 22.9 mg kg-1 13.7 1.67
N37W03T4 UK Er 26.0 mg kg-1
N27W04T1 Spain Er 14.2 mg kg-1 7.94 1.79
N37W03T4 UK Ho 9.16 mg kg-1
N27W04T1 Spain Ho 4.81 mg kg-1 2.73 1.76
N37W03T4 UK Lu 3.21 mg kg-1
N27W04T1 Spain Lu 2.03 mg kg-1 1.06 1.92
N37W03T4 UK Tm 4.03 mg kg-1
N27W04T1 Spain Tm 2.18 mg kg-1 1.12 1.95
N37W03T4 UK Yb 25.0 mg kg-1
N27W04T1 Spain Yb 13.3 mg kg-1 7.11 1.87
N37W03T4 UK Tb 7.01 mg kg-1 3.65 1.92
N37W03T4 UK Y 267 mg kg-1 149 1.79
N19W10T1 Spain Nb 134 mg kg-1
N40E04T4 Norway Nb 132 mg kg-1
N31E06T6 Italy Nb 110 mg kg-1 64.8 1.70
N30E05T4 Italy Ni 2 690 mg kg-1
N26E14T2 Greece Ni 2 070 mg kg-1
N27E12T1 Greece Ni 1 190 mg kg-1
N28E11T1 Albania Ni 1 160 mg kg-1 470 2.47
N26E14T4 Greece Pb 970 mg kg-1
N30E03T5 France Pb 904.00 mg kg-1 487 1.86
N37W03T4 UK U 53.20 mg kg-1
N30E02T2 France U 22.30 mg kg-1 12.6 1.77
N29E01T3 Spain V 537 mg kg-1 322 1.67
N26E12T4 Greece Cd 14.1 mg kg-1
N30E03T2 France Cd 7.11 mg kg-1
N30E04T4 France Cd 6.75 mg kg-1
N35E04T5 Germany Cd 5.38 mg kg-1 3.10 1.74
N29E01T3 Spain Cu 256 mg kg-1
N30E05T1 Italy Cu 215 mg kg-1 120 1.79
N29E01T3 Spain Mo 21.3 mg kg-1 11.5 1.85
N29W01T4 Spain Te 0.93 mg kg-1 0.35 2.66
N30E03T2 France Tl 24 mg kg-1 6.30 3.81
N19W10T1 Spain TiO2 5.45 %
N31E06T4 Italy TiO2 4.17 % 2.67 1.56
N26E12T2 Greece Cr 6 234 mg kg-1
N30E05T4 Italy Cr 3 272 mg kg-1
N30E06T5 Italy Cr 2 666 mg kg-1 1 6160 1.65
N41E05T3 Norway Sn 1060 mg kg-1 65 1.63
N28W01T4 Spain Sr 3 125 mg kg-1 883 3.54
N30E03T5 France Zn 2 904 mg kg-1
N30E04T4 France Zn 2 328 mg kg-1 991 2..35
N29W02T2 Spain Zn 991 mg kg-1
N30E03T3 France Zn 986 mg kg-1 398 2.48

Throughout the text the following notation is used for the correlation coefficients:

  • Very strong correlation: > 0.8;
  • Strong correlation: between 0.6 and 0.8;
  • Good correlation: between 0.4 and 0.6, and
  • Weak correlation: between 0.3 and 0.4.

Because of the large number of samples, even the so-called weak correlations are significant at the 0.01 confidence level.

The use of Pearson’s correlation coefficients rather than Spearman’s rank correlation, and the deletion of outliers, is linked to a factor analysis, performed on the dataset (Batista et al. 2006, Annex 5 in this volume).

When interpreting correlation coefficients of major elements, especially Si, Al and Ca, it should be kept in mind that we work in a closed system where the sum total of all elements (or oxides in the case of majors) add up to approximately 100%. This leads to some form of autocorrelation, as observed in the negative correlation between Si and Ca, and the negative correlation Si-Al. This phenomenon also contributes to the significance of the so-called weak correlations mentioned earlier.

Despite the critical objections that can be made against using linear correlation coefficients in this data set covering very dissimilar regions, the coefficients were found to express general endencies of element associations at the continental scale. These tendencies are similar for soil, stream sediment and floodplain sediment (together forming the solid sample media), and are, therefore, taken to be meaningful indicators of geochemical processes. For this reason, and despite the obvious shortcomings, correlations are mentioned throughout the text. The overall pattern will further be discussed in the factor analysis section by Batista et al. 2006 (Annex 5 in this volume).

The ratios topsoil/subsoil were calculated for all elements (Table 4), averaging all individual ratios between pairs of subsoil and topsoil samples, for the whole dataset including the outliers, provided both subsoil and topsoil were sampled at the same site. For calculation of the correlation coefficient between topsoil and subsoil, likewise, outliers were included (Table 4). The only purpose of this correlation coefficient is to evaluate how systematic an increasing (or decreasing) tendency is from subsoil to topsoil: a high coefficient means that enrichment (ratio >1) or leaching (ratio <1) in topsoil occurs in most of Europe, a low coefficient indicates that enrichment or leaching is more erratic or applies only to certain regions.

Table 4. Ratios of topsoil/subsoil for all published elements.
Element Ratio top/sub soil Correlation REE Toxic / heavy metal
TOC 2.585 0.588
Hg 1.660 0.526 heavy metal
Cd 1.477 0.802 heavy metal
Pb 1.364 0.863 heavy metal
P2O5 1.261 0.653
Sb 1.184 0.612
S 1.140 0.379
Sn 1.134 0.421 heavy metal
Mo 1.105 0.822
Bi 1.100 0.701
Cr 1.096 0.814 heavy metal
Zr 1.088 0.844
Hf 1.080 0.825
MnO 1.075 0.762
Zn 1.075 0.876 heavy metal
Ag 1.053 0.752
TiO2 1.021 0.886
As 1.017 0.479 toxic
SiO2 1.017 0.877
I 1.013 0.561
Nb 1.008 0.931
Ta 0.999 0.894
Te 0.989 0.407
W 0.987 0.313
In 0.982 0.688
Tl 0.982 0.375 heavy metal
Ba 0.976 0.914
Cu 0.967 0.666 heavy metal
Rb 0.965 0.894
V 0.953 0.860
U 0.952 0.666
Na2O 0.951 0.958
Lu 0.947 0.744 heavy REE (Y)
Ga 0.946 0.856
Yb 0.944 0.715 heavy REE (Y)
K2O 0.940 0.891
pH 0.936 0.840
Tm 0.935 0.682 heavy REE (Y)
Er 0.931 0.689 heavy REE (Y)
Al2O3 0.925 0.853
Th 0.925 0.924
Fe2O3 0.923 0.840
Ho 0.923 0.665 heavy REE (Y)
Dy 0.919 0.659 heavy REE (Y)
Co 0.918 0.869
Ni 0.918 0.972 heavy metal
Y 0.914 0.699 Y
Tb 0.910 0.681 heavy REE (Y)
Cs 0.909 0.810
La 0.905 0.832 light REE (Ce)
Be 0.904 0.794
Ce 0.899 0.815 light REE (Ce)
Gd 0.897 0.717 heavy REE (Y)
Pr 0.897 0.806 light REE (Ce)
Nd 0.896 0.799 light REE (Ce)
Sc 0.894 0.842
Sm 0.888 0.759 light REE (Ce)
Eu 0.884 0.829 light REE (Ce)
Sr 0.877 0.880
MgO 0.796 0.831
CaO 0.702 0.819

Acknowledgements

The Geochemical Atlas “Agricultural Soils in northern Europe” by Clemens Reimann et al. (2003) was often consulted for comparison.

PDF version

placeholder text