Modeling Soil Water and Solute Transport—Fast, Simplified Numerical Solutions

Abstract

Modeling of water and solute transport in soils is increasingly important in hydrology and agriculture, but there remain many gaps and unresolved issues. One of these, the speed and robustness of numerical solutions, is important in large‐scale and stochastic modeling. A fast method was developed for solution of the Richards equation for water transport. Brooks–Corey soil hydraulic property descriptions were used with water content as the dependent variable in unsaturated regions and pressure head in saturated regions. Central time weighting was used in unsaturated conditions to improve accuracy and fully implicit weighting in saturated conditions to improve stability. Water fluxes were calculated using matric flux potentials combined with a novel spatial weighting scheme for the gravitational component. Flow across soil property interfaces was calculated by equating fluxes to and from the interface. Results on a test problem involving rainfall, surface ponding, evaporation, and drainage over 400 h showed the new method to be an order of magnitude faster and more robust than an iterative scheme typifying current practice. Execution time for the test problem with <3% error was 0.006 s on a 166‐MHz personal computer. Solutions for solute transport described by the advection–dispersion equation were obtained concurrently without substantial increase in execution time by assuming average water fluxes over several steps of the water transport solution. The methods presented here should also be applicable to two‐ and three‐dimensional problems.