Modelling the concentrations of dissolved contaminants (Cd, Cu, Ni, Pb, Zn) in floodplain soils

Potentially toxic elements in surface soils of the Lower Don floodplain and the Taganrog Bay coast: sources, spatial distribution and pollution assessment

The pollution of floodplain, deltaic and adjacent coastal soils in large fluvial systems, considered an urgent environmental problem, as well as potentially toxic elements in such environments, can negatively affect aquatic ecosystems, as well as pose significant risks to human health. This paper is devoted to the geochemistry of potentially toxic elements in soils of the Lower Don basin, which is one of the largest and most anthropogenically transformed water bodies in Southern Russia, as well as the adjacent areas of the Taganrog Bay coast. The median element concentrations in the soils of the study area were consistent with the world soil average and the contents of elements in background soils. Comparative assessment of the spatial distributions as well as the results of Pearson's correlations, cluster analysis and principal component analysis showed that Cr, Ni, Cu and Zn are predominantly of natural origin; Mn and As are of mixed sources; and Cd and Pb are predominantly of anthropogenic origin. The geochemical anomalies of elements were associated with the impact of local anthropogenic sources. Geochemical background values for Cr, Mn, Ni, Cu, Zn, As, Cd and Pb in the soils of the Lower Don and the Taganrog Bay coast determined using the 'median + 2 median absolute deviations' approach are presented. The highest values of the integrated pollution indices were observed in floodplain soils of small rivers.

NOM-mineral interaction: Significance for speciation of cations and anions

In this study, the nano-scale spatial distribution of natural organic matter (NOM) on the surface of iron (hydr)oxides and its relevance to oxyanion (PO₄^3-) and metal cation (Cd²⁺ and Cu²⁺) adsorption to the assemblage of oxide (goethite) and NOM (humic acids (HA) or fulvic acids (FA)) was investigated with experiments and advanced surface complexation modeling. Both the linear additive Multi-Surface model (MSM) and the more sophisticated Natural Organic Matter-Charge Distribution (NOM-CD) model were used. The MSM model ignores the effects of NOM-mineral interaction on ion adsorption, whereas the NOM-CD model considers this effect. The results showed that with the increase of NOM loading on oxides, deviation between the MSM and NOM-CD model became bigger for PO₄^3-, but smaller for Cd²⁺ and Cu²⁺. Oxyanions bind mainly to oxides and therefore the competitive effect of NOM cannot be neglected, which explains the large difference between these two models for PO₄^3-. On the contrary, at a relatively high NOM loading, a large fraction of NOM extends further away from the surface of oxides. Thus for metal cations that bind mainly to NOM, the influence of NOM-mineral interaction on their adsorption is small and the results of the MSM and NOM-CD model are similar. In top soils, the NOM loading on oxides is often high, therefore the linear additive MSM is applicable for metal cation speciation calculations as reported in many literatures. An approach based on the NOM-CD model was proposed, which can not only calculate the macroscopic solid-solution distribution of both cations and anions, but can also provide information regarding their microscopic surface speciation.

Field observations to establish the impact of fluvial flooding on potentially toxic element (PTE) mobility in floodplain soils

Inundation of river water during flooding deposits contaminated sediments onto floodplain topsoil. Historically, floodplains were considered an important sink for potentially toxic elements (PTEs). With increasing flood frequency and duration, due to climate change and land use change, it is important to understand the impact that further flooding may have on this legacy contamination. In this study a field-based approach was taken, extracting soil pore waters by centrifugation of soils sampled on multiple occasions from multiple locations across a floodplain site, which lies adjacent to the River Loddon in southeast England. Flooding generally decreased pore water PTE concentrations and significantly lower pore water concentrations of Cd, Cu, and Cr were found post-flood compared to pre-flood. The dominant process responsible for this observation was precipitation with sulphides resulting in PTE removal from the pore water post-flood. The changes in pH were found to be associated with the decreased pore water concentration of Cu, which suggests the pH rise may have aided adsorption mechanisms or precipitation with phosphates. The impact of flooding on the release and retention of PTEs in floodplain soils is the net effect of several key processes occurring concurrently. It is important to understand the dominant processes that drive mobility of individual PTEs on specific floodplains so that site-specific predictions can determine the impact of future floods on the environmental fate of legacy contaminants.

Simple models efficiently predict free cadmium Cd2+ in the solutions of low-contaminated agricultural soils

