Modeling Coupled Hydrologic and Chemical Processes: Long-Term Uranium Transport following Phosphorus Fertilization

• Published in: Vadose zone journal Vol. 7; no. 2; pp. 698 - 711
• Main Authors: Jacques, D, Simunek, J, Mallants, D, van Genuchten, M.T
• Format: Journal Article
• Language: English
• Published: Madison Soil Science Society of America 01.05.2008; Soil Science Society
• Subjects: actinides; adsorption; agriculture; alkaline earth metals; Belgium; bioavailability; biochemical processes; boundary conditions; calcium; case studies; chemical fractionation; chemical reactions; coupling; electrochemical properties; Environmental geology; Europe; fertilizer application; fertilizers; Flanders Belgium; geochemistry; hydrochemistry; hydrodynamics; isotopes; layered materials; Lommel Belgium; long term effects; mathematical models; metals; movement; northern Belgium; phosphorus; phosphorus fertilizers; pollutants; pollution; radioactive isotopes; retention; simulation; simulation models; soil depth; soil pH; soil profiles; soil surveys; soil transport processes; soils; solute transport; Spodosols; steady flow; surface layers; surveys; U-238; unsaturated zone; unsteady flow; uranium; vadose zone; Western Europe
• ISSN: 1539-1663; 1539-1663
• DOI: 10.2136/vzj2007.0084

Coupling physical and geochemical processes within one integrated numerical simulator provides a process-based tool for investigating the mobility of contaminants as affected by changing hydrologic regimes and geochemical conditions. We review interactions between physical and biogeochemical processes in the vadose zone, and then present a case study demonstrating these complex interactions. A hypothetical application is presented of the HP1 multicomponent transport simulator to predict the transport of two major elements (Ca and P) and one trace element (U) applied annually for 200 yr to a field soil in the form of an inorganic P fertilizer. Interactions of Ca, P, and U with the solid phase are described using cation exchange and surface complexation reactions. Simulations assuming steady-state or transient flow conditions were analyzed in terms of temporal variations of the linear distribution coefficient, K(d), which depends strongly on pH and the composition of the aqueous phase. If the composition of the aqueous phase is constant, adsorption of Ca and U increases with increasing pH. Due to the annual addition of Ca, P, and U, and competition between P and U for sorption sites, the K(d) of these elements decreased with time near the soil surface. Deeper in the soil, the K(d) of U followed the pH increase because of a lack of competition from P. Because of the combined effects of changing hydrologic and geochemical conditions, the Ca and U distribution coefficients and solute fluxes during the transient simulation exhibited large short-time variations of up to three orders of magnitude.