# Occurrence of Vanadium in Belgian and European Alluvial Soils

Authors: V. Cappuyns; E. Slabbinck
Journal: Applied and Environmental Soil Science
 (2012)

---

## Abstract

Vanadium (V) is a naturally occurring trace element, but total concentrations in soils and sediments are also dependent on the parent material and might be influenced by anthropogenic activities (e.g., steel industry). Despite the fact that threshold values for V in soils and/or sediments exist in various European countries, in Belgium, V is not taken into account when the environmental quality of soils and sediments has to be evaluated, despite the existence of several (diffuse) sources for V. In the first part of the study, the occurrence of V alluvial soils in Belgium was compared with V concentrations in alluvial soils (floodplain soils) across Europe. By analysis of both the Belgian and European data, the relationship between physicochemical soil characteristics and total V concentrations was quantified and some areas polluted with V were detected. A regression equation, in which V concentrations in alluvial soils were expressed as a function of major element composition, was proposed for the Belgian and European data. Additionally, single extractions with CaCl2 (0.01 mol L−1) and ammonium-EDTA (0.05 mol L−1) were used to estimate short- and long-term mobility of V in the alluvial soils.


---

## Body

## 1. Introduction

### 1.1. Vanadium in Soils and Sediments

Vanadium (V) naturally occurs as trace element in soils and sediments. The average concentration in soil is estimated to be around 150μg g−1 [1] and 90 μg g−1 [2]. The concentration of V in soils and sediments also depends on the parent material and the occurrence of V-containing ore minerals in the subsoil [3]. Vanadium concentrations in igneous rocks are higher than in acidic and in siliceous rocks, whereas the V content of metamorphic and sedimentary rocks is intermediate between the content of basic and acidic igneous rocks [3]. In soils, V is mainly associated with Fe(hydr)oxides, clay minerals, and organic matter and it can also occur as discrete mineral phases such as carnotite (K2(UO2)2(VO4)2·3H2O) and vanadinite (Pb5(VO4)3Cl) [2, 4, 5]. As an element of group VB and with atomic number 23 in the periodic table, it can acquire several oxidation states (+III, +IV, +V) [2, 6].Beside its natural occurrence in soils and sediments, several anthropogenic sources can cause an enrichment or even pollution of soils and sediments with V. The main anthropogenic sources of V are (1) the use of V as a catalyst in the metal-, cement-, electronics, and textile industry; (2) the production of certain metals such as iron and uranium, with V as a secondary product [2], vanadium concentrations up to 738 mg/kg [7] and 3505 mg/kg [1] have been reported in mining areas; (3) air emission of V from the combustion of fossil fuels [8]. The majority of the emissions of V into the air are taken up by the soil surface, especially when soils are rich in organic matter [9]. In Flanders (Belgium), the petrochemical and chemical industry are the main anthropogenic sources for V. Other minor contributors to V emissions are the food industry and the non-ferro industry [10].Most food crops and also other plants contain trace amounts of V [6] because small concentrations of V stimulate plant growth and act as a catalyst for nitrogen fixation [11]. When plants and crops are exposed to enhanced V concentrations in soil, they are usually smaller than the same plants grown on noncontaminated soils [1, 2, 12]. V concentrations in porewater above 3 mg L−1 result in important toxic effects for plants [3]. Kasai et al. [13] found that 0.1 μM vanadate had slight effects on seed germination, growth, photosynthesis, respiration, and ATPase activity of wheat and rye grass. The use of certain fertilizers ((NH4)2PO4, Na3PO4) increases the solubility, mobility, and uptake of V by plants [14]. The oxidation state of V is primarily controlled by the Eh (redox potential) and pH of a system. Under oxidizing and more alkaline conditions, vanadium is found in anionic form, as vanadate (HVO42− or H2VO4−, oxidation state +V), whereas vanadyl (VO2+, oxidation state +IV) has cationic properties and occurs under more acid conditions [3]. Complexes with fluoride, sulfate, or oxalate can increase the solubility of V under oxidizing conditions [15]. Acute poisoning of cattle by uptake of V-containing soil particles sticking to grass has been reported by McCrindle et al. [16].Vanadium toxicity for human beings can vary widely, depending on its oxidation state.  Clinical picture of poisoning (starting from 10 mg per day) shows different toxic effects of vanadium on the respiratory, circulatory and central nervous systems, the digestive organs, kidneys, and skin [17]. In the past, V has also been used as a therapeutic for the treatment of diabetes.


### 1.2. Alluvial Soils, Overbank Sediments, Floodplain Sediments

Darnley et al. [18] made a distinction between alluvial soils according to the catchment basin size: floodplain and overbank sediments are alluvial soils of large and small floodplains, respectively, but both consist of fine-grained (silty-clay, clayey-silt) sediments. Floodplain and overbank sediments are deposited during flood events in low energy environments [19]. Overbank sediments are considered a representative sampling medium for geochemical mapping [19]. In areas of historical mining and pollution, only one overbank sediment sample is not always representative for the geochemical characteristics of a large catchment basin [20]. De Vos et al. [4] evaluated the use of overbank sediments as a medium for geochemical mapping in industrialized areas. Samples of lower overbank sediment samples (i.e., samples taken deeper in the overbank sediment profile) reflect the natural, pristine situation, while upper overbank sediments are influenced by anthropogenic chemical contamination. From an environmental point of view, the use of overbank sediments is interesting because background concentrations are an essential reference point to evaluate the pollution status of soils and sediments.


### 1.3. Legal Framework with regard to V Contamination

In the environmental legislation of Flanders and Wallonia (which are the two main regions in Belgium, each with an own environmental legislation), threshold vanadium concentrations are defined for emissions in air, surface water, and groundwater. Although background and intervention values are provided in the Flemish and Walloon Soil Decree for As, Cd, Cr, Cu, Hg, Ni, Pb, Zn, and organic contaminants, no norm values are provided for V. Nevertheless, several European countries established norm values for V in soils (Table1).Table 1
Norm values for vanadium (total concentration in mg kg−1) in soil in European countries [21, 22].
CountryNegligibleWarningUnacceptableUnacceptableriscriscrisc (residential)risc (industrial)The Netherlands42—250—Czech Republic180340450—Slovenia120200500—Finland—100150250Sweden—200——Italy——90250Lithuania——150—Russia——150—The norm values in Table1 differ because of the use of different models, software, and criteria with regard to human toxicology and ecotoxicology [21]. Additionally, natural (background) concentrations of trace elements in soils and sediments are dependent on soil and sediment properties (clay content, organic carbon content, etc.) and on the geological substratum, which makes it difficult to compare these norm values [23]. Moreover, total concentrations of trace elements do not give a good indication of the mobility and toxicity of an element [2].The aim of this study was (1) to evaluate the V-content of alluvial soils (overbank sediments, floodplain sediments and dredged sediment-derived soils) along some Belgian and European rivers and to compare these values with norm values established in other European countries, (2) to establish regression equations that allow to explain the V content in alluvial soils based on major element composition and/or physico-chemical soil properties such as clay- and organic matter content, (3) to assess the actual and potential release (“mobile” and “mobilizable” fractions resp.) of V from these soils by means of standardized European extractions and leaching tests (single extractions with CaCl2 and ammonium-EDTA). It is not the purpose of this study to propose norm values for V in (alluvial) soils, but the regression equations obtained in this study can be useful to establish whether V concentrations that are measured in alluvial soils can be considered as “normal” or not.


## 1.1. Vanadium in Soils and Sediments

Vanadium (V) naturally occurs as trace element in soils and sediments. The average concentration in soil is estimated to be around 150μg g−1 [1] and 90 μg g−1 [2]. The concentration of V in soils and sediments also depends on the parent material and the occurrence of V-containing ore minerals in the subsoil [3]. Vanadium concentrations in igneous rocks are higher than in acidic and in siliceous rocks, whereas the V content of metamorphic and sedimentary rocks is intermediate between the content of basic and acidic igneous rocks [3]. In soils, V is mainly associated with Fe(hydr)oxides, clay minerals, and organic matter and it can also occur as discrete mineral phases such as carnotite (K2(UO2)2(VO4)2·3H2O) and vanadinite (Pb5(VO4)3Cl) [2, 4, 5]. As an element of group VB and with atomic number 23 in the periodic table, it can acquire several oxidation states (+III, +IV, +V) [2, 6].Beside its natural occurrence in soils and sediments, several anthropogenic sources can cause an enrichment or even pollution of soils and sediments with V. The main anthropogenic sources of V are (1) the use of V as a catalyst in the metal-, cement-, electronics, and textile industry; (2) the production of certain metals such as iron and uranium, with V as a secondary product [2], vanadium concentrations up to 738 mg/kg [7] and 3505 mg/kg [1] have been reported in mining areas; (3) air emission of V from the combustion of fossil fuels [8]. The majority of the emissions of V into the air are taken up by the soil surface, especially when soils are rich in organic matter [9]. In Flanders (Belgium), the petrochemical and chemical industry are the main anthropogenic sources for V. Other minor contributors to V emissions are the food industry and the non-ferro industry [10].Most food crops and also other plants contain trace amounts of V [6] because small concentrations of V stimulate plant growth and act as a catalyst for nitrogen fixation [11]. When plants and crops are exposed to enhanced V concentrations in soil, they are usually smaller than the same plants grown on noncontaminated soils [1, 2, 12]. V concentrations in porewater above 3 mg L−1 result in important toxic effects for plants [3]. Kasai et al. [13] found that 0.1 μM vanadate had slight effects on seed germination, growth, photosynthesis, respiration, and ATPase activity of wheat and rye grass. The use of certain fertilizers ((NH4)2PO4, Na3PO4) increases the solubility, mobility, and uptake of V by plants [14]. The oxidation state of V is primarily controlled by the Eh (redox potential) and pH of a system. Under oxidizing and more alkaline conditions, vanadium is found in anionic form, as vanadate (HVO42− or H2VO4−, oxidation state +V), whereas vanadyl (VO2+, oxidation state +IV) has cationic properties and occurs under more acid conditions [3]. Complexes with fluoride, sulfate, or oxalate can increase the solubility of V under oxidizing conditions [15]. Acute poisoning of cattle by uptake of V-containing soil particles sticking to grass has been reported by McCrindle et al. [16].Vanadium toxicity for human beings can vary widely, depending on its oxidation state.  Clinical picture of poisoning (starting from 10 mg per day) shows different toxic effects of vanadium on the respiratory, circulatory and central nervous systems, the digestive organs, kidneys, and skin [17]. In the past, V has also been used as a therapeutic for the treatment of diabetes.


## 1.2. Alluvial Soils, Overbank Sediments, Floodplain Sediments

Darnley et al. [18] made a distinction between alluvial soils according to the catchment basin size: floodplain and overbank sediments are alluvial soils of large and small floodplains, respectively, but both consist of fine-grained (silty-clay, clayey-silt) sediments. Floodplain and overbank sediments are deposited during flood events in low energy environments [19]. Overbank sediments are considered a representative sampling medium for geochemical mapping [19]. In areas of historical mining and pollution, only one overbank sediment sample is not always representative for the geochemical characteristics of a large catchment basin [20]. De Vos et al. [4] evaluated the use of overbank sediments as a medium for geochemical mapping in industrialized areas. Samples of lower overbank sediment samples (i.e., samples taken deeper in the overbank sediment profile) reflect the natural, pristine situation, while upper overbank sediments are influenced by anthropogenic chemical contamination. From an environmental point of view, the use of overbank sediments is interesting because background concentrations are an essential reference point to evaluate the pollution status of soils and sediments.


## 1.3. Legal Framework with regard to V Contamination

