Numerical simulation of heavy metal migration in stabilized soils under multi-field coupling

Abstract Heavy metal contamination of soils has become a global issue. In pursuit of a reference for the management of heavy metal pollution, numerical simulations were used here to investigate the migration pattern of heavy metals in soil under multi-field coupling of rainfall, temperature, and hydraulic conductivity. Models were built using TOUGHREACT, and measured and simulated values were compared with those of existing studies to verify the accuracy of the software in simulating coupled water vapor-thermal transport and rainfall infiltration. Taking Nanjing rainfall data as an example, the transport pattern of zinc in soil under multi-field coupling was explored. Although evaporation and rainwater infiltration were found to increase with the temperature, the effects of evaporation and temperature counteracted each other, resulting in a weak influence of temperature on heavy metal transport. The hydraulic conductivity of stabilized soil was the main factor controlling zinc migration, with the saturation and heavy metal concentration clouds of soil columns tending to stabilize when the hydraulic conductivity of solidified soil was under 1 × 10 -11 m/s.