Speciation of Cd in soil solutions strongly determines the fate of this toxic metal in the environment. Generally, in soil solutions, Cd predominantly binds to the dissolved organic matter (DOM). The determination of the quantity and reactivity of DOM that actually complexes Cd in soil solutions is challenging for operational purposes. Therefore, this study tested whether Cd²⁺ concentration in soil solutions could effectively be predicted by considering complexation with a single mean organic ligand assumed to be a fraction of DOM of unspecified nature or assumed to be purely fulvic acids (FA) with reactivity as described in WHAM VII. The reactivity of the unspecified ligand and the concentration of FA were modelled and fitted to experimental data from 76 agricultural soils with low Cd contents. The optimal reactivity and FA concentration that minimized the relative error (RE) of predictions of the concentration of Cd²⁺ in soil solutions were either considered constant across soils or modelled from soil properties by multiple linear regressions (MLR) or random forests (RF), giving 6 models, the predictive value of which was assessed by 10-folds cross-validation. When the reactivity of the mean ligand and the optimal FA concentration were considered constant across soils, the models were biased and 66.9% of predictions had relative errors below a factor of 2. By contrast, if the reactivity of the mean ligand or the optimal FA concentration were allowed to vary with soil characteristics, these performances increased to 95.5%, soil pH being the main predictor and RF being slightly more efficient than MLR. With more than 95% of the relative errors of prediction below a factor of 2, the models developed in this work could be valuable for assessing Cd speciation in the solution of soils having a low Cd content.

The impact of increased flooding occurrence on the mobility of potentially toxic elements in floodplain soil - A review

The frequency and duration of flooding events is increasing due to land-use changes increasing run-off of precipitation, and climate change causing more intense rainfall events. Floodplain soils situated downstream of urban or industrial catchments, which were traditionally considered a sink of potentially toxic elements (PTEs) arriving from the river reach, may now become a source of legacy pollution to the surrounding environment, if PTEs are mobilised by unprecedented flooding events. When a soil floods, the mobility of PTEs can increase or decrease due to the net effect of five key processes; (i) the soil redox potential decreases which can directly alter the speciation, and hence mobility, of redox sensitive PTEs (e.g. Cr, As), (ii) pH increases which usually decreases the mobility of metal cations (e.g. Cd²⁺, Cu²⁺, Ni²⁺, Pb²⁺, Zn²⁺), (iii) dissolved organic matter (DOM) increases, which chelates and mobilises PTEs, (iv) Fe and Mn hydroxides undergo reductive dissolution, releasing adsorbed and co-precipitated PTEs, and (v) sulphate is reduced and PTEs are immobilised due to precipitation of metal sulphides. These factors may be independent mechanisms, but they interact with one another to affect the mobility of PTEs, meaning the effect of flooding on PTE mobility is not easy to predict. Many of the processes involved in mobilising PTEs are microbially mediated, temperature dependent and the kinetics are poorly understood. Soil mineralogy and texture are properties that change spatially and will affect how the mobility of PTEs in a specific soil may be impacted by flooding. As a result, knowledge based on one river catchment may not be particularly useful for predicting the impacts of flooding at another site. This review provides a critical discussion of the mechanisms controlling the mobility of PTEs in floodplain soils. It summarises current understanding, identifies limitations to existing knowledge, and highlights requirements for further research.

The Effect of Flooding and Drainage Duration on the Release of Trace Elements from Floodplain Soils

Floodplains downstream of urban catchments are sinks for potentially toxic trace elements. An intensification of the hydrological cycle and changing land use will result in floodplains becoming inundated for longer durations in the future. We collected intact soil cores from a floodplain meadow downstream of an urban catchment and subjected them to an inundation/drainage cycle in the laboratory to investigate the effect of flood duration on trace element concentrations in the soil porewater. The porewater concentrations of Ni, Cr, and Zn increased, whereas Cu and Pb decreased with flood duration. All the Cr present in porewaters was identified as Cr(III). Copper concentrations increased after drainage but Pb mobility remained suppressed. Both pH and dissolved organic carbon (DOC) increased with flood duration but were lower in treatments that were drained for the longest duration (which were also the treatments flooded for the shortest duration). The porewater concentrations of Cr and Ni decreased after drainage to levels below those observed before inundation, mirroring the DOC concentrations. We concluded that the duration of floodplain inundation does have an influence on the environmental fate of trace elements but that flooding does not influence all trace elements in the same way. The implications of an intensification of the hydrological cycle over the coming decades are that floodplains may become a source of some trace elements to aquatic and terrestrial ecosystems. Environ Toxicol Chem 2020;39:2124-2135. © 2020 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

The impact of long-term biosolids application (>100 years) on soil metal dynamics