In the environmental legislation of Flanders and Wallonia (which are the two main regions in Belgium, each with an own environmental legislation), threshold vanadium concentrations are defined for emissions in air, surface water, and groundwater. Although background and intervention values are provided in the Flemish and Walloon Soil Decree for As, Cd, Cr, Cu, Hg, Ni, Pb, Zn, and organic contaminants, no norm values are provided for V. Nevertheless, several European countries established norm values for V in soils (Table1).Table 1
Norm values for vanadium (total concentration in mg kg−1) in soil in European countries [21, 22].
CountryNegligibleWarningUnacceptableUnacceptableriscriscrisc (residential)risc (industrial)The Netherlands42—250—Czech Republic180340450—Slovenia120200500—Finland—100150250Sweden—200——Italy——90250Lithuania——150—Russia——150—The norm values in Table1 differ because of the use of different models, software, and criteria with regard to human toxicology and ecotoxicology [21]. Additionally, natural (background) concentrations of trace elements in soils and sediments are dependent on soil and sediment properties (clay content, organic carbon content, etc.) and on the geological substratum, which makes it difficult to compare these norm values [23]. Moreover, total concentrations of trace elements do not give a good indication of the mobility and toxicity of an element [2].The aim of this study was (1) to evaluate the V-content of alluvial soils (overbank sediments, floodplain sediments and dredged sediment-derived soils) along some Belgian and European rivers and to compare these values with norm values established in other European countries, (2) to establish regression equations that allow to explain the V content in alluvial soils based on major element composition and/or physico-chemical soil properties such as clay- and organic matter content, (3) to assess the actual and potential release (“mobile” and “mobilizable” fractions resp.) of V from these soils by means of standardized European extractions and leaching tests (single extractions with CaCl2 and ammonium-EDTA). It is not the purpose of this study to propose norm values for V in (alluvial) soils, but the regression equations obtained in this study can be useful to establish whether V concentrations that are measured in alluvial soils can be considered as “normal” or not.


## 2. Methodology

### 2.1. Sampling Locations and Sampling

#### 2.1.1. Belgian Samples

Between November 2006 and October 2009, alluvial sediments, consisting of overbank sediment- or dredged sediment-derived soils, were sampled along 3 different rivers in Belgium (Figure1).Figure 1
Overview of the Belgian sampling locations.In northern Belgium, overbank sediment- and dredged sediment-derived soils were sampled in the catchment of the Scheldt river (catchment area in Belgium = 13336 km2), which is divided in 11 subcatchments, including the Leie and Demer catchments. The Leie river (north-western part of Belgium) is an important tributary of the river Scheldt and its total catchment area, consisting of Tertiary sands and clays, is about 3675 km2 (in Belgium). The Leie river is characterized by a regular flooding regime. In 1976, a part of the river was straightened and dredged sediments were stored along the river. Heavy metals in the dredged-sediment-derived soils originate from industrial activities and agriculture. Most of the dredged-sediment-derived soils are currently used as meadows or forests. The Grote Beek (Demer catchment, 1272 km2) river (northern Belgium) is underlain by glauconite-rich medium-grained sands (Diestian Formation). Organic- and iron-rich wetland soils have developed along this stream. The river follows a very meandering path and as a consequence several flooding zones occur along the river, which are flooded a few times a year during periods of heavy rainfall. The major source of pollution is a factory that is producing phosphate fertilizers by the refinement of phosphate ores. The majority of the alluvial soils along the Grote Beek consist of unused (i.e., not used for living, industry or agriculture) land, situated in a natural reserve.In eastern Belgium and in the southern part of The Netherlands, overbank sediments were sampled along the Geul river (Figure1). The Geul river, a tributary of the Meuse, flows in sandstones (Upper Devonian) and carbonates and shales (Lower and Middle Carbonian) [24]. The discharge of the Geul river depends largely on the amount of rainfall and is characterized by a flashy regime with a sequence of small floods in autumn and winter and a few major floods in summer [25]. From the Middle Ages until the beginning of the 20th century, extensive Zn-Pb mining and smelting was carried out in Plombières and La Calamine (eastern Belgium). Besides the important amount of waste that is stored in the large mine tailing piles, overbank sediments along the nearby Geul river are severely contaminated with Zn, Pb, and Cd [24]. The Geul drains a predominantly rural catchment with a relatively low population density.At each location, vertical profiles were dug and samples (1 kg of material per sample) were taken at depth every 5 to 20 cm, depending on visual differences in color, organic matter content and/or texture, yielding 3 to 5 samples per horizon. In total, 206 samples were taken. Along the Leie river, samples were taken in dredged sediment-derived soils. Along the Geul, profiles were sampled within overbank sediment-derived soils and along the Grote Beek both overbank sediment- and dredged sediment-derived soils were sampled.


#### 2.1.2. European Samples

“The FOREGS Geochemical Baseline Mapping Programme’s main aim is to provide high-quality, multipurpose environmental geochemical baseline data for Europe. The need for this type of data was justified by the first Working Group on Regional Geochemical Mapping immediately after the Chernobyl accident in 1986, when it was realized that a baseline for radioactive and other polluting elements could not be defined” [26–28]. After several attempts to compile existing regional geochemical databases in Europe, it was clear that the establishment of a harmonized European wide geochemical database was essential. The Forum of European Geological Surveys’ Directors (FOREGS) approved a Geochemical Baseline Mapping Programme in 1996, and in 1997 an agreement on the principles of field and analytical methodologies was obtained and the field methodology was tested, finally resulting in the “FOREGS Geochemical Field Manual,” that was published in 1998 [29]. In the present study, data from floodplain sediments samples in 26 countries that were involved in the FOREGS Geochemical Baseline Mapping Programme are used.From each sampling site, 2 kg of top floodplain sediment (sampling depth 0–25 cm) was collected. The sampling procedure for floodplain soils is described in the field manual [29].


### 2.2. Physicochemical Soil Characterization

#### 2.2.1. Belgian Alluvial Soils

Physicochemical Characterization
All analyses were performed on air-dried soil samples. For the determination of element concentrations and organic carbon content, part of the soil sample was disaggregated in a porcelain mortar and sieved (<1 mm). pH(H2O) was measured in a soil/water suspension (1/2.5). Organic carbon was determined according to the Walkley and Black method [30]. Grain size was determined by laser diffraction analyzis (Malvern Mastersizer S long bed). Total element concentrations (Al, As, Ba, Cd, Cu, Cr, Co, Ca,  Fe, K, Mg, Mn, Ni, Pb, V, and Zn) were determined in all the samples, but the discussion will mainly focus on V. One gram of sample was dissolved in 4 mL HClconc, 2 mL HNO3conc, and 2 mL HFconc in a Teflon beaker. The mixture was gently heated on a hot plate until half dry, subsequently reattacked with the same amount of the three acids and heated until completely dry. The residue was redissolved with 20 mL 2.5 mol L−1 HCl and filtered (Whatman 45). Finally, the solution was diluted to 50 mL with distilled water. This destruction method was chosen because we have a lot of experience with this method in our laboratory and it allows to compare the results with previous studies performed in our laboratory. Concentrations of major elements (Fe, Al, Ca, K, and Mg) were measured with Flame atomic absorption spectrometry (FAAS, Varian AA6), whereas an Induced coupled plasma mass spectrometry (ICP-MS, HP 4500 Series) was used for As, Ba, Cd, Cu, Cr, Co, Mn, Ni, Pb, Vm and Zn. Three certified reference materials, namely, GBW07411 Soil, (a Chinese soil), GBW07411 sediment (a Chinese stream sediment), and SRM 2710 Montana soil (a soil collected from the upper 10 cm of pasture along a creek and contaminated by overbank deposition of contaminated sediments) were analyzed for quality control. Total certified V concentrations in these materials are determined by XRF and/or neutron activation analysis, whereas a dissolution with 3 strong acids (HCl, HNO3, and HF) is used in the present study, followed by a measurement with ICP-MS. Despite the use of hydrofluoric acid, a small residue (usually less than 5% of the dry mass of the sample) is not dissolved. A previous investigation in our laboratory showed that this residue consists of quartz particles and sometimes also phyllosilicates. This possibly explains why measured V concentrations are slightly lower compared to certified values (Table 2). Nevertheless, the use of a dissolution with strong acids and ICP-MS determination allows a better comparison with data from the FOREGS project, where V has been determined by ICP-MS after aqua regia dissolution and by XRF. It is clear that the procedures for the determination of total element concentrations in soils and sediments should be well established, especially in large-scale studies, since this is essential to compare results of different studies and/or laboratories. Additionally, triplicate analysis was performed for some samples. Relative standard deviations on triplicate analysis (3 separate subsamples) were below 5% for all elements, except Ca and Al (below 10%).Table 2
Comparison between certified and measured (average of 3 replications) values for V in 3 certified reference materials.
Certified value (mg kg−1)Measured value (mg kg−1)SRM2710Montana soil76.669.8 ± 4.7GBW07311Soil47 ± 542.6 ± 3.3GBW07411Stream sediment88.5 ± 6.984.1 ± 6.8Single Extractions
Single extractions were performed on 65 samples of alluvial soils from the Leie and Grote Beek rivers. For the CaCl2 and ammonium-EDTA extraction, the protocol of the SMT (Standards Measurement and Testing) program [31] was followed. 20 mL of a 0.01 mol L−1 CaCl2 or 20 mL of a 0.05 mol L−1 ammonium-EDTA solution, respectively, was added to 2 g of air-dried sediment in a centrifuge tube. The suspension was shaken for 3 h or 1 h, respectively, in a reciprocal shaker, centrifuged (3500 rpm, 10 minutes), decanted, and filtered (0.45 μm). After measurement of pH, the CaCl2 extract was acidified with concentrated HNO3 to bring the pH to <2. The ammonium-EDTA extracts were not acidified prior to analysis to prevent precipitation of EDTA salts at very low pH. A reference material (CRM 483) certified for its ammonium-EDTA and CaCl2- extractable content of Cd, Zn, Cu, and Pb was also included. No certified values are provided for V, but the data obtained in this study are presented in Table 3 in order to allow comparison with other studies in the future. Standard deviation between triplicate extractions was less than 5%. The data of the ammonium-EDTA extracts are presented in mg/kg dry matter, since EDTA-extractable element concentrations are often compared to total element concentrations. The V concentrations in the CaCl2 extracts are expressed in μg L−1, to facilitate the comparison with, for example, porewater concentrations of V.Table 3
Values for V measured in the CaCl2 (in μg L−1) and ammonium-EDTA (in mg kg−1) extracts performed on certified reference material CRM 483 according to the BCR extraction protocol for single extractions with CaCl2 0.01 mol L−1 and ammonium-EDTA 0.05 mol L−1.
CaCl2 0.01 mol L−1ammionum-EDTA 0.05 mol L−1Total concentrationμgL-1mg kg−1mg kg−1CRM4835.7 ± 0.94.3 ± 0.164.5 ± 3.7


#### 2.2.2. Alluvial Soils from the FOREGS Database

Details about the analytical methods are provided in Sandström et al. [32]: total concentrations of MgO, P2O5, K2O, CaO, TiO2, V, Cr, MnO, Cs, Ba, La, and Ce were determined using energy dispersive polarized X-ray fluorescence spectrometry. In the floodplain soil samples, V was also determined by ICP-MS after aqua regia destruction. Because the data in both datasets are not always in the same order and data are sometimes missing, a comparison between data obtained with both methods was not performed. In the present study, mainly the data from the XRF analysis will be used.


### 2.3. Statistical Analysis

