S imulation of solute fl uxes in complex terrain typically fails to account for subsurface lateral movement of solutes under transient saturated conditions. This type of solute fl ux could be important in phytoremediation (the use of plants to remove or neutralize contaminants in polluted soil or water), and more generally in evaluating the fate of solutes applied to soil. Solutes could move laterally and remain within the root zone, and thus could be taken up by the vegetation rather than leaching or moving offsite.A number of catchment-scale models have been developed that simulate transient saturated conditions in the vadose zone. The Variable Source Area (VSA) approach recognized that some areas within a catchment contribute to surface and subsurface runoff more than others, and that these areas change with time, expanding and contracting as a function of rain events and antecedent moisture (Hewlett and Hibbert, 1967). Simple assumptions have been made to represent VSA processes in models: declining saturated hydraulic conductivity with depth, steadystate catchment water table response, topographically defi ned water fl ow paths, and linear wetting and drying from the valley bottom upward to the ridge; examples of such models are the Soil Moisture Routing model (SMR) (Frankenberger et al., 1999) and TOPMODEL (Beven and Kirkby, 1979). The SMR represents the vadose zone using a water balance approach, where subsurface lateral fl ow is simulated using a simplifi ed Richards equation. TOPMODEL uses geographic similarity (i.e., contributing area and slope of soil surface) as a basis for modeling saturated zone catchment response; Shaman et al. (2002) added simulation of subsurface lateral fl ow using the same TOPMODEL assumptions. Band et al. (1993) coupled TOPMODEL and the ecosystem model BIOME-BGC to account for the impact of lateral fl uxes on ecological processes. In this coupled model, simulated water and solute fl uxes occur in the saturated zone rather than in the vadose zone. Wigmosta et al. (1994) introduced DHSVM, a semi-three-dimensional hydrology vegetation model that includes canopy interception, evaporation, transpiration, and snow accumulation. Whereas TOPMODEL was not spatially explicit, DHSVM was one of the fi rst spatially explicit models that simulated transient saturated conditions in the vadose zone.A phytoremediation project was established at the U.S. Department of Energy's Savannah River Site (South Carolina) to reduce fl uxes of tritium-contaminated groundwater to surface waters. Contaminated groundwater was collected in a pond and applied by spray irrigation to a catchment of mixed forest on Atlantic Coastal Plain soils. The objectives of this research project were to simulate tritium uptake by the vegetation and to determine if subsurface lateral fl ow at the sand-clay interface impacts tritium uptake by forest vegetation. To simulate water and tritium fl uxes within the catchment, we developed a spatially explicit water and solute transport model. Vertical water fl ow was simulated with...