Biosolids application to arable land is a common, and cost-effective, practice but the impact of prolonged disposal remains uncertain. We evaluated the dynamics of potentially toxic elements (PTEs) at a long-established 'dedicated' sewage treatment farm. Soil metal concentrations exceeded regulations governing application of biosolids to non-dedicated arable land. However, measurement of isotopic exchangeability of Ni, Cu, Zn, Cd and Pb demonstrated support for the 'protection hypothesis' in which biosolids constituents help immobilise potential toxic metals (PTMs). Metal concentrations in a maize crop were strongly, and almost equally, correlated with all 'capacity-based' and 'intensity-based' estimates of soil metal bioavailability. This was attributable to high correlations between soil factors controlling bioavailability (organic matter, phosphate etc.) on a site receiving a single source of PTMs. Isotopic analysis of the maize crop suggested contributions to foliar Pb from soil dust originating from neighbouring fields. There was also clear evidence of metal-specific effects of biosolids on soil metal lability. With increasing metal concentrations there was both decreasing lability of Cd and Pb, due to interaction with increasing phosphate concentrations, and increasing lability of Ni, Cu and Zn due to weaker soil binding. Such different responses to prolonged biosolids disposal to arable soil should be considered when setting regulatory limits.

Soils in lakes: the impact of inundation and storage on surface water quality

The large-scale storage and inundation of contaminated soils and sediments in deep waterlogged former sand pits or in lakes have become a fairly common practice in recent years. Decreasing water depth potentially promotes aquatic biodiversity, but it also poses a risk to water quality as was shown in a previous study on the impact on groundwater. To provide in the urgent need for practical and robust risk indicators for the storage of terrestrial soils in surface waters, the redistribution of metals and nutrients was studied in long-term mesocosm experiments. For a range of surface water turbidity (suspended matter concentrations ranging from 0 to 3000 mg/L), both chemical partitioning and toxicity of pollutants were tested for five distinctly different soils. Increasing turbidity in surface water showed only marginal response on concentrations of heavy metals, phosphorus (P) and nitrogen (N). Toxicity testing with bioluminescent bacteria, and biotic ligand modelling (BLM), indicated no or only minor risk of metals in the aerobic surface water during aerobic mixing under turbid conditions. Subsequent sedimentation of the suspended matter revealed the chemical speciation and transport of heavy metals and nutrients over the aerobic and anaerobic interface. Although negative fluxes occur for Cd and Cu, most soils show release of pollutants from sediment to surface waters. Large differences in fluxes occur for PO₄, SO₄, B, Cr, Fe, Li, Mn and Mo between soils. For an indicator of aerobic chemical availability, dilute nitric acid extraction (0.43 M HNO₃; Aqua nitrosa) performed better than the conventional Aqua regia destruction. Both the equilibrium concentrations in surface waters, and fluxes from sediment, were adequately (r² = 0.81) estimated by a 1 mM CaCl₂ soil extraction procedure. This study has shown that the combination of 0.43 M HNO₃ and 1 mM CaCl₂ extraction procedures can be used to adequately estimate emissions from sediment to surface waters, and assess potential water quality changes, when former sand pits are being filled with soil materials.

Heavy metal behaviour at mineral-organo interfaces: Mechanisms, modelling and influence factors

The mineral-organo composites control the speciation, mobility and bioavailability of heavy metals in soils and sediments by surface adsorption and precipitation. The dynamic changes of soil mineral, organic matter and their associations under redox, aging and microbial activities further complicate the fate of heavy metals. Over the past decades, the wide application of advanced instrumental techniques and modelling has largely extended our understanding on heavy metal behavior within mineral-organo assemblages. In this review, we provide a comprehensive summary of recent progress on heavy metal immobilization by mineral-humic and mineral-microbial composites, with a special focus on the interfacial reaction mechanisms of heavy metal adsorption. The impacts of redox and aging conditions on heavy metal speciations and associations with mineral-organo complexes are discussed. The modelling of heavy metals adsorption and desorption onto synthetic mineral-organo composites and natural soils and sediments are also critically reviewed. Future challenges and prospects in the mineral-organo interface are outlined. More in-depth investigations are warranted, especially on the function and contribution of microorganisms in the immobilization of heavy metals at the complex mineral-organo interface. It has become imperative to use the state-of-the-art methodologies to characterize the interface and develop in situ analytical techniques in future studies.

Health risk assessment of potentially toxic elements in soils along the Central Elbe River, Germany

