Extension of coupled multispecies metal transport and speciation (TRANSPEC) model to soil

Atmospheric deposition of metals emitted from mining operations has raised metal concentrations in the surrounding soils. This repository may be remobilized and act as a source of metals to nearby surface aquatic systems. It is important to understand metal dynamics and the impact of various chemistry and fate parameters on metal movement in the soil environment in order to evaluate risk associated with metals in terrestrial ecosystems and accurately establish critical discharge limits that are protective of aquatic biota. Here we extend our previously developed coupled multispecies metal fate-TRANsport and SPECiation/complexation (TRANSPEC) model, which was applicable to surface aquatic systems. The extended TRANSPEC, termed TRANSPEC-II, estimates the partition coefficient, K(d), between the soil-solid and -soluble phases using site-specific data and a semi-empirical regression model obtained from literature. A geochemical model calculates metal and species fractions in the dissolved and colloidal phases of the soil solution. The multispecies fugacity/aquivalence based fate-transport model then estimates inter-media transport rates such as leaching from soil, soil runoff, and water-sediment exchanges of each metal species. The model is illustratively applied to Ni in the Kelly Lake watershed (Sudbury, Ontario, Canada), where several mining operations are located. The model results suggest that the current atmospheric fallout supplies only 4% of Ni removed from soil through soil runoff and leaching. Soil runoff contributes about 20% of Ni entering into Kelly Lake with the rest coming from other sources. Leaching to groundwater, apart from runoff, is also a major loss process for Ni in the soil. A sensitivity analysis indicates that raising soil pH to above 6 may substantially reduce metal runoff and improve water quality of nearby water bodies that are impacted by runoff.