Modeling the environmental fate of cadmium in a large wastewater irrigation area

The fate of cadmium in soils is governed by spatially heterogeneous processes that proceed from decades to centuries. This study aimed at modeling the fate of Cd within the wastewater irrigation area (WIA) of Braunschweig (Germany). The sandy soils (mainly Dystric Cambisol or Typic Haplumbrept) at this site (28 km2) have received considerable loads of heavy metals by irrigation of municipal wastewater for up to 40 yr. The soils of the WIA are in agricultural use. The main crops are sugar beet (Beta vulgaris L.), potato (Solanum tuberosum L.), and wheat (Triticum aestivum L.). As a result of asparagus (Asparagus officinalis L.) cropping, about 15% of the soils have been converted to Rigosols. In 1996, we measured the vertical distribution (0 to 1.2 m) of soil pH, organic carbon content, and the EDTA-extractable content and the solution phase concentration of Cd at 153 sites. At sites not used for asparagus cultivation, Cd has migrated on average to a depth of about 0.5 m. Due to deep plowing, which accelerates migration, Cd has been displaced on average to about 0.7 m at the Rigosol sites. To model the fate of Cd at the scale of the WIA, we used different parallel soil column approaches. In each column the local model SEFAH was used to simulate both displacement and plant uptake of Cd. The model was fed with measured or randomly generated soil data. The results of retrospective simulations from 1957 to 1996 agreed well with observed Cd profiles. The better the spatial variability of sorption was described, the better the performance. Our simulation results show that Cd pollution of soil at first affects the soil-plant pathway. The breakthrough of Cd to the groundwater is dampened and is delayed for many decades.