Prospective modeling with Hydrus-2D of 50 years Zn and Pb movements in low and moderately metal-contaminated agricultural soils.

Results of detailed modeling of in situ redistribution of heavy metals in pedological horizons of low and moderately metal contaminated soils, considering distinctly different long-term land use, are scarcely reported in literature. We used Hydrus-2D software parameterized with abundant available local soil data to simulate future Zn and Pb movements in soils contaminated by metallurgical fallout in the 20th century. In recent work on comparing different modeling hypotheses, we validated a two-site reactive model set with adjusted chemical kinetic constant values by fitting the 2005 Zn and Pb concentration profiles in soils, with estimated 1901-1963 airborne Zn and Pb loads (Mallmann et al., 2012a). In the present work, we used the same approach to simulate 2005-2055 changes in Zn and Pb depth-distribution and soil-solution concentrations, comparing two hypotheses of chemical equilibrium: i) the validated two-site model (one site at equilibrium and the other involved in kinetic reactions with pore water) set with adjusted kinetic EDTA extraction constants, and ii) a non-linear one-surface site adsorption equilibrium model. Simulated transfers were found generally lower and more realistic when using the two-site model. Simulations showed that consistent Zn redistribution and loss occurred in the moderately contaminated soil until 2055, i.e., more than one century after the main metal deposition, but negligible in low contaminated soils. Transfer of Pb was small in the three soils and under both hypotheses. In 2055, simulated Zn outflow concentrations remained under threshold values for drinking water.

Modeling field-scale vertical movement of zinc and copper in a pig slurry-amended soil in Brazil.

Organic amendments often represent a source of trace metals (TMs) in soils, which may partly leach into the groundwater. The objectives of this study were (1) to validate Hydrus-2D for modeling the transport of Zn and Cu in an Alfisol amended with pig slurry (PS) by comparing numerical simulations and experimental field data, and (2) to model the next 50 years of TM movements under scenarios of suspended or continued PS amendments. First, between 2000 and 2008, we collected detailed Zn and Cu data from a soil profile in Santa Maria, Brazil. Two hypotheses about Zn and Cu reactivity with the solid phase were tested, considering physical, hydraulic, and chemical characteristics of six soil layers. Using a two-site sorption model with a sorption kinetic rate adjusted based on laboratory EDTA extractions, Hydrus simulations of the vertical TM transport were found to satisfactorily describe the soil Zn and Cu concentration profiles. Second, the long-term fate of Zn and Cu in the soil was assessed using the validated parameterized model. Numerical simulations showed that Zn and Cu did not present risks for groundwater pollution. However, future Cu accumulation in the surface soil layer would exceed the Brazilian threshold for agricultural soils.

Using a two site-reactive model for simulating one century changes of Zn and Pb concentration profiles in soils affected by metallurgical fallout.

Predicting the transfer of contaminants in soils is often hampered by lacking validation of mathematical models. Here, we applied Hydrus-2D software to three agricultural soils for simulating the 1900-2005 changes of zinc and lead concentration profiles derived from industrial atmospheric deposition, to validate the tested models with plausible assumptions on past metal inputs to reach the 2005 situation. The models were set with data from previous studies on the geochemical background, estimated temporal metal deposition, and the 2005 metal distributions. Different hypotheses of chemical reactions of metals with the soil solution were examined: 100% equilibrium or partial equilibrium, parameterized following kinetic chemical extractions. Finally, a two-site model with kinetic constant values adjusted at 1% of EDTA extraction parameters satisfactory predicted changes in metal concentration profiles for two arable soils. For a grassland soil however, this model showed limited applicability by ignoring the role of earthworm activity in metal incorporation.