Statistical analysis was performed with the software package SPSS 16.0 for Windows.Descriptive statistics (average, median, minimum, maximum, standard deviation, and variance) were calculated for each variable. The 90th percentile of element concentrations in the upper soil layer of noncontaminated soils is often considered representative for the background value in soils [22, 33]. Because our datasets contained both contaminated and noncontaminated samples, no attempt was made to propose background values. By comparing average and median values, soils with enrichment in V can be detected [33].The normal distribution of the variables was checked by means of the Kolmogorof-Smirnov (K-S) test and correlations between variables were tested by calculating two-tailed Pearson correlation coefficients for the log transformed values. Analysis of variance (ANOVA) was applied to compare averages of the variables between river basins. Multiple linear regression according to the stepwise method was performed to deduce possible causal relationships between the variables. Attention was mainly paid to the possibility of predicting trace element concentrations in sediments based on major element composition, pH, clay, and organic matter content. Different assumptions of the linear regression (normality of the de residues, autocorrelation, quasi-multicollinearity (QMC), and heteroscedasticity) were tested.For the statistical analysis, the recommendations of Webster [34, 35] (reporting mean values with standard errors, performing linear regressions, etc.) were taken into account.


## 2.1. Sampling Locations and Sampling

### 2.1.1. Belgian Samples

Between November 2006 and October 2009, alluvial sediments, consisting of overbank sediment- or dredged sediment-derived soils, were sampled along 3 different rivers in Belgium (Figure1).Figure 1
Overview of the Belgian sampling locations.In northern Belgium, overbank sediment- and dredged sediment-derived soils were sampled in the catchment of the Scheldt river (catchment area in Belgium = 13336 km2), which is divided in 11 subcatchments, including the Leie and Demer catchments. The Leie river (north-western part of Belgium) is an important tributary of the river Scheldt and its total catchment area, consisting of Tertiary sands and clays, is about 3675 km2 (in Belgium). The Leie river is characterized by a regular flooding regime. In 1976, a part of the river was straightened and dredged sediments were stored along the river. Heavy metals in the dredged-sediment-derived soils originate from industrial activities and agriculture. Most of the dredged-sediment-derived soils are currently used as meadows or forests. The Grote Beek (Demer catchment, 1272 km2) river (northern Belgium) is underlain by glauconite-rich medium-grained sands (Diestian Formation). Organic- and iron-rich wetland soils have developed along this stream. The river follows a very meandering path and as a consequence several flooding zones occur along the river, which are flooded a few times a year during periods of heavy rainfall. The major source of pollution is a factory that is producing phosphate fertilizers by the refinement of phosphate ores. The majority of the alluvial soils along the Grote Beek consist of unused (i.e., not used for living, industry or agriculture) land, situated in a natural reserve.In eastern Belgium and in the southern part of The Netherlands, overbank sediments were sampled along the Geul river (Figure1). The Geul river, a tributary of the Meuse, flows in sandstones (Upper Devonian) and carbonates and shales (Lower and Middle Carbonian) [24]. The discharge of the Geul river depends largely on the amount of rainfall and is characterized by a flashy regime with a sequence of small floods in autumn and winter and a few major floods in summer [25]. From the Middle Ages until the beginning of the 20th century, extensive Zn-Pb mining and smelting was carried out in Plombières and La Calamine (eastern Belgium). Besides the important amount of waste that is stored in the large mine tailing piles, overbank sediments along the nearby Geul river are severely contaminated with Zn, Pb, and Cd [24]. The Geul drains a predominantly rural catchment with a relatively low population density.At each location, vertical profiles were dug and samples (1 kg of material per sample) were taken at depth every 5 to 20 cm, depending on visual differences in color, organic matter content and/or texture, yielding 3 to 5 samples per horizon. In total, 206 samples were taken. Along the Leie river, samples were taken in dredged sediment-derived soils. Along the Geul, profiles were sampled within overbank sediment-derived soils and along the Grote Beek both overbank sediment- and dredged sediment-derived soils were sampled.


### 2.1.2. European Samples

“The FOREGS Geochemical Baseline Mapping Programme’s main aim is to provide high-quality, multipurpose environmental geochemical baseline data for Europe. The need for this type of data was justified by the first Working Group on Regional Geochemical Mapping immediately after the Chernobyl accident in 1986, when it was realized that a baseline for radioactive and other polluting elements could not be defined” [26–28]. After several attempts to compile existing regional geochemical databases in Europe, it was clear that the establishment of a harmonized European wide geochemical database was essential. The Forum of European Geological Surveys’ Directors (FOREGS) approved a Geochemical Baseline Mapping Programme in 1996, and in 1997 an agreement on the principles of field and analytical methodologies was obtained and the field methodology was tested, finally resulting in the “FOREGS Geochemical Field Manual,” that was published in 1998 [29]. In the present study, data from floodplain sediments samples in 26 countries that were involved in the FOREGS Geochemical Baseline Mapping Programme are used.From each sampling site, 2 kg of top floodplain sediment (sampling depth 0–25 cm) was collected. The sampling procedure for floodplain soils is described in the field manual [29].


## 2.1.1. Belgian Samples

Between November 2006 and October 2009, alluvial sediments, consisting of overbank sediment- or dredged sediment-derived soils, were sampled along 3 different rivers in Belgium (Figure1).Figure 1
Overview of the Belgian sampling locations.In northern Belgium, overbank sediment- and dredged sediment-derived soils were sampled in the catchment of the Scheldt river (catchment area in Belgium = 13336 km2), which is divided in 11 subcatchments, including the Leie and Demer catchments. The Leie river (north-western part of Belgium) is an important tributary of the river Scheldt and its total catchment area, consisting of Tertiary sands and clays, is about 3675 km2 (in Belgium). The Leie river is characterized by a regular flooding regime. In 1976, a part of the river was straightened and dredged sediments were stored along the river. Heavy metals in the dredged-sediment-derived soils originate from industrial activities and agriculture. Most of the dredged-sediment-derived soils are currently used as meadows or forests. The Grote Beek (Demer catchment, 1272 km2) river (northern Belgium) is underlain by glauconite-rich medium-grained sands (Diestian Formation). Organic- and iron-rich wetland soils have developed along this stream. The river follows a very meandering path and as a consequence several flooding zones occur along the river, which are flooded a few times a year during periods of heavy rainfall. The major source of pollution is a factory that is producing phosphate fertilizers by the refinement of phosphate ores. The majority of the alluvial soils along the Grote Beek consist of unused (i.e., not used for living, industry or agriculture) land, situated in a natural reserve.In eastern Belgium and in the southern part of The Netherlands, overbank sediments were sampled along the Geul river (Figure1). The Geul river, a tributary of the Meuse, flows in sandstones (Upper Devonian) and carbonates and shales (Lower and Middle Carbonian) [24]. The discharge of the Geul river depends largely on the amount of rainfall and is characterized by a flashy regime with a sequence of small floods in autumn and winter and a few major floods in summer [25]. From the Middle Ages until the beginning of the 20th century, extensive Zn-Pb mining and smelting was carried out in Plombières and La Calamine (eastern Belgium). Besides the important amount of waste that is stored in the large mine tailing piles, overbank sediments along the nearby Geul river are severely contaminated with Zn, Pb, and Cd [24]. The Geul drains a predominantly rural catchment with a relatively low population density.At each location, vertical profiles were dug and samples (1 kg of material per sample) were taken at depth every 5 to 20 cm, depending on visual differences in color, organic matter content and/or texture, yielding 3 to 5 samples per horizon. In total, 206 samples were taken. Along the Leie river, samples were taken in dredged sediment-derived soils. Along the Geul, profiles were sampled within overbank sediment-derived soils and along the Grote Beek both overbank sediment- and dredged sediment-derived soils were sampled.


## 2.1.2. European Samples

“The FOREGS Geochemical Baseline Mapping Programme’s main aim is to provide high-quality, multipurpose environmental geochemical baseline data for Europe. The need for this type of data was justified by the first Working Group on Regional Geochemical Mapping immediately after the Chernobyl accident in 1986, when it was realized that a baseline for radioactive and other polluting elements could not be defined” [26–28]. After several attempts to compile existing regional geochemical databases in Europe, it was clear that the establishment of a harmonized European wide geochemical database was essential. The Forum of European Geological Surveys’ Directors (FOREGS) approved a Geochemical Baseline Mapping Programme in 1996, and in 1997 an agreement on the principles of field and analytical methodologies was obtained and the field methodology was tested, finally resulting in the “FOREGS Geochemical Field Manual,” that was published in 1998 [29]. In the present study, data from floodplain sediments samples in 26 countries that were involved in the FOREGS Geochemical Baseline Mapping Programme are used.From each sampling site, 2 kg of top floodplain sediment (sampling depth 0–25 cm) was collected. The sampling procedure for floodplain soils is described in the field manual [29].


## 2.2. Physicochemical Soil Characterization

### 2.2.1. Belgian Alluvial Soils

Physicochemical Characterization
All analyses were performed on air-dried soil samples. For the determination of element concentrations and organic carbon content, part of the soil sample was disaggregated in a porcelain mortar and sieved (<1 mm). pH(H2O) was measured in a soil/water suspension (1/2.5). Organic carbon was determined according to the Walkley and Black method [30]. Grain size was determined by laser diffraction analyzis (Malvern Mastersizer S long bed). Total element concentrations (Al, As, Ba, Cd, Cu, Cr, Co, Ca,  Fe, K, Mg, Mn, Ni, Pb, V, and Zn) were determined in all the samples, but the discussion will mainly focus on V. One gram of sample was dissolved in 4 mL HClconc, 2 mL HNO3conc, and 2 mL HFconc in a Teflon beaker. The mixture was gently heated on a hot plate until half dry, subsequently reattacked with the same amount of the three acids and heated until completely dry. The residue was redissolved with 20 mL 2.5 mol L−1 HCl and filtered (Whatman 45). Finally, the solution was diluted to 50 mL with distilled water. This destruction method was chosen because we have a lot of experience with this method in our laboratory and it allows to compare the results with previous studies performed in our laboratory. Concentrations of major elements (Fe, Al, Ca, K, and Mg) were measured with Flame atomic absorption spectrometry (FAAS, Varian AA6), whereas an Induced coupled plasma mass spectrometry (ICP-MS, HP 4500 Series) was used for As, Ba, Cd, Cu, Cr, Co, Mn, Ni, Pb, Vm and Zn. Three certified reference materials, namely, GBW07411 Soil, (a Chinese soil), GBW07411 sediment (a Chinese stream sediment), and SRM 2710 Montana soil (a soil collected from the upper 10 cm of pasture along a creek and contaminated by overbank deposition of contaminated sediments) were analyzed for quality control. Total certified V concentrations in these materials are determined by XRF and/or neutron activation analysis, whereas a dissolution with 3 strong acids (HCl, HNO3, and HF) is used in the present study, followed by a measurement with ICP-MS. Despite the use of hydrofluoric acid, a small residue (usually less than 5% of the dry mass of the sample) is not dissolved. A previous investigation in our laboratory showed that this residue consists of quartz particles and sometimes also phyllosilicates. This possibly explains why measured V concentrations are slightly lower compared to certified values (Table 2). Nevertheless, the use of a dissolution with strong acids and ICP-MS determination allows a better comparison with data from the FOREGS project, where V has been determined by ICP-MS after aqua regia dissolution and by XRF. It is clear that the procedures for the determination of total element concentrations in soils and sediments should be well established, especially in large-scale studies, since this is essential to compare results of different studies and/or laboratories. Additionally, triplicate analysis was performed for some samples. Relative standard deviations on triplicate analysis (3 separate subsamples) were below 5% for all elements, except Ca and Al (below 10%).Table 2
Comparison between certified and measured (average of 3 replications) values for V in 3 certified reference materials.
Certified value (mg kg−1)Measured value (mg kg−1)SRM2710Montana soil76.669.8 ± 4.7GBW07311Soil47 ± 542.6 ± 3.3GBW07411Stream sediment88.5 ± 6.984.1 ± 6.8Single Extractions
Single extractions were performed on 65 samples of alluvial soils from the Leie and Grote Beek rivers. For the CaCl2 and ammonium-EDTA extraction, the protocol of the SMT (Standards Measurement and Testing) program [31] was followed. 20 mL of a 0.01 mol L−1 CaCl2 or 20 mL of a 0.05 mol L−1 ammonium-EDTA solution, respectively, was added to 2 g of air-dried sediment in a centrifuge tube. The suspension was shaken for 3 h or 1 h, respectively, in a reciprocal shaker, centrifuged (3500 rpm, 10 minutes), decanted, and filtered (0.45 μm). After measurement of pH, the CaCl2 extract was acidified with concentrated HNO3 to bring the pH to <2. The ammonium-EDTA extracts were not acidified prior to analysis to prevent precipitation of EDTA salts at very low pH. A reference material (CRM 483) certified for its ammonium-EDTA and CaCl2- extractable content of Cd, Zn, Cu, and Pb was also included. No certified values are provided for V, but the data obtained in this study are presented in Table 3 in order to allow comparison with other studies in the future. Standard deviation between triplicate extractions was less than 5%. The data of the ammonium-EDTA extracts are presented in mg/kg dry matter, since EDTA-extractable element concentrations are often compared to total element concentrations. The V concentrations in the CaCl2 extracts are expressed in μg L−1, to facilitate the comparison with, for example, porewater concentrations of V.Table 3
Values for V measured in the CaCl2 (in μg L−1) and ammonium-EDTA (in mg kg−1) extracts performed on certified reference material CRM 483 according to the BCR extraction protocol for single extractions with CaCl2 0.01 mol L−1 and ammonium-EDTA 0.05 mol L−1.
CaCl2 0.01 mol L−1ammionum-EDTA 0.05 mol L−1Total concentrationμgL-1mg kg−1mg kg−1CRM4835.7 ± 0.94.3 ± 0.164.5 ± 3.7


