Abstract: A numerical model is developed for the simulation of vertical unsteady-state, two-phase reactive solute-transport through the unsaturated zone. Advective transport is simulated using the method of characteristics, whereas ...