Seven unsaturated-zone solute-transport models were tested with two data sets to select models for use by the Agricultural Chemical Team of the U.S. Geological Survey's National Water-Quality Assessment Program. The data sets were from a bromide tracer test near Merced, California, and an atrazine study in the White River Basin, Indiana. In this study the models are designated either as complex or simple based on the water flux algorithm. The complex models, HYDRUS2D, LEACHP, RZWQM, and VS2DT, use Richards' equation to simulate water flux and are well suited to process understanding. The simple models, CALF, GLEAMS, and PRZM, use a tipping-bucket algorithm and are more amenable to extrapolation because they require fewer input parameters. The purpose of this report is not to endorse a particular model, but to describe useful features, potential capabilities, and possible limitations that emerged from working with the model input data sets. More rigorous assessment of model applicability involves proper calibration, which was beyond the scope of this study.

Uncalibrated ("cold") simulations were run using all seven models to predict the transport of bromide (Merced) and the transport and fate of atrazine and three of its transformation products (White River Basin). Among the complex models, HYDRUS2D successfully predicted both the surface retention and accumulation of bromide at depth at the Merced site, whereas RZWQM and VS2DT predicted only the latter. RZWQM predictions of atrazine were closest to observed values at the White River Basin site, where preferential flow has been observed. LEACHP predicted smaller solute concentrations than observed at both the Merced and White River Basin sites. Among the simple models, CALF predicted the highest values of atrazine and deethylatrazine at the measurement depth of 1.5 meters. CALF includes the Addiscott flow option for preferential flow, and also accepts user-specified dispersivity. PRZM underpredicted solute concentrations, probably because control of dispersion is a problem with this model. GLEAMS has a maximum simulation depth of 1.5 meters, which is limiting for mass-balance purposes because it creates a potential disconnect between unsaturated-zone transport and the water table.

Of the models tested, RZWQM, HYDRUS2D, VS2DT, GLEAMS and PRZM had graphical user interfaces. Extensive documentation was available for RZWQM, HYDRUS2D, and VS2DT. RZWQM can explicitly simulate water and solute flux in macropores, and both HYDRUS2D and VS2DT can simulate water and solute flux in two dimensions. The version of RZWQM tested had a maximum simulation depth of 3 meters. The complex models simulate the formation, transport, and fate of degradates of up to three to five compounds including the parent, with the exception of VS2DT, which simulates the transport and fate of a single compound.