### 2.2.2. Alluvial Soils from the FOREGS Database

Details about the analytical methods are provided in Sandström et al. [32]: total concentrations of MgO, P2O5, K2O, CaO, TiO2, V, Cr, MnO, Cs, Ba, La, and Ce were determined using energy dispersive polarized X-ray fluorescence spectrometry. In the floodplain soil samples, V was also determined by ICP-MS after aqua regia destruction. Because the data in both datasets are not always in the same order and data are sometimes missing, a comparison between data obtained with both methods was not performed. In the present study, mainly the data from the XRF analysis will be used.


## 2.2.1. Belgian Alluvial Soils

Physicochemical Characterization
All analyses were performed on air-dried soil samples. For the determination of element concentrations and organic carbon content, part of the soil sample was disaggregated in a porcelain mortar and sieved (<1 mm). pH(H2O) was measured in a soil/water suspension (1/2.5). Organic carbon was determined according to the Walkley and Black method [30]. Grain size was determined by laser diffraction analyzis (Malvern Mastersizer S long bed). Total element concentrations (Al, As, Ba, Cd, Cu, Cr, Co, Ca,  Fe, K, Mg, Mn, Ni, Pb, V, and Zn) were determined in all the samples, but the discussion will mainly focus on V. One gram of sample was dissolved in 4 mL HClconc, 2 mL HNO3conc, and 2 mL HFconc in a Teflon beaker. The mixture was gently heated on a hot plate until half dry, subsequently reattacked with the same amount of the three acids and heated until completely dry. The residue was redissolved with 20 mL 2.5 mol L−1 HCl and filtered (Whatman 45). Finally, the solution was diluted to 50 mL with distilled water. This destruction method was chosen because we have a lot of experience with this method in our laboratory and it allows to compare the results with previous studies performed in our laboratory. Concentrations of major elements (Fe, Al, Ca, K, and Mg) were measured with Flame atomic absorption spectrometry (FAAS, Varian AA6), whereas an Induced coupled plasma mass spectrometry (ICP-MS, HP 4500 Series) was used for As, Ba, Cd, Cu, Cr, Co, Mn, Ni, Pb, Vm and Zn. Three certified reference materials, namely, GBW07411 Soil, (a Chinese soil), GBW07411 sediment (a Chinese stream sediment), and SRM 2710 Montana soil (a soil collected from the upper 10 cm of pasture along a creek and contaminated by overbank deposition of contaminated sediments) were analyzed for quality control. Total certified V concentrations in these materials are determined by XRF and/or neutron activation analysis, whereas a dissolution with 3 strong acids (HCl, HNO3, and HF) is used in the present study, followed by a measurement with ICP-MS. Despite the use of hydrofluoric acid, a small residue (usually less than 5% of the dry mass of the sample) is not dissolved. A previous investigation in our laboratory showed that this residue consists of quartz particles and sometimes also phyllosilicates. This possibly explains why measured V concentrations are slightly lower compared to certified values (Table 2). Nevertheless, the use of a dissolution with strong acids and ICP-MS determination allows a better comparison with data from the FOREGS project, where V has been determined by ICP-MS after aqua regia dissolution and by XRF. It is clear that the procedures for the determination of total element concentrations in soils and sediments should be well established, especially in large-scale studies, since this is essential to compare results of different studies and/or laboratories. Additionally, triplicate analysis was performed for some samples. Relative standard deviations on triplicate analysis (3 separate subsamples) were below 5% for all elements, except Ca and Al (below 10%).Table 2
Comparison between certified and measured (average of 3 replications) values for V in 3 certified reference materials.
Certified value (mg kg−1)Measured value (mg kg−1)SRM2710Montana soil76.669.8 ± 4.7GBW07311Soil47 ± 542.6 ± 3.3GBW07411Stream sediment88.5 ± 6.984.1 ± 6.8Single Extractions
Single extractions were performed on 65 samples of alluvial soils from the Leie and Grote Beek rivers. For the CaCl2 and ammonium-EDTA extraction, the protocol of the SMT (Standards Measurement and Testing) program [31] was followed. 20 mL of a 0.01 mol L−1 CaCl2 or 20 mL of a 0.05 mol L−1 ammonium-EDTA solution, respectively, was added to 2 g of air-dried sediment in a centrifuge tube. The suspension was shaken for 3 h or 1 h, respectively, in a reciprocal shaker, centrifuged (3500 rpm, 10 minutes), decanted, and filtered (0.45 μm). After measurement of pH, the CaCl2 extract was acidified with concentrated HNO3 to bring the pH to <2. The ammonium-EDTA extracts were not acidified prior to analysis to prevent precipitation of EDTA salts at very low pH. A reference material (CRM 483) certified for its ammonium-EDTA and CaCl2- extractable content of Cd, Zn, Cu, and Pb was also included. No certified values are provided for V, but the data obtained in this study are presented in Table 3 in order to allow comparison with other studies in the future. Standard deviation between triplicate extractions was less than 5%. The data of the ammonium-EDTA extracts are presented in mg/kg dry matter, since EDTA-extractable element concentrations are often compared to total element concentrations. The V concentrations in the CaCl2 extracts are expressed in μg L−1, to facilitate the comparison with, for example, porewater concentrations of V.Table 3
Values for V measured in the CaCl2 (in μg L−1) and ammonium-EDTA (in mg kg−1) extracts performed on certified reference material CRM 483 according to the BCR extraction protocol for single extractions with CaCl2 0.01 mol L−1 and ammonium-EDTA 0.05 mol L−1.
CaCl2 0.01 mol L−1ammionum-EDTA 0.05 mol L−1Total concentrationμgL-1mg kg−1mg kg−1CRM4835.7 ± 0.94.3 ± 0.164.5 ± 3.7


## 2.2.2. Alluvial Soils from the FOREGS Database

Details about the analytical methods are provided in Sandström et al. [32]: total concentrations of MgO, P2O5, K2O, CaO, TiO2, V, Cr, MnO, Cs, Ba, La, and Ce were determined using energy dispersive polarized X-ray fluorescence spectrometry. In the floodplain soil samples, V was also determined by ICP-MS after aqua regia destruction. Because the data in both datasets are not always in the same order and data are sometimes missing, a comparison between data obtained with both methods was not performed. In the present study, mainly the data from the XRF analysis will be used.


## 2.3. Statistical Analysis

Statistical analysis was performed with the software package SPSS 16.0 for Windows.Descriptive statistics (average, median, minimum, maximum, standard deviation, and variance) were calculated for each variable. The 90th percentile of element concentrations in the upper soil layer of noncontaminated soils is often considered representative for the background value in soils [22, 33]. Because our datasets contained both contaminated and noncontaminated samples, no attempt was made to propose background values. By comparing average and median values, soils with enrichment in V can be detected [33].The normal distribution of the variables was checked by means of the Kolmogorof-Smirnov (K-S) test and correlations between variables were tested by calculating two-tailed Pearson correlation coefficients for the log transformed values. Analysis of variance (ANOVA) was applied to compare averages of the variables between river basins. Multiple linear regression according to the stepwise method was performed to deduce possible causal relationships between the variables. Attention was mainly paid to the possibility of predicting trace element concentrations in sediments based on major element composition, pH, clay, and organic matter content. Different assumptions of the linear regression (normality of the de residues, autocorrelation, quasi-multicollinearity (QMC), and heteroscedasticity) were tested.For the statistical analysis, the recommendations of Webster [34, 35] (reporting mean values with standard errors, performing linear regressions, etc.) were taken into account.


## 3. Results and Discussion

### 3.1. Total Concentrations of V in Floodplain Soils

