An approach is developed to simulate leaching of a dissolved chemical constituent in the vadose zone of an aquifer. Specifically, nitrate loading at the water table for different water table depths, for a range of aquifer permeability values, and for different cases of heterogeneity of the aquifer, are considered. Models from the literature are first used to derive soil-water characteristic curves (water retention and hydraulic conductivity) from a grain size distribution curve for unsaturated conditions. Given infiltration from the surface, the initial conditions for the chemical concentration, and the water content profile, leaching of the chemical in the vadose zone is simulated as a function of both time and depth. The methodology is illustrated for a permeable aquifer. Simulations are undertaken using a finite element code for saturated and unsaturated flow. Different scenarios are simulated depending on the heterogeneity of the aquifer and the depth of the water table. Modeling results show that in the example case studied, nitrate concentration loading at the water table does not depend strongly on the position of the water table, but rather on the material properties of the aquifer. The contribution of this endeavor resides in the methodology which allows a prediction of nitrate leaching using only the grain size property of the aquifer. It allows practitioners to obtain a first assessment of leaching with limited data.