Floodplain soils across Central Elbe River (CER), Germany, vary considerably in potentially toxic element (PTE) content. However, there has never been a comprehensive study that links PTE levels with human health risk for children and adults. Our objective was to determine the contamination of 13 PTEs in 94 soil profiles along CER and assess the associated health risk via diverse indices for adults and children. Of 94 soil profiles, we measured soil properties and total content of arsenic, barium, chromium, copper, nickel, lead, rubidium, tin, strontium, vanadium, zinc, and zirconium using x-ray fluorescence spectrometer (XRF). We calculated the Contamination Factor and the Pollution Load Index (PLI), and assessed the health risk for male and female adults as well as for children. Topsoil median contents of Cr (84 mg kg^-1), Cu (42), Ni (33), and Zn (195) exceeded the Precautionary Values for sandy soils according to the German Federal Soil Protection and Contaminated Sites Ordinance, while As, Pb, and V were 32, 73, and 77 mg kg^-1, respectively. Median topsoil PLI was 1.73, indicating elevated multi-element contamination, with 90th percentile and maximum values being 3.20 and 4.31, respectively. All PTE concentrations were higher in top- compared to subsoils. Also at the 50th percentile the most enriched elements were Sn and As, followed by Zr and Rb, while in the 90th percentile Sn and As were followed by Zn, Pb and Cu. Median children's hazard index (HI) was higher than unity (HI = 2.27) and the 90th percentile was 5.53, indicating elevated health risk. Adult median HIs were 0.18 for male and 0.21 for female persons. Arsenic was found to be the primary contributor to total risk, accounting of 57.4% of HI in all three-person groupings, with Cr (17.3%) being the second, and V (10.2%) the third. Children's health is at dramatically higher risk than that of adults; also As, Cr, Pb, and V have a predominant role in contamination-related health risks. The presence of V, a less-expected element, among those of major risk contribution, reveals the necessity of monitoring areas at large scale. Our results demonstrate that our study may serve as a model for similar works studying multi-element-contaminated areas in future.

Potential ecological risks of metal(loid)s in riverine floodplain soils

The quality of soils under different land uses is getting deteriorated throughout the world due to various anthropogenic activities. This deterioration is highly complex in riverine floodplain areas due to contamination by multiple point and non-point sources and change in seasons. Therefore, a study was conducted to analyze seasonal (pre and post-monsoon) variations in physico-chemical characteristics, contents of metal(loid)s (Al, As, Cd, Cr, Co, Cu, Fe, Mn, Mo, Ni, Pb, Sb and Zn) in riverine floodplain soils under three land uses (agricultural, riverbank and roadside) from areas around the rivers Beas and Sutlej in Punjab, India. Further, analysis was done to assess the ecological and genotoxic risks (Allium cepa genotoxicity assay) posed by metal(loid)s in these soils. It was observed that soil samples under the three land uses were slightly alkaline (pre-monsoon) to acidic (post-monsoon) in nature with sandy texture and low soil organic matter. The levels of most metal(loid)s increased in post-monsoon soil samples under the three land uses, which was attributed to increase in soil organic matter, silt and clay contents in post-monsoon samples due to precipitation, flooding and sedimentation. The ecological Risk Index (58.3-104.5) and Modified Risk Index (145.2-178.9) calculated to analyze the level of ecological risks of metal(loid)s revealed that As, Cd and Sb posed moderate to considerable ecological risks in the agricultural and roadside soils in both seasons. Allium cepa genotoxicity assay indicated that the metal(loid)s in studied soils can cause genotoxic effects in biological systems. Therefore, various steps such as reduction in use of agrochemicals, promotion of organic agricultural methods and decontamination of soils using techniques such as phytoremediation etc must be taken to ensure reduction and containment of metal(loid)s in such riverine floodplain areas.

Multi-surface modeling of Ni(II) and Cd(II) partitioning in soils: Effects of salts and solid/liquid ratios

Metal partitioning in soils is a key process controlling metal bioavailability and mobility and is greatly influenced by the solid/liquid ratio. However, metal partitioning is difficult to describe either by a simple partition coefficient or by isotherm adsorption equations. This study investigated the solubility of Ni(II) and Cd(II) in 19 soils as a function of three extraction reagents (water, 0.01 M NaNO₃, and 0.01 M CaCl₂), five solid/liquid ratios (5-400 g/L) and field condition extracted by Rhizon samplers. Thermodynamically based multi-surface models (MSMs) that included generic parameters were used to describe metal partitioning under the studied conditions. The results showed that Ni/Cd solubility depended on the soil type, extraction reagent, and solid/liquid ratio. Soil major background cations (especially Ca²⁺, Mg²⁺, Fe³⁺ and Al³⁺) had a significant effect on the model's prediction ability. The MSM was able to predict the extractable metal in 0.01 M CaCl₂ in various soils at different solid/liquid ratios when soil background cations were included in the calculation; without the background cations, the model was able to predict metal partitioning only at solid/liquid ratios of <100 g/L. In addition, the model failed to predict water-extracted and 0.01 M-NaNO₃-extracted Ni/Cd if background cations were not included, but could reasonably do so if they were included. More importantly, after the background cations were included, MSMs relatively well predicted the Ni/Cd content in soil pore water under 80% field capacity conditions with water as the solution matrix.