Average total V concentration in the 3 Belgian catchments and in the European FOREGS database were in the range 45–87 mg kg−1. The high average V concentrations in the alluvial soils of the Grote Beek catchment is mainly due to 7 outliers with V concentrations between 189 and 301 mg kg−1. Alluvial soils of the Geul river are characterized by significant lower V concentrations compared to the Leie and Grote Beek rivers. V concentrations in floodplain soils are rarely reported in literature. Peh and Miko [36] measured average V concentrations of 70 mg kg−1 in Croatia, which is in line with the values reported here. Ivanov and Kashin [37] recorded (AAS-measurement after aqua regia destruction) average V concentrations of 110 mg kg−1 in alluvial deposits in Transbaikalia (Russia), and Wigginton and Price [38] found (ICP-OES determination after aqua regia destruction) concentrations between 8 and 17 mg kg−1 in floodplain soils of Bayou Creek (UK). On the 4th October 2010, red mud suspension were released from the Ajkai Timfoldgyar Zrt alumina plant. Among other trace elements, mud samples at the source contained important V concentrations (>1000 mg kg−1), and floodplain sediments downstream were significantly affected by the red mud [39].Total V concentrations in floodplain soils from the FOREGS database are represented on Figure2. Since the data presented on Figure 2 are characterized by a right kurtosis and thus not normal distribution, concentrations represented on the map are strongly influenced by locally enhanced V concentrations [40]. Data for the Belgian catchments are not presented on a map, since the data originate from 3 distinct locations, only covering a limited area in Belgium.Figure 2
Total V concentrations in alluvial soils samples in 26 European countries [27].Data observations which lied more than 1.5 interquartile range lower than the first quartile or 1.5 interquartile range higher than the third quartile were considered as an outlier. In the Belgian data for V, 7 outliers were observed, between 189–301 mg kg−1, all coming from samples from the alluvial soil samples of the Grote Beek. De Grote Beek is also characterized by the highest average V concentrations and an important difference between average and median values for V, which indicates an anthropogenic source for V [33]. In this area, phosphate ores are processed by a plant that is emitting its waste water directly into the river. V is known as a secondary product of phosphate production [8]. In accordance to this, the correlation coefficient between V and P in alluvial soils is also very high in the samples from the Grote Beek river.When the V concentrations in these samples are compared with the norm values in Table1, 11 samples of the Grote Beek alluvial plain exceed the value for “unacceptable risk” in residential areas established in Finland, Lithuania, and Russia (150 mg kg−1), whereas 3 samples exceed the Dutch “negligible risk” value of 250 mg kg−1. The samples from the catchment of the Leie river are contaminated with heavy metals and often characterized by an elevated clay- and organic matter content. Only one sample from alluvial plain of the Leie river is characterized by a V content above 150 mg kg−1 (total V concentration is 161 mg kg−1).The samples from the Geul-catchment, finally, are characterized by relatively low total V concentrations (average 45 mg kg−1, Table 4). The average V concentration in the alluvial soils of the Geul river is below the 90th percentile of the European dataset, which could be considered as the background value [33], as well as below the average V concentration in European alluvial soils. The maximal V concentrations is 57 mg kg−1 (Table 4) and thus below all threshold values established in Table 1, except the value for “negligible risk” from The Netherlands (42 mg kg−1). The average V concentration in the soils along the Geul is very close to the median concentration, which indicates that V is not of anthropogenic origin [33]. The subsoil in the Geul catchment consists of sand and shale and other studies also showed that total V concentrations in soils developed on shale are lower compared to more clay-rich substrates [41].Table 4
Summary statistics for total V concentrations (in mg kg−1), based on Belgian (this study) and European [28] data. Standard deviation is reported within brackets. n: number of samples.
datasetareameanmedianminmaxnBelgian samples3 Belgian floodplains69 (43)5318301206floodplain of the Geul45 (6)46305775floodplain of the Grote Beek87 (55)662830194floodplain of the Leie71 (28)691816037FOREGSEurope (XRF)59 (35)561266749Europe (aqua regia)56747In the European FOREGS dataset, “the spatial distribution of V is clearly related to the bedrock geology and mineralization, especially mafic and ultramafic lithology, and also clay-rich soil with high Al2O3 content” [15]. The highest V concentration is found in the mineralized Oslo rift (266 mg kg−1). In England and Germany, elevated V concentrations (224 mg kg−1) are respectively due to industrial pollution and coal and oil combustion [15].


### 3.2. Relation of V with Other Elements

#### 3.2.1. Correlations

Because the majority of the data (trace element and major element concentrations) showed a right kurtosis and were significant on the K-S test, the data were log transformed. Concerning the FOREGS data, the correlation coefficients in Table5 were calculated with the data from the XRF determination, since more elements that are relevant for this paper were measured with this method compared with the aqua regia method. Unfortunately, it was not possible to merge the results of the aqua regia and XRF determination into one table, because the data were not in the same order and there were several missing samples for each of the methods. This inconvenience of the FOREGS database has also been mentioned by [40].Table 5
Two-tailed Pearson correlation coefficients of V with other elements measured in floodplain soils of Belgian and European rivers.
Belgian dataGrote BeekGeulLeieFOREGSMg0.793**0.821**0.391**0.778**0.570**Al0.0300.785**0.686**0.870*0.680**P0.626**0.556**−0.0280.2800.374**K0.681**0.343**—0.605**0.433**Ca0.481**0.088−0.1000.1700.064Cr0.729**0.933**0.1770.631**0.665**Mn−0.0790.516**0.2110.397*0.640**Fe0.668**−0.1180.461**0.827**0.896**Co0.201**0.825**0.399**0.672**0.841**Ni0.424**0.551**0.322**0.717**0.706**Cu0.513**0.651**0.235*0.484**0.649**Zn0.0210.665**0.1520.382*0.556**As0.575**0.411**0.249*0.594**—Se0.468**0.365*0.1840.715**—Sr0.667**0.548**0.601**0.375*—Cd0.428**0.314**0.1210.365*—Ba0.439**0.437**0.561**0.628**0.463**Pb−0.1280.615**0.2110.461**0.347**Ti————0.869**OM———0.657**—Clay———0.713**—Loam———0.054—Sand———−0.511**—S———0.463*—*: significant atα=0.05(two-tailed).**: significant atα=0.01 (two-tailed).In general, Mg, Cr, and Al showed the strongest correlation with V (Table5). Among the trace elements, Cr is characterized by a strongly positive correlation with V in the samples of the Grote Beek (R=0.933). Other studies also demonstrated a high correlation between Cr and V in soils and sediments [23, 42]. Cr and V can substitute for each other in the crystal structure of clay minerals and vanadate and chromate show a similar behavior in soils.In the catchment of the Grote Beek, elevated Fe concentrations were generally found, with a maximal concentration of 7.66% Fe, which could be explained by the occurrence of iron-bearing sands and sandstone as well as glauconite in the subsoil [43].Fe and Al also show a positive correlation, which is often observed in soils and sediments (e.g., [15, 42, 44–46]).Data for the clay and organic carbon content were only available for the floodplain soils of the Leie. Vanadium shows a significant positive correlation with the clay content (0.713), which is in accordance with the positive correlation with Al and other clay-associated elements and with the organic carbon content (0.657). Although the association of V with organic material is also shown in other studies [2, 46], some authors [42] report a significant negative correlation between V and organic matter content. Połedniok and Buhl [5] concluded that the amount of V bound to organic material is lower in industrial areas compared to rural areas.For the XRF-data from the FOREGS database, vanadium in the European floodplain soil samples is strongly positively correlated with Fe and Ti, and to a lesser extent with Al.These positive correlations can be explained by the fact that (hydr)oxides of Fe, Ti and Al are good adsorbers for V ions [47, 48] and because of the occurrence of these elements in clay minerals. Moreover, the atomic radius of the V ion is similar to the radius of Fe and Al [2, 15].


#### 3.2.2. Regression Equations to Predict V Concentrations in Alluvial Soils

Multiple linear regressions were performed with Fe, Mg, Ca, K, Al, P, and the clay (<2μm fraction), and organic matter content as independent variables. Clay and organic matter were included because they are already used in the standardization of background values for heavy metals in soils and sediments [49]. Fe is a major component of Fe(hydr)oxides, and is also a constituent of sheet silicates. In soils and sediments, Ca and Mg are dominantly found in sheet silicates when the parent material is Mg-rich (e.g., ultrabasic rocks) and in carbonate minerals (e.g., CaCO3) in nonacidic soils.Regression equations were constructed according to the stepwise method [50]. In order to only include the most significant independent variables. With exception of the soils from the Grote Beek, all regression equations in Table 6 are characterized by one or more outliers. Omission of the outliers from the dataset resulted in a better fit for the regression equations (Table 6). Additionally, the QMC, homoscedasticity, normality of the residues in independent residues improved in most of the cases.Table 6
Regression equations for the different areas with V (logtransformed) as dependent variable.
AreaRegression equationR²3 Belgian floodplainsalog V = −3.172 + 0.433 log Al + 0.367 log Mg + 0.399 log Fe + 0.058 log P**0.831blog V = −2.956 + 0.424 log Al + 0.341 log Mg + 0.312 log Fe + 0.068 log P**0.920floodplain of the Geulalog V = −1.424 + 0.507 log Al + 0.174 log Fe**0.606blog V = −1.342 + 0.503 log Al + 0.160 log Fe**0.634floodplain of the Grote Beekalog V = −2.820 + 0.761 log Al + 0.196 log Mg +0.207 log P **0.809floodplain of the Leiealog V = −4.244 + 0.945 log Al + 0.368 log Mg + 0.315 log clay**0.824blog V = −3.870 + 0.775 log Al + 0.472 log Mg + 0.330 log clay**0.900Europealog V = 1.041 + 1.32 log Fe + 0.413 log Al + 0.102 log Mg**0.836blog V = 1.108 + 1.33 log Fe + 0.279 log Al + 0.086 log Mg**0.880a In the first regression equation, outliers are included.b The second regression equation is without outliers. **Significant at 0.01 level.In general, total V concentrations could be predicted very well by the independent variables Al, Fe, and Mg (Table6). An example of a good fit between measured and predicted values is presented in Figure 3 for the alluvial soils of the Leie.Figure 3
Predicted versus measured V concentrations in the floodplain soils of the Leie river.For the Belgian dataset, Al was an important explaining variable, together with Fe, Mg, and/or P (Table3). Particularly for the Grote Beek alluvial soils, P is an explaining variable in the regression equation for V, which can be related to the emissions from the phosphate ore processing plant.Despite the fact that major elements such as Fe, and Mg are also a component of sheet silicates, they are also major components of resp. Fe-(hydr)oxides and carbonates. Moreover, the attribution of the term “clay fraction” to the <2μm fraction is sometimes misleading. Beside clay minerals, Fe-rich minerals such as Fe-oxi/hydroxides and phyllosilicates are often concentrated in the clay fraction (<2 μm) [51]. A geochemical approach using the total concentrations of conservative elements may be used to overcome these difficulties (e.g., [52]). Before processing the dataset, the potential anthropogenic input of the element and diagenetic processes that may alter its concentrations in sediments should be checked [53]. For the European data, Fe, Al, and Mg were the most significant explaining variables (Table 6).


### 3.3. Single Extractions

#### 3.3.1. Prediction of “Mobile” Metal Concentrations

The composition of soil porewater is important from an environmental point of view because it gives an indication of the “actual mobility” of heavy metals and because the uptake of trace elements by plants occurs via the porewater. Moreover, porewater is also the carrier for elements to the groundwater.The composition of the CaCl2 extract is often considered representative for porewater composition. Contaminants in the porewater can move to deeper soil layers or to the groundwater by infiltrating water. In the present study, average CaCl2-extractable V concentrations were 2 μg L−1 for the samples from the Leie and 4 μg L−1 for the samples from the Grote Beek (Table 7), which is far below the value of 3 mg L−1, which is proposed by Edwards et al. [3] and would cause significant toxic effects to plants.Table 7
Summary statistics for V concentrations in the CaCl2 (in μg L−1) and ammonium-EDTA (in mg kg−1) extracts. Standard deviation is reported within brackets. n: number of samples.
meanmedianminmaxnLeie riverCaCl2 (μg L−1)2.0 (1.7)1.8<0.16.739ammonium-EDTA (mg kg−1)2.4 (2.2)2.0<0.19.839Grote Beek riverCaCl2 (μg L−1)4.4 (4.5)2.0<0.112.227ammonium-EDTA (mg kg−1)11.2 (12.1)5.01.335.327Metal partitioning in soils can be quantified by models in which metal concentrations in the porewater are described as a function of the metal binding solid phases such as Fe- and Al-(hydr)oxides, organic matter, and clay and as a function of soil characteristics that influence heavy metal partitioning, such as pH. The ratio between total metal content bound to a soil relative to its concentration in the soil solution is often represented byKd coefficients. However, such a model assumes that the sorption capacity of a material is independent of the soil properties (organic matter content, pH, clay content, etc.) and therefore, single Kd values are not appropriate to predict metal solubility in soil. Therefore, several authors indicated that metal solubility could be predicted from soil properties. For example, Houba et al. [54] observed a significant relationship between the distribution coefficient (Kd) of Cd in soils and the pH of the CaCl2 extract. According to McBride et al. [55] most of the variability of metal solubility in soil is explained by pH, organic matter content and total metal concentrations.A semimechanistic approach was developed by Sauvé et al. [56] which is based on the assumption that exchangeable metals and protons compete for adsorption on soil-exchange sites. In this model, pH, total metal concentrations, and soil organic matter are used to predict dissolved metal concentrations.In the present study, a modified version of this competitive adsorption model was applied, since ammonium-EDTA-extractable metal concentrations were used instead of total metal concentrations. EDTA extractions are often used to estimate the potentially “available pool” of metals (i.e., the pool that can deliver metals from the solid phase of the soil to the soil solution in a relatively short time period). Stepwise multiple linear regressions was performed with pH(CaCl2), organic carbon content, clay content, CaCl2-extractable V in P concentrations, dissolved organic carbon (DOC) content in the CaCl2-extracts, and ammonium-EDTA-extractable V in P concentrations. Phosphorus was included because phosphate can influence the release of V from soils [7, 13, 57].The distribution of an element between the liquid and solid phase depends strongly on the speciation of the selected element. Because vanadium is known to occur in multiple oxidation states, the actual distribution of the oxidation states was first calculated using the speciation code VisualMinteQ. According to these calculations, the main vanadium species encountered in the CaCl2 extracts is HVO42−, which means that V occurs as an anion that will have the tendency to be desorbed when pH rises.For the Grote Beek floodplain soils, only the pH of the CaCl2 extract and ammonium-EDTA-extractable V concentrations were retained in the final regression model, as the other variables did not contribute to an improvement of the regression equation. For the Grote Beek river catchment, “mobile” V concentrations in floodplain soils could be predicted using the following fitting equation:  (1)log⁡[V]s=5.356-1.041pH+1.338log⁡[V]aR2=0.840
with [V]s = the V concentration in the CaCl2 extract (mg L−1), [V]a = the “available pool” of V (mg kg−1).For the floodplain soils along the Leie, the pH of the CaCl2 extract, EDTA-extractable V-concentrations and organic carbon content were retained in the final regression model, resulting in the following regression equation: (2)log⁡⁡[V]s=1.313-1.969log⁡clay+0.378log⁡[P]a+1.150log⁡OCR2=0.590
with [V]s = the V concentration in the CaCl2 extract (mg L−1), [P]a = the “available pool” of P (mg kg−1), OC = the organic carbon content (%), clay = the clay content (in %).


