Cadmium availability to wheat grain in soils treated with sewage sludge or metal salts

Grain Cd concentrations were determined in wheat (Triticum aestivum L.) grown in 1999, 2001 and 2003, at six sludge cake field experiments. Three of these sites also had comparisons with Cd availability from metal amended liquid sludge and metal salts. Grain Cd concentrations in all years and at all sites were significantly linearly correlated with NH4NO3 extractable Cd and soil total Cd (P<0.001). Soil extractability was greater in the liquid sludge and metal salt experiments than in the cake experiments, as were grain Cd concentrations. Across all the sites, NH4NO3 extractable soil Cd was no better at predicting grain Cd than soil total Cd. Stepwise multiple linear regression analysis showed that soil total Cd, pH and organic carbon were the only significant (P<0.001) variables influencing wheat grain Cd concentrations, explaining 78% of the variance across all field experiments (1408 plots). This regression predicted that the current UK soil total Cd limit of 3 mg kg(-1) was not sufficiently protective against producing grain above the European Union (EU) grain Cd Maximum Permissible Concentration (MPC) of 0.235 mg Cd kg(-1) dry weight, unless the soil pH was > 6.8. Our predictions show that grain would be below the MPC with > 95% confidence with the proposed new EU draft regulations permitting maximum total Cd concentrations in soils receiving sludge of 0.5 mg kg(-1) for soils of pH 5-6, 1 mg kg(-1) for soils of pH 6-7, and 1.5 mg kg(-1) for soils of pH > or = 7.     

A reaction-based paradigm to model reactive chemical transport in groundwater with general kinetic and equilibrium reactions

This paper presents a reaction-based water quality transport model in subsurface flow systems. Transport of chemical species with a variety of chemical and physical processes is mathematically described by M partial differential equations (PDEs). Decomposition via Gauss-Jordan column reduction of the reaction network transforms M species reactive transport equations into two sets of equations: a set of thermodynamic equilibrium equations representing N(E) equilibrium reactions and a set of reactive transport equations of M-N(E) kinetic-variables involving no equilibrium reactions (a kinetic-variable is a linear combination of species). The elimination of equilibrium reactions from reactive transport equations allows robust and efficient numerical integration. The model solves the PDEs of kinetic-variables rather than individual chemical species, which reduces the number of reactive transport equations and simplifies the reaction terms in the equations. A variety of numerical methods are investigated for solving the coupled transport and reaction equations. Simulation comparisons with exact solutions were performed to verify numerical accuracy and assess the effectiveness of various numerical strategies to deal with different application circumstances. Two validation examples involving simulations of uranium transport in soil columns are presented to evaluate the ability of the model to simulate reactive transport with complex reaction networks involving both kinetic and equilibrium reactions.