A multi-surface model to predict Cd phytoavailability to wheat (Triticum aestivum L.)

The prediction of metal bioavailability in soils is critical to metal risk assessments and the amount of dissolved metal in soils is a key factor determining bioavailability. Because the recently developed geochemical multi-surface models (MSMs) offer a promising tool for the determination of metal partitioning in soils, in this study, a MSM based on generic parameters was used to assess the bioavailability of Cd in wheat (Triticum aestivum L.) growing in 12 soils with a wide range of properties. The amount of MSM-calculated dissolved Cd correlated strongly with the amount of Cd uptake by wheat (R²=0.873 for roots and R²=0.837 for shoots), and the model's performance was better than that of chemical extraction methods (0.01M CaCl₂, 0.43M HNO₃ and soil total Cd). The reactive fraction of soil organic matter, the soil/solution ratio, and the inclusion/exclusion of background cations influenced the calculation results. The best calculation condition was optimized. The application of the MSM was also examined in 84 wheat-soil samples from the field. The amount of Cd in wheat seeds had a stronger correlation with the amount of MSM-predicted Cd than with the amount of Cd obtained using chemical extraction methods. Our results suggested that MSM-calculated Cd is an effective indicator of the bioavailability of Cd in soils and demonstrated the utility of the method as a tool to assess the risk of Cd contamination in soils.

Multisurface modeling of Ni bioavailability to wheat (Triticum aestivum L.) in various soils

Continual efforts have been made to determine a simple and universal method of estimating heavy metal phytoavailability in terrestrial systems. In the present study, a mechanism-based multi-surface model (MSM) was developed to predict the partition of Ni(II) in soil-solution phases and its bioaccumulation in wheat (Triticum aestivum L.) in 19 Chinese soils with a wide range of soil properties. MSM successfully predicted the Ni(II) dissolution in 0.01 M CaCl₂ extracting solution (R² = 0.875). The two-site model for clay fraction improved the prediction, particularly for alkaline soils, because of the additional consideration of edge sites. More crucially, the calculated dissolved Ni(II) was highly correlated with the metal accumulation in wheat (R² = 0.820 for roots and 0.817 for shoots). The correlation coefficients for the MSM and various chemical extraction methods have the following order: soil pore water > MSM ≈ diffuse gradient technique (DGT) > soil total Ni > 0.43 M HNO₃ > 0.01 M CaCl₂. The results suggested that the dissolved Ni(II) calculated using MSM can serve as an effective indicator of the bioavailability of Ni(II) in various soils; hence, MSM can be used as an supplement for metal risk prediction and assessment besides chemical extraction techniques.

Modelling the potential mobility of Cd, Cu, Ni, Pb and Zn in Mollic Fluvisols

European floodplain soils are frequently contaminated with potentially toxic inorganic substances. We used a multi-surface model to estimate the aqueous concentrations of Cd, Cu, Ni, Pb and Zn in three Mollic Fluvisols from the Central Elbe River (Germany). The model considered complexation in solution and interactions with soil organic matter (SOM), a clay mineral and hydrous Al, Fe and Mn oxides. The amounts of reactive metals were derived from extraction with 0.43 M HNO₃. Modelling was carried out as a function of pH (soil pH ± 1.4) because it varies in floodplain soils owing to redox processes that consume or release protons. The fraction of reactive metals, which were dissolved according to the modelling, was predominantly <1%. Depending on soil properties, especially pH and contents of SOM and minerals of the clay fraction, the modelled concentrations partially exceeded the trigger values for the soil-groundwater pathway of the German soil legislation. This differentiation by soil properties was given for Ni, Pb and Zn. On the other hand, Cd was more mobile, i.e., the trigger values were mostly exceeded. Copper represented the opposite, as the modelling did not predict exceeding the trigger values in any horizon. Except for Pb and partially Zn (where oxides were more important), SOM was the most important adsorbent for metals. However, given the special composition and dynamics of SOM in mollic horizons, we suggest further quantitative and qualitative investigations on SOM and on its interaction with metals to improve the prediction of contaminant dynamics.