#### 3.3.2. “Potential Availability” of V

Ammonium-EDTA was also used to estimate the “mobilizable” metal concentrations. EDTA is a nonselective reagent that exhibits a strong capacity to complex metals. EDTA was shown to dissolve carbonates, thereby mobilizing occluded elements [58]. Borggaard [59] showed that EDTA extracts amorphous Fe-oxides, but this dissolution is very slow in the presence of other metal-chelate complexes [60]. It is also able to form organometal complexes, which compete with organic matter in soil. Additionally, complexes of EDTA with heavy metals such as Cd, Cu, Pb, Ni, and Zn are more stable than complexes between V and EDTA, which results in a lower ammonium-EDTA-extractable fraction of V when the heavy metals mentioned before are present in high concentrations [61]. The average EDTA-extractable V concentration in the floodplain soils of the Leie was 2.46 mg kg−1 (Table 7) which is comparable with average EDTA-extractable V concentrations (3.40 mg kg−1) reported by Gäbler et al. [61]. In the samples from the Grote Beek river, the average ammonium-EDTA-extractable V concentration was 11.2 mg/kg. Floodplain soil samples along the Grote Beek are contaminated with V and display elevated V concentrations, which can also explain for the higher ‘potentially mobile V-content in these soils. In relative concentrations, however, only 3.5% and 9.8% of the total V content in the floodplain soil samples of, respectively, the Leie and the Grote Beek were extracted with ammonium-EDTA, pointing to the fact that the majority of the V in the samples is characterized by a very low mobility. In the highly contaminated samples of the Grote Beek, a more detailed study of V mobility may nevertheless be useful to elucidate the potential release of V under changing environmental conditions (changes in pH, Eh, etc.).


## 3.1. Total Concentrations of V in Floodplain Soils

Average total V concentration in the 3 Belgian catchments and in the European FOREGS database were in the range 45–87 mg kg−1. The high average V concentrations in the alluvial soils of the Grote Beek catchment is mainly due to 7 outliers with V concentrations between 189 and 301 mg kg−1. Alluvial soils of the Geul river are characterized by significant lower V concentrations compared to the Leie and Grote Beek rivers. V concentrations in floodplain soils are rarely reported in literature. Peh and Miko [36] measured average V concentrations of 70 mg kg−1 in Croatia, which is in line with the values reported here. Ivanov and Kashin [37] recorded (AAS-measurement after aqua regia destruction) average V concentrations of 110 mg kg−1 in alluvial deposits in Transbaikalia (Russia), and Wigginton and Price [38] found (ICP-OES determination after aqua regia destruction) concentrations between 8 and 17 mg kg−1 in floodplain soils of Bayou Creek (UK). On the 4th October 2010, red mud suspension were released from the Ajkai Timfoldgyar Zrt alumina plant. Among other trace elements, mud samples at the source contained important V concentrations (>1000 mg kg−1), and floodplain sediments downstream were significantly affected by the red mud [39].Total V concentrations in floodplain soils from the FOREGS database are represented on Figure2. Since the data presented on Figure 2 are characterized by a right kurtosis and thus not normal distribution, concentrations represented on the map are strongly influenced by locally enhanced V concentrations [40]. Data for the Belgian catchments are not presented on a map, since the data originate from 3 distinct locations, only covering a limited area in Belgium.Figure 2
Total V concentrations in alluvial soils samples in 26 European countries [27].Data observations which lied more than 1.5 interquartile range lower than the first quartile or 1.5 interquartile range higher than the third quartile were considered as an outlier. In the Belgian data for V, 7 outliers were observed, between 189–301 mg kg−1, all coming from samples from the alluvial soil samples of the Grote Beek. De Grote Beek is also characterized by the highest average V concentrations and an important difference between average and median values for V, which indicates an anthropogenic source for V [33]. In this area, phosphate ores are processed by a plant that is emitting its waste water directly into the river. V is known as a secondary product of phosphate production [8]. In accordance to this, the correlation coefficient between V and P in alluvial soils is also very high in the samples from the Grote Beek river.When the V concentrations in these samples are compared with the norm values in Table1, 11 samples of the Grote Beek alluvial plain exceed the value for “unacceptable risk” in residential areas established in Finland, Lithuania, and Russia (150 mg kg−1), whereas 3 samples exceed the Dutch “negligible risk” value of 250 mg kg−1. The samples from the catchment of the Leie river are contaminated with heavy metals and often characterized by an elevated clay- and organic matter content. Only one sample from alluvial plain of the Leie river is characterized by a V content above 150 mg kg−1 (total V concentration is 161 mg kg−1).The samples from the Geul-catchment, finally, are characterized by relatively low total V concentrations (average 45 mg kg−1, Table 4). The average V concentration in the alluvial soils of the Geul river is below the 90th percentile of the European dataset, which could be considered as the background value [33], as well as below the average V concentration in European alluvial soils. The maximal V concentrations is 57 mg kg−1 (Table 4) and thus below all threshold values established in Table 1, except the value for “negligible risk” from The Netherlands (42 mg kg−1). The average V concentration in the soils along the Geul is very close to the median concentration, which indicates that V is not of anthropogenic origin [33]. The subsoil in the Geul catchment consists of sand and shale and other studies also showed that total V concentrations in soils developed on shale are lower compared to more clay-rich substrates [41].Table 4
Summary statistics for total V concentrations (in mg kg−1), based on Belgian (this study) and European [28] data. Standard deviation is reported within brackets. n: number of samples.
datasetareameanmedianminmaxnBelgian samples3 Belgian floodplains69 (43)5318301206floodplain of the Geul45 (6)46305775floodplain of the Grote Beek87 (55)662830194floodplain of the Leie71 (28)691816037FOREGSEurope (XRF)59 (35)561266749Europe (aqua regia)56747In the European FOREGS dataset, “the spatial distribution of V is clearly related to the bedrock geology and mineralization, especially mafic and ultramafic lithology, and also clay-rich soil with high Al2O3 content” [15]. The highest V concentration is found in the mineralized Oslo rift (266 mg kg−1). In England and Germany, elevated V concentrations (224 mg kg−1) are respectively due to industrial pollution and coal and oil combustion [15].


## 3.2. Relation of V with Other Elements

### 3.2.1. Correlations

Because the majority of the data (trace element and major element concentrations) showed a right kurtosis and were significant on the K-S test, the data were log transformed. Concerning the FOREGS data, the correlation coefficients in Table5 were calculated with the data from the XRF determination, since more elements that are relevant for this paper were measured with this method compared with the aqua regia method. Unfortunately, it was not possible to merge the results of the aqua regia and XRF determination into one table, because the data were not in the same order and there were several missing samples for each of the methods. This inconvenience of the FOREGS database has also been mentioned by [40].Table 5
Two-tailed Pearson correlation coefficients of V with other elements measured in floodplain soils of Belgian and European rivers.
Belgian dataGrote BeekGeulLeieFOREGSMg0.793**0.821**0.391**0.778**0.570**Al0.0300.785**0.686**0.870*0.680**P0.626**0.556**−0.0280.2800.374**K0.681**0.343**—0.605**0.433**Ca0.481**0.088−0.1000.1700.064Cr0.729**0.933**0.1770.631**0.665**Mn−0.0790.516**0.2110.397*0.640**Fe0.668**−0.1180.461**0.827**0.896**Co0.201**0.825**0.399**0.672**0.841**Ni0.424**0.551**0.322**0.717**0.706**Cu0.513**0.651**0.235*0.484**0.649**Zn0.0210.665**0.1520.382*0.556**As0.575**0.411**0.249*0.594**—Se0.468**0.365*0.1840.715**—Sr0.667**0.548**0.601**0.375*—Cd0.428**0.314**0.1210.365*—Ba0.439**0.437**0.561**0.628**0.463**Pb−0.1280.615**0.2110.461**0.347**Ti————0.869**OM———0.657**—Clay———0.713**—Loam———0.054—Sand———−0.511**—S———0.463*—*: significant atα=0.05(two-tailed).**: significant atα=0.01 (two-tailed).In general, Mg, Cr, and Al showed the strongest correlation with V (Table5). Among the trace elements, Cr is characterized by a strongly positive correlation with V in the samples of the Grote Beek (R=0.933). Other studies also demonstrated a high correlation between Cr and V in soils and sediments [23, 42]. Cr and V can substitute for each other in the crystal structure of clay minerals and vanadate and chromate show a similar behavior in soils.In the catchment of the Grote Beek, elevated Fe concentrations were generally found, with a maximal concentration of 7.66% Fe, which could be explained by the occurrence of iron-bearing sands and sandstone as well as glauconite in the subsoil [43].Fe and Al also show a positive correlation, which is often observed in soils and sediments (e.g., [15, 42, 44–46]).Data for the clay and organic carbon content were only available for the floodplain soils of the Leie. Vanadium shows a significant positive correlation with the clay content (0.713), which is in accordance with the positive correlation with Al and other clay-associated elements and with the organic carbon content (0.657). Although the association of V with organic material is also shown in other studies [2, 46], some authors [42] report a significant negative correlation between V and organic matter content. Połedniok and Buhl [5] concluded that the amount of V bound to organic material is lower in industrial areas compared to rural areas.For the XRF-data from the FOREGS database, vanadium in the European floodplain soil samples is strongly positively correlated with Fe and Ti, and to a lesser extent with Al.These positive correlations can be explained by the fact that (hydr)oxides of Fe, Ti and Al are good adsorbers for V ions [47, 48] and because of the occurrence of these elements in clay minerals. Moreover, the atomic radius of the V ion is similar to the radius of Fe and Al [2, 15].


### 3.2.2. Regression Equations to Predict V Concentrations in Alluvial Soils

Multiple linear regressions were performed with Fe, Mg, Ca, K, Al, P, and the clay (<2μm fraction), and organic matter content as independent variables. Clay and organic matter were included because they are already used in the standardization of background values for heavy metals in soils and sediments [49]. Fe is a major component of Fe(hydr)oxides, and is also a constituent of sheet silicates. In soils and sediments, Ca and Mg are dominantly found in sheet silicates when the parent material is Mg-rich (e.g., ultrabasic rocks) and in carbonate minerals (e.g., CaCO3) in nonacidic soils.Regression equations were constructed according to the stepwise method [50]. In order to only include the most significant independent variables. With exception of the soils from the Grote Beek, all regression equations in Table 6 are characterized by one or more outliers. Omission of the outliers from the dataset resulted in a better fit for the regression equations (Table 6). Additionally, the QMC, homoscedasticity, normality of the residues in independent residues improved in most of the cases.Table 6
Regression equations for the different areas with V (logtransformed) as dependent variable.
AreaRegression equationR²3 Belgian floodplainsalog V = −3.172 + 0.433 log Al + 0.367 log Mg + 0.399 log Fe + 0.058 log P**0.831blog V = −2.956 + 0.424 log Al + 0.341 log Mg + 0.312 log Fe + 0.068 log P**0.920floodplain of the Geulalog V = −1.424 + 0.507 log Al + 0.174 log Fe**0.606blog V = −1.342 + 0.503 log Al + 0.160 log Fe**0.634floodplain of the Grote Beekalog V = −2.820 + 0.761 log Al + 0.196 log Mg +0.207 log P **0.809floodplain of the Leiealog V = −4.244 + 0.945 log Al + 0.368 log Mg + 0.315 log clay**0.824blog V = −3.870 + 0.775 log Al + 0.472 log Mg + 0.330 log clay**0.900Europealog V = 1.041 + 1.32 log Fe + 0.413 log Al + 0.102 log Mg**0.836blog V = 1.108 + 1.33 log Fe + 0.279 log Al + 0.086 log Mg**0.880a In the first regression equation, outliers are included.b The second regression equation is without outliers. **Significant at 0.01 level.In general, total V concentrations could be predicted very well by the independent variables Al, Fe, and Mg (Table6). An example of a good fit between measured and predicted values is presented in Figure 3 for the alluvial soils of the Leie.Figure 3
Predicted versus measured V concentrations in the floodplain soils of the Leie river.For the Belgian dataset, Al was an important explaining variable, together with Fe, Mg, and/or P (Table3). Particularly for the Grote Beek alluvial soils, P is an explaining variable in the regression equation for V, which can be related to the emissions from the phosphate ore processing plant.Despite the fact that major elements such as Fe, and Mg are also a component of sheet silicates, they are also major components of resp. Fe-(hydr)oxides and carbonates. Moreover, the attribution of the term “clay fraction” to the <2μm fraction is sometimes misleading. Beside clay minerals, Fe-rich minerals such as Fe-oxi/hydroxides and phyllosilicates are often concentrated in the clay fraction (<2 μm) [51]. A geochemical approach using the total concentrations of conservative elements may be used to overcome these difficulties (e.g., [52]). Before processing the dataset, the potential anthropogenic input of the element and diagenetic processes that may alter its concentrations in sediments should be checked [53]. For the European data, Fe, Al, and Mg were the most significant explaining variables (Table 6).


## 3.2.1. Correlations

Because the majority of the data (trace element and major element concentrations) showed a right kurtosis and were significant on the K-S test, the data were log transformed. Concerning the FOREGS data, the correlation coefficients in Table5 were calculated with the data from the XRF determination, since more elements that are relevant for this paper were measured with this method compared with the aqua regia method. Unfortunately, it was not possible to merge the results of the aqua regia and XRF determination into one table, because the data were not in the same order and there were several missing samples for each of the methods. This inconvenience of the FOREGS database has also been mentioned by [40].Table 5
Two-tailed Pearson correlation coefficients of V with other elements measured in floodplain soils of Belgian and European rivers.
Belgian dataGrote BeekGeulLeieFOREGSMg0.793**0.821**0.391**0.778**0.570**Al0.0300.785**0.686**0.870*0.680**P0.626**0.556**−0.0280.2800.374**K0.681**0.343**—0.605**0.433**Ca0.481**0.088−0.1000.1700.064Cr0.729**0.933**0.1770.631**0.665**Mn−0.0790.516**0.2110.397*0.640**Fe0.668**−0.1180.461**0.827**0.896**Co0.201**0.825**0.399**0.672**0.841**Ni0.424**0.551**0.322**0.717**0.706**Cu0.513**0.651**0.235*0.484**0.649**Zn0.0210.665**0.1520.382*0.556**As0.575**0.411**0.249*0.594**—Se0.468**0.365*0.1840.715**—Sr0.667**0.548**0.601**0.375*—Cd0.428**0.314**0.1210.365*—Ba0.439**0.437**0.561**0.628**0.463**Pb−0.1280.615**0.2110.461**0.347**Ti————0.869**OM———0.657**—Clay———0.713**—Loam———0.054—Sand———−0.511**—S———0.463*—*: significant atα=0.05(two-tailed).**: significant atα=0.01 (two-tailed).In general, Mg, Cr, and Al showed the strongest correlation with V (Table5). Among the trace elements, Cr is characterized by a strongly positive correlation with V in the samples of the Grote Beek (R=0.933). Other studies also demonstrated a high correlation between Cr and V in soils and sediments [23, 42]. Cr and V can substitute for each other in the crystal structure of clay minerals and vanadate and chromate show a similar behavior in soils.In the catchment of the Grote Beek, elevated Fe concentrations were generally found, with a maximal concentration of 7.66% Fe, which could be explained by the occurrence of iron-bearing sands and sandstone as well as glauconite in the subsoil [43].Fe and Al also show a positive correlation, which is often observed in soils and sediments (e.g., [15, 42, 44–46]).Data for the clay and organic carbon content were only available for the floodplain soils of the Leie. Vanadium shows a significant positive correlation with the clay content (0.713), which is in accordance with the positive correlation with Al and other clay-associated elements and with the organic carbon content (0.657). Although the association of V with organic material is also shown in other studies [2, 46], some authors [42] report a significant negative correlation between V and organic matter content. Połedniok and Buhl [5] concluded that the amount of V bound to organic material is lower in industrial areas compared to rural areas.For the XRF-data from the FOREGS database, vanadium in the European floodplain soil samples is strongly positively correlated with Fe and Ti, and to a lesser extent with Al.These positive correlations can be explained by the fact that (hydr)oxides of Fe, Ti and Al are good adsorbers for V ions [47, 48] and because of the occurrence of these elements in clay minerals. Moreover, the atomic radius of the V ion is similar to the radius of Fe and Al [2, 15].


## 3.2.2. Regression Equations to Predict V Concentrations in Alluvial Soils

Multiple linear regressions were performed with Fe, Mg, Ca, K, Al, P, and the clay (<2μm fraction), and organic matter content as independent variables. Clay and organic matter were included because they are already used in the standardization of background values for heavy metals in soils and sediments [49]. Fe is a major component of Fe(hydr)oxides, and is also a constituent of sheet silicates. In soils and sediments, Ca and Mg are dominantly found in sheet silicates when the parent material is Mg-rich (e.g., ultrabasic rocks) and in carbonate minerals (e.g., CaCO3) in nonacidic soils.Regression equations were constructed according to the stepwise method [50]. In order to only include the most significant independent variables. With exception of the soils from the Grote Beek, all regression equations in Table 6 are characterized by one or more outliers. Omission of the outliers from the dataset resulted in a better fit for the regression equations (Table 6). Additionally, the QMC, homoscedasticity, normality of the residues in independent residues improved in most of the cases.Table 6
Regression equations for the different areas with V (logtransformed) as dependent variable.
AreaRegression equationR²3 Belgian floodplainsalog V = −3.172 + 0.433 log Al + 0.367 log Mg + 0.399 log Fe + 0.058 log P**0.831blog V = −2.956 + 0.424 log Al + 0.341 log Mg + 0.312 log Fe + 0.068 log P**0.920floodplain of the Geulalog V = −1.424 + 0.507 log Al + 0.174 log Fe**0.606blog V = −1.342 + 0.503 log Al + 0.160 log Fe**0.634floodplain of the Grote Beekalog V = −2.820 + 0.761 log Al + 0.196 log Mg +0.207 log P **0.809floodplain of the Leiealog V = −4.244 + 0.945 log Al + 0.368 log Mg + 0.315 log clay**0.824blog V = −3.870 + 0.775 log Al + 0.472 log Mg + 0.330 log clay**0.900Europealog V = 1.041 + 1.32 log Fe + 0.413 log Al + 0.102 log Mg**0.836blog V = 1.108 + 1.33 log Fe + 0.279 log Al + 0.086 log Mg**0.880a In the first regression equation, outliers are included.b The second regression equation is without outliers. **Significant at 0.01 level.In general, total V concentrations could be predicted very well by the independent variables Al, Fe, and Mg (Table6). An example of a good fit between measured and predicted values is presented in Figure 3 for the alluvial soils of the Leie.Figure 3
Predicted versus measured V concentrations in the floodplain soils of the Leie river.For the Belgian dataset, Al was an important explaining variable, together with Fe, Mg, and/or P (Table3). Particularly for the Grote Beek alluvial soils, P is an explaining variable in the regression equation for V, which can be related to the emissions from the phosphate ore processing plant.Despite the fact that major elements such as Fe, and Mg are also a component of sheet silicates, they are also major components of resp. Fe-(hydr)oxides and carbonates. Moreover, the attribution of the term “clay fraction” to the <2μm fraction is sometimes misleading. Beside clay minerals, Fe-rich minerals such as Fe-oxi/hydroxides and phyllosilicates are often concentrated in the clay fraction (<2 μm) [51]. A geochemical approach using the total concentrations of conservative elements may be used to overcome these difficulties (e.g., [52]). Before processing the dataset, the potential anthropogenic input of the element and diagenetic processes that may alter its concentrations in sediments should be checked [53]. For the European data, Fe, Al, and Mg were the most significant explaining variables (Table 6).


## 3.3. Single Extractions

### 3.3.1. Prediction of “Mobile” Metal Concentrations

The composition of soil porewater is important from an environmental point of view because it gives an indication of the “actual mobility” of heavy metals and because the uptake of trace elements by plants occurs via the porewater. Moreover, porewater is also the carrier for elements to the groundwater.The composition of the CaCl2 extract is often considered representative for porewater composition. Contaminants in the porewater can move to deeper soil layers or to the groundwater by infiltrating water. In the present study, average CaCl2-extractable V concentrations were 2 μg L−1 for the samples from the Leie and 4 μg L−1 for the samples from the Grote Beek (Table 7), which is far below the value of 3 mg L−1, which is proposed by Edwards et al. [3] and would cause significant toxic effects to plants.Table 7
Summary statistics for V concentrations in the CaCl2 (in μg L−1) and ammonium-EDTA (in mg kg−1) extracts. Standard deviation is reported within brackets. n: number of samples.
meanmedianminmaxnLeie riverCaCl2 (μg L−1)2.0 (1.7)1.8<0.16.739ammonium-EDTA (mg kg−1)2.4 (2.2)2.0<0.19.839Grote Beek riverCaCl2 (μg L−1)4.4 (4.5)2.0<0.112.227ammonium-EDTA (mg kg−1)11.2 (12.1)5.01.335.327Metal partitioning in soils can be quantified by models in which metal concentrations in the porewater are described as a function of the metal binding solid phases such as Fe- and Al-(hydr)oxides, organic matter, and clay and as a function of soil characteristics that influence heavy metal partitioning, such as pH. The ratio between total metal content bound to a soil relative to its concentration in the soil solution is often represented byKd coefficients. However, such a model assumes that the sorption capacity of a material is independent of the soil properties (organic matter content, pH, clay content, etc.) and therefore, single Kd values are not appropriate to predict metal solubility in soil. Therefore, several authors indicated that metal solubility could be predicted from soil properties. For example, Houba et al. [54] observed a significant relationship between the distribution coefficient (Kd) of Cd in soils and the pH of the CaCl2 extract. According to McBride et al. [55] most of the variability of metal solubility in soil is explained by pH, organic matter content and total metal concentrations.A semimechanistic approach was developed by Sauvé et al. [56] which is based on the assumption that exchangeable metals and protons compete for adsorption on soil-exchange sites. In this model, pH, total metal concentrations, and soil organic matter are used to predict dissolved metal concentrations.In the present study, a modified version of this competitive adsorption model was applied, since ammonium-EDTA-extractable metal concentrations were used instead of total metal concentrations. EDTA extractions are often used to estimate the potentially “available pool” of metals (i.e., the pool that can deliver metals from the solid phase of the soil to the soil solution in a relatively short time period). Stepwise multiple linear regressions was performed with pH(CaCl2), organic carbon content, clay content, CaCl2-extractable V in P concentrations, dissolved organic carbon (DOC) content in the CaCl2-extracts, and ammonium-EDTA-extractable V in P concentrations. Phosphorus was included because phosphate can influence the release of V from soils [7, 13, 57].The distribution of an element between the liquid and solid phase depends strongly on the speciation of the selected element. Because vanadium is known to occur in multiple oxidation states, the actual distribution of the oxidation states was first calculated using the speciation code VisualMinteQ. According to these calculations, the main vanadium species encountered in the CaCl2 extracts is HVO42−, which means that V occurs as an anion that will have the tendency to be desorbed when pH rises.For the Grote Beek floodplain soils, only the pH of the CaCl2 extract and ammonium-EDTA-extractable V concentrations were retained in the final regression model, as the other variables did not contribute to an improvement of the regression equation. For the Grote Beek river catchment, “mobile” V concentrations in floodplain soils could be predicted using the following fitting equation:  (1)log⁡[V]s=5.356-1.041pH+1.338log⁡[V]aR2=0.840
with [V]s = the V concentration in the CaCl2 extract (mg L−1), [V]a = the “available pool” of V (mg kg−1).For the floodplain soils along the Leie, the pH of the CaCl2 extract, EDTA-extractable V-concentrations and organic carbon content were retained in the final regression model, resulting in the following regression equation: (2)log⁡⁡[V]s=1.313-1.969log⁡clay+0.378log⁡[P]a+1.150log⁡OCR2=0.590
with [V]s = the V concentration in the CaCl2 extract (mg L−1), [P]a = the “available pool” of P (mg kg−1), OC = the organic carbon content (%), clay = the clay content (in %).


### 3.3.2. “Potential Availability” of V

Ammonium-EDTA was also used to estimate the “mobilizable” metal concentrations. EDTA is a nonselective reagent that exhibits a strong capacity to complex metals. EDTA was shown to dissolve carbonates, thereby mobilizing occluded elements [58]. Borggaard [59] showed that EDTA extracts amorphous Fe-oxides, but this dissolution is very slow in the presence of other metal-chelate complexes [60]. It is also able to form organometal complexes, which compete with organic matter in soil. Additionally, complexes of EDTA with heavy metals such as Cd, Cu, Pb, Ni, and Zn are more stable than complexes between V and EDTA, which results in a lower ammonium-EDTA-extractable fraction of V when the heavy metals mentioned before are present in high concentrations [61]. The average EDTA-extractable V concentration in the floodplain soils of the Leie was 2.46 mg kg−1 (Table 7) which is comparable with average EDTA-extractable V concentrations (3.40 mg kg−1) reported by Gäbler et al. [61]. In the samples from the Grote Beek river, the average ammonium-EDTA-extractable V concentration was 11.2 mg/kg. Floodplain soil samples along the Grote Beek are contaminated with V and display elevated V concentrations, which can also explain for the higher ‘potentially mobile V-content in these soils. In relative concentrations, however, only 3.5% and 9.8% of the total V content in the floodplain soil samples of, respectively, the Leie and the Grote Beek were extracted with ammonium-EDTA, pointing to the fact that the majority of the V in the samples is characterized by a very low mobility. In the highly contaminated samples of the Grote Beek, a more detailed study of V mobility may nevertheless be useful to elucidate the potential release of V under changing environmental conditions (changes in pH, Eh, etc.).


## 3.3.1. Prediction of “Mobile” Metal Concentrations

The composition of soil porewater is important from an environmental point of view because it gives an indication of the “actual mobility” of heavy metals and because the uptake of trace elements by plants occurs via the porewater. Moreover, porewater is also the carrier for elements to the groundwater.The composition of the CaCl2 extract is often considered representative for porewater composition. Contaminants in the porewater can move to deeper soil layers or to the groundwater by infiltrating water. In the present study, average CaCl2-extractable V concentrations were 2 μg L−1 for the samples from the Leie and 4 μg L−1 for the samples from the Grote Beek (Table 7), which is far below the value of 3 mg L−1, which is proposed by Edwards et al. [3] and would cause significant toxic effects to plants.Table 7
Summary statistics for V concentrations in the CaCl2 (in μg L−1) and ammonium-EDTA (in mg kg−1) extracts. Standard deviation is reported within brackets. n: number of samples.
meanmedianminmaxnLeie riverCaCl2 (μg L−1)2.0 (1.7)1.8<0.16.739ammonium-EDTA (mg kg−1)2.4 (2.2)2.0<0.19.839Grote Beek riverCaCl2 (μg L−1)4.4 (4.5)2.0<0.112.227ammonium-EDTA (mg kg−1)11.2 (12.1)5.01.335.327Metal partitioning in soils can be quantified by models in which metal concentrations in the porewater are described as a function of the metal binding solid phases such as Fe- and Al-(hydr)oxides, organic matter, and clay and as a function of soil characteristics that influence heavy metal partitioning, such as pH. The ratio between total metal content bound to a soil relative to its concentration in the soil solution is often represented byKd coefficients. However, such a model assumes that the sorption capacity of a material is independent of the soil properties (organic matter content, pH, clay content, etc.) and therefore, single Kd values are not appropriate to predict metal solubility in soil. Therefore, several authors indicated that metal solubility could be predicted from soil properties. For example, Houba et al. [54] observed a significant relationship between the distribution coefficient (Kd) of Cd in soils and the pH of the CaCl2 extract. According to McBride et al. [55] most of the variability of metal solubility in soil is explained by pH, organic matter content and total metal concentrations.A semimechanistic approach was developed by Sauvé et al. [56] which is based on the assumption that exchangeable metals and protons compete for adsorption on soil-exchange sites. In this model, pH, total metal concentrations, and soil organic matter are used to predict dissolved metal concentrations.In the present study, a modified version of this competitive adsorption model was applied, since ammonium-EDTA-extractable metal concentrations were used instead of total metal concentrations. EDTA extractions are often used to estimate the potentially “available pool” of metals (i.e., the pool that can deliver metals from the solid phase of the soil to the soil solution in a relatively short time period). Stepwise multiple linear regressions was performed with pH(CaCl2), organic carbon content, clay content, CaCl2-extractable V in P concentrations, dissolved organic carbon (DOC) content in the CaCl2-extracts, and ammonium-EDTA-extractable V in P concentrations. Phosphorus was included because phosphate can influence the release of V from soils [7, 13, 57].The distribution of an element between the liquid and solid phase depends strongly on the speciation of the selected element. Because vanadium is known to occur in multiple oxidation states, the actual distribution of the oxidation states was first calculated using the speciation code VisualMinteQ. According to these calculations, the main vanadium species encountered in the CaCl2 extracts is HVO42−, which means that V occurs as an anion that will have the tendency to be desorbed when pH rises.For the Grote Beek floodplain soils, only the pH of the CaCl2 extract and ammonium-EDTA-extractable V concentrations were retained in the final regression model, as the other variables did not contribute to an improvement of the regression equation. For the Grote Beek river catchment, “mobile” V concentrations in floodplain soils could be predicted using the following fitting equation:  (1)log⁡[V]s=5.356-1.041pH+1.338log⁡[V]aR2=0.840
with [V]s = the V concentration in the CaCl2 extract (mg L−1), [V]a = the “available pool” of V (mg kg−1).For the floodplain soils along the Leie, the pH of the CaCl2 extract, EDTA-extractable V-concentrations and organic carbon content were retained in the final regression model, resulting in the following regression equation: (2)log⁡⁡[V]s=1.313-1.969log⁡clay+0.378log⁡[P]a+1.150log⁡OCR2=0.590
with [V]s = the V concentration in the CaCl2 extract (mg L−1), [P]a = the “available pool” of P (mg kg−1), OC = the organic carbon content (%), clay = the clay content (in %).


## 3.3.2. “Potential Availability” of V

Ammonium-EDTA was also used to estimate the “mobilizable” metal concentrations. EDTA is a nonselective reagent that exhibits a strong capacity to complex metals. EDTA was shown to dissolve carbonates, thereby mobilizing occluded elements [58]. Borggaard [59] showed that EDTA extracts amorphous Fe-oxides, but this dissolution is very slow in the presence of other metal-chelate complexes [60]. It is also able to form organometal complexes, which compete with organic matter in soil. Additionally, complexes of EDTA with heavy metals such as Cd, Cu, Pb, Ni, and Zn are more stable than complexes between V and EDTA, which results in a lower ammonium-EDTA-extractable fraction of V when the heavy metals mentioned before are present in high concentrations [61]. The average EDTA-extractable V concentration in the floodplain soils of the Leie was 2.46 mg kg−1 (Table 7) which is comparable with average EDTA-extractable V concentrations (3.40 mg kg−1) reported by Gäbler et al. [61]. In the samples from the Grote Beek river, the average ammonium-EDTA-extractable V concentration was 11.2 mg/kg. Floodplain soil samples along the Grote Beek are contaminated with V and display elevated V concentrations, which can also explain for the higher ‘potentially mobile V-content in these soils. In relative concentrations, however, only 3.5% and 9.8% of the total V content in the floodplain soil samples of, respectively, the Leie and the Grote Beek were extracted with ammonium-EDTA, pointing to the fact that the majority of the V in the samples is characterized by a very low mobility. In the highly contaminated samples of the Grote Beek, a more detailed study of V mobility may nevertheless be useful to elucidate the potential release of V under changing environmental conditions (changes in pH, Eh, etc.).


## 4. Conclusion

Average V concentrations in Belgian alluvial soils (69 mg kg−1) are in general comparable with the average V content in European alluvial soils (56–59 mg kg−1). In one of the river catchments analyzed in this study, a contamination of alluvial soils, originating from a phosphate ore processing plant, was detected. By analysis of both the Belgian and European data, the relationship between concentrations of major elements and V was quantified. Fe, Al, and Mg were the most significant variables that could predict the total V content in the alluvial soils and the obtained regression equation could be used to predict V concentrations in alluvial soils. For the alluvial soils contaminated by the emissions from the phosphate ore treatment plant, P was also an explaining variable. “Mobile” V concentrations, as estimated by the amount of V released by a single extraction with CaCl2 0.01 mol L−1, were low, even in the most contaminated soil samples. Despite the low actual mobility of V, it might nevertheless be useful to study the effect of changing environmental conditions such as soil acidification and fluctuating redox conditions (e.g., by inundation of the floodplain soils). This will be the subject of a followup study, in order to obtain a better insight in the geochemistry of V in soils and sediments.Finally, it is clear that the procedures for the determination of total element concentrations in soils and sediments should be well established, especially in large-scale studies, since this is essential to compare results of different studies and/or laboratories.



---

*Source: 979501-2012-04-10.xml*